Changeset 17 for trunk/src

Timestamp:

12/17/10 11:00:48 (13 years ago)

Author:

pinsard

Message:

oaflux_mask_30N30S.pro no more hard coded directories

File:

• trunk/src/oaflux_mask_30N30S.pro (modified) (6 diffs)

trunk/src/oaflux_mask_30N30S.pro

@@ -7 +7 @@
 ; ======================
 ;
+; Produce a NetCDF file with  OAFLUX mask over global tropical oceans 
+;
 ;     .. graphviz::
 ;
@@ -13 +15 @@
 ;           rankdir="LR",
 ;           ]
-;           file_oaflux [shape=ellipse,fontname=Courier,label="/Volumes/vialardj/OAFLUX/lh_oaflux_2004.nc"];
-;           mask [shape=ellipse,fontname=Courier,label="/Users/jv/data/OAFLUX/mask_oaflux_30N30S.nc"];
+;           file_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lh_oaflux_2004.nc"];
+;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/mask_oaflux_30N30S.nc"];
 ;
 ;           oaflux_mask_30N30S [shape=box,
@@ -31 +33 @@
 ; :ref:`guide data OAFLUX <data_in_oaflux>`
 ;
+; :ref:`tropflux_profile.sh`
+;
+; :func:`report`
+; :func:`initncdf`
+; :func:`ncdf_lec`
+; :ref:`ncdf_quickwrite`
+; :func:`isadirectory`
+;
 ; EXAMPLES
 ; ========
@@ -41 +51 @@
 ; ====
 ;
-; hard coded directory - usage of ${TROPFLUX_ID}
+; SAXO new (keep compatibility false)
 ;
 ; coding rules
 ;
+; pro -> function
+;
+; NetCDF CF (may be area_type) 
+; 
+; global attributes corrections :
+;  - written now : OAFLUX mask over the Indian Ocean: 40E-120E, 30S-30N
+;  - while data latitude=[-29.5,+29.5] longitude=[30.5,379.5]
+;
+; KNOWN ISSUES
+; ============
+;
+; test of existence of fullfilename_in not very efficient because
+; MUST_EXIST keyword of :func:`isafile` not yet implemented
+;
 ; EVOLUTIONS
 ; ==========
+;
+; - fplod 20101217T081915Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * remove hard coded directory - usage of ${TROPFLUX_ID} and ${TROPFLUX_OD}
+;   * add IO test
 ;
 ; - fplod 20101216T141137Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -59 +88 @@
 ;
 PRO oaflux_mask_30N30S
-
+;
 @common
-
 ;
-;; Read oaflux grid and mask
+; check for input file
 ;
-dir='/Volumes/vialardj/OAFLUX/'
-dir='/Users/jv/data/OAFLUX/'
-
-yy=2004 & cy=string(yy,format='(i4.4)')
-fi=dir+'lh_oaflux_'+cy+'.nc'
-initncdf, fi
-lh=ncdf_lec(fi,var='lhtfl',count=[jpi,jpj,1]) & lh=float(lh)
+; 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 input filename
+yyyy=2004
+cy=string(yyyy,format='(i4.4)')
+filename_in='lh_oaflux_'+cy+'.nc'
+;
+; check if this file exists
+fullfilename_in = isafile(iodirin + filename_in, NEW=0, /MUST_EXIST)
+IF fullfilename_in[0] EQ '' THEN BEGIN
+   msg = 'eee : the file ' + fullfilename_in + ' 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
+;
+fullfilename_out=iodirout+'mask_oaflux_30N30S.nc'
+; 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
+;
+; Read oaflux grid and mask
+initncdf, fullfilename_in
+lh=ncdf_lec(fullfilename_in,var='lhtfl',count=[jpi,jpj,1]) & lh=float(lh)
 valmask=1.e20
 ind=where(lh ge 32000,compl=nind) & lh(ind)=valmask & lh(nind)=-0.1*lh(nind)
@@ -88 +180 @@
 
 
-ncfile='!'+dir+'mask_oaflux_30N30S.nc'
+ncfile='!'+fullfilename_out
 lon_attr={units:'degrees_east',long_name:'Longitude'}
 lat_attr={units:'degrees_north',long_name:'Latitude'}