- 积分
- 310
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-6
- 最后登录
- 1970-1-1
![[l一媞l] 粉丝数:7 微博数:233 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 一一一一 于 2018-3-15 11:49 编辑
有变量pds的数据,目标是计算每天有多少个满足pds<=-0.04的格点数n,将与时间一一对应的数据n输出
写数据的gs如下:
'reinit'
'open f:\data\pds.ctl'
'set gxout fwrite'
'set fwrite f:\data\test.grd'
'set lev 850'
t2=1
while(t2<=100)
'set t 't2''
m=0
a=105.0
while(a<=125.0)
'set lon 'a''
b=22.0
while(b<=31.0)
'set lat 'b''
if(pds<=-0.04)
m=m+1
endif
b=b+0.5
endwhile
a=a+0.5
endwhile
'define n='m''
'd n'
t2=t2+1
endwhile
'disable fwrite'
*
写出来的ctl如下:
dset f:\data\test.grd
title n
undef -9.99e+08
xdef 1 levels 1
ydef 1 levels 1
zdef 1 levels 1
tdef 100 linear 00z01may1979 1440mn
vars 1
n 0 99 n
endvars
|
|