| 1 | pro turb_section_period |
|---|
| 2 | ; |
|---|
| 3 | ; |
|---|
| 4 | ; |
|---|
| 5 | what_I_plot=0. |
|---|
| 6 | overcontour=0. |
|---|
| 7 | SPAWN, '\rm param_plot.idl ; cp turb_section_period_inc.pro param_plot.idl' |
|---|
| 8 | ; |
|---|
| 9 | ; |
|---|
| 10 | ; |
|---|
| 11 | filefile = "./MERIDIANI_tau050_wind10_section_plume.nc" |
|---|
| 12 | getcdf, file=filefile, charvar="W", invar=zew |
|---|
| 13 | getcdf, file=filefile, charvar="PHTOT", invar=zevar |
|---|
| 14 | zew = reform(zew) |
|---|
| 15 | zevar = reform(zevar) |
|---|
| 16 | |
|---|
| 17 | nx = n_elements(zevar(*,0)) |
|---|
| 18 | nz = n_elements(zevar(0,*)) |
|---|
| 19 | |
|---|
| 20 | xx = findgen(nx) * 100. / 1000. |
|---|
| 21 | zz = reform(zevar(0,*)) / 3.72 / 1000. + ( 1473.35 / 1000. ) |
|---|
| 22 | |
|---|
| 23 | zew = zew[*,0:nz-2] |
|---|
| 24 | zz = zz[0:nz-2] |
|---|
| 25 | |
|---|
| 26 | |
|---|
| 27 | PS_Start, filename='sectionLES.ps' |
|---|
| 28 | !P.Charsize = 1.2 |
|---|
| 29 | !p.charthick = 2.0 |
|---|
| 30 | !p.thick = 2.0 |
|---|
| 31 | !x.thick = 2.0 |
|---|
| 32 | !y.thick = 2.0 |
|---|
| 33 | |
|---|
| 34 | what_I_plot = zew |
|---|
| 35 | lon = xx |
|---|
| 36 | lat = zz |
|---|
| 37 | |
|---|
| 38 | map_latlon, $ |
|---|
| 39 | what_I_plot, $ ; 2D field |
|---|
| 40 | lon, $ ; 1D latitude OR 2D |
|---|
| 41 | lat, $ ; 1D longitude OR 2D |
|---|
| 42 | minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 43 | maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 44 | overcontour=overcontour, $ ; another 2D field to overplot with contour lines (=0: no) |
|---|
| 45 | overvector_x=overvector_x, $ ; wind vector - x component (=0: no) |
|---|
| 46 | overvector_y=overvector_y, $ ; wind vector - y component (=0: no) |
|---|
| 47 | ct=pal, $ ; color table (33-rainbow is default) |
|---|
| 48 | colors=colors, $ ; number of colors/levels (32 is default) |
|---|
| 49 | title=title_user, $ ; title of the plot ('' is default) |
|---|
| 50 | format=format ; format of colorbar annotations ('(F6.2)' is default) |
|---|
| 51 | |
|---|
| 52 | PS_End, /PNG |
|---|
| 53 | end |
|---|
