| 1 | |
|---|
| 2 | |
|---|
| 3 | PS_Start, FILENAME='velmax.ps' |
|---|
| 4 | !P.Charsize = 1.2 |
|---|
| 5 | |
|---|
| 6 | |
|---|
| 7 | history_interval_s = 400. |
|---|
| 8 | history_interval_s = 100. |
|---|
| 9 | smoothampl=3700/history_interval_s |
|---|
| 10 | ;smoothampl=0. ;; no smooth |
|---|
| 11 | ;smoothampl=20. |
|---|
| 12 | smoothampl=5. |
|---|
| 13 | |
|---|
| 14 | lev=1 |
|---|
| 15 | ;lev=0 qu'il faut faire ??? |
|---|
| 16 | |
|---|
| 17 | |
|---|
| 18 | @report_tasi.inc |
|---|
| 19 | |
|---|
| 20 | |
|---|
| 21 | |
|---|
| 22 | |
|---|
| 23 | ;restore, filename='REAL_50m_145_145_201_12km/case_A/addturb.dat' |
|---|
| 24 | restore, filename=saveplot+cases(0)+'/addturb.dat' |
|---|
| 25 | restore, filename=saveplot+cases(0)+'/getturb.dat' |
|---|
| 26 | print, h(lev) |
|---|
| 27 | velmaxye = reform(velmax(lev,*)) |
|---|
| 28 | velmaxye = SMOOTH(TEMPORARY(velmaxye), [smoothampl], /EDGE_TRUNCATE) |
|---|
| 29 | plot, localtime, velmaxye, $ |
|---|
| 30 | xrange=[8.,17.], xtickinterval=1., ytickinterval=2.0, $ |
|---|
| 31 | yrange=[0.,14.], ytitle='Horizontal velocity (m s!U-1!N)', xtitle='Local time (h)' |
|---|
| 32 | |
|---|
| 33 | ;restore, filename='REAL_50m_145_145_201_12km/case_B/addturb.dat' |
|---|
| 34 | restore, filename=saveplot+cases(1)+'/addturb.dat' |
|---|
| 35 | restore, filename=saveplot+cases(1)+'/getturb.dat' |
|---|
| 36 | print, h(lev) |
|---|
| 37 | velmaxye = reform(velmax(lev,*)) |
|---|
| 38 | velmaxye = SMOOTH(TEMPORARY(velmaxye), [smoothampl], /EDGE_TRUNCATE) |
|---|
| 39 | oplot, localtime, velmaxye, linestyle=1 |
|---|
| 40 | |
|---|
| 41 | ;restore, filename='REAL_50m_145_145_201_12km/case_C/addturb.dat' |
|---|
| 42 | restore, filename=saveplot+cases(2)+'/addturb.dat' |
|---|
| 43 | restore, filename=saveplot+cases(2)+'/getturb.dat' |
|---|
| 44 | print, h(lev) |
|---|
| 45 | velmaxye = reform(velmax(lev,*)) |
|---|
| 46 | velmaxye = SMOOTH(TEMPORARY(velmaxye), [smoothampl], /EDGE_TRUNCATE) |
|---|
| 47 | oplot, localtime, velmaxye, linestyle=2 |
|---|
| 48 | |
|---|
| 49 | ;restore, filename='REAL_50m_145_145_201_12km/case_I/addturb.dat' |
|---|
| 50 | restore, filename=saveplot+cases(3)+'/addturb.dat' |
|---|
| 51 | restore, filename=saveplot+cases(3)+'/getturb.dat' |
|---|
| 52 | print, h(lev) |
|---|
| 53 | velmaxye = reform(velmax(lev,*)) |
|---|
| 54 | velmaxye = SMOOTH(TEMPORARY(velmaxye), [smoothampl], /EDGE_TRUNCATE) |
|---|
| 55 | oplot, localtime, velmaxye, linestyle=3 |
|---|
| 56 | |
|---|
| 57 | PS_END, /PNG |
|---|
