script available?

[Rodney Teck]

Fank

I've attached a script which I came across sometime back it may not do
exactly as you want but see how you get on.

Rod Teck+


On Wed, 2009-02-11 at 18:05, Frank Colby wrote:
> Hi Everyone,
>
> Does anyone have a script to compute frontogenesis?
>
> Thanks,
>
> Frank Colby
> UMass Lowell

--
Dr.Rodney Teck
ICARUS Climatic Research Unit
St.Ann's
NUI Maynooth
Co Kildare
Ireland

tel: ++353 (1)7083769
email: rodney.teck at nuim.ie
http://icarus.nuim.ie

-------------- next part --------------
'reinit'
'open lipr8710.ctl'
'open lisf8710.ctl'
'set lon 280 360'
'set lat 20 72.5'
'set t 5'
'set map 1'

'set gxout shaded'
'rhmax=max(rh, lev=850, lev=700)'
'n=100-(100-rhmax)*2'
defcol_clouds_l()
'set csmooth linear'
'd n'
'run cbarc 1.8 7.6 0'
'set lev 700'
'set ccolor 16'
'set gxout contour'
'set cstyle 1'
'set cint 20'
'set clab on'
'set clskip 3'
'd z/10'
'set ccolor 2'
'set gxout contour'
'set cstyle 5'
'set cint 20'
'set clab on'
'set clskip 5'
'd 0.1*(z(lev=500)-z(lev=1000))'
'set csmooth linear'
'set gxout contour'
'set lev 850'
'dtx = cdiff(tmp,x) '
'dty = cdiff(tmp,y) '
'dt = sqrt(dtx*dtx+dty*dty)'
'dttx = cdiff(dt,x) '
'dtty = cdiff(dt,y) '
'dtt = sqrt(dttx*dttx+dtty*dtty)'
'rhmax=max(rh, lev=850, lev=500)'
'fr = dt*0.15*rhmax'
'dx = cdiff(lon,x)*3.1416/180 '
'dy = cdiff(lat,y)*3.1416/180 '
'adv = -1*((u*dtx)/(cos(lat*3.1416/180)*dx)+v*dty/dy)/6.37e6'
'zadv = abs(adv)/adv'
defcol_fr_warm()
'd fr*zadv'
defcol_fr_cold()
'd fr*zadv'
pull dummy

'set clab forced'
'set cthick 6'
'set cstyle 3'
'set clskip 3'
'd mslp.2/100'

*disp_pres()
'draw title Frontal Zones(AT-850), 1.5-3 km Clouds[%],H700,H500/1000'
*'printim at112_8.png x800 y500 black'
*'enable print f:/front110105.gmf'
*'print'
*'disable print'

function defcol_clouds_l()
'set rgb 81 10 10 10'
'set rgb 82 20 20 20'
'set rgb 83 30 30 30'
'set rgb 84 40 40 40'
'set rgb 85 50 50 50'
'set rgb 86 60 60 60'
'set rgb 87 70 70 70'
'set rgb 88 80 80 80'
'set rgb 89 90 90 90'
'set rgb 90 100 100 100'
'set rgb 91 110 110 110'
'set rgb 92 120 120 120'
'set rgb 93 130 130 130'
'set rgb 94 140 140 140'
'set rgb 95 150 150 150'
'set rgb 96 160 160 160'
'set rgb 97 170 170 170'
'set rgb 98 180 180 180'
'set rgb 99 190 190 190'
'set clevs 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100'
'set ccols 0 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99'
return

function defcol_fr_warm()
'run define_colors'
'set clevs  32'
'set ccols  28'
'set cthick 10'
'set clab off'
return


function defcol_fr_cold()
'set clevs  -29'
'set ccols  47'
'set clab off'
return