110.f90

e额外e额外 饿繁荣如同热 饿繁荣如同热e额外 饿繁荣如同热e额外 饿繁荣如同热

SUBROUTINE OUT (q11,Ccl)
REAL Ccl(9,24),q11(25)
REAL ts(24),q1_add(24),q1_aver(24)
REAL CLr2(72),CLr(72),CLR1(24,72),Ccl1(18,24)
REAL add(2,24),CLR_1(3,24),over(24)
INTEGER b(24),num(24)
COMMON K,TIME
!********************************************************************
! 以下循环计算未知参数的值
!********************************************************************
DO I0=1,9
DO I1=1,TIME
Ccl1(2*I0,I1)=Ccl(I0,I1)
ENDDO
ENDDO
DO I0=1,8
DO I1=1,TIME
Ccl1(2*I0+1,I1)=(Ccl1(2*(I0+1),I1)+Ccl1(2*I0,I1))/2
IF(K.NE.1) Ccl1(1,I1)=Ccl1(2,I1)
IF(K.EQ.1) Ccl1(1,I1)=(Ccl1(2,I1))/2
ENDDO
ENDDO
WRITE(*,24)((Ccl1(I0,I1),I0=1,18),I1=1,TIME)
24 FORMAT(18F4.3)
!********************************************************************
! 初始化参数
!********************************************************************
N=0
R=0
J1=0
T=0
W=0
q11_add=0
!!!!****************************************************************
DO I=1,TIME
num(I)=0
IF(q11(I).EQ.0.0) THEN
N=N+1
b(N)=I
ENDIF
ENDDO
!********************************************************************
!以下四个循环计算不同特征的数列的分段数目以及每段元素的和
!********************************************************************

IF(q11(1).EQ.0.0.AND.q11(TIME).EQ.0.0) THEN
DO J=2,N
X=0
J1=J1+1
ts(J1)=b(J)-b(J-1)-1
DO L=b(J-1),b(J)
X=q11(L)+X
ENDDO
R=R+1
q1_add(R)=X
num(R)=b(J-1)+1
PRINT *,num(7)
ENDDO
ENDIF
!!****************************************************************
IF(q11(1).EQ.0.0.AND.q11(TIME).NE.0.0) THEN
DO J=2,N+1
b(N+1)=TIME+1
J1=J1+1
ts(J1)=b(J)-b(J-1)-1
ENDDO
b(N+1)=TIME
DO K1=2,N+1
X=0
DO L=b(K1-1),b(K1)
X=q11(L)+X
ENDDO
R=R+1
q1_add(R)=X
num(R)=b(K1-1)+1
ENDDO 
ENDIF
!!****************************************************************	
IF(q11(1).NE.0.0.AND.q11(TIME).EQ.0.0) THEN
DO J=1,N
X=0
IF(J.EQ.1.0) THEN
S=0
ELSE
S=b(J-1)
ENDIF
ts(J)=b(J)-S-1
ENDDO
DO 208 M=N+1,2,-1
     b(M)=b(M-1)
208 CONTINUE
     b(1)=1
DO 209,K1=1,N
X=0
DO 210,L=b(K1),b(K1+1)
X=q11(L)+X
210 CONTINUE
R=R+1
q1_add(R)=X
num(R)=b(K1)+1
209 CONTINUE
ENDIF
!!!****************************************************************
IF(q11(1).NE.0.0.AND.q11(TIME).NE.0.0) THEN
DO 211,J=1,N+1
X=0
IF(J.EQ.1.0) THEN
S=0
ELSE
S=b(J-1)
ENDIF
b(N+1)=TIME+1
ts(J)=b(J)-S-1
211 CONTINUE
 DO 212 M=N+1,2,-1
     b(M)=b(M-1)
212 CONTINUE
     b(1)=1
	 b(N+2)=TIME
DO 213,K1=1,N+1
X=0
DO 214,L=b(K1),b(K1+1)
X=q11(L)+X
214 CONTINUE
R=R+1
q1_add(R)=X
num(R)=b(K1)+1
213 CONTINUE
ENDIF

!!****************************************************************
DO 1200,K1=1,TIME
q11_add=q11_add+q11(K1)
1200 CONTINUE

