Changeset 72 for trunk/src/lwr_correction_ncdf.pro

Timestamp:

08/08/11 17:08:15 (13 years ago)

Author:

pinsard

Message:

start to work on correction tools

File:

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

trunk/src/lwr_correction_ncdf.pro

@@ -6 +6 @@
 ; lwr_correction_ncdf.pro
 ; =======================
+;
+; Correction of lwr on OAFLUX grid
+;
+; :file:`${PROJECT_OD}/erai_lwr_19890101_20091231_oafluxgrid.nc` have been
+; produced by :ref:`interp_erai_lwr_1989_2009.pro`.
+;
+; Corrected lwr on OAFLUX grid is written in
+; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc`
+; if this file not already exists.
+;
+; This file will be used by :ref:`TropFlux_19890101_20091231.pro`.
 ;
 ;     .. graphviz::
@@ -14 +25 @@
 ;           ]
 ;
-;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_lwr_19890102_20091231_oafluxgrid.nc"];
-;
-;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_cor/full_cor/TropFlux_lwr_19890101_20091231_v1.nc"];
+;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_lwr_19890101_20091231_oafluxgrid.nc"];
+;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc"];
+;
 ;           lwr_correction_ncdf [shape=box,
 ;           fontname=Courier,
@@ -23 +34 @@
 ;           label="${PROJECT}/src/lwr_correction_ncdf.pro"];
 ;
-;           {file_in} -> {lwr_correction_ncdf} -> {ncfile}
+;           {file_in} -> {lwr_correction_ncdf} -> {file_out}
 ;
 ;        }
@@ -29 +40 @@
 ; SEE ALSO
 ; ========
+;
+; :ref:`project_profile.sh`
+;
+; :ref:`mooring_corrections`
+;
+; :ref:`interp_erai_lwr_1989_2009.pro`
 ;
 ; :func:`initncdf <saxo:initncdf>`
@@ -44 +61 @@
 ;  IDL> lwr_correction_ncdf
 ;
-;
 ; TODO
 ; ====
 ;
-; hard coded directory - usage of ${PROJECT_ID}
-;
 ; coding rules
+;
+; lwr long name 'Air Temperature at 2m" : ???
+;
+; 19890101 is not provided in ERA-Intermim file str : a possible workaround is
+; to copy 19890102
+; dataset to 19890101 using nco before processing. see
+; ref:`compute_erai_daily_region_2d.sh`
+;
+; here I (fp) replace hard coded date 19890102 by 19890101. must check
+; consistency with this module and interp_erai_lwr_1989_2009.pro.
+;
+; KNOWN ISSUES
+; ============
+;
+; test of existence of fullfilename_in not very efficient because
+; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented
 ;
 ; EVOLUTIONS
 ; ==========
 ;
+; - fplod 20110808T120234Z cratos (Linux)
+;
+;   * usage of ${PROJECT_OD}
+;   * complete description
+;   * replace erai_lwr_19890102_20091231_oafluxgrid.nc by erai_lwr_19890101_20091231
+;   * remove v1 in output filename
+;
 ; - fplod 20101215T113707Z aedon.locean-ipsl.upmc.fr (Darwin)
 ;
@@ -69 +106 @@
 ;-
 pro lwr_correction_ncdf
-@common
+;
+@cm_4cal
+@cm_4data
+@cm_4mesh
+@cm_4data
+@cm_project
+;
+; test if ${PROJECT_OD} defined
+CASE project_od_env OF
+  '' : BEGIN
+         msg = 'eee : ${PROJECT_OD} is not defined'
+         ras = report(msg)
+       STOP
+       END
+  ELSE: BEGIN
+          msg = 'iii : ${PROJECT_OD} is ' + project_od_env
+          ras = report(msg)
+        END
+ ENDCASE
+;
+; check if output data will be possible
+iodirout = isadirectory(project_od_env)
+;
+; existence and protection for reading
+IF (FILE_TEST(iodirout, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN
+   msg = 'eee : the directory' + iodirout  + ' is not accessible.'
+   ras = report(msg)
+   STOP
+ENDIF
+;
+; existence and protection for writing
+IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN
+    msg = 'eee : the directory' + iodirout  + ' was not found.'
+    ras = report(msg)
+    STOP
+ENDIF
+;
 da1=19880101 & da2=20091231
-
-file='/Volumes/PRAVEEN/TropFlux/input_uncor/erai_lwr_19890102_20091231_oafluxgrid.nc'
-initncdf, file
-lwr=read_ncdf('lwr',da1,da2,file=file,/nostr)
+;
+; build data filename
+filename='erai_lwr_19890101_20091231_oafluxgrid.nc'
+;
+; check if this file exists
+fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST)
+IF fullfilename[0] EQ '' THEN BEGIN
+   msg = 'eee : the file ' + fullfilename + ' was not found.'
+   ras = report(msg)
+   STOP
+ENDIF
+;
+; build output filename
+filename_out = 'TropFlux_lwr_19890101_20091231.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
+;
+initncdf, fullfilename
+lwr=read_ncdf('lwr',da1,da2,file=fullfilename,/nostr)
 help, lwr
-
+;
 lwr_mean=grossemoyenne(lwr,'t',/nan)
 help,lwr_mean
@@ -116 +209 @@
 time=time-julday(1,1,1950) & jpt=n_elements(time)
 
-ncfile='!/Volumes/PRAVEEN/TropFlux/input_cor/full_cor/TropFlux_lwr_19890101_20091231_v1.nc'
+ncfile='!' + fullfilename_out
 lon_attr={units:'degrees_east',long_name:'Longitude'}
 lat_attr={units:'degrees_north',long_name:'Latitude'}
@@ -123 +216 @@
 globattr={source:'Basic data obtained from ERAI.  Mean correction for  bias and correction for variability are applied',timerange:cda0+' - '+cda1}
 
-
 ncfields = 'lwr[longitude,latitude,time]=lwr_new:lwr_attr; ' $
                       + 'longitude[]=lon:lon_attr; ' $