micaps多时次单变量批处理及其后grads画图

平流层的萝卜

Joucly 发表于 2013-5-18 10:30

                               
登录/注册后可看大图

楼主·· 我运行fortran 出现这个问题是怎么回事啊 ？
谢谢了

你的sample读的有问题，把你的程序段和sample贴出来看看吧

Joucly

fortran程序就是用的你那个···我只是想试试··以后好用 结果就出错了
只是把parameter(nt=32)变了

这个是CMD命令处理后的

平流层的萝卜

Joucly 发表于 2013-5-18 11:19

                               
登录/注册后可看大图

fortran程序就是用的你那个···我只是想试试··以后好用 结果就出错了
只是把parameter(nt=32)变了

...

CMD命令用的有点问题，到该目录下，dir /b *.000>filename.txt

平流层的萝卜

Joucly 发表于 2013-5-18 11:19

                               
登录/注册后可看大图

fortran程序就是用的你那个···我只是想试试··以后好用 结果就出错了
只是把parameter(nt=32)变了

...

你的filename应该也没正确生成吧？

Joucly

改了··
flename 是这样

Joucly

运行出来还是这样

平流层的萝卜

Joucly 发表于 2013-5-18 15:14

                               
登录/注册后可看大图

运行出来还是这样

看看你读的程序

Joucly

parameter(nt=32)
integer::n,num
real,allocatable::lon(:),lat(:),vi(:)        !动态数组分别用来储存站点的经纬度和能见度值
real:: temp(14)
character*8,allocatable::sta(:)                !动态数组，储存站号，站号必须是8个字符，5个的话会stnmap出不出来
character*8 ss(112) !112为华北112个站点，用于储存站点号         
character*12 filename(nt)        !用于储存248个时次的文件名
open(11,file='sample.txt')  ！这个sample.txt里边按行储存了华北站点站号。
read(11,*)
do i=1,112
        read(11,*) ss(i)
enddo
close(11)
open(11,file='filename.txt')
do i=1,nt
        read(11,*) filename(i)
        print*,filename(i)
enddo
close(11)               
open(12,file='vis.grd',form='binary')        !储存248个时次的能见度场
do k=1,nt
        print*,'文件数',k
        open(11,file=filename(k))
        read(11,*)
        read(11,*) n,n,n,n,n        ！将该时次的站点数赋值于n
        allocate(lat(n))
        allocate(lon(n))
        allocate(vi(n))
        allocate(sta(n))
        !print*,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,14),vi(i)
                if(vi(i)==9999) then !若为缺测站点，直接pass
                continue
                else        !不是缺测站点，接下来判断是否属于华北地区
                num=0
                do ii=1,112
                        if(sta(i)==ss(ii)) then
                        num=1
                        !print*,k,'   ',sta(i),'    ',ss(ii),'    ',num
                        endif
                enddo
                !print*,vi(i)
               
                if(num==1) then   ！如果该站点属于华北并未缺测，则写入grd中。
                        tim=0.0;nlev=1;nflag=1
                        write(12)  sta(i),lat(i),lon(i),tim,nlev,nflag,vi(i)
                endif
                endif
        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
end

Joucly

这个就是我的程序···  和你一样  只是读取的文件没有那么多

做个霸气的木头

楼主辛苦了~