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

*calc wind speed at every grid point 
*magnitude function is mag(u,v)
define wind = mag(u,v)

*display x-section of magnitude 
set grads off
set lat 45
set cint 1
d wind 
draw title Vert Cross Section of Wind Speed (45N) - CNTROLPI DJF 
draw ylab sigma level
draw xlab longitude
print
disable print