Climatology of Seasonal mean using Daily data in GrADS

Prompt 'Enter the folder name: '
pull fld
'!mkdir ../'fld
ans=y
while(ans=y)

prompt 'Enter the starting year:'
pull y1
prompt 'Enter the ending year:'
pull yn
'reinit'
prompt ' enter the starting month as number: '
pull m1
n=m1
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmm1=mmm
prompt 'Enter the starting day: '
pull d1
prompt ' enter the ending month as number: '
pull mn
n=mn
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmmn=mmm
prompt 'Enter the ending day: '
pull dn
prompt 'Enter the Key word (Active/Break) : '
pull key


'open sst6098.ctl'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1''y1' 'dn''mmmn''yn''
'define t=ave(sst,time='d1''mmm1',time='dn''mmmn')'
'set time 'dn''mmmn''yn''
'define ct=ave(t,time='d1''mmm1''y1',time='dn''mmmn''yn')'
'undefine t'
'set grads off'
'set gxout shaded'
'set clevs 298.5 299 299.5 300 300.5 301 301.5 302 302.5 303 303.5 304'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'd ct'
'cbarn'
'set gxout contour'
*'set ccolor 0'
*'d ct'
'undefine ct'
'close 1'
'draw title sst climatology 'key''
'printim ../'fld'/'key'_'sst'.jpg white'

'clear'
'sdfopen slp.nc'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1''y1' 'dn''mmmn''yn''
'define p=ave(slp,time='d1''mmm1',time='dn''mmmn')'
'set time 'dn''mmmn''yn''
'define cp=ave(p,time='d1''mmm1''y1',time='dn''mmmn''yn')'
'undefine p'
'set grads off'
'set gxout shaded'
*'set clevs 1006.0 1007.0 1008.0 1009.0 1010.0 1011.0 1012.0 1013.0 1014.0 1015.0 1016.0'
*'set ccols 9 14 4 11 5 13 3 10 7 12 8 2'
'd cp'
'cbarn'
'set gxout contour'
'set ccolor 0'
'd cp'
'undefine cp'
'close 1'

'draw title slp climatology 'key''
'printim ../'fld'/'key'_'slp'.jpg white'

'clear'
'open u850.ctl'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1''y1' 'dn''mmmn''yn''
'define uw=ave(u,time='d1''mmm1',time='dn''mmmn')'
'set time 'dn''mmmn''yn''
'define cu=ave(uw,time='d1''mmm1''y1',time='dn''mmmn''yn')'
'undefine uw'
'set grads off'
'set gxout shaded'
'd cu'
'set gxout contour'
'd cu'
'cbarn'
'draw title u wind climatology 'key''
'printim ../'fld'/'key'_'uwind'.jpg white'
'undefine cu'
'close 1'

prompt 'Do you have another one: '
pull ans
endwhile

Code for Climatology of Selected time Average using daily data in GrADS

Prompt 'Enter the folder name: '
pull fld
'!mkdir ../'fld
num=1
'reinit'
'sdfopen slp.nc'
'open sst6098.ctl'
'open u850.ctl'
'set time 31dec1960'
'set lat -25 25'
'set lon 50 180'
'define sp=0'
'define st=0'
'define su=0'
'close 3'
'close 2'
'close 1'
prompt 'Enter the Key word (Active/Break) : '
pull key

prompt 'Enter the year:'
pull yyyy

while(yyyy>1950)
prompt ' enter the starting month as number: '
pull m1
n=m1
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmm1=mmm
prompt 'Enter the starting day: '
pull d1
prompt ' enter the ending month as number: '
pull mn
n=mn
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmmn=mmm
prompt 'Enter the ending day: '
pull dn

