sdfopen /grove3/iesuw/ccm1/CCM1-21KEXP-SEAE2025.nc
enable print xy.lat.prt
set grid off

*define meridional average over all 5 seasons
*meridional avg calc for lats 35N-55N 
set t 1 5 
define aveu21k = ave(U,lat=35,lat=55)

*draw xy-plot
*set y to 1 because holding y dim constant
set y 1
set grads off

*set time to DJF (first time in file) and draw avg for that month
*set cmark to open circle 
set t 1
set cmark 2
d aveu21k

*set time to JJA (third time in file)
*set cmark to closed circle 
set t 3
set cmark 3
d aveu21k


*set time to ANN (fifth time in file)
*set cmark to closed square
set t 5
set cmark 5 
d aveu21k

*label plot
draw ylab U (m/s)
draw xlab latitude 
draw title Meridional Average (35N-55N) of U (DJF,JJA,ANN) 

*draw key for lines
*first draw mark symbols
*draw mark mark_number x_pos y_pos size
draw mark 2 9 7.5 .1
draw mark 3 9 7 .1
draw mark 5 9 6.5 .1

*draw strings next to marks
*draw string x_pos y_pos
draw string 9.25 7.45 DJF
draw string 9.25 6.95 JJA
draw string 9.25 6.45 ANN


print
disable print