Changeset 97 for trunk/src/interp_erai_dewt_1989_2009.pro
- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/interp_erai_dewt_1989_2009.pro
r94 r97 235 235 initncdf, fullfilename 236 236 domdef 237 latin=reform(gphit(0,*)) 237 latin=reform(gphit(0,*)) 238 238 lonin=reform(glamt(*,0)) 239 239 print, 'lat grid ',min(latin),max(latin),latin(1)-latin(0) 240 240 print, 'lon grid ',min(lonin),max(lonin),lonin(1)-lonin(0) 241 241 d2min=read_ncdf("d2",19881201,20101231,file=fullfilename,/nostr) 242 242 ; 243 243 timein=time 244 244 jptin=jpt 245 245 tab=d2min(*,*,0) 246 246 mskin=glamt*0.+1. 247 247 ; 248 248 initncdf, fullfilename_msk 249 249 domdef … … 253 253 print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 254 254 mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 255 255 ; 256 256 help, d2min,lonin,latin,mskin,lonout,latout,mskout 257 257 ; 258 258 si=size(d2min) 259 259 d2mout=fltarr(jpi,jpj,jptin) … … 266 266 d2mout(*,*,jt)=d2mout(*,*,jt)*mskout+(1.-mskout)*1.e20 267 267 endfor 268 268 ; 269 269 timein=timein 270 270 jptin=jpt 271 271 ; 272 272 initncdf, fullfilename_msk 273 273 cda0=string(jul2date(timein(0)),format='(i8.8)') 274 274 cda1=string(jul2date(timein(jpt-1)),format='(i8.8)') 275 276 time=timegen(7670, units='days', start=julday(1,1,1989)) 275 ; 276 time=timegen(7670, units='days', start=julday(1,1,1989)) 277 277 jpt=n_elements(time) 278 278 cda0=string(jul2date(time(0)),format='(i8.8)') 279 279 cda1=string(jul2date(time(jpt-1)),format='(i8.8)') 280 280 timein=time-julday(1,1,1950,00,00) 281 281 ; 282 282 lat=latout 283 283 lon=lonout … … 288 288 d2m_attr={units:'degK',missing_value:1e20,long_name:'Dew Point Temperature at 2m',short_name:'d2m',axis:'TYX'} 289 289 globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 290 290 ; 291 291 ncfields = 'd2m[longitude,latitude,time]=d2mout:d2m_attr; ' $ 292 292 + 'longitude[]=lon:lon_attr; ' $ 293 293 + 'latitude[]=lat:lat_attr; ' $ 294 294 + 'time[*time]=timein:time_attr ' $ 295 295 + ' @ globattr' 296 296 ; 297 297 @ncdf_quickwrite 298 298 ; 299 299 end
Note: See TracChangeset
for help on using the changeset viewer.