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

[Roland Schweitzer <Roland.Schweitzer@xxxxxxxx>]

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>
>  
>  
>  
>  
>  
>  
>  
>