Changeset 170

Timestamp:

03/20/12 12:03:36 (12 years ago)

Author:

pinsard

Message:

interp_erai_sst is now a function

Location:

trunk

Files:

• docs/docs_dev/source/guides/data_content.rst (modified) (1 diff)
• docs/docs_dev/source/guides/generate_tropflux.rst (modified) (1 diff)
• src/interp_erai_sst.pro (moved) (moved from trunk/src/interp_erai_sst_1989_2009.pro) (18 diffs)
• src/oaflux_mask_30n30s.pro (modified) (1 diff)
• src/sst_correction_ncdf.pro (modified) (2 diffs)
• src/tropflux.sh (modified) (3 diffs)

trunk/docs/docs_dev/source/guides/data_content.rst

@@ -137 +137 @@
  - :func:`interp_erai_dewt`
  - :func:`interp_erai_lwr`
- - :ref:`interp_erai_sst_1989_2009.pro`
+ - :func:`interp_erai_sst`
  - :func:`interp_erai_t2m`
  - :ref:`interp_erai_ws_1989_2009.pro`

trunk/docs/docs_dev/source/guides/generate_tropflux.rst

@@ -102 +102 @@
  - :func:`interp_erai_dewt`
  - :func:`interp_erai_lwr`
- - :ref:`interp_erai_sst_1989_2009.pro`
+ - :func:`interp_erai_sst`
  - :func:`interp_erai_msl.pro`
  - :func:`interp_erai_t2m`

trunk/src/interp_erai_sst.pro

@@ -1 +1 @@
 ;+
 ;
-; .. _interp_erai_sst_1989_2009.pro:
-;
-; =============================
-; interp_erai_sst_1989_2009.pro
-; =============================
+;
+; ===================
+; interp_erai_sst.pro
+; ===================
+;
+; .. function:: interp_erai_sst(yyyymmddb,yyyymmdde)
 ;
 ; DESCRIPTION
@@ -25 +26 @@
 ; Interpolated sst
 ; is written in
-; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc`
+; :file:`${PROJECT_OD}/erai_sst_{yyyymmdd}_{yyyymmdd}_oafluxgrid.nc`
 ; if this file not already exists.
 ;
 ; This output file
-; :file:`${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc`
+; :file:`${PROJECT_OD}/erai_sst_{yyyymmdd}_{yyyymmdd}_oafluxgrid.nc`
 ; must be processed after by
 ; :ref:`sst_correction_ncdf.pro`.
@@ -41 +42 @@
 ;     .. graphviz::
 ;
-;        digraph interp_erai_sst_1989_2009 {
+;        digraph interp_erai_sst {
 ;
 ;           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"];
-;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc"];
-;
-;           interp_erai_sst_1989_2009 [shape=box,
+;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_sst_{yyyymmdd}_{yyyymmdd}_oafluxgrid.nc"];
+;
+;           interp_erai_sst [shape=box,
 ;           fontname=Courier,
 ;           color=blue,
-;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/interp_erai_sst_1989_2009.pro",
-;           label="${PROJECT}/src/interp_erai_sst_1989_2009.pro"];
-;
-;           {file_in mask} -> {interp_erai_sst_1989_2009} -> {file_out}
+;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/interp_erai_sst.pro",
+;           label="${PROJECT}/src/interp_erai_sst.pro"];
+;
+;           {file_in mask} -> {interp_erai_sst} -> {file_out}
 ;
 ;        }
@@ -61 +62 @@
 ;
 ; :ref:`interpolate_data`
+;
+; Used by :ref:`tropflux.sh`
 ;
 ; :ref:`project_profile.sh`
@@ -81 +84 @@
 ;
 ;  IDL> .compile file_interp
-;  IDL> interp_erai_sst_1989_2009
+;  IDL> yyyymmddb = 19890101
+;  IDL> yyyymmdde = 20091231
+;  IDL> result = interp_erai_sst(yyyymmddb, yyymmdde)
+;  IDL> print, result
+;
 ;
 ; TODO
@@ -94 +101 @@
 ; use real output of :ref:`compute_erai_daily_region_2d.sh`.
 ;
+; scientific validation (strange look of data with ncview now)
+;
+; check [yyyymmddb,yyyymmdde] validity
+;
+; check if [yyyymmddb,yyyymmdde] is included in ERA-I file
+;
+; what happen if yyyymmdde > 20091231 : need to read an other ERA-I file
+;
+; add an example with only one month (ie 200001)
+;
 ; KNOWN ISSUES
 ; ============
@@ -105 +122 @@
 ; $Id$
 ;
-; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_sst_1989_2009.pro $
+; $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/tropflux/svn/trunk/src/interp_erai_sst.pro $
+;
+; - fplod 20120320
+;
+;   * taking project_overwite into account
+;   * try to add compile_opt seems to be incompatible with ncdf_quickwrite
+;   * pro -> function
 ;
 ; - pinsard 2011-08-23T09:09:16Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -136 +159 @@
 ;
 ;-
-pro interp_erai_sst_1989_2009
+function interp_erai_sst $
+         , yyyymmddb $
+         , yyyymmdde
+;
+;++compile_opt idl2, strictarrsubs, logical_predicate
 ;
 @cm_4cal
@@ -144 +171 @@
 @cm_project
 ;
+; Return to caller if errors
+ON_ERROR, 2
+;
+result = -1
+;
+usage = 'result = interp_erai_lwr(yyyymmddb, yyyymmdde)
+nparam = N_PARAMS()
+IF (nparam NE 2) THEN BEGIN
+   ras = report(['Incorrect number of arguments.' $
+   + '!C' $
+   + 'Usage : ' + usage])
+   return, result
+ENDIF
+;
 ; check for input directory
 ;
@@ -151 +192 @@
      msg = 'eee : ${PROJECT_ID} is not defined'
      ras = report(msg)
-     STOP
+     return, result
            END
  ELSE: BEGIN
@@ -165 +206 @@
    msg = 'eee : the directory' + iodirin  + ' is not accessible.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -178 +219 @@
    msg = 'eee : the file ' + fullfilename_msk + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -191 +232 @@
    msg = 'eee : the file ' + fullfilename + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -199 +240 @@
          msg = 'eee : ${PROJECT_OD} is not defined'
          ras = report(msg)
-       STOP
+       return, result
        END
   ELSE: BEGIN
@@ -214 +255 @@
     msg = 'eee : the directory' + iodirout  + ' was not found.'
     ras = report(msg)
-    STOP
+    return, result
 ENDIF
 ;
 ; build output filename
-filename_out = 'erai_sst_19890101_20091231_oafluxgrid.nc'
+filename_out = 'erai_sst_' + string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '_oafluxgrid.nc'
 fullfilename_out = iodirout + filename_out
 ; in order to avoid unexpected overwritten
-IF (FILE_TEST(fullfilename_out) EQ 1) THEN BEGIN
+IF ((FILE_TEST(fullfilename_out) EQ 1)  AND (project_overwrite EQ 0)) THEN BEGIN
    msg = 'eee : the file ' + fullfilename_out  + ' already exists.'
    ras = report(msg)
-   STOP
+   return, result
 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)
+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=fullfilename,/nostr)
 ;
 timein=time
 jptin=jpt
-tab=sstkin(*,*,0)
+tab=sstkin[*,*,0]
 mskin=glamt*0.+1.
 ind=where((tab eq 1.e9) or finite(tab,/nan))
-mskin(ind)=0.
+mskin[ind]=0.
 ;
 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)
+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=fullfilename_msk,/nostr)
 ;
@@ -255 +296 @@
 for jt=0,jptin-1 do begin
   ; ++print, 'Interpolation jt=',jt,' / ',jptin-1
-  tab=reform(sstkin(*,*,jt))
-  sstkout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+  tab=reform(sstkin[*,*,jt])
+  sstkout[*,*,jt]=call_interp2d(tab,lonin,latin,mskin $
       , lonout,latout,method='bilinear' $
       , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
-  sstkout(*,*,jt)=sstkout(*,*,jt)*mskout+(1.-mskout)*1.e20
+  sstkout[*,*,jt]=sstkout[*,*,jt]*mskout+(1.-mskout)*1.e20
 endfor
 ;
@@ -265 +306 @@
 time=timegen(7670, start=julday(1,1,1989,0), units='days')
 jpt=n_elements(time)
-cda0=string(jul2date(time(0)),format='(i8.8)')
-cda1=string(jul2date(time(jpt-1)),format='(i8.8)')
+cda0=string(jul2date(time[0]),format='(i8.8)')
+cda1=string(jul2date(time[jpt-1]),format='(i8.8)')
 timein=timein-julday(1,1,1950,00,00)
 ;
@@ -286 +327 @@
 @ncdf_quickwrite
 ;
+result = 0
+return, result
+;
+
 end

trunk/src/oaflux_mask_30n30s.pro

@@ -65 +65 @@
 ; :func:`interp_erai_dewt`,
 ; :func:`interp_erai_lwr`,
-; :ref:`interp_erai_sst_1989_2009.pro`,
+; :func:`interp_erai_sst`,
 ; :func:`interp_erai_t2m`,
 ; :ref:`interp_erai_ws_1989_2009.pro`,

trunk/src/sst_correction_ncdf.pro

@@ -13 +13 @@
 ; ++
 ; has been produced by
-; :ref:`interp_erai_sst_1989_2009.pro`.
+; :func:`interp_erai_sst`.
 ;
 ; Corrected sst on OAFLUX grid
@@ -61 +61 @@
 ; :ref:`mooring_corrections`
 ;
-; :ref:`interp_erai_sst_1989_2009.pro`
+; :func:`interp_erai_sst`
 ;
 ; :func:`initncdf <saxo:initncdf>`

trunk/src/tropflux.sh

@@ -91 +91 @@
 # :func:`oaflux_mask_30n30s`, :func:`interp_erai_t2m`, 
 # :func:`interp_erai_msl`, :func:`interp_erai_dewt`, :func:`interp_erai_lwr`,
-# etc.
+# :func:`interp_erai_sst`, etc.
 #
 # Next step : +validation ++ comparison
@@ -112 +112 @@
 #
 # $URL$
+#
+# - pinsard 20120320
+#
+#   * handling error of inter_erai_sst (now function)
 #
 # - pinsard 20120319
@@ -317 +321 @@
    exit
 ENDIF
-;interp_erai_sst_1989_2009
+result = interp_erai_sst(${yyyymmddb}, ${yyyymmdde})
+IF result < 0 THEN BEGIN
+   msg = 'eee : pb after interp_erai_t2m' + string(result)
+   err = report(msg)
+   exit
+ENDIF
 result = interp_erai_t2m(${yyyymmddb}, ${yyyymmdde})
 IF result < 0 THEN BEGIN