<div dir="ltr"><div>Dear Sir,<br></div> Thanks. It worked. I also searched similar issues in grads archives and found one solution with if else statement by Serge Claudio and it worked well for me. But Still I could not include the second if condition which is if(W>=6&W<=26);'define c=(0.6+0.07*W)/1000';. generally this will be included in the else statement except only during certain conditions such as cyclones. But Still any way to include this too in the script..?. The following is the currently working script with the output figure attached. <br>
<div class="gmail_extra"><br>'reinit'
<br>'sdfopen <a href="http://uwnd.2009.nc">uwnd.2009.nc</a>'
<br>'sdfopen <a href="http://vwnd.2009.nc">vwnd.2009.nc</a>'
<br>'set lev 1000'<br>'set lat -10 45'<br>'set lon 30 110'
<br>'define W=mag(uwnd,vwnd.2)'<br>if (W>= 3&W<= 6)<br>'define cd=(0.29+3.1/W+7.7/W**2)/1000'<br>else<br>'define cd=(0.6+0.07*W)/1000'<br>endif<br>'define taux=1.3*cd*W*uwnd'
<br>'define tauy=1.3*cd*W*vwnd.2'
<br>'d cd'
<br>'printim cdrag2.jpg white'<br><br><br>-- <br>Sincerely<br>***********************************************<br>Jayakrishnan.P.R<br>Teaching Associate (Atmospheric Science)<br>Academy of Climate Change Education and Research<br>
Kerala Agricultural University<br>Vellanikkara-680 656<br><br>and also<br><br>Research Fellow (Part Time)<br>Department of Atmospheric Sciences<br>CUSAT, Cochin-682 016<br>Kerala, India.<br>Mob: 09496026547<br><br>***********************************************
</div></div>