[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[IDV #YIX-913245]: Error loading 3D isosurface

> There appears to still be an issue.  It now plots the sigma based
> extinction contours, however they seem to be computed incorrectly.
> Towards the surface levels, contours seem to be inverted.  Comparison
> images are attached - one using the sigma coordinates, one with
> arbitrary altitude coords.
> 
> Let me know if the problem may be on my end.
> 
> Walter

Walter,
       I understand your results, and I would say both images give the correct 
displays. The reason you see the inverted contours is that in the sigma coord. 
display the sigma coordinate being transferred to pressure coordinate, and 
resulting a 3D surface displaying in the IDV, but the arbitrary coords is a 2D 
surface, you can tilt the view point to get the feeling of it. Otherwise, the 
color shaded plan view may be better option for you. The two displays from 
different coordinates will have a little difference. Let me know if you see it 
differently, and remember to attach a bundle with your email so I can reproduce 
your display locally.


Yuan
> 
> On Fri, Jun 11, 2010 at 1:34 PM, Unidata IDV Support
> <address@hidden> wrote:
> >> > I've uploaded the file - let me know if you find the problem.  Thanks
> >> > for your time!
> >> >
> >> > Walter
> >>
> >> Walter,
> >> I think it is a bug in the IDV when transforming the vertical coordinate, 
> >> I am working on it.
> >>
> >>
> >
> > Walter,
> >
> > I fixed the problem on the IDV and the local test works fine, you can 
> > install a nightly release tomorrow, let me know if
> > you see any issue again.
> >
> >
> >
> > Yuan
> >>
> >> Yuan
> >> >
> >> > On Thu, Jun 10, 2010 at 1:05 PM, Unidata IDV Support
> >> > <address@hidden> wrote:
> >> > >> Full Name: W. Session
> >> > >> Email Address: address@hidden
> >> > >> Organization: FSU Dept. of Meteorology
> >> > >> Package Version: 2.8u1 build date:2010-02-02 22:30 UTC
> >> > >> Operating System: Linux
> >> > >> Hardware: Java: home: /users/sessions/prgms/IDV_2.8u1/jre version: 
> >> > >> 1.6.0_17 j3d:1.5.2 fcs (build4)
> >> > >> Description of problem: When I attempt to create a 3D isosurface, I 
> >> > >> get the following error:
> >> > >>
> >> > >> Couldn't get data
> >> > >> SampledSet.init_samples: sample#2 length 1620000 doesn't match 
> >> > >> expected length 1624500
> >> > >>
> >> > >> The netcdf file is a lat/lon grid (-179.5-179.5 W/E, -89.5-89.5 S/N   
> >> > >> 1x1 deg) with 25 vertical sigma levels.  The 3D fields I'm interested 
> >> > >> in are pressure and extinction (sigma x lat x lon).  The file opens 
> >> > >> fine in NCVIEW, NCL, and IDL, so this may be a CF standards issue.
> >> > >>
> >> > >> IDV 2.8u1 on RH Linux.
> >> > >>
> >> > > Hi Walter,
> >> > >   If the dataset follows the CF conventions, the IDV should be able to 
> >> > > display it. Could you upload a sample file here:
> >> > > http://motherlode.ucar.edu/repository/alias/idvupload
> >> > >
> >> > > I would like to trace the problem locally.
> >> > >
> >> > >
> >> > > Yuan
> >> > >> Here is the ncdump header information for the file, if this helps:
> >> > >>
> >> > >> netcdf farop {
> >> > >> dimensions:
> >> > >> lon = 360 ;
> >> > >> lat = 180 ;
> >> > >> sigma = 25 ;
> >> > >> variables:
> >> > >> float longitude(lon) ;
> >> > >> longitude:long_name = "longitude" ;
> >> > >> longitude:standard_name = "longitude" ;
> >> > >> longitude:units = "degrees_east" ;
> >> > >> float latitude(lat) ;
> >> > >> latitude:long_name = "latitude" ;
> >> > >> latitude:standard_name = "latitude" ;
> >> > >> latitude:units = "degrees_north" ;
> >> > >> float pressure(sigma, lat, lon) ;
> >> > >> pressure:long_name = "pressure" ;
> >> > >> pressure:units = "hPa" ;
> >> > >> pressure:positive = "down" ;
> >> > >> pressure:standard_name = "air_pressure" ;
> >> > >> float surface_pressure(lat, lon) ;
> >> > >> surface_pressure:long_name = "surface pressure" ;
> >> > >> surface_pressure:units = "hPa" ;
> >> > >> surface_pressure:coordinates = "longitude latitude" ;
> >> > >> float aerosol_extinction(sigma, lat, lon) ;
> >> > >> aerosol_extinction:long_name = "aerosol extinction" ;
> >> > >> aerosol_extinction:units = "1/m" ;
> >> > >> aerosol_extinction:coordinates = "pressure longitude latitude" ;
> >> > >> float ptop ;
> >> > >> ptop:long_name = "pressure at top of model" ;
> >> > >> ptop:units = "hPa" ;
> >> > >> float sigma(sigma) ;
> >> > >> sigma:long_name = "sigma at layer midpoints" ;
> >> > >> sigma:standard_name = "atmosphere_sigma_coordinate" ;
> >> > >> sigma:formula_terms = "sigma: sigma ps: surface_pressure ptop: ptop" ;
> >> > >>
> >> > >> // global attributes:
> >> > >> :Conventions = "CF-1.0" ;
> >> > >> }
> >> > >>
> >> > >>
> >> > >> ******************
> >> > >> Stack trace:
> >> > >> visad.SetException: SampledSet.init_samples: sample#2 length 1620000 
> >> > >> doesn't match expected length 1624500
> >> > >> at visad.SampledSet.init_samples(SampledSet.java:232)
> >> > >> at visad.GriddedSet.<init>(GriddedSet.java:78)
> >> > >> at visad.Gridded3DSet.<init>(Gridded3DSet.java:83)
> >> > >> at visad.GriddedSet.create(GriddedSet.java:218)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.makeDomainWithVerticalTransform(GeoGridAdapter.java:1018)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.getSpatialDomainSetInner(GeoGridAdapter.java:441)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.getSpatialDomainSet(GeoGridAdapter.java:356)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.getFlatField(GeoGridAdapter.java:1143)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.makeSequence(GeoGridAdapter.java:1332)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridAdapter.getSequence(GeoGridAdapter.java:1087)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridDataSource.makeFieldImpl(GeoGridDataSource.java:1708)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GeoGridDataSource.getDataInner(GeoGridDataSource.java:1280)
> >> > >> at ucar.unidata.data.DataSourceImpl.getData(DataSourceImpl.java:2235)
> >> > >> at 
> >> > >> ucar.unidata.data.DirectDataChoice.getData(DirectDataChoice.java:332)
> >> > >> at ucar.unidata.data.DataChoice.getData(DataChoice.java:637)
> >> > >> at ucar.unidata.data.DataInstance.getData(DataInstance.java:239)
> >> > >> at ucar.unidata.data.DataInstance.getData(DataInstance.java:203)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GridDataInstance.init(GridDataInstance.java:202)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GridDataInstance.<init>(GridDataInstance.java:159)
> >> > >> at 
> >> > >> ucar.unidata.data.grid.GridDataInstance.<init>(GridDataInstance.java:140)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.GridDisplayControl.doMakeDataInstance(GridDisplayControl.java:299)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.DisplayControlImpl.initializeDataInstance(DisplayControlImpl.java:3040)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.DisplayControlImpl.setData(DisplayControlImpl.java:3022)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.ThreeDSurfaceControl.setData(ThreeDSurfaceControl.java:249)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.ThreeDSurfaceControl.init(ThreeDSurfaceControl.java:155)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.DisplayControlImpl.init(DisplayControlImpl.java:1302)
> >> > >> at 
> >> > >> ucar.unidata.idv.control.DisplayControlImpl.init(DisplayControlImpl.java:1003)
> >> > >> at 
> >> > >> ucar.unidata.idv.ControlDescriptor.initControl(ControlDescriptor.java:987)
> >> > >> at 
> >> > >> ucar.unidata.idv.ControlDescriptor$1.run(ControlDescriptor.java:912)
> >> > >> at ucar.unidata.util.Misc$3.run(Misc.java:1096)
> >> > >>
> >> > >>
> >> > >>
> >> > >
> >> > >
> >> > > Ticket Details
> >> > > ===================
> >> > > Ticket ID: YIX-913245
> >> > > Department: Support IDV
> >> > > Priority: Normal
> >> > > Status: Open
> >> > >
> >> > >
> >> >
> >> >
> >>
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: YIX-913245
> > Department: Support IDV
> > Priority: Normal
> > Status: Closed
> >
> >
> 
> 


Ticket Details
===================
Ticket ID: YIX-913245
Department: Support IDV
Priority: Normal
Status: Open