<div dir="ltr">Hey <div><br></div><div>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 </div><div><br></div><div>
<div>'reinit'</div><div>'open wrfout.ctl'</div><div><br></div><div>'set lat 20 40'</div><div>'set lon 60 80'</div><div>'set lev 1'</div><div>'set gxout shaded'</div><div>'set mpdset hires'</div>
<div>'set t 1 696'</div><div>*============================ aod calculation ==============================</div><div>'do i=1,696'</div><div>xx=0.</div><div>yy=0.</div><div> 'do j=0,29'</div><div><br>
</div><div> 'x=f->TAUAER2(i,j,,)'</div><div> ' y=f->TAUAER3(i,j,,)'</div><div> 'xx=xx+x'</div><div> ' yy=yy+y'</div><div><br></div><div>'AE550nm=-log(xx/yy)/log(400.0/600.0)'</div>
<div>'beta=yy/(600.0^(-1.0*AE550nm))'</div><div>'AOD500nm=beta*(550.0^(-1.0*AE550nm))'</div><div><br></div><div><br></div><div><br></div><div>'print AOD550nm'</div><div><br></div><div><br></div><div>
*=====================================================================</div><div><br></div><div>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</div>
<div><br></div>-- <br><div dir="ltr"><div><span style="background-color:rgb(255,255,255);font-family:arial,sans-serif;font-size:13px"><i>Kind Regards</i></span></div><div><span style="background-color:rgb(255,255,255);font-family:arial,sans-serif;font-size:13px"><i><br>
</i></span></div><div><b style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Muhammad Zeeshaan Shahid</b></div><div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
PhD Research Fellow</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Institute of Atmospheric Physics (IAP)<br>Chinese Academy of Sciences (CAS)<br>P.O.Box 9804, Beijing 100029<br>
P.R. China<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"><br></div></div><div><br></div></div>
</div></div>