计算积温的fortran源代码~~希望对大家有所帮助~

miss_timo

SUBROUTINE HDJW(TONEYEAR,IDAYS)
REAL TONEYEAR(366),XAVE(IDAYS),JIWEN
INTEGER XSHI,XZHONG
  XAVE=0.
  XAVE(1)=(TONEYEAR(1)+TONEYEAR(2)+TONEYEAR(3))/FLOAT(3)
  XAVE(2)=(TONEYEAR(1)+TONEYEAR(2)+TONEYEAR(3)+TONEYEAR(4))/FLOAT(4)
  XAVE(IDAYS)=(TONEYEAR(IDAYS-2)+TONEYEAR(IDAYS-1)+TONEYEAR(IDAYS))/FLOAT(3)
  XAVE(IDAYS-1)=(TONEYEAR(IDAYS-3)+TONEYEAR(IDAYS-2)+TONEYEAR(IDAYS-1)+TONEYEAR(IDAYS))/FLOAT(4)
  DO K=3,IDAYS-2
    XAVE(K)=(TONEYEAR(K-2)+TONEYEAR(K-1)+TONEYEAR(K)+TONEYEAR(K+1)+TONEYEAR(K+2))/FLOAT(5)
  ENDDO
  XSHI=0
  XZHONG=0
loop1:  DO K=1,IDAYS,1
    IF(XAVE(K).GE.FLOAT(10))THEN
     DO M=-2,2,1
    IF(TONEYEAR(K+M).GE.FLOAT(10))THEN
       XSHI=K+M
    EXIT loop1
    ENDIF
   ENDDO
  ENDIF
   ENDDO loop1   
loop2: DO K=IDAYS,XSHI,-1
    IF(XAVE(K).GE.FLOAT(10))THEN
     DO M=2,-2,-1
    IF(TONEYEAR(K+M).GE.FLOAT(10))THEN
       XZHONG=K+M
    EXIT loop2
    ENDIF
  ENDDO
  ENDIF
   ENDDO loop2
  JIWEN=0.
  DO K=XSHI,XZHONG
  JIWEN=JIWEN+TONEYEAR(K)
  ENDDO
  OPEN(2,FILE='E:\result.txt')
    WRITE(2,100)'起点温度:',TONEYEAR(XSHI), '   终点温度:',TONEYEAR(XZHONG),'   积温:',JIWEN
  100 format(1x,3(a,f8.1))
  END

baiduwenku

有对算法的描述吗？

baiduwenku

没有注释，看着有些费劲。