Changeset 19 for trunk/src/interp_olr_30n30s_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_olr_30n30s_1989_2009.pro (modified) (5 diffs)

trunk/src/interp_olr_30n30s_1989_2009.pro

@@ -15 +15 @@
 ;           file_in [shape=ellipse,fontname=Courier,label="/Users/pkb/data/heat_budget/NOAA_OLR/new/olr.day.mean_new.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="/Users/pkb/data/heat_budget/NOAA_OLR/newolr_oafluxgrid_30n30s_19890101_20091231.nc"];
@@ -32 +32 @@
 ; ========
 ;
+; :ref:`tropflux_profile.sh`
+;
 ; EXAMPLES
 ; ========
@@ -46 +48 @@
 ; 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 20101215T113118Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -64 +76 @@
 pro interp_olr_30n30s_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
 
 filein='/Users/pkb/data/heat_budget/NOAA_OLR/new/olr.day.mean_new.nc'
-gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc'
 
 initncdf, filein
@@ -79 +126 @@
 timein=time & jptin=jpt
 
-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, olrin,lonin,latin,mskin,lonout,latout,mskout