Changeset 156 for trunk/src

Timestamp:

03/06/12 14:54:23 (12 years ago)

Author:

pinsard

Message:

try to add compile_opt but seems to be incompatible with ncdf_quickwrite

Location:

trunk/src

Files:

• interp_erai_t2m.pro (modified) (5 diffs)
• oaflux_mask_30n30s.pro (modified) (3 diffs)

trunk/src/interp_erai_t2m.pro

@@ -168 +168 @@
 ; - fplod 20120306
 ;
+;   * try to add compile_opt seems to be incompatible with ncdf_quickwrite
 ;   * pro -> function
 ;
@@ -234 +235 @@
          , yyyymmddb $
          , yyyymmdde
+;
+;++compile_opt idl2, strictarrsubs, logical_predicate
 ;
 @cm_4cal
@@ -341 +344 @@
 initncdf, fullfilename_t2
 domdef
-latin=reform(gphit(0,*))
-lonin=reform(glamt(*,0))
-print, 'lat grid from data',min(latin),max(latin),latin(1)-latin(0)
-print, 'lon grid from data',min(lonin),max(lonin),lonin(1)-lonin(0)
+latin=reform(gphit[0,*])
+lonin=reform(glamt[*,0])
+print, 'lat grid from data',min(latin),max(latin),latin[1]-latin[0]
+print, 'lon grid from data',min(lonin),max(lonin),lonin[1]-lonin[0]
 ;
 ; get time in t2 file
@@ -375 +378 @@
 initncdf, fullfilename_msk
 domdef
-latout=reform(gphit(0,*))
-lonout=reform(glamt(*,0))
-print, 'lat grid from mask ',min(latout),max(latout),latout(1)-latout(0)
-print, 'lon grid from mask ',min(lonout),max(lonout),lonout(1)-lonout(0)
+latout=reform(gphit[0,*])
+lonout=reform(glamt[*,0])
+print, 'lat grid from mask ',min(latout),max(latout),latout[1]-latout[0]
+print, 'lon grid from mask ',min(lonout),max(lonout),lonout[1]-lonout[0]
 mskout=read_ncdf("msk", file=fullfilename_msk,/nostr)
 ;
@@ -386 +389 @@
 t2mout=fltarr(jpi,jpj,jptin)
 for jt=0,jptin-1 do begin
-  tab=reform(t2min(*,*,jt))
-  t2mout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
+  tab=reform(t2min[*,*,jt])
+  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
+  t2mout[*,*,jt]=t2mout[*,*,jt]*mskout+(1.-mskout)*1.e20
 endfor
 ;

trunk/src/oaflux_mask_30n30s.pro

@@ -105 +105 @@
 ; $URL$
 ;
+; - fplod 20120306
+;
+;   * try to add compile_opt
+;
 ; - fplod 20120305
 ;
@@ -130 +134 @@
 ;
 FUNCTION oaflux_mask_30n30s
+;
+;++ compile_opt idl2, strictarrsubs, logical_predicate
 ;
 @cm_4cal
@@ -179 +185 @@
 valmask=1.e20
 ind=where(lh ge 32000,compl=nind)
-lh(ind)=valmask
-lh(nind)=-0.1*lh(nind)
+lh[ind]=valmask
+lh[nind]=-0.1*lh[nind]
 mask_out=(lh ne valmask)
 domdef, 0.,360.,-30.,30.
 
-lon=reform(glamt(*,0))
+lon=reform(glamt[*,0])
 ind1=where(lon le 20.)
 ind2=where(lon ge 30.)
 
-lat_out=reform(gphit(0,firstyt:lastyt))
-
-lon_out=[reform(glamt(ind2,0)),reform(glamt(ind1,0))+360]
-
+lat_out=reform(gphit[0,firstyt:lastyt])
+
+lon_out=[reform(glamt[ind2,0]),reform(glamt[ind1,0])+360]
+
+; ++ this line is incompatible with 
+; compile_opt idl2, strictarrsubs, logical_predicate
 mask_out=mask_out([ind2,ind1],firstyt:lastyt)*1.