'clear'
'sdfopen slp.nc'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1'1960'
'define p=ave(slp,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'set time 31dec1960'
'define sp=sp+p'
'undefine p'
'define cp=sp/'num''
'set grads off'
'set gxout shaded'
*'set clevs 1006.0 1007.0 1008.0 1009.0 1010.0 1011.0 1012.0 1013.0 1014.0 1015.0 1016.0'
*'set ccols 9 14 4 11 5 13 3 10 7 12 8 2'
'd cp'
'cbarn'
'set gxout contour'
'd cp'
'undefine cp'
'draw title slp climatology 'key''
'printim ../'fld'/'key'_'slp'.jpg white'
'close 1'

'clear'
if (yyyy<=1998)
'open sst6098.ctl'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1''1960''
'define t=ave(sst,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'define st=st+t'
'undefine t'
'define ct=st/'num''
'set grads off'
'set gxout shaded'
'set clevs 298 298.5 299 299.5 300 300.5 301 301.5 302 302.5 303 303.5 304'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'd ct'
'cbarn'
'set gxout contour'
'd ct'
'draw title sst climatology 'key''
'printim ../'fld'/'key'_'sst'.jpg white'
'undefine ct'
'close 1'
endif

'clear'
'open u850.ctl'
'set lat -25 25'
'set lon 50 180'
*'set mproj scaled'
'set time 'd1''mmm1''1960''
'define uwnd=ave(u,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'define su=su+uwnd'
'undefine uwnd'
'define cu=su/'num''
'set grads off'
'set gxout shaded'
'd cu'
'set gxout contour'
'd cu'
'draw title u wind climatology 'key''
'printim ../'fld'/'key'_'uwind'.jpg white'
'undefine cu'
'close 1'

prompt 'Enter the next year:'
pull yyyy

num=num+1
endwhile

Spatial distribution

Prompt 'Enter the folder name: '
pull fld
'!mkdir ../'fld

prompt 'Enter the starting year:'
pull y1
prompt 'Enter the starting year:'
pull yn
yyyy=y1
While (yyyy<=yn)

ans='yes'
while(ans='yes')
'reinit'
prompt 'enter the onset month as number'
pull m1
n=m1
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmm1=mmm
prompt mmm1' '
prompt 'Enter the onset day: '
pull d1
prompt 'enter the withdrowal month as number'
pull mn
n=mn
if(n=1)
mmm='jan'
endif
if(n=2)
mmm='feb'
endif
if(n=3)
mmm='mar'
endif
if(n=4)
mmm='apr'
endif
if(n=5)
mmm='may'
endif
if(n=6)
mmm='jun'
endif
if(n=7)
mmm='jul'
endif
if(n=8)
mmm='aug'
endif
if(n=9)
mmm='sep'
endif
if(n=10)
mmm='oct'
endif
if(n=11)
mmm='nov'
endif
if(n=12)
mmm='dec'
endif
mmmn=mmm
prompt mmmn' '
prompt 'Enter the withdraw day: '
pull dn
'open rf_5103.ctl'
'set lat 5 25'
'set lon 60 90'
'set time 'dn''mmm''yyyy''
'define smr=ave(rf,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'set grads off'
'set gxout shaded'
'set clev 10 20 30 40 50 60 70 80 90 100'
'd r'
'undefine r'
'close 1'

if (yyyy<=1998)
'open sst6098.ctl'
'set lat 5 25'
'set lon 60 90'
'set time 'dd''mmm''yyyy''
'define t=ave(sst,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'set grads off'
'set gxout shaded'
'set clev 298 298.5 299 299.5 300 300.5 301 301.5 302'
'd t'
'undefine t'
'close 1'
endif

'sdfopen slp.nc'
'set lat 5 25'
'set lon 60 90'
'set time 'dd''mmm''yyyy''
'define p=ave(slp,time='d1''mmm1''yyyy',time='dn''mmmn''yyyy')'
'set grads off'
'set gxout contour'
'set ccolor 15'
'd p'
'undefine p'
'close 1'

'open u850.ctl'
'set lat 5 25'
'set lon 60 90'
'set time 'dd''mmm''yyyy''
'define uwnd=ave(u,t-2,t+2)'
'set grads off'
'set gxout contour'
'set ccolor 2'
*'d uwnd'
'undefine uwnd'
'close 1'

