[gradsusr] How to calculate AOD in GrADS

[Zeeshaan Shahid]

Hey

I am new to GrADS, and I would like to calculate Aerosol Optical Depth at
550nm from WRF-Chem output result, I tried to make a script by myself  as
following

'reinit'
'open wrfout.ctl'

'set lat 20 40'
'set lon 60 80'
'set lev 1'
'set gxout shaded'
'set mpdset hires'
'set t 1 696'
*============================ aod calculation ==============================
'do i=1,696'
xx=0.
yy=0.
 'do j=0,29'

  'x=f->TAUAER2(i,j,,)'
 ' y=f->TAUAER3(i,j,,)'
 'xx=xx+x'
 ' yy=yy+y'

'AE550nm=-log(xx/yy)/log(400.0/600.0)'
'beta=yy/(600.0^(-1.0*AE550nm))'
'AOD500nm=beta*(550.0^(-1.0*AE550nm))'



'print AOD550nm'


*=====================================================================

But its not working, I would like to request to you that please guide me
how i can calculate AOD, I shall be very thankful to you for your kind help

-- 
*Kind Regards*

*Muhammad Zeeshaan Shahid*
PhD Research Fellow
Institute of Atmospheric Physics (IAP)
Chinese Academy of Sciences (CAS)
P.O.Box 9804, Beijing 100029
P.R. China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131129/ad4e9d8d/attachment.html