;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;; W 
colors        = 32
minfield_init = -6.
maxfield_init =  6. 
format        = '(F4.1)'
pal           = 4
title_user    = 'Vertical velocity (m s!U-1!N)'
;;;;;;;; tk
minfield_init =  110.
maxfield_init =  210.
format        = '(I0)'
pal	       = 33
title_user    = 'Temperature (K)'
;;;;;;; tk - tsat
minfield_init =  -4. ;-8.
maxfield_init =  4. ;+8.
ndiv	      =  8	
format        = '(F4.1)'
pal           =  0
title_user    = 'T!Datm!N - T!Dsat!N (K)'
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
missing_value=1.e30
lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max
lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
flag_cb = 'true'
windowx = [5.,55.]
windowy = [0.,100.]
windowy = [75.,85.]
windowx = [5., 35.] 
intervalx = 5.
intervaly = 5.
windowy = [0.,130.]
windowx = [10., 110.] 
intervalx = 10.
intervaly = 10.
title_axis = ['horizontal coordinate (x5km)','altitude (km)']
poscb=0.55
;isotropic='false'
lev = 10.*findgen(100)
lev = [-6.,-4.,-2.,2.,4.,6.]
lev = [-8.,-5.,-1.,1.,5.,8.]
lev = [-6.,-5.,-4.,-3.,-2.,-1.,0.]
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;