这程序把我折磨疯了，为什么就是不能出现正确的图呢？

TWiligHT

你们看看这图怎么了？

实习内容计算1948-2013年1月500hPa高度场平均图，2008年1月500hPa高度距平图，2008年1月500hPa高度纬偏图，程序是这样的，如果把三个变量写入一个grd文件就会出错，如果单个写就会出正确的图，程序如下：
program EX1
real a(144,73,12,65),ave1(144,73),ab(144,73),wp(144,73),n(73)
open(2,file='e:\duanqishixi\hgt500.grd',form='binary') !补充正确路径

!补充数据输出路径
      do it=1,65
do imo=1,12
do j=1,73
do i=1,144
read(2)a(i,j,imo,it)
enddo;enddo;enddo;enddo
    do j=1,73
do i=1,144
do it=1,65
ave1(i,j)=ave1(i,j)+a(i,j,1,it)      !65年1月的500hPa高度场累jia
end do
end do
end do
do j=1,73
do i=1,144
      ave1(i,j)=ave1(i,j)/65.0
end do;end do
      do j=1,73
do i=1,144
      ab(i,j)=a(i,j,1,61)-ave1(i,j) !2008年1月的高度场距平
      end do
end do
      do j=1,73
do i=1,144
      n(j)=n(j)+a(i,j,1,61)
    end do
n(j)=n(j)/144.0
      end do
   do j=1,73
   do i=1,144
    wp(i,j)=n(j)
end do
end do
do j=1,73
do i=1,144
    wp(i,j)=a(i,j,1,61)-wp(i,j)!2008年1月500hPa高度纬偏场
end do;end do
    open(8,file='e:\duanqishixi\swx\new1.grd',form='binary')
    do j=1,73
do i=1,144
write(8)  ave1(i,j)
    write(8)  ab(i,j)
write(8)   wp(i,j)
end do
end do
close(8)
end
我新编的ctl你们看看有错没？
dset e:\duanqishixi\swx\new1.grd
undef -9.99E+33
title  NCEP Reanalysis PROJECT
xdef 144 linear 0.000 2.500
ydef 73 linear -90.000 2.500
zdef 1 levels 500
tdef 780 linear jan1948 1mo
vars 3
ave1 1 0
ab   1 0
wp   1 0
endvars
;
这是画出来的平均图很奇葩的准备好你们的眼睛就是上面的那个图。
快来救救我吧1

lqouc

变量要在外层循环，好好看看grads的数据格式要求吧。

TWiligHT

lqouc 发表于 2013-10-30 15:34

                               
登录/注册后可看大图

变量要在外层循环，好好看看grads的数据格式要求吧。

大神，说的具体点呗？变量在外层循环说的是fortran?

lqouc

在fortran输出的时候要把变量也循环，是放在空间格点的外层循环，如果有时间的话时间在最外层。
你做的这些用grads直接计算很简单的，用不着这么费事还要写个fortran

TWiligHT

lqouc 发表于 2013-10-30 15:57

                               
登录/注册后可看大图

在fortran输出的时候要把变量也循环，是放在空间格点的外层循环，如果有时间的话时间在最外层。
你做的这些 ...

我去试试，没办法啊，老师非得让我们编个frotran

才仔

什么都不会{:5_275:}

TWiligHT

TWiligHT 发表于 2013-10-30 16:01

                               
登录/注册后可看大图

我去试试，没办法啊，老师非得让我们编个frotran

大神我膜拜你，我给改了图出来啦

TWiligHT

才仔 发表于 2013-10-30 16:03

                               
登录/注册后可看大图

什么都不会

孩子不要灰心，上面的大神会帮助你的。

river

这个实习要求用fortran是为了提高fortran的能力吧，其实直接用grads很简单就出图了。楼主之所以用fortran出不来最主要的问题是没有理解grads所需资料的排列方式，看下手册就知道了，数据排放顺序为先经度、纬度、高度，然后是物理变量，最后是时次变化。一个x、y数据场构成一个记录。简单说就是（x，y）,z,VAR,t

还有就是你计算的距平场 平均场 纬偏场都只有一个时次，你的ctl怎么能写780个呢

TWiligHT

river 发表于 2013-10-30 18:13

                               
登录/注册后可看大图

这个实习要求用fortran是为了提高fortran的能力吧，其实直接用grads很简单就出图了。楼主之所以用fortran出 ...

谢谢啦，我改啦，都改啦，出图啦，撒花。