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

*calc div at every grid point 
*divergence function is hdivg(u,v)
*from grads guide: "The numerical stability of calc horizontal divergence is very low.
*Please use the function with caution"
define div = hdivg(u,v)

*display map of divergence at sigma=0.245
set grads off
set lev 0.245 
d div
draw title Divergence at sigma=0.245 - CNTROLPI DJF 
print
disable print