[Thread Prev][Thread Next][Index]

[las_users] wrong plots of data with curvilinear grids




LAS 7.0.2    tomcat 6     ferret 6.1

problem: A netcdf file with curvilinear coordinates is properly plotted by Ferret, but wrongly by LAS.

LAS plots only the lower half of the (lat/lon) data, leaving the upper half blank, as if there were missing values.

Probably I wrote a wrong xml file, but I did my best. If the xml is right, then may be there are problems in LAS ferret/scripts dealing with curvilinear transformations (may be LAS_transforms or transform_curvi_xy ?).
I'm not able to solve the problem on my own, thus I ask for your help.

You can access LAS at
              https://clima.casaccia.enea.it:8443/las
and then select the file "ernesto"

I prefere not to attach any file, thus please find at
              ftp://clima.casaccia.enea.it/emanuele.lombardi
the following files:
ferret.gif = the output of ferret shade/k=1/l=1 votemper,nav_lon,nav_lat
              LAS.gif        =  the outut of LAS
              ernesto.nc   =  the netcdf file
              ernesto.xml =  the xml I use in LAS
              ernesto.cdl  =  the "ncdump -h" of the netcdf file

Thanks for any help,
Emanuele

Here it is the "ncdump -h" of the netcdf file

netcdf ernesto {
dimensions:
        deptht = 2 ;
        y = 376 ;
        x = 360 ;
        time_counter = UNLIMITED ; // (2 currently)
variables:
        float deptht(deptht) ;
                deptht:long_name = "Vertical levels" ;
                deptht:units = "m" ;
                deptht:title = "deptht" ;
        float nav_lat(y, x) ;
                nav_lat:units = "degrees_north" ;
                nav_lat:long_name = "Latitude" ;
                nav_lat:nav_model = "Default grid" ;
        float nav_lon(y, x) ;
                nav_lon:units = "degrees_east" ;
                nav_lon:long_name = "Longitude" ;
                nav_lon:nav_model = "Default grid" ;
        float time_counter(time_counter) ;
time_counter:units = "days since 2009-01-20 00:00:00" ;
                time_counter:calendar = "gregorian" ;
                time_counter:title = "Time" ;
                time_counter:long_name = "Time axis" ;
                time_counter:time_origin = "2009-01-20 00:00:00" ;
        float votemper(time_counter, deptht, y, x) ;
                votemper:units = "C" ;
                votemper:missing_value = 0.f ;
                votemper:valid_min = 10.f ;
                votemper:valid_max = 35.f ;
                votemper:long_name = "Temperat" ;
                votemper:short_name = "votemper" ;
                votemper:online_operation = "ave(only(x))" ;
                votemper:axis = "TZYX" ;
                votemper:interval_operation = 600.f ;
                votemper:interval_write = 3600.f ;
votemper:associate = "time_counter deptht nav_lat nav_lon" ;

// global attributes:
                :Conventions = "GDT 1.3" ;
                :file_name = "    " ;
                :production = "ENEA CLIMMOD" ;
                :TimeStamp = "    " ;
                :associate_file = "    " ;
:history = "Fri Jan 23 22:24:52 2009: ncks -v votemper,nav_lon,nav_lat -d time_counter,0,1 -d deptht,0,1 TIRRfpE090120_01_grid_T.nc ernesto.nc" ;
}


And this is the xml I wrote for the above netcdf file:

<datasets>
  <ciccio name="ernesto1" url="/home/bancadat/ernesto.nc">
  <properties>
   <ferret>
       <curvi_coord_lon>nav_lon</curvi_coord_lon>
       <curvi_coord_lat>nav_lat</curvi_coord_lat>
   </ferret>
  </properties>
  <variables>
   <votemper name="Potential Temperature" units="C">
     <link match="/lasdata/grids/grid_xyzt-ciccio"/>
   </votemper>
  </variables>
 </ciccio>
</datasets>
<grids>
 <grid_xyzt-ciccio>
  <link match="/lasdata/axes/x-ciccio"/>
  <link match="/lasdata/axes/y-ciccio"/>
  <link match="/lasdata/axes/z-ciccio"/>
  <link match="/lasdata/axes/t-ciccio"/>
 </grid_xyzt-ciccio>
</grids>
<axes>
 <x-ciccio type="x" units="degrees_east">
  <arange start="8.8125" size="360" step="0.0217" />
 </x-ciccio>
 <y-ciccio type="y" units="degrees_north">
  <arange start="36.688" size="376" step="0.0208" />
 </y-ciccio>
 <z-ciccio type="z" units="m">
   <v>1.47</v>
   <v>4.55</v>
 </z-ciccio>
 <t-ciccio type="t" units="t">
  <v>27-jan-2009:00:00:00</v>
  <v>27-jan-2009:00:59:59</v>
 </t-ciccio>
</axes>








[Thread Prev][Thread Next][Index]


Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement