Changeset 97 for trunk/src/interp_erai_lwr_1989_2009.pro

Timestamp:

08/30/11 14:38:29 (13 years ago)

Author:

pinsard

Message:

suppress blank lines trailing blank

File:

• trunk/src/interp_erai_lwr_1989_2009.pro (modified) (7 diffs)

trunk/src/interp_erai_lwr_1989_2009.pro

@@ -12 +12 @@
 ; Interpolation of str from ERA-I grid to OAFLUX grid
 ;
-; :file:`${PROJECT_ID}/20c3m_erai_str_TROP_1989_2009.nc` 
-; containing str from ERA-I 
+; :file:`${PROJECT_ID}/20c3m_erai_str_TROP_1989_2009.nc`
+; containing str from ERA-I
 ; has been produced
 ; by :ref:`compute_erai_daily_region_2d.sh`.
 ;
-; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` 
+; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc`
 ; containing OAFLUX mask
-; has been produced by 
+; has been produced by
 ; :ref:`oaflux_mask_30N30S.pro`.
 ;
 ; Interpolated str is written in
-; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc` 
+; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc`
 ; if this file not already exists.
 ;
-; This output file 
-; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc` 
+; This output file
+; :file:`${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc`
 ; must be processed after by :ref:`add_198900101.sh`.
 ;
@@ -207 +207 @@
    STOP
 ENDIF
-
-st=19890101 
+;
+st=19890101
 en=20090801
-
+;
 ; build data filename
 filename='20c3m_erai_str_TROP_1989_2009.nc'
@@ -260 +260 @@
 lwrin=read_ncdf("str",0,7668,/timestep,file=fullfilename,/nostr)
 lwrin=-1*lwrin
-
+;
 initncdf, fullfilename
 domdef
-latin=reform(gphit(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)
-
-timein=time 
+;
+timein=time
 jptin=jpt
 tab=lwrin(*,*,0)
 mskin=glamt*0.+1.
-
+;
 initncdf, fullfilename_msk
 domdef
@@ -279 +279 @@
 print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0)
 mskout=read_ncdf("msk", file=fullfilename_msk,/nostr)
-
+;
 si=size(lwrin)
 lwrout=fltarr(jpi,jpj,jptin)
-
+;
 for jt=0,jptin-1 do begin
   ; ++print, 'Interpolation jt=',jt,' / ',jptin-1
-
+;
   tab=reform(lwrin(*,*,jt))
   lwrout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $
@@ -291 +291 @@
       , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
   lwrout(*,*,jt)=lwrout(*,*,jt)*mskout+(1.-mskout)*1.e20
-
+;
 endfor
-
-time=timegen(7669, start=julday(1,2,1989,0), units='days') 
+;
+time=timegen(7669, start=julday(1,2,1989,0), units='days')
 jpt=n_elements(time)
-tt=time 
+tt=time
 jptin=jpt
 cda0=string(jul2date(tt(0)),format='(i8.8)')
 cda1=string(jul2date(tt(jpt-1)),format='(i8.8)')
 tt=tt-julday(1,1,1950,00,00,00)
-xlon=reform(glamt(*,0) ) 
+xlon=reform(glamt(*,0) )
 ylat=reform(gphit(0,*))
-
+;
 initncdf, fullfilename_msk
 valmask=1.e20
-
+;
 ylat=latout
 xlon=lonout
@@ -313 +313 @@
 lat_attr={units:'degrees_north',long_name:'Latitude'}
 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'}
-
+;
 lwr_attr={units:'W/m2',missing_value:valmask,long_name:'Surface net longwave radiation ',short_name:'lwr',axis:'TYX'}
 time_attr={units:'days since 1950-01-01 00:00:00',long_name:'Time axis',time_origin:' 1950-JAN-01 00:00:00'}
 globattr={source:'Tropical ocean fluxes obtained from ERA Interim',timerange:cda0+' - '+cda1}
-
+;
 ncfields = 'lwr[longitude,latitude,time]=lwrout:lwr_attr; ' $
                       + 'longitude[]=xlon:lon_attr; ' $
@@ -323 +323 @@
                       + 'tt[*time]=tt:time_attr ' $
                       + ' @ globattr'
-
+;
 @ncdf_quickwrite
-
+;
 end