Changeset 28 for trunk/src

Timestamp:

01/04/11 12:34:16 (13 years ago)

Author:

pinsard

Message:

start to work on correction tools

File:

• trunk/src/t2m_correction_ncdf.pro (modified) (8 diffs)

trunk/src/t2m_correction_ncdf.pro

@@ -6 +6 @@
 ; t2m_correction_ncdf.pro
 ; =======================
+;
+; Mean correction for air temperature bias and correction for variability are
+; applied.
+;
+; :file:`${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc` have been produced by :ref:`interp_erai_t2m_1989_2009.pro`.
+; It contains air temperature at 2 m height from ERA-I interpolated on OAFLUX grid.
+;
+; Corrected air temperature at 2 m height is written in
+; :file:`${TROPFLUX_OD}/TropFlux_t2m_19890101_20091231_v50.nc`.
 ;
 ;     .. graphviz::
@@ -14 +23 @@
 ;           ]
 ;
-;           file_t2m [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc"];
-;
-;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_t2m_19890101_20091231_v20.nc"];
+;           file_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"];
+;
+;           ncfile [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/TropFlux_t2m_19890101_20091231_v50.nc"];
 ;
 ;           t2m_correction_ncdf [shape=box,
@@ -30 +39 @@
 ; SEE ALSO
 ; ========
+;
+; :ref:`tropflux_profile.sh`
+;
+; :ref:`mooring_corrections`
+;
+; :ref:`interp_erai_t2m_1989_2009.pro`
 ;
 ; :func:`initncdf <saxo:initncdf>`
@@ -38 +53 @@
 ; :func:`jul2date <saxo:jul2date>`
 ; :func:`ncdf_quickwrite <saxo:ncdf_quickwrite>`
+; :func:`ncdf_getatt <saxo:ncdf_getatt>`
 ;
 ; EXAMPLES
@@ -49 +65 @@
 ; ====
 ;
-; hard coded directory - usage of ${TROPFLUX_ID}
+; make it work !!
 ;
 ; coding rules
+;
+; understand usage of jtt
+;
+; check time values
+;
+; hard coded time in module name and in output filename
+;
+; why v50 in filename output
+;
+; hard coded correction values
+;
+; check side effect of replacement of read_ncdf by ncdf_lec ::
+;
+;    Probleme d'adequation entre les tailles du domaine nx*ny*jpt 350*60*1 et du tableau 350*60*7670
+;
+;
+; hard coded attributes
+;
+; CF conventions
+;
+; KNOWN ISSUES
+; ============
+;
+; test of existence of fullfilename_msk not very efficient because
+; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented
 ;
 ; EVOLUTIONS
 ; ==========
 ;
+; - fplod 20110104T093758Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * complete header
+;   * replace /Volumes/Iomega_HDD/TropFlux/input_uncor/ by ${TROPFLUX_ID}
+;   * replace /Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/ by
+;     ${TROPFLUX_OD}
+;   * use :func:`ncdf_getatt <saxo:ncdf_getatt>` for attributes of lat and long
+;     variables
+;   * same problem of time axis and memory like in interp_erai_t2m_1989_2009
+;     using read_ncdf.
+;     replace read_ncdf by netcdf_lec
+;
 ; - fplod 20101215T114503Z aedon.locean-ipsl.upmc.fr (Darwin)
 ;
@@ -70 +123 @@
 ;-
 pro t2m_correction_ncdf
+;
 @common
-da1=19890101 & da2=20091231
-
-file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc'
-initncdf, file
-t2m=read_ncdf('t2m',da1,da2,file=file,/nostr) & t2m=t2m-273.15
+; check for input directory
+;
+; test if ${TROPFLUX_ID} defined
+tropflux_id_env=GETENV('TROPFLUX_ID')
+CASE tropflux_id_env OF
+    ''  :  BEGIN
+     msg = 'eee : ${TROPFLUX_ID} is not defined'
+     ras = report(msg)
+     STOP
+           END
+ ELSE: BEGIN
+     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env
+     ras = report(msg)
+       END
+ENDCASE
+;
+iodirin = isadirectory(tropflux_id_env)
+;
+; existence and protection of ${TROPFLUX_ID}
+IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN
+   msg = 'eee : the directory' + iodirin  + ' is not accessible.'
+   ras = report(msg)
+   STOP
+ENDIF
+;
+; build uncorrected t2 filename
+da1=19890101
+da2=20091231
+filename_t2_uncor='erai_t2m_'+ string(da1,format='(i8.8)')+'_'+ string(da2,format='(i8.8)')+'_oafluxgrid.nc'
+;
+; check if this file exists
+fullfilename_t2_uncor = isafile(iodirin + filename_t2_uncor, NEW=0, /MUST_EXIST)
+IF fullfilename_t2_uncor[0] EQ '' THEN BEGIN
+   msg = 'eee : the file ' + fullfilename_t2_uncor + ' was not found.'
+   ras = report(msg)
+   STOP
+ENDIF
+; test if ${TROPFLUX_OD} defined
+tropflux_od_env=GETENV('TROPFLUX_OD')
+CASE tropflux_od_env OF
+  '' : BEGIN
+         msg = 'eee : ${TROPFLUX_OD} is not defined'
+         ras = report(msg)
+       STOP
+       END
+  ELSE: BEGIN
+          msg = 'iii : ${TROPFLUX_OD} is ' + tropflux_od_env
+          ras = report(msg)
+        END
+ ENDCASE
+;
+; check if output data will be possible
+iodirout = isadirectory(tropflux_od_env)
+;
+; existence and protection
+IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN
+    msg = 'eee : the directory' + iodirout  + ' was not found.'
+    ras = report(msg)
+    STOP
+ENDIF
+;
+; build output filename
+filename_out='TropFlux_t2m_'+ string(da1,format='(i8.8)')+'_'+ string(da2,format='(i8.8)')+'_v50.nc'
+fullfilename_out = iodirout + filename_out
+; in order to avoid unexpected overwritten
+IF (FILE_TEST(fullfilename_out) EQ 1) THEN BEGIN
+   msg = 'eee : the file ' + fullfilename_out  + ' already exists.'
+   ras = report(msg)
+   STOP
+ENDIF
+;
+; define grid parameters with t2 file
+initncdf, fullfilename_t2_uncor
+;
+; get time in t2 file
+timein=ncdf_lec(fullfilename_t2_uncor,var='time')
+jptin=n_elements(timein)
+print, 'time steps from data ', jptin
+print, 'The first 10 time values (variable timein) = ', timein[0:9]
+;
+; find first and last dates yyyymmdd
+; they will be written in global attributes of output file
+da=jul2date(julday(01, 01, 1957,timein[0]))
+cda0=string(da,format='(i8.8)')
+da=jul2date(julday(01, 01, 1957,timein[jptin-1]))
+cda1=string(da,format='(i8.8)')
+print, 'first date ', cda0
+print, 'last date ' , cda1
+;
+; read t2 data
+t2m=ncdf_lec(fullfilename_t2_uncor,var='t2m')
+;
+t2m=t2m-273.15
 help, t2m
-
+;
+jpt=jptin
 t2m_mean=grossemoyenne(t2m,'t',/nan)
-help, w_mean
-
-tt=time & jpt=n_elements(time)
-caldat, time,mon,day,yea
+help, t2m_mean
+
 t2m_m=t2m*0.
 bias_cor=t2m*0.
@@ -89 +230 @@
 ;; line fit ->  -0.0741607x + 1.67601   ;; (2000-2009)
 
-for jt=0,jpt-1 do begin
-  jtt=(time(jt)-julday(1,1,yea(jt))) < 364
+for jt=0,jptin-1 do begin
+  caldat, julday(01, 01, 1957,timein[jt]),mon,day,yea
+  ;++print,day
+  jtt=(julday(01, 01, 1957,timein[jt])-julday(1,1,yea)) < 364
   t=reform(t2m_mean(*,*))
   t2m_m(*,*,jt)=t
@@ -116 +259 @@
 lat=reform(gphit(0,0:jpj-1))
 lon=reform(glamt(0:jpi-1,0))
-time=timegen(7670, units='days', start=julday(1,1,1989)) & jpt=n_elements(time)
-
-cda0=string(jul2date(time(0)),format='(i8.8)')
-cda1=string(jul2date(time(jpt-1)),format='(i8.8)')
-
-time=time-julday(1,1,1950) & jpt=n_elements(time)
-
-;ncfile='!/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_t2m_19890101_20091231_v20.nc'
-ncfile='!/Users/pkb/data/TropFlux/TropFlux_t2m_19890101_20091231_v50.nc'
-lon_attr={units:'degrees_east',long_name:'Longitude'}
-lat_attr={units:'degrees_north',long_name:'Latitude'}
+cda0=string(da1)
+cda1=string(da2)
+
+time=time-julday(1,1,1950)
+jpt=n_elements(time)
+
+ncfile='!' + fullfilename_out
+
+ncdf_getatt, fullfilename_t2_uncor, 'longitude', units=units
+ncdf_getatt, fullfilename_t2_uncor, 'longitude', long_name=long_name
+lon_attr={units:units, long_name:long_name}
+ncdf_getatt, fullfilename_t2_uncor, 'latitude', units=units
+ncdf_getatt, fullfilename_t2_uncor, 'latitude', long_name=long_name
+lat_attr={units:units, long_name:long_name}
+
 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:'1950-JAN-01 00:00:00'}
-t2m_attr={units:'degK',missing_value:1.e20,long_name:'Air Temperature at 2m',short_name:'t2m',axis:'TYX'}
+
+ncdf_getatt, fullfilename_t2_uncor, 't2m', units=units
+ncdf_getatt, fullfilename_t2_uncor, 't2m', long_name=long_name
+t2m_attr={units:units,missing_value:1.e20,long_name:long_name,short_name:'t2m',axis:'TYX'}
 globattr={source:'Basic data obtained from ERAI.  Mean correction for air temperautre bias and correction for variability are applied',timerange:cda0+' - '+cda1}