Index: trunk/MESOSCALE_DEV/PLOT/SPEC/LES/getturb.pro =================================================================== --- trunk/MESOSCALE_DEV/PLOT/SPEC/LES/getturb.pro (revision 324) +++ trunk/MESOSCALE_DEV/PLOT/SPEC/LES/getturb.pro (revision 1181) @@ -20,4 +20,7 @@ ;p0=610. & t0=220. & r_cp=1.0/3.9 & grav = 3.72 & R=191. ;print, 'ATTENTION ATTENTION R/cp !!!!', r_cp + +; INTERCOMP INTERCOMP +r_cp = 192./770. ; @@ -32,5 +35,6 @@ endif else begin ;PREF_SET, 'IDL_PATH', '/home/spiga/Save/SOURCES/IDL/fsc_psconfig:', /COMMIT - PREF_SET, 'IDL_PATH', '/home/spiga/SVN/trunk/mesoscale/PLOT/MINIMAL:', /COMMIT + ;PREF_SET, 'IDL_PATH', '/home/spiga/SVN/trunk/mesoscale/PLOT/MINIMAL:', /COMMIT + PREF_SET, 'IDL_PATH', '/home/spiga/MODELES/MESOSCALE_DEV/PLOT/MINIMAL:', /COMMIT endelse @@ -172,9 +176,6 @@ ;@include_ustar.pro - ;@include_hfx.pro - ;@include_w.pro - ;@include_tprime.pro @@ -374,25 +375,25 @@ ;;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = localtime -;zey = pbl_height -;set_name = 'HEIGHT.ps' -;set_title = 'Boundary layer height (km)' -;set_titlex = 'Local Time (h)' -;set_titley = 'Boundary layer height (km)' -;set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' -;set_xrange = [11.,17.] ;[8.,17.] -;set_yrange = [0.,8.] -;set_tickx = 1. -;set_ticky = 1. -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin -;;oplot, zefield, zey, psym=5 -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; + +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +zefield = localtime +zey = pbl_height +set_name = 'HEIGHT.ps' +set_title = 'Boundary layer height (km)' +set_titlex = 'Local Time (h)' +set_titley = 'Boundary layer height (km)' +set_subtitle = '' ;'Criterion is : static stability > 1.5 K.m!U-1!N' +set_xrange = [7.,19.] ;[11.,17.] ;[8.,17.] +set_yrange = [0.,9.] +set_tickx = 1. +set_ticky = 1. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +if (saveps eq 'true') then PS_Start, FILENAME=set_name +!P.Charsize = 1.2 +plot, zefield, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, linestyle=altlin +;oplot, zefield, zey, psym=5 +if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;zefield = localtime @@ -437,39 +438,39 @@ ;;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = transpose(tke) -;zex = localtime -;zey = h -;set_name = 'TKE.ps' -;set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[++] -;set_titlex = 'Local Time (h)' -;set_titley = 'Altitude above surface (km)' -;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [11.,17.] ;[8.,17.] -;set_yrange = [0.,8.] ;; [0.,200.] -;set_tickx = 1. -;set_ticky = 1. ;; 50. -;minval = 0. -;maxval = 20. -;nlev = maxval-minval -;pal = 22 -;rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;;; 0. levels -;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;;; 1. background -;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -; ;; 2. color field -; loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) -; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;;; 3. contour field -;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;;; 4. choose output -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; + +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +zefield = transpose(tke) +zex = localtime +zey = h +set_name = 'TKE.ps' +set_title = "Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[++] +set_titlex = 'Local Time (h)' +set_titley = 'Altitude above surface (km)' +set_subtitle = '' ;'Mean over the simulation domain' +set_xrange = [7.,19.] ;[8.,17.] +set_yrange = [0.,9.] ;; [0.,200.] +set_tickx = 1. +set_ticky = 1. ;; 50. +minval = 0. +maxval = 20. +nlev = maxval-minval +pal = 22 +rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +if (saveps eq 'true') then PS_Start, FILENAME=set_name +!P.Charsize = 1.2 +;; 0. levels +lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;; 1. background +loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 + ;; 2. color field + loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) + contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;; 3. contour field +loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;; 4. choose output +if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;zefield = 2.*transpose(ztke) @@ -478,38 +479,38 @@ ;;maxval = 12. ;;pal = 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = transpose(ztke) -;zex = localtime -;zey = h -;set_name = 'zTKE.ps' -;set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] -;set_titlex = 'Local Time (h)' -;set_titley = 'Altitude above surface (km)' -;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [8.,17.] -;set_yrange = [0.,7.] -;set_tickx = 1. -;set_ticky = 1. -;minval = 0. -;maxval = 8. -;nlev = maxval-minval -;pal = 22 -;rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;;; 0. levels -;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;;; 1. background -;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;;; 2. color field -;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) -; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;;; 3. contour field -;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;;; 4. choose output -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +zefield = transpose(ztke) +zex = localtime +zey = h +set_name = 'zTKE.ps' +set_title = "Vertical Turbulent Kinetic Energy (m!U2!N.s!U-2!N)" ;0.5[] +set_titlex = 'Local Time (h)' +set_titley = 'Altitude above surface (km)' +set_subtitle = '' ;'Mean over the simulation domain' +set_xrange = [7.,19.] +set_yrange = [0.,9.] +set_tickx = 1. +set_ticky = 1. +minval = 0. +maxval = 8. +nlev = maxval-minval +pal = 22 +rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +if (saveps eq 'true') then PS_Start, FILENAME=set_name +!P.Charsize = 1.2 +;; 0. levels +lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;; 1. background +loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;; 2. color field +loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) + contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;; 3. contour field +loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;; 4. choose output +if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + ;;restore, filename='addturb.dat' ;;modvar = SMOOTH(TEMPORARY(modvar), [0,smoothampl], /EDGE_TRUNCATE) @@ -547,50 +548,50 @@ ;;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -; -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = transpose(wt) -;zex = localtime -;zey = h -;set_name = 'HF.ps' -;set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux " ; -;set_titlex = 'Local Time (h)' -;set_titley = 'Altitude above surface (km)' -;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [8.,17.] -;set_yrange = [0.,8.] ;; [0.,200.] -;set_tickx = 1. -;set_ticky = 1. ;; 50. -;minval = -1.5 ;-2. -;maxval = 2.5 ;2. -;nlev = floor(maxval-minval)*10 -;pal = 33 -;rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;;minval = -0.8 -;;maxval = 1.2 -;;pal = 0 -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;;; 0. levels -;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;;; 1. background -;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;;; 2. color field -;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) -; ;;;-------------------------------------------------------------------------------------------------------------------------------- -; ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone -; ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5 -; nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) -; TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b -; ;;;-------------------------------------------------------------------------------------------------------------------------------- -; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;;; 3. contour field -;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;;; 4. choose output -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; + + +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +zefield = transpose(wt) +zex = localtime +zey = h +set_name = 'HF.ps' +set_title = "Turbulent Heat Flux (K.m.s!U-1!N)" ;"Vertical Eddy Heat Flux " ; +set_titlex = 'Local Time (h)' +set_titley = 'Altitude above surface (km)' +set_subtitle = '' ;'Mean over the simulation domain' +set_xrange = [7.,19.] +set_yrange = [0.,9.] ;; [0.,200.] +set_tickx = 1. +set_ticky = 1. ;; 50. +minval = -1.5 ;-2. +maxval = 2.5 ;2. +nlev = floor(maxval-minval)*10 +pal = 33 +rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;minval = -0.8 +;maxval = 1.2 +;pal = 0 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +if (saveps eq 'true') then PS_Start, FILENAME=set_name +!P.Charsize = 1.2 +;; 0. levels +lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;; 1. background +loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;; 2. color field +loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) + ;;;-------------------------------------------------------------------------------------------------------------------------------- + ;;; WHITE ZONE - 1. get location of interval in the CT - 2. change the CT to have a white zone + ulim=0.09 & dlim=-0.09 & w=where(lev le dlim) & n1=w[n_elements(w)-1] & w=where(lev ge ulim) & n2=w[0] & yy=BYTSCL(lev) & nd=yy[n1] & nu=yy[n2]-5 + nu = nd + (nu-nd)/2 ;; otherwise the interval is too large (because we removed 0) + TVLCT, r, g, b, /Get & r[nd:nu]=255 & g[nd:nu]=255 & b[nd:nu]=255 & TVLCT, r, g, b + ;;;-------------------------------------------------------------------------------------------------------------------------------- + contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;; 3. contour field +loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;; 4. choose output +if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + ;;;restore, filename='tpot_profB' ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -750,6 +751,6 @@ ;if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -;if (n_elements(wmax) eq 0) then stop + +if (n_elements(wmax) eq 0) then stop ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; @@ -822,43 +823,43 @@ if (saveps eq 'true') then PS_End, /PNG ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -; -;restore, filename='addturb.dat' -;velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;zefield = transpose(velmax) -;zex = localtime -;zey = h -;set_name = 'velmax.ps' -;set_title = "Maximum horizontal wind speed (m.s!U-1!N)" -;set_titlex = 'Local Time (h)' -;set_titley = 'Altitude above surface (km)' -;set_subtitle = '' ;'Mean over the simulation domain' -;set_xrange = [11.,17.] ;[8.,17.] -;set_yrange = [0.,4.] -;set_tickx = 1. -;set_ticky = 1. -;minval = 0. -;maxval = 30. ;12. -;nlev = maxval-minval -;pal = 22 -;rrr = 'no' -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -;if (saveps eq 'true') then PS_Start, FILENAME=set_name -;!P.Charsize = 1.2 -;;; 0. levels -;lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] -;;;; 1. background -;loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 -;;;; 2. color field -;loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) -; contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill -;;; 3. contour field -;loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) -;;; 4. choose output -;if (saveps eq 'true') then PS_End, /PNG -;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -; -; + + +restore, filename='addturb.dat' +velmax = SMOOTH(TEMPORARY(velmax), [0,smoothampl], /EDGE_TRUNCATE) +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +zefield = transpose(velmax) +zex = localtime +zey = h +set_name = 'velmax.ps' +set_title = "Maximum horizontal wind speed (m.s!U-1!N)" +set_titlex = 'Local Time (h)' +set_titley = 'Altitude above surface (km)' +set_subtitle = '' ;'Mean over the simulation domain' +set_xrange = [07.,19.] ;[8.,17.] +set_yrange = [0.,4.] +set_tickx = 1. +set_ticky = 1. +minval = 0. +maxval = 30. ;12. +nlev = maxval-minval +pal = 22 +rrr = 'no' +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +if (saveps eq 'true') then PS_Start, FILENAME=set_name +!P.Charsize = 1.2 +;; 0. levels +lev = minval + (maxval-minval)*findgen(nlev+1)/float(nlev) & if (minval ne 0.) then lev = lev[where(lev ne 0.)] +;;; 1. background +loadct, 0 & contour, /NODATA, zefield, zex, zey, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0 +;;; 2. color field +loadct, pal & if (rrr eq 'yes') then TVLCT, r, g, b, /Get & if (rrr eq 'yes') then TVLCT, Reverse(r), Reverse(g), Reverse(b) + contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /overplot, /cell_fill +;; 3. contour field +loadct, 0 & contour, zefield, zex, zey, levels=lev, c_labels=findgen(n_elements(lev))*0.+1., /noerase, xtitle=set_titlex, xrange=set_xrange, xtickinterval=set_tickx, ytitle=set_titley, yrange=set_yrange, ytickinterval=set_ticky, title=set_title, subtitle=set_subtitle, color=0, C_LINESTYLE = (lev LT 0.0) +;; 4. choose output +if (saveps eq 'true') then PS_End, /PNG +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; + + ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;zefield = localtime