求助下，FORTRAN程序在此，如何编GRADS的CTL文件？

Zery冬日的激情

subroutine cgw(ua,va,za,rm,f,d,m,n)
dimension ua(m,n),va(m,n),za(m,n),rm(m,n),f(m,n)
          do i=1,m
          do j=1,n
      if(j==1) then
      va(i,j)=rm(i,j)/f(i,j)*((za(i,j+1)-za(i,j))/d)*9.8
      else if(j==n) then
      va(i,j)=rm(i,j)/f(i,j)*((za(i,j)-za(i,j-1))/d)*9.8
      else
      va(i,j)=rm(i,j)/f(i,j)*((za(i,j+1)-za(i,j-1))/(2*d))*9.8
      end if
          end do
          end do

          do i=1,m
          do j=1,n
      if(i==1) then
      ua(i,j)=-1*rm(i,j)/f(i,j)*((za(i+1,j)-za(i,j))/d)*9.8
      else if(i==m) then
      ua(i,j)=-1*rm(i,j)/f(i,j)*((za(i,j)-za(i-1,j))/d)*9.8
      else
      ua(i,j)=-1*rm(i,j)/f(i,j)*((za(i+1,j)-za(i-1,j))/(2*d))*9.8
      end if
          end do
          end do
          return
      end

Zery冬日的激情

    subroutine ssip(a,w,s,m,n,k,l)  dimension a(m,n),w(m,n)
          m1=m-2
          n1=n-2
          do i=3,m1
          do j=3,n1
          w(i,j)=a(i,j)+s*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4.0*a(i,j))/4.0
          end do
          end do
          a(3:m1,3:n1)=w(3:m1,3:n1)
          
          if(1/=1) then
          do i=3,m1
          do j=3,n1
          w(i,j)=a(i,j)-s*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4.0*a(i,j))/4.0
          end do
            end do
      a(3:m1,3:n1)=w(3:m1,3:n1)
      end if
      return
                  end

lqouc

哎呀，好高端的子程序，完全看不懂啊。ctl什么的最难得到了，貌似要拿一血再超神以后才可以得到，大气所的手册上真是一点也没讲啊，马红云和朱禾的书上也找不到。我也是小白，画图不行只好帮顶贴了。

liutaoz21

我会说完全不知道楼主说的什么么？

Zery冬日的激情

liutaoz21 发表于 2013-5-28 19:16

                               
登录/注册后可看大图

我会说完全不知道楼主说的什么么？

就是这个FORTRAN怎么样编个GS文件，我要用GRADS画图

Zery冬日的激情

lqouc 发表于 2013-5-28 18:23

                               
登录/注册后可看大图

哎呀，好高端的子程序，完全看不懂啊。ctl什么的最难得到了，貌似要拿一血再超神以后才可以得到，大气所的手 ...

这么艰难啊，这是个数值预报的子秩序，我们要最后用grads画图，怎么样编个gs文件啊

river

楼主先去看GrADS的基础去吧，好歹你有个雏形我们可以帮忙检查错误，你直接问怎么得到ctl，我只能说去看书吧

天道酬勤

你得有数据，才能编gs文件呀。楼主根据程序算出数据再说吧。你这问题真心不好回答呀。

wygohappy

已知条件不够

善人/jw

LZ你最基础的都没搞清吧
ctl是什么，怎么来的，gs又是什么，怎么来的
写看看手册，了解最基础的吧