Changeset 169


Ignore:
Timestamp:
03/20/12 10:39:23 (12 years ago)
Author:
jv
Message:

COrrection inversion NS t2m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/interp_erai_t2m.pro

    r166 r169  
    356356print, 'lon grid from data',min(lonin),max(lonin),lonin[1]-lonin[0] 
    357357; 
    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] 
     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] 
    363363; 
    364364; find first and last dates yyyymmdd 
    365 ; 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 
     365;; 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 
    372372; 
    373373; read t2 data 
    374374;++ pb memory t2min=read_ncdf("t2",0,jptin-1,/timestep,file=fullfilename_t2,/nostr) 
    375375; 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_ncdf 
     376;t2min_read_ncdf=read_ncdf("t2",0,10,/timestep,file=fullfilename_t2,/nostr) 
     377;help, t2min_read_ncdf 
    378378;++print, 'The first 10 time values (variable time) = ', time[0:9] 
    379379;++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 
     382t2min=read_ncdf('t2',yyyymmddb-.5d,yyyymmdde,file=fullfilename_t2,/nost) 
     383;; time:units = "hours since 1957-01-01 00:00:0.0" ; 
     384timein=24.d*(time-julday(1,1,1957,0,0,0)) & jptin=n_elements(timein) 
     385da=jul2date(time(0)) & cda0=string(da,format='(i8.8)') 
     386da=jul2date(time(jpt-1)) & cda1=string(da,format='(i8.8)') 
     387print, 'first date ', cda0 
     388print, 'last date ' , cda1 
     389help, t2min 
     390 
     391 
    381392; 
    382393mskin=glamt*0.+1. 
Note: See TracChangeset for help on using the changeset viewer.