source: trunk/UTIL/IDLplot/SPEC/LES/tsurf.pro @ 1306

Last change on this file since 1306 was 85, checked in by aslmd, 15 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
File size: 894 bytes
Line 
1
2domain='d01' & filesWRF = FindFile('wrfout_'+domain+'_????-??-??_??:??:??') & nf=n_elements(filesWRF)
3id=ncdf_open(filesWRF(0))
4NCDF_DIMINQ, id, NCDF_DIMID(id, 'west_east'    ), toto, nx & NCDF_DIMINQ, id, NCDF_DIMID(id, 'south_north'  ), toto, ny
5NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time'         ), toto, nt
6NCDF_CLOSE, id
7id=ncdf_open(filesWRF(nf-1))  ;; for interrupted runs
8NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time'         ), toto, ntlast
9NCDF_CLOSE, id
10yeye = 0 & nttot = (nf-1)*nt + ntlast
11
12history_interval_s = 100.
13OPENR, 22, 'input_coord' & READF, 22, lonu & READF, 22, latu & READF, 22, lsu & READF, 22, lctu & CLOSE, 22
14
15localtime = lctu + history_interval_s*findgen(nttot)/3700.
16
17
18cdfid = ncdf_open('surf.nc')
19varid=ncdf_varid(cdfid,'TSURF')
20ncdf_varget, cdfid, varid, tsurf
21
22;print, tsurf(20,20,*)
23
24yeyey = reform(tsurf(20,20,*))
25
26w = where(yeyey eq max(yeyey))
27print, localtime[w]
28
Note: See TracBrowser for help on using the repository browser.