'draw title contour 'yyyy' 'mmm' 'dd' '
'printim ../'fld'/contour'yyyy''n''dd'.jpg white'


prompt 'Do you want to plot another pulse: '
pull ans
endwhile
yyyy=yyyy+1
endwhile

GrADS code to plot Rainfall slp and zonal wind of Northeast monsoon over southern peninsular India

Prompt 'Enter the folder name: '
pull fld
'!mkdir ../'fld

yyyy=1951
while(yyyy<=2003)

'reinit'
'set grads off'

'set vpage 0 8 0 11 '
'set parea .5 8 8.5 10'
'set mproj scaled'
'sdfopen slp.nc'
'set lat 10'
'set lon 75'
'set time 01Oct'yyyy' 31Dec'yyyy
'define ap=aave(slp,lon=65,lon=75,lat=5,lat=15)'
'define pentp=ave(ap,t-2,t+2)'
'define meanp=ave(ap,time=01Oct'yyyy',time=31Dec'yyyy')'
'd pentp'
'd meanp'
'draw title slp lon=65,lon=75,lat=5,lat=15 'yyyy' '
'close 1'
'undefine ap'
'undefine pentp'
'undefine meanp'

'set vpage 0 8 0 11 '
'set parea .5 8 6.5 8'
'set mproj scaled'
'open u850.ctl'
'set lat 10'
'set lon 75'
'set time 01Oct'yyyy' 31Dec'yyyy
'define au=aave(u,lon=70,lon=80,lat=13,lat=15)'
'define pentu=ave(au,t-2,t+2)'
'define meanu=ave(au,time=01Oct'yyyy',time=31Dec'yyyy')'
'd pentu'
'd meanu'
'draw title u850 lon=70,lon=80,lat=13,lat=15 'yyyy' '
'close 1'
'undefine au'
'undefine pentu'
'undefine meanu'


'set vpage 0 8 0 11'
'set parea .5 8 4.5 6'
'set mproj scaled'
'open rf_5103.ctl'
'set lat 10'
'set lon 75'
'set time 01Oct'yyyy' 31Dec'yyyy
'define r=aave(rf,lon=70,lon=80,lat=8,lat=13)'
'define pentr=ave(r,t-2,t+2)'
'define meanr=ave(r,time=01Oct'yyyy',time=31Dec'yyyy')'
'd pentr'
'd meanr'
'draw title Rain lon=70,lon=80,lat=8,lat=13 'yyyy' '
'close 1'
'undefine r'
'undefine pentr'
'undefine meanr'

'set vpage 0 8 0 11 '
'set parea .5 8 2.5 4'
'set mproj scaled'
'open u850.ctl'
'set lat 10'
'set lon 75'
'set time 01Oct'yyyy' 31Dec'yyyy
'define au=aave(u,lon=70,lon=80,lat=8,lat=13)'
'define pentu=ave(au,t-2,t+2)'
'define meanu=ave(au,time=01Oct'yyyy',time=31Dec'yyyy')'
'd pentu'
'd meanu'
'draw title u850 lon=70,lon=80,lat=8,lat=13 'yyyy' '
'close 1'
'undefine au'
'undefine pentu'
'undefine meanu'

'set vpage 0 8 0 11 '
'set parea .5 8 0.5 2'
'set mproj scaled'
'sdfopen slp.nc'
'set lat 10'
'set lon 75'
'set time 01Oct'yyyy' 31Dec'yyyy
'define ap=aave(slp,lon=80,lon=90,lat=5,lat=15)'
'define pentp=ave(ap,t-2,t+2)'
'define meanp=ave(ap,time=01Oct'yyyy',time=31Dec'yyyy')'
'd pentp'
'd meanp'
'draw title slp lon=80,lon=90,lat=5,lat=15 'yyyy' '
'close 1'
'undefine ap'
'undefine pentp'
'undefine meanp'

'printim ../'fld'/wind'yyyy'.jpg white'

yyyy=yyyy+1
endwhile