grads网格赋值问题

wly

我想通过网格的下垫面类型来赋值饱和度，比如下垫面类型值为1的地方，n赋值为0，下垫面类型值为2的地方，n赋值为1，值为3和4的地方，n赋值为0.53，我这样写完，画出n值就是0.53，是不是要写个循环啊，求大神赐教下怎么写？


if(cnp=1)
  'define n=0'
else
if(cnp=2)
  'define n=1'
else
  'define n=0.53'
endif
endif

gwk

把你gs文件都发出来看看，这一小段看不出问题

wly

'open 0418.ctl'
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.5 10.0 2 8.0'

'set xaxis 0 1800 200'
'set yaxis 1 1260 200'

'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'

'set grads off'
'set grid off'
'set mpdraw off'

'set cint 2'
'set csmooth on'

'set font 4'
'set string 1 bc 3 0'
'set strsiz 0.15 0.2'
'draw string 6 1.2 X(m)'
'set string 1 l 3 90'
'set strsiz 0.15 0.2'
'draw string 0.5 4.5 Y(m)'
'set string 1 br 3 0'

'set strsiz 0.1 0.15'
*'draw string 10.2 0.38 W/m`a2`n'
'draw string 10.2 0.6 W/m`a2`n'
'set gxout shaded'

'define c=7.63*(tsfc-273.15)/(241.9+tsfc-273.15)'
'define qq=pow(10,c)'
'define qq=qq*6.11*100.0'
'define qtg=0.622*qq/101900'

if(cnp=1)
  'define n=0'
else
if(cnp=2)
  'define n=1'
else
  'define n=0.53'
endif
endif

'define qa=q(lev=60)'
'define wspd=mag(u(lev=60),v(lev=60))'

'define he=n*1.29*2500000*0.0025*wspd*(qtg-qa)'

'd he'

'cbarn 1 0'

wly

wly 发表于 2018-9-13 19:15
'open 0418.ctl'
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.5 10.0 2 8.0'

n=0的时候，he就为0，但是这样写画出的图n为0的地方he并不为0

gwk

wly 发表于 2018-9-13 19:15
'open 0418.ctl'
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.5 10.0 2 8.0'

cnp是什么？

wly

gwk 发表于 2018-9-13 19:44
cnp是什么？

是地表利用类型

wly

gwk 发表于 2018-9-13 19:44
cnp是什么？

是地表利用类型

gwk

wly 发表于 2018-9-13 20:12
是地表利用类型

它从哪儿来的？数据读入了？

wly

gwk 发表于 2018-9-13 20:55
它从哪儿来的？数据读入了？

是模拟的结果

wly

gwk 发表于 2018-9-13 20:55
它从哪儿来的？数据读入了？

是模拟的结果