Re: Monthly time element in a generated xml by addXml.pl

[Joe McLean <joe.mclean@xxxxxxxx>]

Matsutani-san,

Two comments:
1.  Adding the attribute:
     modulo = " " ;
     to the time coordinate variable in the netcdf file will also solve your problem.

2.  the XML tag <salt_anal1degmon_nc_time type="t" units="days"> should be
     <salt_anal1degmon_nc_time type="t" units="month">

Joe
----------------------------------------


Roland Schweitzer wrote:

> Hirohisa,
>
> Try this.
>
> Tell Ferret that the time axis is a climatological axis.
>
> Change the variable declaration to:
> <salt name="Ocean Salinity, analyzed mean, 1-deg grid, Monthly"
> units="g/kg">
> <link match="/lasdata/grids/salt_anal1degmon_nc_time_level_lat_lon_grid"/>
> <properties>
>      <ferret>
>            <init_script>init_modulo_t</init_script>
>       </ferret>
>   </properties>
> </salt>
>
> The time axis arange should look like:
>
> <arange start="0001-01-01 00" step="1" size="12" />
>
> Roland
>
> Matsutani Hirohisa wrote:
> > Hello,
> >
> > I am in a difficulty that
> > I can't display the correct monthly data in Las.
> > Please give me a good advise.
> >
> >
> > My data is as follows
> >
> > ・瘢雹ncdump -h <file>
> > ---------------------------------------------------------------
> >
> > netcdf salt.anal1degmon {
> > dimensions:
> >  lon = 360 ;
> >  lat = 180 ;
> >  level = 24 ;
> >  time = UNLIMITED ; // (12 currently)
> > variables:
> >  float lat(lat) ;
> >   lat:units = "degrees_north" ;
> >   lat:long_name = "Latitude" ;
> >   lat:actual_range = 89.5f, -89.5f ;
> >  float lon(lon) ;
> >   lon:units = "degrees_east" ;
> >   lon:long_name = "Longitude" ;
> >   lon:actual_range = 0.5f, 359.5f ;
> >  float level(level) ;
> >   level:units = "meters" ;
> >   level:positive = "down" ;
> >   level:long_name = "Level" ;
> >   level:actual_range = 0.f, 1500.f ;
> >  double time(time) ;
> >   time:units = "days since 1-1-1 00:00:0.0" ;
> >   time:long_name = "Time" ;
> >   time:actual_range = 0., 334. ;
> >   time:delta_t = "0000-01-00 00:00:00" ;
> >   time:avg_period = "0097-00-00 00:00:00" ;
> >   time:prev_avg_period = "0000-01-00 00:00:00" ;
> >   time:ltm_range = 693597., 729026. ;
> >  float salt(time, level, lat, lon) ;
> >   salt:long_name = "Ocean Salinity, analyzed mean, 1-deg grid, Monthly" ;
> >   salt:valid_range = 0.f, 100.f ;
> >   salt:actual_range = 3.5243f, 41.4294f ;
> >   salt:units = "g/kg" ;
> >   salt:add_offset = 0.f ;
> >   salt:scale_factor = 1.f ;
> >   salt:missing_value = -9.96921e+36f ;
> >   salt:var_desc = "Ocean Salinity, analyzed mean\n",
> >     "BS" ;
> >   salt:dataset = "NODC World Ocean Atlas 1998\n",
> >     "AU" ;
> >   salt:level_desc = "Multiple Levels\n",
> >     "F" ;
> >   salt:statistic = "Analyzed Mean\n",
> >     "AM" ;
> >   salt:parent_stat = "Mean\n",
> >     "M" ;
> >
> > // global attributes:
> >   :title = "NODC World Ocean Atlas 1998" ;
> >   :history = "Mon Jul 17 16:39:02 2000: ncrcat -c -d level,0,23 -v
> > lat,lon,level,time,salt salt.anal1deg.nc /data/hoop/salt.anal1deg.nc\n",
> >     "created February 2000 by Hoop" ;
> >   :platform = "Marine Analyses" ;
> >   :Conventions = "COARDS" ;
> > }
> > ------------------------------------------------------------------------------
> >
> > ・瘢雹the generated xml(by　addXml.pl)
> > -----------------------------------------------------------------------------------------------------
> > <datasets>
> >  <salt_anal1degmon_nc name="NODC World Ocean Atlas 1998"
> > url="http://localhost/cgi-bin/nph-nc/data/salt.anal1degmon.nc"; doc="">
> >   <variables>
> >    <salt name="Ocean Salinity, analyzed mean, 1-deg grid, Monthly"
> > units="g/kg">
> >     <link
> > match="/lasdata/grids/salt_anal1degmon_nc_time_level_lat_lon_grid"/>
> >    </salt>
> >   </variables>
> >  </salt_anal1degmon_nc>
> > </datasets>
> > <grids>
> >  <salt_anal1degmon_nc_time_level_lat_lon_grid>
> >   <link match="/lasdata/axes/salt_anal1degmon_nc_time"/>
> >   <link match="/lasdata/axes/salt_anal1degmon_nc_level"/>
> >   <link match="/lasdata/axes/salt_anal1degmon_nc_lat"/>
> >   <link match="/lasdata/axes/salt_anal1degmon_nc_lon"/>
> >  </salt_anal1degmon_nc_time_level_lat_lon_grid>
> > </grids>
> > <axes>
> >  <salt_anal1degmon_nc_time type="t" units="days">
> >      <!-- My try setting arange start="1-1-16" step="1" size="12"/  my
> > try!!!  -->
> >   <v>Jan</v>
> >   <v>Feb</v>
> >   <v>Mar</v>
> >   <v>Apr</v>
> >   <v>May</v>
> >   <v>Jun</v>
> >   <v>Jul</v>
> >   <v>Aug</v>
> >   <v>Sep</v>
> >   <v>Oct</v>
> >   <v>Nov</v>
> >   <v>Dec</v>
> >  </salt_anal1degmon_nc_time>
> >
> >  <salt_anal1degmon_nc_level type="z" units="meters">
> >   <v>0
> >   </v>
> >   <v>10
> >   </v>
> > .........
> >   <v>1500
> >   </v>
> >  </salt_anal1degmon_nc_level>
> >  <salt_anal1degmon_nc_lat type="y" units="degrees_north">
> >   <arange start="89.5" step="-1" size="180"/>
> >  </salt_anal1degmon_nc_lat>
> >  <salt_anal1degmon_nc_lon type="x" units="degrees_east">
> >   <arange start="0.5" step="1" size="360"/>
> >  </salt_anal1degmon_nc_lon>
> > </axes>
> > -------------------------------------------------------------------------
> >
> > In Las's operation ,I did
> >
> > Select time : 1-Jan and next key
> > others set are  default
> >
> > ,then the following error is send
> >
> > -------------------------------------------------
> > Illeagal Limits
> > There are nodata available at 01-JAN  00:00
> > The data beign at 16-DEC 12:00 and continued at 16-DEC 00:00
> > ----------------------------------------------------
> >
> > Why is  such the data range set?
> >
> > Please teach me how to change the xml or...??? Help me!
> >
> > ,where
> > I found the second best plan that
> > set units ="index" and <v>1</v>...<v>12</v>
> > in the above xml.
> > The result is ok in such a way but  I can't recognize the monthly name.
> >
> > regars.
> >
> > Hirohisa Matsutani
> > e-mail: hmatsutani@nifty.com <mailto:hmatsutani@nifty.com>
> >
> >
> >
> >
> >
> >
> >
> >