micaps 地面观测资料画6小时降水

kyo1981

反复学习了。。。。。。。。。。。。。

i_gotta_go

已收藏！

流水2012

厉害哦，赞一个

梦之蓝

好高大上呀

chenhongfengye

这里的rain.map从哪里来的呀？

311wen

您好，想请教一下你参考这个帖子试着写了16个时次的降水
提示如下：想请教一下原因


program p
parameter(nt=16)
integer::n,nlev,nflag
real,allocatable::lon(:),lat(:),pp(:)     !动态数组存储经纬度和p
real:: temp(9)
character*8,allocatable::sta(:)         !动态数组存储站号，8个字符
character*12 filename(nt)        !用于储存16个时次的文件名

open(11,file='filename.txt')       ！
do i=1,nt
        read(11,*) filename(i)
        print*,filename(i)
enddo
close(11)  

open(12,file='p.grd',form='binary')        !储存248个时次的能见度场
tim=0.0;nlev=1;nflag=1
do k=1,nt
        open(11,file=filename(k))
        read(11,*)
        read(11,*) n,n,n,n,n        !将该时次的站点数赋值于n
        allocate(lat(n))
        allocate(lon(n))
        allocate(pp(n))
        allocate(sta(n))
        close(11)
        open(11,file=filename(k))   ！获知n后，重新读取。
        read(11,*)
        read(11,*)
        do i=1,n        !只读降水资料
                read(11,*) sta(i),lon(i),lat(i),(temp(j),j=1,9),pp(i)
                write(12)  sta(i),lat(i),lon(i),tim,nlev,nflag,pp(i)
        enddo
    nlev=0
    write(12) sta(n),lat(n),lon(n),tim,nlev,nflag         !一个时次的能见度场输入完毕
        deallocate(lat)
        deallocate(lon)
        deallocate(vi)
        deallocate(sta)
        close(11)
enddo
close(12)
end

Yezi

有用！学习了。。

nunu18

好东西谢谢楼主分享

moji珏丫头

mark......

holton

太谢谢了啊