Changeset 324 for trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro

Timestamp:

Oct 21, 2011, 11:31:28 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE : Preparatory commit for the ultimate option mars=42 which

would allow mesoscale modeling with photochemistry.

[see r76 method to add 'mars' options]
Modified module_lmd_driver.F and Registry.EM and runmeso
Modified readmeteo.F90 and introduced an option -DPHOTOCHEM
...
Transparent to the casual user
Option mars=42 not yet finished though -- so do not use!

File:

• trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro (modified) (5 diffs)

trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro

@@ -58 +58 @@
 nts=1  & nte=100  & nxs=50   & nxe=50
 nts=1  & nte=100  & nxs=60   & nxe=60
-nts=1  & nte=100  & nxs=20   & nxe=20
+;nts=1  & nte=100  & nxs=20   & nxe=20
+nts=29 & nte=29   & nxs=60   & nxe=60
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;
@@ -150 +151 @@
      alt = [40.,80.]
      alt = [0.,120.]
+     alt = [10.,110.]
      ;;;;;;;;;;;
      ;;;;;;;;;;;
@@ -156 +158 @@
      minfield_init = 85.  ;100.
      maxfield_init = 185. ;250.
+     minfield_init = 90.
+     maxfield_init = 230.
      title_axis = ['Temperature (K)','Altitude (km)']
      overplot = mcd_tpot
@@ -161 +165 @@
      SPAWN, "echo 'intervaly=10.' >> param_plot.idl"
      SPAWN, "echo 'intervalx=5.' >> param_plot.idl"
+     SPAWN, "echo 'intervalx=10.' >> param_plot.idl"
      ;;;;;;;;;;;
 
@@ -216 +221 @@
         ;;oplot, yeye, vert/1000., linestyle=1
 
+
+loadct, 4
+for i=16,37 do begin  ;10-37
+for j=0,40 do begin
+dec = j - 20  
+overoverplot = reform(invar(nx+dec,middle+dec,*,i))
+oplot, overoverplot, column, color=220  ;psym=3
+endfor
+endfor
+
+!p.thick = 4.0
+oplot, overplot, overplot_column, linestyle=2
+!p.thick = 2.0
+
+!p.thick = 4.0
+oplot, reform(invar(nx,middle,*,28)), column
+!p.thick = 2.0
+
 @tempcond.inc
 overplot = reform(overplot(nx,middle,*,nt))
+w = where(overplot le 0.) & overplot[w] = !VALUES.F_NAN
 oplot, overplot, overplot_column, linestyle=1