Question

我在GrADS中遇到此代码时遇到问题。错误一直说我的两个if语句无效。我恢复了一些明星，但它给了我一个全新的错误。

 "open NAM12Z0709.ctl"
 "set display color white"
 "clear"
 "set mpdset hires"
 "set grad off"
 "set lev 1000" 
 'run rgbset.gs'
 'set lat 36 38'
 'set lon -125 -120'
  'define lcl=((TMPsfc-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 expo=(exp((2257*wr)/(1004*TMPprs)))'
 'define numer=(-2.06327*TMPprs/(6516))'
 'define denom=(pow((RHprs/100),0.46))'
 'define thetae=(expo*numer/denom)'
 'define b=((-6.67*thetae-0.667)*(-0.1*VVELprs))'
  if (HGTprs <= lcl)
    'b=0'
  endif
  if (VVELprs >= 0)
    'b=0'
  endif
  'set clevs 0.0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1.0'
  'set ccols 1 4 13 3 10 7 12 8 2'
  'set gxout shaded'
  'set grads off'
  'd b'

  "printim testing12Zcloud.png png"
  'clear'
  'quit'

Answer 1

你在Grads的变量/字段和Grads Script的变量之间感到困惑。具体而言，您不能将HGTprs和VVELprs用作Grads脚本变量。

GrADS代码错误

1 个答案: