fortran循环问题

jueshizhiyao

代码是这样的，求线性系数
program main
    implicit none
    integer,parameter ::a=407,b=32
    integer :: i,j,year(a,b),yr(b)
    real :: lat(a),lon(a),raw(b),norm(b),wind(a,b),relation
    character*5 station(a)
    real sfenzi,sfenmu1,sfenmu2,s1,s2  
    real amean,bmean  
     
   open(10, file="C:\Users\Lsz\Desktop\datasets\Observation\407stations.txt",status="old")
    do i=1,a
        read(10,*)station(i),lat(i),lon(i)
    end do
    close(10)
   
    do i=1,a
        open(11,file="C:\Users\Lsz\Desktop\datasets\Average(new)\summer\yearly\300hPa\"//station(i)//".txt",status="old")
        do j=1,b
            read(11,*)year(i,j),wind(i,j) !stid(i,j),
        end do
        close(11)
    end do
   
    open(12, file="C:\Users\Lsz\Desktop\datasets\Average(new)\summer\EASM.txt",status="old")
    do j=1,b
        read(12,*)yr(j),raw(j),norm(j)
    end do
    close(12)
    s2=0.00
    bmean=0
    do j=1,b
         s2=s2+raw(j)
         bmean=s2/b
    enddo
   
    open(13,file="C:\Users\Lsz\Desktop\300!!!!!!!.txt",status="unknown")
   
    do i=1,a
         s1=0.00
         
         sfenzi=0.00
         sfenmu1=0.00
         sfenmu2=0.00
         amean=0
         do j=1,32
         s1=s1+wind(i,j)
         amean=s1/32         
         
         sfenzi=sfenzi+(wind(i,j)-amean)*(raw(j)-bmean)
         sfenmu1=sfenmu1+(wind(i,j)-amean)**2
         sfenmu2=sfenmu2+(raw(j)-bmean)**2
   
         end do
    relation=sfenzi/sqrt(sfenmu1*sfenmu2)
    write(13,100)relation !station(i),lat(i),lon(i),
    enddo
100 format(f10.5) !i5,f8.2,f8.2,
   
    end program main

是这样的，所求的sfenzi,sfenmu1,sfenmu2中所用的amean应该是32次循环后最终的amean，而现在是每一层的amean都用上了，但是循环体内部又不能同时出现两个循环，我该怎么办？主要牵扯到二维数组数据文件的提取，所以挺麻烦的，求大神解答下！

扫眉才子

你把s1=s1+wind(i,j)放在do =1,32循环外面 不就行了？？是我对你问题的描述理解有误吗

jueshizhiyao

扫眉才子 发表于 2017-4-17 08:55
你把s1=s1+wind(i,j)放在do =1,32循环外面 不就行了？？是我对你问题的描述理解有误吗

可能是我没说清楚，不可以放外面，因为j循环中还要做wind(i,1)到wind(i,32)的累加求平均

stcya

不知道有无理解错，楼主试试将红色那块程序换成，看看能不能行

do j=1,32
    s1=sum(wind(i,:))
    amean=s1/32         
    sfenzi=sfenzi+(wind(i,j)-amean)*(raw(j)-bmean)
    sfenmu1=sfenmu1+(wind(i,j)-amean)**2
    sfenmu2=sfenmu2+(raw(j)-bmean)**2
end do
或者——————————————————————————

do j=1,32
    s1=s1+wind(i,j)
    amean=s1/32   
enddo

do j=1,32
    sfenzi=sfenzi+(wind(i,j)-amean)*(raw(j)-bmean)
    sfenmu1=sfenmu1+(wind(i,j)-amean)**2
    sfenmu2=sfenmu2+(raw(j)-bmean)**2
end do

金色牛神

哈哈，这倒没那么复杂。既然已经知道不能交叉循环，那就把amean的计算单独提出一个循环算成结果就行了，循环逻辑写错了。
do j=1,32
         s1=s1+wind(i,j)
         amean=s1/32   
上面代码改一下：
DO K = 1,32
    s1=s1+wind(i,K)     
END DO
amean=s1/32
do j=1,32
  继续下面计算sfenzi...之类的

jueshizhiyao

stcya 发表于 2017-4-17 11:09
不知道有无理解错，楼主试试将红色那块程序换成，看看能不能行

do j=1,32

可以的，就是这个意思。不过循环是二层循环，改一下就可以了！多谢！

jueshizhiyao

金色牛神 发表于 2017-4-17 11:21
哈哈，这倒没那么复杂。既然已经知道不能交叉循环，那就把amean的计算单独提出一个循环算成结果就行了，循 ...

嗯嗯可以的，多谢了大神