!********************************************************************
!以下循环计算数列中没有0点的分段数目以及每段元素的和
!********************************************************************
IF(N.EQ.0) THEN
DO 100 K5=1,TIME
ts(K5)=1
q1_add(K5)=q11(K5)
num(K5)=K5-1
100 CONTINUE 
ENDIF
!********************************************************************
! 初始化参数
!********************************************************************
DO 101 ,K8=1,24
q1_aver(K8)=0
101 CONTINUE
DO 12,K0=1,72
CLr2(K0)=0
CLR(K0)=0
12 CONTINUE
DO 11,K0=1,72
DO 10,K1=1,24
CLR1(K1,K0)=0
10 CONTINUE
11 CONTINUE
!********************************************************************
!以下循环计算数列中每段元素的平均数
!********************************************************************
DO 1000,L=1,TIME
W=W+1
IF(q1_add(L).EQ.0.OR.ts(L).EQ.0) THEN
q1_aver(W)=0
ELSE
q1_aver(W)=q1_add(L)/ts(L)
WRITE(*,*) ts(L),L  
ENDIF
WRITE(*,*) q1_add(L),L
!********************************************************************
!以下循环查参数计算q1*Ccl,并且将数列的列下标偏移8+num(L)
!********************************************************************
IF(q11(1).NE.0) THEN
num(1)=0
ENDIF
DO 999,L1=1,TIME
IF(ts(L).EQ.0.OR.q11_add.EQ.0)THEN
CLR1(L,L1)=0
ELSE
CLR1(L,L1+8+num(L))=Ccl1(ts(L),L1)*q1_aver(W)
ENDIF
999 CONTINUE
1000 CONTINUE
!********************************************************************
!以下循环将数列的列下标相同的项相加，使得二维数组变成一维
!********************************************************************
DO 1009, K2=1,3*TIME
DO 1010, K1=1,TIME 
CLr2(K2)=CLr2(K2)+CLR1(K1,K2)
1010 CONTINUE
1009 CONTINUE
WRITE(*,*)(CLr2(J),J=1,3*TIME)  
!********************************************************************
!以下循环将数列的列下标界限设置为TIME，使得叠加后的一维数组变成二维
!********************************************************************
DO 1011 K1=1,3
DO 1012 K2=1,TIME
CLR_1(K1,K2)=CLr2((K1-1)*TIME+K2)
1012 CONTINUE
1011 CONTINUE
!********************************************************************
!以下循环将数列的列下标相同的项相加，使得二维数组CLR_1变成一维CLr，元素个数TIME
!********************************************************************
DO 1013 K1=1,3
DO 1014 K2=1,TIME
CLr(K2)=CLR_1(K1,K2)+CLr(K2)
1014 CONTINUE
1013 CONTINUE
!********************************************************************
!以下代码实现将两种情况的输出数值存在二维数组add中的功能
!********************************************************************
N5=N5+1
IF(N5.GT.2.0) N5=0 
DO 888 K0=1,TIME
add(N5,K0)=CLr(K0)    
WRITE(*,*),K0,CLr(K0)  
888 CONTINUE
!********************************************************************
!以下代码实现当两种情况计算完毕后，实现两种情况数值的叠加
!********************************************************************
IF(N5.EQ.2.0) THEN
DO 889 K1=1,2
DO 887 K2=1,TIME
 over(K2)=over(K2)+add(K1,K2)
887 CONTINUE
889 CONTINUE
!********************************************************************
!以下代码实现参数输入和结果输出功能
!********************************************************************
WRITE (*,99)
99 FORMAT(//)
IF(K.EQ.1) PRINT *,'人体散热引起的各时刻冷负荷为(W):'
IF(K.EQ.2) PRINT *,'设备和用具散热引起的各时刻冷负荷为(W)：'
IF(K.EQ.3) PRINT *,'照明散热引起的各时刻冷负荷为(W):'
DO 250,K0=1,TIME
IF(K.EQ.1) THEN
WRITE(*,70)K0,over(K0)
70 FORMAT(//1X, I2':00'F12.4)
ELSEIF(K.EQ.2) THEN
WRITE(*,71)K0,over(K0)
71 FORMAT(//1X, I2':00'F12.4)
ELSEIF(K.EQ.3) THEN
WRITE(*,72)K0-1,over(K0)
72 FORMAT(//1X, I2':00'F12.4)
ENDIF
250 CONTINUE
ENDIF
RETURN
END
!!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
!!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PROGRAM MAIN
REAL Ccl_r(9,24)
REAL q11_y(24),q11_w(24)
REAL Ccl_sz(9,24),Ccl_sw(9,24)
REAL Ccl_dm_24(9,24),Ccl_dm_16(9,16),Ccl_dm_12(9,12)
REAL Ccl_dan_24(9,24),Ccl_dan_16(9,16),Ccl_dan_12(9,12)
COMMON K,TIME
WRITE(*,90)
90 FORMAT(1X,'请输入所要计算冷负荷的类型:'//1X,&
             '[1]-人体散热引起的冷负荷'//1X,&
			 '[2]-设备和用具显热散热引起的冷负荷'//1X,&
             '[3]-照明散热引起的冷负荷'//)			 
89 PRINT *,'请输入1,2,3'
READ(*,*) K
IF(K.NE.1.0.AND.K.NE.2.0.AND.K.NE.3.0) GOTO 89
!!****************************************************************
IF(K.EQ.1.0) THEN
TIME=24
WRITE(*,111)
111 FORMAT(//1X,'输入24小时的显热(W)：'//)
READ(*,*)(q11_y(I),I=1,TIME)
CALL OUT(q11_w,Ccl_r)
CALL OUT(q11_y,Ccl_r)
ENDIF
!!****************************************************************
IF(K.EQ.2.0) THEN
TIME=24
WRITE(*,112)
112 FORMAT(//1X,'输入有罩设备24小时的显热(W)：'/1X,&
                '(没有该类型设备输入0)')