Changeset 169
- Timestamp:
- 03/20/12 10:39:23 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/interp_erai_t2m.pro
r166 r169 356 356 print, 'lon grid from data',min(lonin),max(lonin),lonin[1]-lonin[0] 357 357 ; 358 ; get time in t2 file359 timein=ncdf_lec(fullfilename_t2,var='time')360 jptin=n_elements(timein)361 print, 'time steps from data ', jptin362 print, 'The first 10 time values (variable timein) = ', timein[0:9]358 ;; get time in t2 file 359 ;timein=ncdf_lec(fullfilename_t2,var='time') 360 ;jptin=n_elements(timein) 361 ;print, 'time steps from data ', jptin 362 ;print, 'The first 10 time values (variable timein) = ', timein[0:9] 363 363 ; 364 364 ; find first and last dates yyyymmdd 365 ; they will be written in global attributes of output file366 da=jul2date(julday(01, 01, 1957,timein[0]))367 cda0=string(da,format='(i8.8)')368 da=jul2date(julday(01, 01, 1957,timein[jptin-1]))369 cda1=string(da,format='(i8.8)')370 print, 'first date ', cda0371 print, 'last date ' , cda1365 ;; they will be written in global attributes of output file 366 ;da=jul2date(julday(01, 01, 1957,timein[0])) 367 ;cda0=string(da,format='(i8.8)') 368 ;da=jul2date(julday(01, 01, 1957,timein[jptin-1])) 369 ;cda1=string(da,format='(i8.8)') 370 ;print, 'first date ', cda0 371 ;print, 'last date ' , cda1 372 372 ; 373 373 ; read t2 data 374 374 ;++ pb memory t2min=read_ncdf("t2",0,jptin-1,/timestep,file=fullfilename_t2,/nostr) 375 375 ; the following line is here just to prepare replacement of read_ncdf by ncdf_lec 376 t2min_read_ncdf=read_ncdf("t2",0,10,/timestep,file=fullfilename_t2,/nostr)377 help, t2min_read_ncdf376 ;t2min_read_ncdf=read_ncdf("t2",0,10,/timestep,file=fullfilename_t2,/nostr) 377 ;help, t2min_read_ncdf 378 378 ;++print, 'The first 10 time values (variable time) = ', time[0:9] 379 379 ;++print, 'time steps after read_ncdf (variable jpt) ', jpt 380 t2min=ncdf_lec(fullfilename_t2,var='t2') 380 ;t2min=ncdf_lec(fullfilename_t2,var='t2') 381 382 t2min=read_ncdf('t2',yyyymmddb-.5d,yyyymmdde,file=fullfilename_t2,/nost) 383 ;; time:units = "hours since 1957-01-01 00:00:0.0" ; 384 timein=24.d*(time-julday(1,1,1957,0,0,0)) & jptin=n_elements(timein) 385 da=jul2date(time(0)) & cda0=string(da,format='(i8.8)') 386 da=jul2date(time(jpt-1)) & cda1=string(da,format='(i8.8)') 387 print, 'first date ', cda0 388 print, 'last date ' , cda1 389 help, t2min 390 391 381 392 ; 382 393 mskin=glamt*0.+1.
Note: See TracChangeset
for help on using the changeset viewer.