| 1 | pro profile_uvt |
|---|
| 2 | ; |
|---|
| 3 | ; |
|---|
| 4 | ; |
|---|
| 5 | what_I_plot=0. |
|---|
| 6 | overcontour=0 |
|---|
| 7 | mention='' |
|---|
| 8 | @map_uvt_inc.pro |
|---|
| 9 | SPAWN, '\rm param_plot.idl ; cp map_uvt_inc.pro param_plot.idl ; cp -f map_uvt_inc.pro '+save_ps+'.map_uvt_inc.pro' |
|---|
| 10 | if (n_elements(coord2d) eq 0) then coord2d='false' |
|---|
| 11 | ; |
|---|
| 12 | ; |
|---|
| 13 | ; |
|---|
| 14 | zefile=save_ps |
|---|
| 15 | PS_Start, filename=zefile+'.ps' |
|---|
| 16 | print, zefile+'.ps' |
|---|
| 17 | !P.Charsize = 1.2 |
|---|
| 18 | !p.charthick = 2.0 |
|---|
| 19 | !p.thick = 2.0 |
|---|
| 20 | !x.thick = 2.0 |
|---|
| 21 | !y.thick = 2.0 |
|---|
| 22 | !p.multi=[0,2,2] |
|---|
| 23 | !P.Charsize = 1.0 |
|---|
| 24 | ; |
|---|
| 25 | ; |
|---|
| 26 | ; |
|---|
| 27 | if (n_elements(field1) ne 0) then getcdf, file=filename, charvar=field1, invar=cfield1 |
|---|
| 28 | u = getget(filename, 'Um', count=[1,1,0,1], offset=[zex, zey, 0,ntime]) |
|---|
| 29 | v = getget(filename, 'Vm', count=[1,1,0,1], offset=[zex, zey, 0,ntime]) |
|---|
| 30 | w = getget(filename, 'WAVE', count=[1,1,0,1], offset=[zex, zey, 0,ntime]) |
|---|
| 31 | u2 = getget(filename, 'Um', count=[1,1,0,1], offset=[zex2,zey2,0,ntime]) |
|---|
| 32 | v2 = getget(filename, 'Vm', count=[1,1,0,1], offset=[zex2,zey2,0,ntime]) |
|---|
| 33 | w2 = getget(filename, 'WAVE', count=[1,1,0,1], offset=[zex2,zey2,0,ntime]) |
|---|
| 34 | tk = getget(filename, 'tk', count=[1,1,0,1], offset=[zex, zey, 0,ntime]) |
|---|
| 35 | tk2 = getget(filename, 'tk', count=[1,1,0,1], offset=[zex2, zey2, 0,ntime]) |
|---|
| 36 | tpot = getget(filename, 'tpot', count=[1,1,0,1], offset=[zex, zey, 0,ntime]) |
|---|
| 37 | tpot2 = getget(filename, 'tpot', count=[1,1,0,1], offset=[zex2, zey2, 0,ntime]) |
|---|
| 38 | getcdf, file=filename, charvar='XLONG', invar=longi |
|---|
| 39 | getcdf, file=filename, charvar='XLAT', invar=lati |
|---|
| 40 | getcdf, file=filename, charvar='HGT', invar=hgt |
|---|
| 41 | getcdf, file=filename, charvar='vert', invar=vert |
|---|
| 42 | getcdf, file=filename, charvar='TSURF', invar=tsurf |
|---|
| 43 | getcdf, file=filename, charvar='USTM', invar=ustar |
|---|
| 44 | print, longi(zex ,zey ,0), lati(zex ,zey ,0) |
|---|
| 45 | print, longi(zex2,zey2,0), lati(zex2,zey2,0) |
|---|
| 46 | ; |
|---|
| 47 | ; |
|---|
| 48 | ; |
|---|
| 49 | loadct, 0 |
|---|
| 50 | contour, reform(hgt(*,*,ntime)), $ |
|---|
| 51 | reform(longi(*,*,ntime)), $ |
|---|
| 52 | reform(lati(*,*,ntime)), $ |
|---|
| 53 | nlevels=20, $ |
|---|
| 54 | xtitle='Longitude', $ |
|---|
| 55 | ytitle='Latitude', $ |
|---|
| 56 | ;xrange=[-146.,-126.], $ |
|---|
| 57 | ;yrange=[11.,27.], $ |
|---|
| 58 | xrange=[-141.,-133.], $ |
|---|
| 59 | yrange=[18.,26.], $ |
|---|
| 60 | xtickinterval=1., $ |
|---|
| 61 | ytickinterval=1., $ |
|---|
| 62 | max_value=22000., $ |
|---|
| 63 | min_value=-4000., $ |
|---|
| 64 | /cell_fill |
|---|
| 65 | xyouts, longi(zex ,zey ,0), lati(zex ,zey ,0), '+ full', color=255, charsize=1.5 |
|---|
| 66 | xyouts, longi(zex2 ,zey2 ,0), lati(zex2 ,zey2 ,0), '+ dashed', color=255, charsize=1.5 |
|---|
| 67 | ; |
|---|
| 68 | ; |
|---|
| 69 | ; |
|---|
| 70 | ;plot, tsurf(zex ,zey ,*), yrange=[150, 300] |
|---|
| 71 | ;oplot, tsurf(zex2 ,zey2 ,*), linestyle=1 |
|---|
| 72 | ; |
|---|
| 73 | ; |
|---|
| 74 | ; |
|---|
| 75 | what_I_plot = tk & column = vert |
|---|
| 76 | ;yeye = tsurf(zex,zey,ntime) & what_I_plot = [yeye,reform(what_I_plot)] & column = [0.01, vert] |
|---|
| 77 | overplot = tk2 & overplot_column = vert |
|---|
| 78 | ;yeye = tsurf(zex2,zey2,ntime) & overplot = [yeye,reform(overplot)] & overplot_column = [0.01, vert] |
|---|
| 79 | print, min(what_I_plot), max(what_I_plot) |
|---|
| 80 | print, min(overplot), max(overplot) |
|---|
| 81 | mention='three' |
|---|
| 82 | ; |
|---|
| 83 | ; |
|---|
| 84 | ; |
|---|
| 85 | profile, $ |
|---|
| 86 | what_I_plot, $ ; 1D vertical profile |
|---|
| 87 | column, $ ; altitudes |
|---|
| 88 | alt=alt, $ ; altitude range [altmin, altmax] |
|---|
| 89 | minfield=minfield_init2, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 90 | maxfield=maxfield_init2, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 91 | inprofile=overplot, $ ; another vertical profile to overplot |
|---|
| 92 | incolumn=overplot_column, $ ; altitudes of the other vertical profile (in case /= column) |
|---|
| 93 | discrete=discrete, $ ; show the profile points (= type of points in !psym) |
|---|
| 94 | title_plot=title_user, $ ; title of the plot ('Profile' is default) |
|---|
| 95 | title_axis=title_axis2, $ ; title of the [x,y] axis (['Field','Altitude'] is default) |
|---|
| 96 | mention=mention ; add text precision within the plot window (default is nothing or '') |
|---|
| 97 | xyouts, 210.,20.,'T!Ds!N = '+string(tsurf(zex,zey,ntime),'(I0)')+' K' |
|---|
| 98 | xyouts, 163.,20.,'T!Ds!N = '+string(tsurf(zex2,zey2,ntime),'(I0)')+' K' |
|---|
| 99 | ; |
|---|
| 100 | ; |
|---|
| 101 | ; |
|---|
| 102 | what_I_plot = w & column = vert |
|---|
| 103 | overplot = w2 & overplot_column = vert |
|---|
| 104 | ; |
|---|
| 105 | ; |
|---|
| 106 | ; |
|---|
| 107 | profile, $ |
|---|
| 108 | what_I_plot, $ ; 1D vertical profile |
|---|
| 109 | column, $ ; altitudes |
|---|
| 110 | alt=alt, $ ; altitude range [altmin, altmax] |
|---|
| 111 | minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 112 | maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 113 | inprofile=overplot, $ ; another vertical profile to overplot |
|---|
| 114 | incolumn=overplot_column, $ ; altitudes of the other vertical profile (in case /= column) |
|---|
| 115 | discrete=discrete, $ ; show the profile points (= type of points in !psym) |
|---|
| 116 | title_plot=title_user, $ ; title of the plot ('Profile' is default) |
|---|
| 117 | title_axis=title_axis, $ ; title of the [x,y] axis (['Field','Altitude'] is default) |
|---|
| 118 | mention=mention ; add text precision within the plot window (default is nothing or '') |
|---|
| 119 | ; |
|---|
| 120 | ; |
|---|
| 121 | ; |
|---|
| 122 | ;plot, tpot, vert |
|---|
| 123 | ;oplot, tpot2, vert, linestyle=2 |
|---|
| 124 | ;PS_End, /PNG |
|---|
| 125 | ;stop |
|---|
| 126 | what_I_plot = u*u + v*v & what_I_plot = sqrt(what_I_plot) & column = vert |
|---|
| 127 | overplot = u2*u2 + v2*v2 & overplot = sqrt(overplot) & overplot_column = vert |
|---|
| 128 | print, min(what_I_plot), max(what_I_plot) |
|---|
| 129 | print, min(overplot), max(overplot) |
|---|
| 130 | mention='two' |
|---|
| 131 | ; |
|---|
| 132 | ; |
|---|
| 133 | ; |
|---|
| 134 | profile, $ |
|---|
| 135 | what_I_plot, $ ; 1D vertical profile |
|---|
| 136 | column, $ ; altitudes |
|---|
| 137 | alt=alt, $ ; altitude range [altmin, altmax] |
|---|
| 138 | minfield=minfield_init2, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 139 | maxfield=maxfield_init2, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 140 | inprofile=overplot, $ ; another vertical profile to overplot |
|---|
| 141 | incolumn=overplot_column, $ ; altitudes of the other vertical profile (in case /= column) |
|---|
| 142 | discrete=discrete, $ ; show the profile points (= type of points in !psym) |
|---|
| 143 | title_plot=title_user, $ ; title of the plot ('Profile' is default) |
|---|
| 144 | title_axis=title_axis2, $ ; title of the [x,y] axis (['Field','Altitude'] is default) |
|---|
| 145 | mention=mention ; add text precision within the plot window (default is nothing or '') |
|---|
| 146 | xyouts, 23.5,20.,'u!D*!N = '+string(ustar(zex,zey,ntime),'(F4.1)')+' m s!U-1!N' |
|---|
| 147 | xyouts, 5.,20.,'u!D*!N = '+string(ustar(zex2,zey2,ntime),'(F4.1)')+' m s!U-1!N' |
|---|
| 148 | ; |
|---|
| 149 | ; |
|---|
| 150 | ; |
|---|
| 151 | PS_End, /PNG |
|---|
| 152 | end |
|---|
