E-P通量剖面的经度问题

酒国

利用http://bbs.06climate.com/forum.php?mod=viewthread&tid=28456中的fortran程序，算的某年E-P通量。

有个问题是。经度值为1个，是全球范围的一个平均吗？全球的和某一区域的肯定是不相同的吧？

如果要算某个范围，比如60-150E，应该如何实现呢？是不是提取数据的时候只提取60-150E这个范围的就可以了？


下面程序来自上面那个链接中的fortran。

do kt=1,nt
          do l=1,nz
          read(1) ((u(i,j,l,kt),i=1,nx),j=1,ny)
          read(2) ((v(i,j,l,kt),i=1,nx),j=1,ny)
          read(3) ((t(i,j,l,kt),i=1,nx),j=1,ny)
      enddo;enddo

      do kt=1,nt
      do l=1,nz
          do j=1,ny
          avev(j,l,kt)=0.0
          aveu(j,l,kt)=0.0
          do i=1,nx
      avev(j,l,kt)=avev(j,l,kt)+v(i,j,l,kt)/nx
      aveu(j,l,kt)=aveu(j,l,kt)+u(i,j,l,kt)/nx
      enddo;enddo;enddo;enddo

      do kt=1,nt
      do l=1,nz
      do j=1,ny
      zz=0.0
      do i=1,nx
      avevu(j,l,kt)=(u(i,j,l,kt)-aveu(j,l,kt))*(v(i,j,l,kt)-avev(j,l,kt))+zz
      zz=avevu(j,l,kt)
      enddo;enddo;enddo;enddo

      do l=1,nz
      do j=1,ny
          ep1(j,l)=0.0
          do kt=1,nt
          ep1(j,l)=ep1(j,l)-a*cos((-90+(j-1)*2.5)/180*3.14)*(avevu(j,l,kt)/nx)/nt
      enddo;enddo;enddo

第一张图为E-P通量，第二张为纬向风距平分布。

weinihou

这代码是全部的纬向平均啊

如果你想搞区域的，你把求数值的那里 i 的循环不要做到 nx

酒国

weinihou 发表于 2018-1-23 15:31
这代码是全部的纬向平均啊

如果你想搞区域的，你把求数值的那里 i 的循环不要做到 nx

谢谢您的解答。是这个意思吗？parameter(nx=17,ny=37,nz=17,nt=55)
我这里nx是对应某一区域。您说的不要做到nx是不是就是这个意思？

weinihou

酒国 发表于 2018-1-23 15:48
谢谢您的解答。是这个意思吗？parameter(nx=17,ny=37,nz=17,nt=55)
我这里nx是对应某一区域。您说的不要 ...

恩，是这个意思

要么是取用数据的时候只取对应区域的，要么就计算的时候只取对应范围的，这样纬向平均就是一定范围的了，其他看样子都不需要变

酒国

weinihou 发表于 2018-1-24 09:43
恩，是这个意思

要么是取用数据的时候只取对应区域的，要么就计算的时候只取对应范围的，这样纬向平均 ...

谢谢。我按照这个思路做了一下。在这个链接中。但E-P通量正负分布与纬向风距平的正负不对应。你看有可能是什么原因呢？
http://bbs.06climate.com/forum.php?mod=viewthread&tid=58637

weinihou

酒国 发表于 2018-1-24 10:34
谢谢。我按照这个思路做了一下。在这个链接中。但E-P通量正负分布与纬向风距平的正负不对应。你看有可能 ...

我具体不做这方面工作，我就简单提点儿自己的理解了，在那边吧