Changeset 542

Timestamp:

04/30/12 12:49:45 (12 years ago)

Author:

pinsard

Message:

fix some time issue in netcdf writing tool

File:

• trunk/src/idl_amsu_netcdf.pro (modified) (11 diffs)

trunk/src/idl_amsu_netcdf.pro

@@ -43 +43 @@
 ; ========
 ;
-; To produce :file:`${PROJECT_OD}/AMSU/a5_20060801_20060812_01000w2000n_00100w3000n.nc`::
+; To produce
+; :file:`${PROJECT_OD}/AMSU/a5_20060801_20060812_01000w2000n_00100w3000n.nc`::
 ;
 ;    numch = 'a5'
@@ -52 +53 @@
 ;    idl_amsu_netcdf, numch , datain, lats, lons, timesin
 ;
+; TIPS
+; ====
+;
+; File produced by can be read using :func:`read_ncdf <saxo:read_ncdf>` ::
+;
+;  tb=read_ncdf('moyenne_tb', 20060801, 20060802, file='/homedata/pinsard/varamma_d/AMSU/a5mrt0.5rh1.0_20060801_20060804_02400w0400s_02400e1900n.nc')
+;
 ; SEE ALSO
 ; ========
@@ -65 +73 @@
 ;
 ;
+; ++ test si toutes les valeurs sont des valeurs manquantes/aberantes
+; si ou ne rien écrire
+;
+; improve tip
+;
 ; check arg : time in julian days. check for noon !!
 ;
@@ -85 +98 @@
 ;
 ; $URL$
+;
+; - fplod 20120430
+;
+;   * time unlimited
+;   * time origine is a ttribute of time not of tb
+;   * fix double vs fltarr
+;   * fix missing value
+;
+;   nb : now cdo showdate and showtime work
 ;
 ; - lelod 20120410
@@ -134 +156 @@
 ;
 ;-
-pro idl_amsu_netcdf, numch, temp,rtemp,rhdeg, datain, lats, lons, timesin
+pro idl_amsu_netcdf, numch, temp, rtemp, rhdeg, datain, lats, lons, timesin
 ;
 compile_opt idl2, strictarrsubs
@@ -149 +171 @@
 rt=string(rtemp,format='(f3.1)')
 rh=string(rhdeg,format='(f3.1)')
-prefix = numch + temp + 'rt'+rt+'rh'+rh +'_'
-
+prefix = numch + temp + 'rt' + rt + 'rh' + rh + '_'
+;
 look = 'filename'
 scale = 100.
@@ -166 +188 @@
         + '.nc'
 ;
-; definition des tableaux qui vont contenir les valeurs des fichiers ascii
+; definition des tableaux qui vont contenir les valeurs à écrire
+; dans les fichiers NetCDF
 moyenne_tb = fltarr(nblon, nblat, nbtime)
 xlon = fltarr(nblon)
@@ -184 +207 @@
 print,"nb valeurs manquantes", nnbad
 ;
+; ++ test si toutes les valeurs sont des valeurs manquantes/aberantes
+; si ou ne rien écrire
+;
 ;
 ;---------------------------------------
@@ -197 +223 @@
 print, 'ouverture pour ecriture de ', cdfout
 id = ncdf_create(cdfout,/clobber)
-;++ncdf_control, id, /nofill
+ncdf_control, id, /nofill
+;
 ; dimension
 dimidx = ncdf_dimdef(id, 'xlon', nblon)
 dimidy = ncdf_dimdef(id, 'ylat', nblat)
-dimidt = ncdf_dimdef(id, 'time', nbtime)
+dimidt = ncdf_dimdef(id, 'time', /UNLIMITED)
 ;
 ; attributs globaux
@@ -209 +236 @@
 creation_date = string(jul2date(SYSTIME(/JULIAN, /UTC )), format='(I8)')
 ncdf_attput, id, 'Creation_date', creation_date, /GLOBAL
-;print,creation_date
+;print, creation_date
 ;
 ; declaration des variables
 varid = lonarr(4)
-varid[0] = ncdf_vardef(id, 'xlon', [dimidx], /double)
-varid[1] = ncdf_vardef(id, 'ylat', [dimidy], /double)
-varid[2] = ncdf_vardef(id, 'time', [dimidt], /double)
+varid[0] = ncdf_vardef(id, 'xlon', [dimidx], /float)
+varid[1] = ncdf_vardef(id, 'ylat', [dimidy], /float)
+varid[2] = ncdf_vardef(id, 'time', [dimidt], /float)
 varid[3] = ncdf_vardef(id, 'moyenne_tb',[dimidx,dimidy,dimidt], /float)
 ;
@@ -221 +248 @@
 ncdf_attput, id, varid[0], 'units', 'degrees_east'
 ncdf_attput, id, varid[0], 'title', 'Longitude'
+;
 ncdf_attput, id, varid[1], 'units', 'degrees_north'
 ncdf_attput, id, varid[1], 'title', 'Latitude'
+;
 ncdf_attput, id, varid[2], 'units', 'days since 2000-01-01 00:00:00'
-ncdf_attput, id, varid[2], 'long_name', 'Time axis'
-ncdf_attput, id, varid[3], 'time_origin', '2000-01-01 00:00:00'
+ncdf_attput, id, varid[2], 'long_name', 'time'
+ncdf_attput, id, varid[2], 'short_name', 'time'
+ncdf_attput, id, varid[2], 'axis', 'T'
+ncdf_attput, id, varid[2], 'time_origin', '2000-01-01 00:00:00'
+;
 ncdf_attput, id, varid[3], 'short_name', 'amsu_tb'
-ncdf_attput, id, varid[3], 'missing_value', '-9999.00'
+ncdf_attput, id, varid[3], 'missing_value', -9999.00
 ncdf_attput, id, varid[3], 'unit', 'degree_kelvin'
 ;