pro gravitwave2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
file='../EN_COURS/gw/ok_hires_wrfout_d01_9999-09-09_09:00:00_z'
file='../EN_COURS/gw/_wrfout_d01_9999-09-09_09:00:00_z'
file='../LMD_MM_MARS/TESTGW/_wrfout_d01_9999-09-09_09:00:00_z'
;file='../LMD_MM_MARS/TESTGW/wind15_wrfout_d01_9999-09-09_09:00:00_z'
file='../TESTGW/wind20c_wrfout_d01_9999-09-09_09:00:00_z'
charvar='W' & charvarc='W'
;charvar='tk' & charvarc='W'
charvar='tk' & charvarc='W' & cond=1
charvar='tk' & charvarc='tk' & cond=1
;charvar='tpot' & charvarc='tpot'
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;
;
;
what_I_plot=0. & overcontour=0.
SPAWN, '\rm param_plot.idl ; cp gravitwave_inc.pro param_plot.idl'
;
;
;
getcdf, $
        file=file, $
        charvar=charvar, $
        invar=invar
getcdf, $
        file=file, $
        charvar=charvarc, $
        invar=invarc
getcdf, $
        file=file, $
        charvar='vert', $
        invar=vert
getcdf, $
        file=file, $
        charvar='PTOT', $
        invar=columnp
;
;
;
for nt=1,15 do begin
zefile='gravitwave_'+charvar+'_'+charvarc+'_'+string(nt+100,'(I0)')
if (n_elements(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
;
;
;
what_I_plot = reform(invar(*,1,*,nt))
overcontour = reform(invarc(*,1,*,nt))
   if (n_elements(cond) ne 0) then begin
     column = columnp
     ;column = reform(columnp(*,1,*,1))
     @tempcond.inc
     yeye = reform(overplot(*,1,*,nt))
     w = where(yeye le 0.) & yeye[w] = 0. 
     if (charvar eq 'tk') then what_I_plot = what_I_plot - yeye 
     if (charvarc eq 'tk') then overcontour = overcontour - yeye
   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
;
;
;
map_latlon, $
        what_I_plot, $                          ; 2D field
        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)
        overcontour=overcontour, $              ; another 2D field to overplot with contour lines (=0: no)
;        overvector_x=overvector_x, $            ; wind vector - x component (=0: no)
;        overvector_y=overvector_y, $            ; wind vector - y component (=0: no)
;        ct=pal, $                               ; color table (33-rainbow is default)
;        colors=colors, $                        ; number of colors/levels (32 is default)
;        title=title_user, $                     ; title of the plot ('' is default)
        format=format                           ; format of colorbar annotations ('(F6.2)' is default)
;
;
;
PS_End, /PNG
endfor
end