爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 633|回复: 1

[求助] Fortran

[复制链接]

新浪微博达人勋

发表于 2023-10-3 10:16:32 | 显示全部楼层 |阅读模式

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

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

x
     求助,结果出错,不知道怎么改。
结果:forrtl: severe (161): Program Exception -array bounds exceeded
Source
Routine
h500kk.exe
Image PC
0040151F Unknown Line
Unknown Unknown
Incrementally linked image——PC correlation disabled.
      program main
      parameter(ii=37,jj=17,mo=12,yr=4)
        integer i,j,m,y
      real var(ii,jj,mo,yr)
      real ano(ii,jj,mo,yr)
      real ave(ii,jj,mo,yr)
      real std(ii,jj,mo,yr)
      open(1,file='H500.dat')


cccccccccccccccc 读数据资料 h500  cccccccccccccccccccccccc
       do  y=1,4
       do  m=1,12
       read(1,*) !空读一行
       read(1,*) ((var(i,j,m,y),i=1,ii),j=1,jj)
       enddo
       enddo
10    format(37f8.1)
20    format(37f7.2)
       close(1)


ccccccccccccccccc以下自己编程计算ccccccccccccccccccccccccc  

        do i=1,37,1
        do j=1,17,1
        ave(i,j,2,1)=0.0
        do y=1,4
        ave(i,j,2,y)=ave(i,j,2,y)+var(i,j,2,y)
        end do
      ave(i,j,2,y)=ave(i,j,2,y)/4
        end do
        end do

        do i=1,37,1
        do j=1,17,1
        std(i,j,2,1)=0.0
        do y=1,4
        ano(i,j,2,y)=var(i,j,2,y)-ave(i,j,2,y)
        std(i,j,2,y)=std(i,j,2,y)+ano(i,j,2,y)*ano(i,j,2,y)
        end do
      std(i,j,2,y)=sqrt(std(i,j,2,y)/4)
        end do
        end do
       do  iy=1,4
       do  m=1,12
        write(10) var(1,1,m,iy)
        end do
        end do

        open(10,file='h500.grd',form="binary")
        write(10) ((ave(i,j,2,y),i=1,37),j=1,17)
        close(10)
        open(11,file='ano-h500.grd',form="binary")
        write(11) ((ano(i,j,2,y),i=1,37),j=1,17)
        close(11)
        open(12,file='std-h500.grd',form="binary")
        write(12) ((ave(i,j,2,y),i=1,37),j=1,17)
        close(12)
      end

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-10-19 15:22:46 | 显示全部楼层
open(10,file='h500.grd',form="binary")
        write(10) ((ave(i,j,2,y),i=1,37),j=1,17)
        close(10)
        open(11,file='ano-h500.grd',form="binary")
        write(11) ((ano(i,j,2,y),i=1,37),j=1,17)
        close(11)
        open(12,file='std-h500.grd',form="binary")
        write(12) ((ave(i,j,2,y),i=1,37),j=1,17)

里面的y需要给出具体值,比如 1 ,2 等。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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