Changeset 70 for trunk/src

Timestamp:

08/08/11 12:41:36 (13 years ago)

Author:

pinsard

Message:

fix call_interp2d call

Location:

trunk/src

Files:

• interp_erai_dewt_1989_2009.pro (modified) (4 diffs, 1 prop)
• interp_erai_lwr_1989_2009.pro (modified) (6 diffs, 1 prop)
• interp_erai_sst_1989_2009.pro (modified) (3 diffs, 1 prop)
• interp_erai_t2m_1989_2009.pro (modified) (9 diffs, 1 prop)
• interp_erai_ws_1989_2009.pro (modified) (3 diffs, 1 prop)

trunk/src/interp_erai_dewt_1989_2009.pro

@@ -4 +4 @@
 ;
 ; ==============================
-; interp_erai_dewt_1989_2009.pro 
+; interp_erai_dewt_1989_2009.pro
 ; ==============================
 ;
@@ -73 +73 @@
 ;   Variable is undefined: OUTMASK_IND.
 ;
+; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call
+;
+; use real output of :ref:`compute_erai_daily_region_2d.sh`.
+;
 ; coding rules
 ;
@@ -87 +91 @@
 ;
 ; $URL$
+;
+; - fplod 20110808T094156Z cratos (Linux)
+;
+;   * add OUTMASK_IND and SET_OUTMSKVAL to call_interp2d call
 ;
 ; - pinsard 2011-07-05T07:33:36Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -220 +228 @@
   print, 'Interpolation jt=',jt,' / ',jptin-1
   tab=reform(d2min(*,*,jt))
-  d2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  d2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+      , lonout,latout,method='bilinear' $
+      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
   d2mout(*,*,jt)=d2mout(*,*,jt)*mskout+(1.-mskout)*1.e20
 endfor

trunk/src/interp_erai_lwr_1989_2009.pro

@@ -75 +75 @@
 ;   Variable is undefined: OUTMASK_IND.
 ;
+; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call
+;
+; use real output of :ref:`compute_erai_daily_region_2d.sh`.
+;
+; 19890101 is not provided in ERA-Intermim file str : a possible workaround is to copy 19890102
+; dataset to 19890101 using nco before processing. see :ref:`compute_erai_daily_region_2d.sh`
+;
 ; coding rules
 ;
@@ -89 +96 @@
 ;
 ; $URL$
+;
+; - fplod 20110808T094156Z cratos (Linux)
+;
+;   * add OUTMASK_IND and SET_OUTMSKVAL to call_interp2d call
+;   * fix typo
+;     .. note::
+;
+;        IDL do not warn when asking size of an undefined array::
+;
+;           IDL> must_crash=size(ginette)
 ;
 ; - pinsard 2011-07-05T07:33:36Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -203 +220 @@
 lwrin=read_ncdf("str",0,7668,/timestep,file=fullfilename,/nostr)
 lwrin=-1*lwrin
-help, lwrin
 
 initncdf, fullfilename
@@ -222 +238 @@
 mskout=read_ncdf("msk", file=fullfilename_msk,/nostr)
 
-help, lwrin
-
-si=size(swrin)
+
+si=size(lwrin)
 lwrout=fltarr(jpi,jpj,jptin)
 
@@ -231 +246 @@
 
   tab=reform(lwrin(*,*,jt))
-  lwrout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  lwrout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+     , lonout,latout,method='bilinear' $
+     , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
+  help, OUTMASK_IND
+  help, SET_OUTMSKVAL
   lwrout(*,*,jt)=lwrout(*,*,jt)*mskout+(1.-mskout)*1.e20
 
@@ -257 +276 @@
 globattr={source:'Tropical ocean fluxes obtained from ERA Interim',timerange:cda0+' - '+cda1}
 
-help,lwrout
-
 ncfields = 'lwr[longitude,latitude,time]=lwrout:lwr_attr; ' $
                       + 'longitude[]=xlon:lon_attr; ' $

trunk/src/interp_erai_sst_1989_2009.pro

@@ -75 +75 @@
 ;   Variable is undefined: OUTMASK_IND.
 ;
+; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call
+;
 ; coding rules
 ;
@@ -91 +93 @@
 ;
 ; $URL$
+;
+; - fplod 20110808T094156Z cratos (Linux)
+;
+;   * add OUTMASK_IND and SET_OUTMSKVAL to call_interp2d call
 ;
 ; - pinsard 2011-07-04T15:45:35Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -226 +232 @@
   print, 'Interpolation jt=',jt,' / ',jptin-1
   tab=reform(sstkin(*,*,jt))
