[gradsusr] Coding Problem

[JP Kalb]

Raymond, it is in the script. Although, part of my script got cut off when I was writing the email."open NAM12Ztest.ctl""set display color white""clear""set mpdset hires""set grads off""set lev 925"t=1while (t <= 2)'set t 't'query dim'diminfo=result'run rgbset.gs' 'set lat 36 39' 'set lon -125 -120''set clevs 0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1.0''set ccols 0 4 13 3 10 7 12 8 2''set gxout shaded''define lcl=((TMP2m-DPT2m)*400)''define tc=(TMPprs-273.15)''define es=(6.11*exp((17.625*tc)/(tc+243.05)))''define e=(es*RHprs/100)''define wr=(622*e/(925-e))''define ht=(HGTprs*3.28)''define expo=(exp((2257*wr)/(1004*TMPprs)))''define numer=(-2.06327*TMPprs/(6516))''define denom=(pow((RHprs/100),0.46))''define dthetae=(expo*numer/denom)''define b=((-6.67*dthetae-0.667)*(-0.1*VVELprs))''define c=(0)''d b'        if (ht <= lcl)        ' d c'        endif        if (VVELprs >= 0)        ' d c'        endif'set grads off''d b'
"printim testing12Zcloud.png png"'clear'endwhile'quit'
                                                                         JP Kalb 

    On Monday, February 29, 2016 3:53 PM, raymond ordinario <raymond1723 at gmail.com> wrote:
 

 Hi JP,
it seems that your script omits increment and endwhile line
t=t+1endwhile 
on the end of your process.

Thanks,
Raymond
On Tue, Mar 1, 2016 at 6:27 AM, JP Kalb <jpkalb1990 at yahoo.com> wrote:

Hi there,I am new to learning GrADS. I am having trouble with my if loops for my code."open NAM12Ztest.ctl""set display color white""clear""set mpdset hires""set grads off""set lev 925"t=1while (t <= 2)'set t 't'query dim'diminfo=result'run rgbset.gs' 'set lat 36 39' 'set lon -125 -120''set clevs 0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1.0''set ccols 0 4 13 3 10 7 12 8 2''set gxout shaded''define lcl=((TMP2m-DPT2m)*400)''define tc=(TMPprs-273.15)''define es=(6.11*exp((17.625*tc)/(tc+243.05)))''define e=(es*RHprs/100)''define wr=(622*e/(925-e))''define ht=(HGTprs*3.28)''define expo=(exp((2257*wr)/(1004*TMPprs)))''define numer=(-2.06327*TMPprs/(6516))''define denom=(pow((RHprs/100),0.46))''define dthetae=(expo*numer/denom)''define b=((-6.67*dthetae-0.667)*(-0.1*VVELprs))''define c=(0)''d b'        if (ht <= lcl)        ' d c'        endif        if (VVELprs >= 0)        ' d c'        endif'set grads off'
Any help would be appreciated.Thanks, JP Kalb
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr





-- 
-------------------------------------------------------------------------------------------------------------------------------Raymond C. Ordinario
MS Atmospheric Science
Ateneo De Manila University
Quezon City, Philippines

Weather Specialist I
Technique Application and Meteorological Satellite Section
Weather DIvision
Philippine Atmospheric, Geophysical and Astronomical Services Admistration
Quezon City, Philippines

  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160229/63822c59/attachment-0001.html