Ignore:
Timestamp:
08/30/11 14:38:29 (13 years ago)
Author:
pinsard
Message:

suppress blank lines trailing blank

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/interp_erai_dewt_1989_2009.pro

    r94 r97  
    235235initncdf, fullfilename 
    236236domdef 
    237 latin=reform(gphit(0,*))  
     237latin=reform(gphit(0,*)) 
    238238lonin=reform(glamt(*,0)) 
    239239print, 'lat grid ',min(latin),max(latin),latin(1)-latin(0) 
    240240print, 'lon grid ',min(lonin),max(lonin),lonin(1)-lonin(0) 
    241241d2min=read_ncdf("d2",19881201,20101231,file=fullfilename,/nostr) 
    242  
     242; 
    243243timein=time 
    244244jptin=jpt 
    245245tab=d2min(*,*,0) 
    246246mskin=glamt*0.+1. 
    247  
     247; 
    248248initncdf, fullfilename_msk 
    249249domdef 
     
    253253print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    254254mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    255  
     255; 
    256256help, d2min,lonin,latin,mskin,lonout,latout,mskout 
    257  
     257; 
    258258si=size(d2min) 
    259259d2mout=fltarr(jpi,jpj,jptin) 
     
    266266  d2mout(*,*,jt)=d2mout(*,*,jt)*mskout+(1.-mskout)*1.e20 
    267267endfor 
    268  
     268; 
    269269timein=timein 
    270270jptin=jpt 
    271  
     271; 
    272272initncdf, fullfilename_msk 
    273273cda0=string(jul2date(timein(0)),format='(i8.8)') 
    274274cda1=string(jul2date(timein(jpt-1)),format='(i8.8)') 
    275  
    276 time=timegen(7670, units='days', start=julday(1,1,1989))  
     275; 
     276time=timegen(7670, units='days', start=julday(1,1,1989)) 
    277277jpt=n_elements(time) 
    278278cda0=string(jul2date(time(0)),format='(i8.8)') 
    279279cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 
    280280timein=time-julday(1,1,1950,00,00) 
    281  
     281; 
    282282lat=latout 
    283283lon=lonout 
     
    288288d2m_attr={units:'degK',missing_value:1e20,long_name:'Dew Point Temperature at 2m',short_name:'d2m',axis:'TYX'} 
    289289globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    290  
     290; 
    291291ncfields = 'd2m[longitude,latitude,time]=d2mout:d2m_attr; ' $ 
    292292                      + 'longitude[]=lon:lon_attr; ' $ 
    293293                      + 'latitude[]=lat:lat_attr; ' $ 
    294                       + 'time[*time]=timein:time_attr ' $ 
     294                      + 'time[*time]=timein:time_attr ' $ 
    295295                      + ' @ globattr' 
    296  
     296; 
    297297@ncdf_quickwrite 
    298  
     298; 
    299299end 
Note: See TracChangeset for help on using the changeset viewer.