source: trunk/MESOSCALE_DEV/PLOT/SPEC/LES/include_hfx.pro @ 211

Last change on this file since 211 was 85, checked in by aslmd, 14 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
File size: 1.9 KB
Line 
1        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
2        OPENR, 22, 'input_coord' & READF, 22, lonu & READF, 22, latu & READF, 22, lsu & READF, 22, lctu & CLOSE, 22
3        OPENR, 23, 'input_more' & READF, 23, hgtu, tsurfu & CLOSE, 23
4        domain='d01' & filesWRF = FindFile('wrfout_'+domain+'_????-??-??_??:??:??') & nf=n_elements(filesWRF)
5        id=ncdf_open(filesWRF(0))
6                NCDF_DIMINQ, id, NCDF_DIMID(id, 'west_east'    ), toto, nx & NCDF_DIMINQ, id, NCDF_DIMID(id, 'south_north'  ), toto, ny
7                NCDF_DIMINQ, id, NCDF_DIMID(id, 'bottom_top'   ), toto, nz & NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time'         ), toto, nt
8                NCDF_CLOSE, id
9        id=ncdf_open(filesWRF(nf-1))  ;; for interrupted runs
10                NCDF_DIMINQ, id, NCDF_DIMID(id, 'Time'         ), toto, ntlast
11                NCDF_CLOSE, id
12        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
13        yeye = 0 & nttot = (nf-1)*nt + ntlast & localtime = lctu + history_interval_s*findgen(nttot)/3700.
14        velmax = fltarr(nz,nttot) & velmin = fltarr(nz,nttot)
15        ustarm = fltarr(nttot) & ustarmax = fltarr(nttot)
16        for loop  = 0, nf-1 do begin
17                                                  timetime = SYSTIME(1)
18          if (loop ne nf-1) then nloop2=nt else nloop2=ntlast
19        for loop2 = 0, nloop2-1 do begin
20!QUIET=1
21ustar = getget(filesWRF(loop), 'HFX', count=[0,0,0,1], offset=[0,0,0,loop2])
22ustarmax(yeye) = max(reform(ustar(*,*,0,0)))
23ustarm(yeye) = reform(TOTAL(TOTAL(TEMPORARY(ustar),1),1) / float(nx) / float(ny))
24print, ustarm(yeye), ustarmax(yeye)
25!QUIET=0
26        yeye = TEMPORARY(yeye) + 1
27        endfor
28        endfor
29        save, ustarm, localtime, filename='addturb.dat' ;; pour l'instant ecrase systematiquement
30        stop
31        ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
32
Note: See TracBrowser for help on using the repository browser.