source: trunk/mesoscale/PLOT/SPEC/LES/comp_ustarmax.pro @ 129

Last change on this file since 129 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.5 KB
RevLine 
[85]1
2PS_Start, FILENAME='ustarmax.ps'
3!P.Charsize = 1.2
4
5
6
7history_interval_s = 400.
8history_interval_s = 100.
9smoothampl=3700/history_interval_s
10smoothampl=0.  ;; no smooth
11;smoothampl=20.
12
13lev=1
14
15
16;restore, filename='REAL_50m_145_145_201_12km/case_A/addturb.dat'
17;ustarmaxye = SMOOTH(TEMPORARY(ustarmax), [smoothampl], /EDGE_TRUNCATE)
18;ustarmye = SMOOTH(TEMPORARY(ustarm), [smoothampl], /EDGE_TRUNCATE)
19;plot, localtime, ustarmaxye, yrange=[0.0,1.0], xtickinterval=1., ytickinterval=0.1, ytitle='m/s', xtitle='Local time (h)'
20;oplot, localtime, ustarmye
21
22restore, filename='REAL_50m_145_145_201_12km/case_B/addturb.dat'
23ustarmaxye = SMOOTH(TEMPORARY(ustarmax), [smoothampl], /EDGE_TRUNCATE)
24ustarmye = SMOOTH(TEMPORARY(ustarm), [smoothampl], /EDGE_TRUNCATE)
25plot, localtime, ustarmaxye, yrange=[0.0,1.2], xtickinterval=1., ytickinterval=0.1, ytitle='m/s', xtitle='Local time (h)'
26oplot, localtime, ustarmye
27
28restore, filename='REAL_50m_145_145_201_12km/case_C/addturb.dat'
29ustarmaxye = SMOOTH(TEMPORARY(ustarmax), [smoothampl], /EDGE_TRUNCATE)
30ustarmye = SMOOTH(TEMPORARY(ustarm), [smoothampl], /EDGE_TRUNCATE)
31oplot, localtime, ustarmaxye, linestyle=2
32oplot, localtime, ustarmye, linestyle=2
33
34;restore, filename='REAL_50m_145_145_201_12km/case_I/addturb.dat'
35;ustarmaxye = SMOOTH(TEMPORARY(ustarmax), [smoothampl], /EDGE_TRUNCATE)
36;oplot, localtime, ustarmaxye, linestyle=3
37
38;restore, filename='REAL_50m_145_145_201_12km/case_HIGH/addturb.dat'
39;ustarmaxye = SMOOTH(TEMPORARY(ustarmax), [smoothampl], /EDGE_TRUNCATE)
40;oplot, localtime, ustarmaxye, linestyle=4
41
42
43
44
45PS_END, /PNG
Note: See TracBrowser for help on using the repository browser.