| [85] | 1 | pro similarity |
|---|
| 2 | |
|---|
| 3 | !p.multi = [0,2,1] |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | lim=0.30 & zelim=(-3.85/alog(lim)+0.07*alog(lim))*exp(-4.61*lim) |
|---|
| 7 | up=1. & & zeup= - 1.55 * up + 1.27 |
|---|
| 8 | slope = (zelim - zeup) / (lim - up) & ord = zeup - slope * up |
|---|
| 9 | t1 = findgen(1001.)*lim/1000. & param1 = (-3.85/alog(t1)+0.07*alog(t1))*exp(-4.61*t1) |
|---|
| 10 | t2 = lim + findgen(1001.)*(1.0 - lim)/1000. & param2 = slope * t2 + ord |
|---|
| 11 | ;;print, slope, ' slope' |
|---|
| 12 | ;;print, ord, ' ord' |
|---|
| 13 | ;; -1.51562 slope 1.23562 ord |
|---|
| 14 | |
|---|
| 15 | |
|---|
| 16 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 17 | set_title = "Similarity functions: vertical eddy heat flux" |
|---|
| 18 | set_titlex = "Dimensionless vertical heat flux <w'T'>/<w'T'>!Dmax!N" |
|---|
| 19 | set_titley = 'Dimensionless height z/z!Di!N' |
|---|
| 20 | set_subtitle = '' |
|---|
| 21 | set_xrange = [-0.5,1.] |
|---|
| 22 | set_yrange = [0.0,1.0] |
|---|
| 23 | set_tickx = 0.25 |
|---|
| 24 | set_ticky = 0.1 |
|---|
| 25 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 26 | ;zesim = 1 - 1.2 * a_h |
|---|
| 27 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 28 | PS_Start, FILENAME='similarity.ps' |
|---|
| 29 | !P.Charsize = 1.2 |
|---|
| 30 | plot, param1, t1, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle;, /ISOTROPIC |
|---|
| 31 | oplot, param2, t2 |
|---|
| 32 | oplot, [0.,0.], set_yrange, linestyle=1 |
|---|
| 33 | oplot, 1 - 1.2 * t1, t1, linestyle=2 |
|---|
| 34 | oplot, 1 - 1.2 * t2, t2, linestyle=2 |
|---|
| 35 | xyouts, 0.6, 0.47, 'MARS' |
|---|
| 36 | xyouts, 0.25, 0.3, 'EARTH' |
|---|
| 37 | ;PS_End, /PNG |
|---|
| 38 | |
|---|
| 39 | |
|---|
| 40 | t1 = findgen(2001.)*1./2000. |
|---|
| 41 | param1 = 2.05 * t1^(2./3.) * ( 1. - 0.64 * t1 )^2 |
|---|
| 42 | param2 = 1.80 * t1^(2./3.) * ( 1. - 0.80 * t1 )^2 |
|---|
| 43 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 44 | set_title = "Similarity functions: vert. vel. variance" |
|---|
| 45 | set_titlex = "Dimensionless vert. vel. variance <w'!U2!N>/w!D*!N!U2!N" |
|---|
| 46 | set_titley = 'Dimensionless height z/z!Di!N' |
|---|
| 47 | set_subtitle = '' |
|---|
| 48 | set_xrange = [0.,1.] |
|---|
| 49 | set_yrange = [0.0,1.0] |
|---|
| 50 | set_tickx = 0.25 |
|---|
| 51 | set_ticky = 0.1 |
|---|
| 52 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 53 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 54 | ;PS_Start, FILENAME='similarity.ps' |
|---|
| 55 | !P.Charsize = 1.2 |
|---|
| 56 | plot, param1, t1, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle;, /ISOTROPIC |
|---|
| 57 | oplot, param2, t1, linestyle=2 |
|---|
| 58 | xyouts, 0.2, 0.3, 'EARTH' |
|---|
| 59 | xyouts, 0.67, 0.4, 'MARS' |
|---|
| 60 | PS_End, /PNG |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | |
|---|
| 64 | end |
|---|
