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

平流层的萝卜

landu_cainingze 发表于 2018-12-12 13:52
改了程序的几个地方，出的数据没什么问题，可是还是生成map文件的时候出错
萝卜可以帮忙看一下吗

你的图没发上来

landu_cainingze

landu_cainingze 发表于 2018-12-12 13:52
改了程序的几个地方，出的数据没什么问题，可是还是生成map文件的时候出错
萝卜可以帮忙看一下吗

改了程序的几个地方，出的数据没什么问题，可是还是生成map文件的时候出错
萝卜可以帮忙看一下吗

landu_cainingze

平流层的萝卜 发表于 2018-12-12 13:53
你的图没发上来

第一次发图，捣鼓了半天才传上去，不好意思
这是fortran程序：

parameter(nt=92)

integer::n,num

real,allocatable::lon(:),lat(:),altitude(:),tmax(:) !动态数组分别用来储存站点的经纬度和日最高温度

character*8,allocatable::sta(:) !动态数组，储存站号，站号必须是8个字符，5个的话会stnmap不出来

character*12 filename(nt) !用于储存92个时次的文件名6.1-8.31

open(11,file='filename.txt')

do i=1,nt

read(11,*) filename(i)

!print*,filename(i)

enddo

close(11)

open(12,file='tmax.grd',form='binary') !储存80个时次的最高气温场

do k=1,nt

       open(13,file=filename(k))

       read(13,*)

       read(13,*)

       read(13,*)

       read(13,*)

       read(13,*) n,n !将该时次的站点数赋值于 n

       allocate(lat(n))

       allocate(lon(n))

       allocate(tmax(n))

       allocate(altitude(n))

       allocate(sta(n))

       !print*,n

              do j=1,n !此循环为该时次的站点循环

                     read(13,*)sta(j),lon(j),lat(j),altitude(j),tmax(j)

                     !print*,tmax(j)

                     tim=0.0;nlev=1;nflag=1

                     write(12)sta(j),lat(j),lon(j),tim,nlev,nflag,tmax(j)               

              enddo

       nlev=0

       write(12)sta(n),lat(n),lon(n),tim,nlev,nflag

       !一个时次的高温场输入完毕

       deallocate(lat)

       deallocate(lon)

       deallocate(tmax)

       deallocate(altitude)

       deallocate(sta)

       close(13)

       close(12)

!      pause

enddo

end

landu_cainingze

这是ctl文件
dset  D:\18wang\learnGrads\batchDrawStation\tmax.grd
dtype    station
stnmap   D:\18wang\learnGrads\batchDrawStation\tmax.map
undef    9999
title    daily tmax
tdef 92 linear  02z01Jun2018 1dy
vars 1
tmax  0 99  daily tmax
endvars

平流层的萝卜

landu_cainingze 发表于 2018-12-12 14:03
改了程序的几个地方，出的数据没什么问题，可是还是生成map文件的时候出错
萝卜可以帮忙看一下吗

你在应用stnmap命令的时候有问题吧

平流层的萝卜

landu_cainingze 发表于 2018-12-12 14:07
第一次发图，捣鼓了半天才传上去，不好意思
这是fortran程序：
parameter(nt=92)integer::n,numreal,al ...

把close(12)放在大循环的外边

landu_cainingze

平流层的萝卜 发表于 2018-12-12 19:20
把close(12)放在大循环的外边

太感谢萝卜了！！是循环的问题！！解决了！！！

平流层的萝卜

landu_cainingze 发表于 2018-12-13 12:59
太感谢萝卜了！！是循环的问题！！解决了！！！

解决就好

landu_cainingze

平流层的萝卜 发表于 2018-12-13 19:07
解决就好

{:5_238:}

fengzaojianche

cmd能给个详细过程吗，我这种菜鸟做不出来