sdfopen /grove3/iesuw/ccm1/CCM1-CNTROLPI-MONE2025.nc
enable print vorticity.prt
set grid off
set t 1 
set z 1 12

*calc vort at every grid point 
*vorticity function is hcurl(u,v)
define vort = hcurl(u,v)

*display a vertical cross-section of vorticity at a single longitude
set grads off
set lon 40
set cint .25e-05
d vort
draw title Vert Cross Section of Vorticity (40E) - CNTROLPI DJF 
draw ylab sigma level
draw xlab latitude
print
disable print