- 积分
- 1662
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
原本的数据格式是nc文件,想提取某个格点925hPa的相对湿度,首先是将nc转化成grd格式,gs如下:
'reinit'
'sdfopen E:\keti\wumai\ncep\Pressurelevel\rhum.2016.nc'
'set gxout fwrite'
'set fwrite E:\keti\wumai\grd\925rh_2016.grd'
it=1
while(it<=728)
'set lon 122.5'
'set lat 40'
'set lev 925'
'set t 'it''
'd rhum'
it=it+1
endwhile
'disable fwrite'
;
由编写了ctl文件,并画图对比,是两个一样的
然后又将grd格式转化成txt,fortran程序如下
PROGRAM EX01 !grd文件转化成txt文件
integer,parameter::year=3,num=1460
integer n(year),rh(1,1,1,1460),iy,it
character name(year)*16
data n/1460,1460,728/
data name/"925rh_2014.grd","925rh_2015.grd", "925rh_2016.grd"/
open(2,file='E:\keti\wumai\rh.txt')
do iy=1,year
open(1,file='E:\keti\wumai\grd\'//trim(name(iy)),form='binary')
! open(1,file='E:\keti\wumai\grd\925rh_2016.grd',form='binary')
do it=1,n(iy)
read(1) rh(1,1,1,it)
enddo
close(1)
do it=1,n(iy)
write(2,*)rh(1,1,1,it)
enddo
enddo
close(2)
print*,rh(1,1,1,1)
end
但是输出来的相对湿度明显数据不对,如下:
1099956224
1108869120
1098907648
1104150528
1108869120
1112801284
1109131264
1115160580
1109655552
1101529088
1100480512
。。。。。相对湿度不应该是100以内的数吗,请教大神们帮忙!!!
|
|