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

TIANYU

您好，请问做出单时次的图像后，如何求月平均

冰枫

好东西，赞一个

TIANYU

我在尝试求月平均过程中，发现累加日观测数据时，缺测区域会不断扩大，请问您这该如何解决呢

double001

好帖，{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}{:eb302:}

小萌萌

萝卜大大你好，所用的数据格式和你的比较我看了是一样的，我想从众多的站里面挑出单个的一个站点的要素值，问题是直接用这个fortran程序阔以么？{:5_213:}

平流层的萝卜

小萌萌 发表于 2016-3-7 15:46
萝卜大大你好，所用的数据格式和你的比较我看了是一样的，我想从众多的站里面挑出单个的一个站点的要素值， ...

那应该是不难的，只要加一个判断语句，把站点号判断一下。再对每个文件都这么循环，就可以了。

小萌萌

平流层的萝卜 发表于 2016-3-7 15:56
那应该是不难的，只要加一个判断语句，把站点号判断一下。再对每个文件都这么循环，就可以了。

昂好的，那个我后续再加上，刚刚运行了下，感觉有点错。我把程序和结果放上来大大看看~
代码：
program surface
          parameter(nt=2126)
          integer::n,num,n1,n2,n3,n4
           real,allocatable::lon(:),lat(:),ws(:)    !动态数组分别用来储存站点的经纬度和风速
           real:: temp(9)!?
          character*8,allocatable::sta(:)      !动态数组，储存站号，站号必须是8个字符，5个的话会stnmap出不出来   
                  !character*8 ss(112) !112为华北112个站点，用于储存站点号
          character*12 filename(nt)        !用于储存2126个时次的文件名
          open(1,file='F:\plot\filename.txt') !这个.txt里边按行储存了青岛站点站号
           
                   !read(11,*)
           !do i=1,112
           !read(11,*) ss(i)
           !enddo
           !close(11)

                   do i=1,nt
               read(1,*) filename(i) !?read(11,*) ss(i)
               print*,'Filename:',filename(i)
            end do
            close(1)      
                       
                       
                                
        open(2,file='F:\plot\2015surface\windspeed.grd',form='binary')  !储存2126个时次的风速
        do k=1,nt
        open(3,file=filename(k))
        read(3,*)
        read(3,*) n1,n2,n2,n4,n    !将该时次的站点数赋值于n
        allocate(lat(n)) !分配动态数组
        allocate(lon(n))
        allocate(ws(n))
        allocate(sta(n))
        close(3)
        open(4,file=filename(k))   !获知n后，重新读取数据。
        read(4,*)
        read(4,*)
          do i=1,n        !只读风速，别的要素暂略
             read(4,*) sta(i),lon(i),lat(i),(temp(j),j=1,9),ws(i)  
              tim=0.0;nlev=1;nflag=1
              write(2)  sta(i),lat(i),lon(i),tim,nlev,nflag,ws(i)  !一个时次的降水场输入完毕
          enddo
           nlev=0
            write(2) sta(n),lat(n),lon(n),tim,nlev,nflag
        !一个时次的降水场输入完毕，告诉Grads 该时次的数据结束。这是一个特殊标记。
        deallocate(lat)
        deallocate(lon)
        deallocate(ws)
        deallocate(sta)
!释放动态数组
        close(4)
        enddo
         close(2)
      end

小萌萌

平流层的萝卜 发表于 2016-3-7 15:56
那应该是不难的，只要加一个判断语句，把站点号判断一下。再对每个文件都这么循环，就可以了。

运行结果

平流层的萝卜

小萌萌 发表于 2016-3-7 16:04
运行结果

检查一下你的出错的文件，有可能是文件的数据不全

小萌萌

平流层的萝卜 发表于 2016-3-7 18:17
检查一下你的出错的文件，有可能是文件的数据不全

好的，一共有2126个文件，但是只出来了1/3/的量，不知道哪里出问题了，麻烦大大看一下消息，详细的问题我在那里提问了，感谢。