Index: trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro
===================================================================
--- trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro	(revision 117)
+++ trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro	(revision 131)
@@ -13,14 +13,20 @@
 folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_2D.morepoints.160818/'
 folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544/'
-folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440/'
-folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468/'
+folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints.161542/'
+;folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440/'
+;folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468/'
 file=folder+'wrfout_d01_9999-09-09_09:00:00_z'
 charvar='W' & charvarc='W'
-charvar='tk' & charvarc='W'
+charvar='W' & charvarc='tk' & cond=0
+;charvar='tk' & charvarc='W'
 ;charvar='tk' & charvarc='W' & cond=1
-charvar='tk' & charvarc='tk' & cond=1
+;charvar='tk' & charvarc='tk' & cond=1
+;charvar='tk' & charvarc='tk' & cond=0
 ;;charvar='tpot' & charvarc='tpot'
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+charvar='HR_NLTE' & charvarc='tk' & cond=0
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+	;folder='/home/aslmd/GRAVITWAVE/'
+	;file=folder+'GW_MARS_highwind_3D_loweropacity_morepoints_diff.nc'
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;
 ;
@@ -39,8 +45,12 @@
         charvar=charvarc, $
         invar=invarc
+;getcdf, $
+;        file=file, $
+;        charvar='vert', $
+;        invar=vert
 getcdf, $
-        file=file, $
-        charvar='vert', $
-        invar=vert
+         file='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D.157077/wrfout_d01_9999-09-09_09:00:00_z', $
+         charvar='vert', $
+         invar=vert
 getcdf, $
         file=file, $
@@ -53,23 +63,20 @@
 ;for nt=27,27 do begin
 for nt=1,40 do begin
+;for nt=9,9 do begin
+;for nt=0,20 do begin
 zefile='gravitwave_'+charvar+'_'+charvarc+'_'+string(nt+100,'(I0)')
-if (n_elements(cond) ne 0) then zefile = 'cond_' + zefile
+if (n_elements(cond) ne 0) then if (cond ne 0) then zefile = 'cond_' + zefile
 ;
 ;
 ;
-  PS_Start, filename=zefile+'.ps'
-  print, zefile+'.ps'
-  !P.Charsize = 1.2
-  !p.charthick = 2.0
-  !p.thick = 2.0
-  !x.thick = 2.0
-  !y.thick = 2.0
+  PS_Start, filename=zefile+'.ps' & print, zefile+'.ps'
 ;
 ;
 ;
 s = size(invar) & middle = floor(s[2]/2) & print, 'PLOT at y subs ', middle
-what_I_plot = reform(invar(*,middle,*,nt))
+what_I_plot = reform(invar(*,middle,*,nt)) & print, min(what_I_plot), max(what_I_plot)
 overcontour = reform(invarc(*,middle,*,nt))
    if (n_elements(cond) ne 0) then begin
+   if (cond ne 0) then begin
      column = columnp
      ;column = reform(columnp(*,1,*,1))
@@ -80,16 +87,26 @@
      if (charvarc eq 'tk') then overcontour = overcontour - yeye 
    endif
+   endif
 xx = findgen(n_elements(what_I_plot(*,0)))
 zz = vert / 1000.
-
-	;w = where(what_I_plot lt 0.)
-	w = where(what_I_plot lt -1.)
-	;w = where(what_I_plot lt -3.)
-	if (w(0) ne -1) then begin
-	for i=0,n_elements(w)-1 do begin
-	 ind = ARRAY_INDICES(what_I_plot, w(i))
-	 print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
-	endfor
-	endif
+;
+	;;w = where(what_I_plot lt 2.)
+	;;w = where(what_I_plot lt 0.)
+	;w = where(what_I_plot lt -1.)
+	;;w = where(what_I_plot lt -3.)
+	;if (w(0) ne -1) then begin
+	;for i=0,n_elements(w)-1 do begin
+	; ind = ARRAY_INDICES(what_I_plot, w(i))
+	; print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
+	;endfor
+	;endif
+;
+;
+;
+if (charvar eq 'HR_NLTE') then begin
+	what_I_plot = what_I_plot * 3700. ;; en K/hour
+	;w = where((zz le 60.) and (zz ge 50.)) ;; pour enlever la ligne d'altitude ou le NLTE est active
+        ;if (w(0) ne -1) then what_I_plot(*,w) = 0.
+endif
 ;
 ;
@@ -99,6 +116,6 @@
         xx, $                                   ; 1D latitude
         zz, $                                   ; 1D longitude
-;        minfield=minfield_init, $               ; minimum value of plotted field (=0: calculate)
-;        maxfield=maxfield_init, $               ; maximum value of plotted field (=0: calculate)
+;        minfield_init, $               ; minimum value of plotted field (=0: calculate)
+;        maxfield_init, $               ; maximum value of plotted field (=0: calculate)
         overcontour=overcontour, $              ; another 2D field to overplot with contour lines (=0: no)
 ;        overvector_x=overvector_x, $            ; wind vector - x component (=0: no)
