Changeset 19 for trunk/src/interp_erai_ws_1989_2009.pro

Timestamp:

12/17/10 16:18:08 (14 years ago)

Author:

pinsard

Message:

remove hard coded directory for mask_oaflux_30N30S.nc

File:

• trunk/src/interp_erai_ws_1989_2009.pro (modified) (7 diffs)

trunk/src/interp_erai_ws_1989_2009.pro

@@ -16 +16 @@
 ;           file_v10 [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_v10_TROP_1989_2009.nc"];
 ;
-;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"];
+;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"];
 ;
 ;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc"];
@@ -33 +33 @@
 ; ========
 ;
+; :ref:`tropflux_profile.sh`
+;
 ; EXAMPLES
 ; ========
@@ -48 +50 @@
 ; coding rules
 ;
+; 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 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * remove hard coded directory for mask_oaflux_30N30S.nc
 ;
 ; - fplod 20101215T112657Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -66 +78 @@
 pro interp_erai_ws_1989_2009
 @common
+;
+; 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 mask filename
+filename_msk='mask_oaflux_30N30S.nc'
+;
+; check if this file exists
+fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST)
+IF fullfilename_msk[0] EQ '' THEN BEGIN
+   msg = 'eee : the file ' + fullfilename_msk + ' was not found.'
+   ras = report(msg)
+   STOP
+ENDIF
 
 st=19880101 & en=20100930
@@ -77 +125 @@
 v10in=read_ncdf("v10",st, en,file=filein,/nostr)
 
-gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc'
-
 initncdf, filein
 domdef
@@ -89 +135 @@
 mskin=glamt*0.+1.
 
-initncdf, gridout
+initncdf, fullfilename_msk
 domdef
 latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0))
 print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0)
 print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0)
-mskout=read_ncdf("msk", file=gridout,/nostr)
+mskout=read_ncdf("msk", file=fullfilename_msk,/nostr)
 
 help, u10in, v10in,lonin,latin,mskin,lonout,latout,mskout
@@ -121 +167 @@
 xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*))
 
-initncdf, gridout
+initncdf, fullfilename_msk
 valmask=1.e20