请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 99|回复: 0

[求助] Fortran用冬季1978~2007年共30年热带太平洋海温场资料求气候场并画图

[复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
fortran 程序
program sx01
implicit none
integer,parameter::nx=32,ny=7,nz=1,nt=30
real w(nx,ny,nz,nt),qhc(nx,ny,nz)
integer ii,ij,iz,it
w(:,:,:,:)=-9.99e33
open(11,file='G:\16qxtj\NCEP_TPSST_30y_Wt.dat',form='binary')
do it=1,nt
read(11) (((w(ii,ij,iz,it),ii=1,nx),ij=1,ny),iz=1,nz)
enddo
do iz=1,nz
do ij=1,ny
do ii=1,nx
qhc(ii,ij,iz)=0.0
do it=1,nt
qhc(ii,ij,iz)=qhc(ii,ij,iz)+w(ii,ij,iz,it)
enddo
qhc(ii,ij,iz)=qhc(ii,ij,iz)/30.
enddo
enddo
enddo
open(12,file='G:\16qxtj\qhc.grd',form='binary')
do iz=1,nz
write(12) ((qhc(ii,ij,iz),ii=1,nx),ij=1,ny)
enddo
end
ctl文件
dset G:\16qxtj\qhc.grd
undef  -9.99e-33
title  qhc
xdef  32  linear  121.88  5.650
ydef   7  levels     -18.0950  -12.3808   -6.6666   -0.9524    4.7618   10.4761   16.1902
zdef   1  levels   1000
tdef  1  linear  0z01Jan1979  1yr
vars   1
qhc   0   99  TPSST
endvars
画出来的图很奇怪,但是不知道错在哪里了,请大神指导,谢谢~
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表