怎么编写计算逆温层厚度的程序？

颂与梵唱

program main mie
implicit none
integer,parameter::x=85,y=69,z=37,t=124
real tem(x,y,z,t),inv(x,y,t),inve(x,y,t),geo(x,y,z,t),inver(x,y,t)
open(1,file='d:/dachuang/datas/2018.01.grd/erapressurelevel2018.01.t.grd',form='binary')
  do it=1,t
    do iz=1,z
        read(1)((tem(i,j,iz,it),i=1,x),j=1,y)
    enddo
  enddo
close(1)
open(2,file='d:/dachuang/datas/2018.01.grd/erapressurelevel2018.01.z.grd',form='binary')
  do it=1,t
    do iz=1,z
        read(2)((geo(i,j,iz,it),i=1,x),j=1,y)
    enddo
  enddo
close(2)
do j=1,y
    do i=1,x
        do it=1,t
            do iz=1,z-1
                izz=iz+1
                if(tem(i,j,iz,it)<tem(i,j,izz,it))
                do k=izz,37
                    if(tem(i,j,k,it)>=tem(i,j,k+1,it))
                    break
                    endif
                enddo
            break
            enddo
            inv(i,j,it)=tem(i,j,izz,it)-tem(i,j,iz,it)
            inve(i,j,it)=geo(i,j,izz,it)-geo(i,j,iz,it)
            inver(i,j,it)=tem(i,j,iz,it)
        enddo
    enddo
enddo
open(3,file='d:/dachuang/datas/inv.grd',form='binary')
open(4,file='d:/dachuang/datas/inve.grd',form='binary')
open(5,file='d:/dachuang/datas/inver.grd',form='binary')
do it=1,t
    write(3)((inv(i,j,it),i=1,x),j=1,y)
    write(4)((inve(i,j,it),i=1,x),j=i,y)
    write(5)((inver(i,j,it),i=1,x),j=i,y)
enddo
close(3)
close(4)
close(5)
end
我的思路就是利用逆温层上层温度大于下层温度，判断出来哪一层开始产生逆温，然后在这个基础上再找出开始出现上层温度小于下层温度的层次，这时候就是逆温层的顶，然后写了判断，但是出不来grd数据，有没有大神帮忙看依稀

gwk

do j=1,y
    do i=1,x
        do it=1,t
            do iz=1,z-1
                izz=iz+1
                if(tem(i,j,iz,it)<tem(i,j,izz,it))
                do k=izz,37
                    if(tem(i,j,k,it)>=tem(i,j,k+1,it))
                    inv(i,j,it)=tem(i,j,k,it)-tem(i,j,iz,it)
                    inve(i,j,it)=geo(i,j,k,it)-geo(i,j,iz,it)
                    inver(i,j,it)=tem(i,j,iz,it)
                    endif
                enddo                endif
            enddo
        enddo
    enddo
enddo
循环这边修改成这样试一试。

颂与梵唱

gwk 发表于 2018-7-29 09:39
do j=1,y
    do i=1,x
        do it=1,t

我后面把break改成了exit就有用了，我把语法串了
十分感谢

wjkhaida

你好  请问你的资料是用的什么资料?

颂与梵唱

wjkhaida 发表于 2019-3-6 14:41
你好  请问你的资料是用的什么资料?

EC的数据 0.75*0.75

wjkhaida

颂与梵唱 发表于 2019-3-6 20:50
EC的数据 0.75*0.75

谢谢  您能说的具体一点吗？  EC的哪个数据？就是哪个数据集里面有这个垂直温度的数据？

咯叽咯叽

您好啊，想问您几个问题
1.您用的具体是EC的哪个数据呢
2.程序运行后，怎么看厚度呢，以及和经度纬度的关系

颂与梵唱

wjkhaida 发表于 2019-3-7 09:29
谢谢  您能说的具体一点吗？  EC的哪个数据？就是哪个数据集里面有这个垂直温度的数据？

EC-interim

颂与梵唱

咯叽咯叽 发表于 2019-3-31 19:51
您好啊，想问您几个问题
1.您用的具体是EC的哪个数据呢
2.程序运行后，怎么看厚度呢，以及和经度纬度的关 ...

1.EC-interim  
2.逆温层厚度是指位势高度场的厚度以及逆温层上下两层的温度差  
   逆温层厚度是一个数组 和经纬度的关系根据格点数就可以算出来了