Changeset 66 for trunk/src/interp_erai_sst_1989_2009.pro

Timestamp:

07/05/11 13:58:48 (13 years ago)

Author:

pinsard

Message:

no more hard coded path in iterp_erai*.pro

File:

• trunk/src/interp_erai_sst_1989_2009.pro (modified) (10 diffs)

trunk/src/interp_erai_sst_1989_2009.pro

@@ -3 +3 @@
 ; .. _interp_erai_sst_1989_2009.pro:
 ;
-; =============================
-; interp_erai_sst_1989_2009.pro
-; =============================
+; ==============================================================================
+; interp_erai_sst_1989_2009.pro - Interpolation of sst from ERA-I grid to OAFLUX
+; ==============================================================================
+;
+; Interpolation of sst from ERA-I grid to OAFLUX grid
+;
+; :file:`${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc` containing sst from ERA-I have been produced
+; by :ref:`compute_erai_daily_region_2d.sh`.
+;
+; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been p roduced by :ref:`oaflux_mask_30N30S.pro`.
+;
+; Interpolated sst is written in
+; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` if this file not already exists.
+;
+; This output file :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc` must be processed after by :ref:`sst_correction_ncdf.pro`.
+;
 ;
 ;     .. graphviz::
@@ -13 +26 @@
 ;           rankdir="LR",
 ;           ]
-;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_sstk_TROP_1989_2009.nc"];
+;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc"];
 ;
 ;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"];
 ;
-;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc"];
+;           ncfile [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc"];
 ;
 ;           interp_erai_sst_1989_2009 [shape=box,
@@ -42 +55 @@
 ; :func:`ncdf_quickwrite <saxo:ncdf_quickwrite>`
 ;
+; :ref:`sst_correction_ncdf`
+;
 ; EXAMPLES
 ; ========
@@ -47 +62 @@
 ; ::
 ;
+;  IDL> .compile file_interp
 ;  IDL> interp_erai_sst_1989_2009
 ;
@@ -52 +68 @@
 ; ====
 ;
-; hard coded directory - usage of ${PROJECT_ID}
+; make it work : pb on loholt1::
+;
+;   Variable is undefined: OUTMASK_IND.
 ;
 ; coding rules
+;
+; use real output of :ref:`compute_erai_daily_region_2d.sh`.
 ;
 ; KNOWN ISSUES
@@ -65 +85 @@
 ; ==========
 ;
+; $Id$
+;
+; $URL$
+;
+; - pinsard 2011-07-04T15:45:35Z loholt1.ipsl.polytechnique.fr (Linux)
+;
+;   * usage of ${PROJECT_ID} and $PROJECT_OD
+;
 ; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin)
 ;
@@ -78 +106 @@
 ;
 ; - pbk 2008
-;
 ;
 ;   * creation
@@ -84 +111 @@
 ;-
 pro interp_erai_sst_1989_2009
-@common
+;
+@cm_4cal
+@cm_4data
+@cm_4mesh
+@cm_4data
 @cm_project
 ;
@@ -121 +152 @@
    STOP
 ENDIF
-
-filein='/Volumes/PRAVEEN/ERAI_global/20c3m_erai_sstk_TROP_1989_2009.nc'
-
-initncdf, filein
+;
+; build data filename
+filename='20c3m_erai_sstk_TROP_1989_2009.nc'
+;
+; check if this file exists
+fullfilename = isafile(iodirin + filename, 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
+;
+; 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
+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 = 'erai_sst_19890101_20091231_oafluxgrid.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
 domdef
 latin=reform(gphit(0,*)) & lonin=reform(glamt(*,0))
 print, 'lat grid ',min(latin),max(latin),latin(1)-latin(0)
 print, 'lon grid ',min(lonin),max(lonin),lonin(1)-lonin(0)
-sstkin=read_ncdf("sstk",19871201,20100930,file=filein,/nostr)
+sstkin=read_ncdf("sstk",19871201,20100930,file=fullfilename,/nostr)
 
 timein=time & jptin=jpt
@@ -166 +239 @@
 lat=latout
 lon=lonout
-ncfile='/Volumes/PRAVEEN/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc'
+ncfile='!' + fullfilename_out
 lon_attr={units:'degrees_east',long_name:'Longitude'}
 lat_attr={units:'degrees_north',long_name:'Latitude'}