Changeset 131 for trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro

Timestamp:

May 25, 2011, 4:19:36 PM (14 years ago)

Author:

aslmd

Message:

mars

M 130 mars/libf/phymars/callradite.F
M 130 mars/README
no more need to modify callradite.F prior to compilation
[but still dimradmars.h must be modified]

--> in callradite.F we now have

-- DEFAULT name_iaer(1) is "dust_conrath"
-- IF (doubleq.AND.active) name_iaer(1) = "dust_doubleq"
-- IF (water.AND.activice) name_iaer(2) = "h2o_ice"

M 130 000-USERS
small meeting about SVN for "planeto" team

M 130 mars/libf/phymars/meso_physiq.F
purely cosmetic changes

M 130 mesoscale/PLOT/MINIMAL/map_latlon.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwave2.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwave_inc.pro
A 0 mesoscale/PLOT/SPEC/MAP/defs/THARSIS_newphys_bugnoRAC_RACgcm_TAUTES.map_uvt_inc.pro
M 130 mesoscale/PLOT/SPEC/MAP/map_uvt.pro
graphic stuff for GW studies and H2O cycle

File:

• trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro (modified) (5 diffs)

trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro

@@ -13 +13 @@
 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 +45 @@
         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 +63 @@
 ;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 +87 @@
      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 +116 @@
         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)