clear
reset graphics

set parea off

*open f:/muam/Muam_Mar_eQBO/0/resEPa.ctl
*open f:/muam/Muam_Mar_eQBO/1/resEPa.ctl
open f:/muam/Muam_Mar_eQBO/resid/resEPa_wqbo_ave.ctl
open f:/muam/Muam_Mar_eQBO/resid/resEPa_eqbo_ave.ctl
open f:/muam/Muam_Mar_eQBO/resid/resEPa_disp_ave.ctl

set lat 0 90
set lev 0 100
set lon 0
set t 1
set font 1

set parea 1 6.5 1 6.5

set gxout shaded
set frame on
set annot 1 15
set strsiz 0.2

set vpage 0.5 8. 1. 7.65
*set parea 0.75 6.25 1 6.5

set grid off
set grads off

set rgb 90 180 220 255
set rgb 91 150 200 245
set rgb 92 150 170 200
set rgb 93 110 175 100
set rgb 94 70 175 80
set rgb 95 205 100 0
set rgb 96 195 90 0
set rgb 97 170 30 155
set rgb 98 160 0 160
set rgb 99 150 0 0
set rgb 75 0 0 0

set rbcols 90 91 92 93 94 95 96 97 98 99
*set clevs -200 -150 -50 0 100 150 250 300 400 
set rgb 75 0 0 0

set clevs -8 -6 -4 -2 0 2 4 7 10 


set digsize 0.29 
set xlopts 1 1 0.19
set ylopts 1 1 0.19

set rbrange -8 10
set cint 2
*set rbrange -0.8 2.4
*set cint 0.4

d smth9(gradF.1-gradF.2)  / 100
*d smth9(U.1-U.2)*100
*d smth9(U.1)*1
cbarn 0.8 1 6.8 3.7


set strmden 3
*set arrowhead 2.22
set cthick 15

set gxout contour
set clevs 0
set cstyle 1
set ccols 1
set cthick 20
set clab off
set cterp on

set annot 1 10
draw ylab Altitude, km
draw xlab Latitude, deg

d smth9(gradF.1-gradF.2) 

set cstyle 1
*set clevs 1.996
set clevs 1.996
*set clevs 1
set ccols 75
set cthick 20
set clopts 1 9 0.39
set clab off
set cterp on
*d smth9(smth9(abs(Wres.1-Wres.2)))*60/Wres.3/1
*d smth9(smth9(abs(U.1-U.2)))*960/D.3/1

enable print f:/muam/Muam_Mar_eQBO/resid/Fdiv_w-e_mar_ave.prn
print
disable print

close 3
close 2
close 1