大家帮忙看看这个程序哪里出错了

╯①路霺笑

program ldwd
!此程序求露点温度
implicit none
!rh为相对湿度，t为绝对温度，td为露点绝对温度，p为气压值。
real,parameter::a=17.2693882,b=35.86
integer i,j,irec
integer,parameter::nx=4,ny=3
real rh(nx,ny),t(nx,ny),es(nx,ny),e(nx,ny),td(nx,ny)
open (1,file='e:/rh.grd',form='binary')
open (2,file='e:/tmp.grd',form='binary')
print*,'ok1'
do i=1,nx
do j=1,ny
read(1,rec=irec)(rh(i,j))
       irec=irec+1     
           continue

read(2,rec=irec)(td(i,j))
       irec=irec+1     
           continue
es(i,j)=6.1078*exp(a*(t(i,j)-273.16)/(t(i,j)-b))
e(i,j)=es(i,j)*rh(i,j)/100
do while(e(i,j)<es(i,j))
t(i,j)=t(i,j)-0.05
es(i,j)=6.1078*exp(a*(t(i,j)-273.16)/(t(i,j)-b))
end do
if(e(i,j)>=es(i,j))then
td(i,j)=t(i,j)
end if
end do
end do

do j=1,ny
do i=1,nx
print*,"Td=",td(i,j)
end do
end do
open (3,file='e:/td.grd',form='binary')
do j=1,ny
do i=1,nx
write(3) td(i,j)
end do
end do
end


我是通过相对湿度和温度的grd文件来求取露点温度，grd文件资料时间和空间只有一层，
运行没错，但最后显示的是


请大家帮我看看哪里出错了，谢谢

言深深

irec没有赋初值，默认为0，溢出了···
在程序前面添加irec=1看看

╯①路霺笑

言深深 发表于 2014-5-5 21:05
irec没有赋初值，默认为0，溢出了···
在程序前面添加irec=1看看

最后做出来所有td都是0……

言深深

╯①路霺笑 发表于 2014-5-5 21:12
最后做出来所有td都是0……

这就没法保证了···自己逐步检查吧···

╯①路霺笑

言深深 发表于 2014-5-5 21:45
这就没法保证了···自己逐步检查吧···

恩，谢谢，我已经弄出来了

quentin011099

请问下LZ中间的do while 语句里-0.05的判断依据是什么啊，还有这样露点温度算出来是对的么？谢谢啦

xmc21

quentin011099 发表于 2014-9-19 00:41
请问下LZ中间的do while 语句里-0.05的判断依据是什么啊，还有这样露点温度算出来是对的么？谢谢啦

应该是逐步逼近的思路，t以0.05度为间隔减少，直到水汽压e达到饱和水汽压es

quentin011099

xmc21 发表于 2014-10-21 14:54
应该是逐步逼近的思路，t以0.05度为间隔减少，直到水汽压e达到饱和水汽压es

嗯嗯，后面理解这个程序意思了。谢谢啦

wanghao

十分感谢
十分感谢
十分感谢
十分感谢
十分感谢

micolian

请问下，在GRIB2数据中，如何计算露点温度呢