-  sstkout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  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
 endfor

trunk/src/interp_erai_t2m_1989_2009.pro

@@ -76 +76 @@
 ;  IDL> interp_erai_t2m_1989_2009
 ;
-;
 ; TODO
 ; ====
@@ -97 +96 @@
 ; why two "initncdf, fullfilename_msk"
 ;
-; hard coded attributes for t2m (missing value, short name, axis) and time (origin) : use ncdf_getatt 
+; hard coded attributes for t2m (missing value, short name, axis) and time (origin) : use ncdf_getatt
 ;
 ; CF conventions
+;
+; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call
+;
+; remove ncdf_quickwrite exercices
 ;
 ; KNOWN ISSUES
@@ -109 +112 @@
 ; EVOLUTIONS
 ; ==========
+;
+; $Id$
+;
+; $URL$
+;
+; - fplod 20110808T094156Z cratos (Linux)
+;
+;   * add OUTMASK_IND and SET_OUTMSKVAL to call_interp2d call
 ;
 ; - fplod 20110103T153734Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -251 +262 @@
 initncdf, fullfilename_t2
 domdef
-latin=reform(gphit(0,*)) 
+latin=reform(gphit(0,*))
 lonin=reform(glamt(*,0))
 print, 'lat grid from data',min(latin),max(latin),latin(1)-latin(0)
@@ -285 +296 @@
 initncdf, fullfilename_msk
 domdef
-latout=reform(gphit(0,*)) 
+latout=reform(gphit(0,*))
 lonout=reform(glamt(*,0))
 print, 'lat grid from mask ',min(latout),max(latout),latout(1)-latout(0)
@@ -297 +308 @@
 for jt=0,jptin-1 do begin
   tab=reform(t2min(*,*,jt))
-  t2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  t2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+      , lonout,latout,method='bilinear' $
+      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
   t2mout(*,*,jt)=t2mout(*,*,jt)*mskout+(1.-mskout)*1.e20
 endfor
@@ -324 +337 @@
 help, t2mout
 help, timein
-;ncfields = 't2m[longitude,latitude,*time]=t2mout:t2m_attr; ' $
-;                      + 'longitude[]=lon:lon_attr; ' $
-;                      + 'latitude[]=lat:lat_attr; ' $
-;                     + 'time[]=timein:time_attr ' $
-;                      + ' @globattr'
+ncfields = 't2m[longitude,latitude,*time]=t2mout:t2m_attr; ' $
+                      + 'longitude[]=lon:lon_attr; ' $
+                      + 'latitude[]=lat:lat_attr; ' $
+                      + 'time[]=timein:time_attr ' $
+                      + ' @globattr'
+@ncdf_quickwrite
 
 ; ok ncfields = 'longitude[]=lon:lon_attr;latitude[]=lat:lat_attr'
@@ -354 +368 @@
 help, lat
 help, lon
-help, t2m 
+help, t2m
 help, time
 help, timein
@@ -368 +382 @@
 ;ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'
 
-timepb=findgen(7670)
-ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'
-@ncdf_quickwrite
+;++timepb=findgen(7670)
+;++ncfields = 't2m[longitude,latitude,*time];longitude[]=lon:lon_attr;latitude[]=lat:lat_attr;time[]=timepb:time_attr'
+;++@ncdf_quickwrite
 
 

trunk/src/interp_erai_ws_1989_2009.pro

@@ -78 +78 @@
 ;   Variable is undefined: OUTMASK_IND.
 ;
+; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call
+;
+; use real output of :ref:`compute_erai_daily_region_2d.sh`.
+;
 ; coding rules
 ;
@@ -92 +96 @@
 ;
 ; $URL$
+;
+; - fplod 20110808T094156Z cratos (Linux)
+;
+;   * add OUTMASK_IND and SET_OUTMSKVAL to call_interp2d call
 ;
 ; - pinsard 2011-07-05T07:33:36Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -245 +253 @@
   print, 'Interpolation jt=',jt,' / ',jptin-1
   tab=reform(u10in(*,*,jt))
-  u10out(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  u10out(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+      , lonout,latout,method='bilinear' $
+      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
   u10out(*,*,jt)=u10out(*,*,jt)*mskout+(1.-mskout)*1.e20
 
   tab=reform(v10in(*,*,jt))
-  v10out(*,*,jt)=call_interp2d(tab,lonin,latin,mskin,lonout,latout,method='bilinear')
+  v10out(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+      , lonout,latout,method='bilinear' $
+      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
   v10out(*,*,jt)=v10out(*,*,jt)*mskout+(1.-mskout)*1.e20