Index: trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro
===================================================================
--- trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro	(revision 207)
+++ trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro	(revision 324)
@@ -58,5 +58,6 @@
 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,4 +151,5 @@
      alt = [40.,80.]
      alt = [0.,120.]
+     alt = [10.,110.]
      ;;;;;;;;;;;
      ;;;;;;;;;;;
@@ -156,4 +158,6 @@
      minfield_init = 85.  ;100.
      maxfield_init = 185. ;250.
+     minfield_init = 90.
+     maxfield_init = 230.
      title_axis = ['Temperature (K)','Altitude (km)']
      overplot = mcd_tpot
@@ -161,4 +165,5 @@
      SPAWN, "echo 'intervaly=10.' >> param_plot.idl"
      SPAWN, "echo 'intervalx=5.' >> param_plot.idl"
+     SPAWN, "echo 'intervalx=10.' >> param_plot.idl"
      ;;;;;;;;;;;
 
@@ -216,6 +221,25 @@
 	;;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