From c852e41e3ced4965552ffbb9425a747088ca1819 Mon Sep 17 00:00:00 2001
From: zzx <784670282@qq.com>
Date: Fri, 9 May 2025 17:55:58 +0800
Subject: [PATCH] =?UTF-8?q?=E6=9B=B4=E6=96=B0=20Fortran/SUB=5FQZ.f90?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 Fortran/SUB_QZ.f90 | 44 +++++++++++++++++++++++++++-----------------
 1 file changed, 27 insertions(+), 17 deletions(-)

diff --git a/Fortran/SUB_QZ.f90 b/Fortran/SUB_QZ.f90
index 45b7c1c..f91e24e 100644
--- a/Fortran/SUB_QZ.f90
+++ b/Fortran/SUB_QZ.f90
@@ -1,5 +1,5 @@
 !	==================================================================================
-!	 用四点隐格式计算未知时层水位Z和流量Q的子程序
+!	 鐢ㄥ洓鐐归殣鏍煎紡璁＄畻鏈煡鏃跺眰姘翠綅Z鍜屾祦閲廞鐨勫瓙绋嬪簭
 !	==================================================================================
 	subroutine sub_QZ(	NRIVER				,&
 						NSECT				,&
@@ -62,14 +62,18 @@
     INTEGER::Ns(nriver),Pc(nriver),Nrc(krc,nriver),Lc(krc,nriver)
     REAL::Dric(krc,nriver),Qj(ndata,krc,nriver),Asave(krc,nriver)
     ! section
-	REAL::ds(mriver,nriver),bd(nsect,nriver),zd(nsect,nriver),sm(nsect,nriver), rough(nsect,nriver)
+	REAL::ds(mriver,nriver),bd(nsect,nriver),zd(nsect,nriver),&
+		sm(nsect,nriver), rough(nsect,nriver)
     ! boundary
-	INTEGER::UB1(nriver),UB2(nriver),DB1(nriver),DB2(nriver),NUB(2,nriver),NDB(2,nriver)
-    REAL::UBV(ndata,nriver),Aphi(2,nriver),DBV(ndata,nriver),Gate(4,nriver),ql,Zctr
+	INTEGER::UB1(nriver),UB2(nriver),DB1(nriver),DB2(nriver),&
+		NUB(2,nriver),NDB(2,nriver)
+    REAL::UBV(ndata,nriver),Aphi(2,nriver),DBV(ndata,nriver),&
+		Gate(4,nriver),ql,Zctr
     ! calcu
 	REAL::dt,sita,Bsor1,Bsor2
     ! zzqq
-    REAL::Z0(nsect,nriver),Q0(nsect,nriver),Z(nsect,nriver),Q(nsect,nriver),V(nsect,nriver)
+    REAL::Z0(nsect,nriver),Q0(nsect,nriver),Z(nsect,nriver),&
+		Q(nsect,nriver),V(nsect,nriver)
     ! r_bv
 	REAL::condu,condd(3)
     ! BARC
@@ -83,7 +87,7 @@
 	if(UB1(river).eq.2)Q(1,river)=condu
 
 !	----------------------------------------------------------------------------------
-!	计算递推公式系数L、M、P、R
+!	璁＄畻閫掓帹鍏紡绯绘暟L銆丮銆丳銆丷
 !	----------------------------------------------------------------------------------
 	P(1)=condu
 	R(1)=0.0
@@ -137,9 +141,9 @@
 
 	if(Ik.ne.1)then		  
 !	------------------------------------------------------
-!    非汊点河段计算
+!    闈炴眾鐐规渤娈佃绠�
 !
-!	计算河段M点处的水面宽B、面积A、水力半径R、流量Q和水位Z
+!	璁＄畻娌虫M鐐瑰鐨勬按闈㈠B銆侀潰绉疉銆佹按鍔涘崐寰凴銆佹祦閲廞鍜屾按浣峑
 !	------------------------------------------------------
 	call sub_sect(			NRIVER				,&
 							NSECT				,&
@@ -187,11 +191,12 @@
 	Zm=f(Z0(k,river),Z(k,river),Z0(k+1,river),Z(k+1,river))
 
 !	------------------------------------------------------
-!	计算河段方程中的系数
+!	璁＄畻娌虫鏂圭▼涓殑绯绘暟
 !	------------------------------------------------------
       a1=1.0
 	c1=2.0*sita*dt/(ds(k,river)*Bm)
-	e1=Z0(k,river)+Z0(k+1,river)+(1.0-sita)/sita*c1*(Q0(k,river)-Q0(k+1,river))+2.0*dt*ql/Bm
+	e1=Z0(k,river)+Z0(k+1,river)+(1.0-sita)/sita*c1*(Q0(k,river)&
+		-Q0(k+1,river))+2.0*dt*ql/Bm
 	a2=2.0*sita*dt/ds(k,river)*((Qm/Am)**2.0*Bm-9.8*Am)
 	c2=1.0-4.0*sita*dt/ds(k,river)*Qm/Am
 	d2=1.0+4.0*sita*dt/ds(k,river)*Qm/Am
@@ -227,11 +232,14 @@
 							Rs					,&
 							Cs					)
 	AZm2=As
-	e2=(1.0-sita)/sita*a2*(Z0(k+1,river)-Z0(k,river))+(1.0-4.0*(1.0-sita)*dt/ds(k,river)*Qm/Am)*Q0(k+1,river)+(1.0+4.0*(1.0-sita)*dt/ds(k,river)*Qm/Am)*Q0(k,river)+2.0*dt*(Qm/Am)**2.0*(AZm2-AZm1)/ds(k,river)-2.0*dt*9.8*abs(Qm)*Qm/(Am*Cm*Cm*Rm)
+	e2=(1.0-sita)/sita*a2*(Z0(k+1,river)-Z0(k,river))+(1.0-4.0*&
+		(1.0-sita)*dt/ds(k,river)*Qm/Am)*Q0(k+1,river)+(1.0+4.0*&
+		(1.0-sita)*dt/ds(k,river)*Qm/Am)*Q0(k,river)+2.0*dt*&
+		(Qm/Am)**2.0*(AZm2-AZm1)/ds(k,river)-2.0*dt*9.8*abs(Qm)*Qm/(Am*Cm*Cm*Rm)
 	else
 
 !	--------------------------------------------------------
-!	汊点计算处理（普通汊点Lc=1，集中入流Lc=2，调蓄水库Lc=3）
+!	姹婄偣璁＄畻澶勭悊锛堟櫘閫氭眾鐐筁c=1锛岄泦涓叆娴丩c=2锛岃皟钃勬按搴揕c=3锛�
 !	--------------------------------------------------------
 	    if(Lc(Ic,river).eq.1)then
 		do 35 j=1,nriver
@@ -279,7 +287,7 @@
 	end if
 
 !	------------------------------------------------------
-!	求上游边界条件为水位过程线的递推系数
+!	姹備笂娓歌竟鐣屾潯浠朵负姘翠綅杩囩▼绾跨殑閫掓帹绯绘暟
 !	------------------------------------------------------
 	if(UB1(river).eq.1)then
 	Y1=a1*R(k)-c1
@@ -293,7 +301,7 @@
 	R(k+1)=(d2*Y1-c1*Y2)/Y5
 	end if
 !	------------------------------------------------------
-!	求上游边界条件为流量过程线的递推系数
+!	姹備笂娓歌竟鐣屾潯浠朵负娴侀噺杩囩▼绾跨殑閫掓帹绯绘暟
 !	------------------------------------------------------
 	if(UB1(river).eq.2)then
 	Y1=a1-c1*R(k)
@@ -317,10 +325,11 @@
 	CC1=CC2
 100	continue
 !	----------------------------------------------------------------------------------
-!	回代求出流量Q和水位Z
+!	鍥炰唬姹傚嚭娴侀噺Q鍜屾按浣峑
 !	----------------------------------------------------------------------------------
 	if(UB1(river).eq.1)then
-	Q(Ns(river),river)=(condd(3)-condd(1)*P(Ns(river)))/(condd(1)*R(Ns(river))+condd(2))
+	Q(Ns(river),river)=(condd(3)-condd(1)*P(Ns(river)))/&
+		(condd(1)*R(Ns(river))+condd(2))
 	do 120 i=Ns(river),2,-1
 	Z(i,river)=P(i)+R(i)*Q(i,river)
 	Q(i-1,river)=L(i)+M(i)*Q(i,river)
@@ -328,7 +337,8 @@
 	end if
 
 	if(UB1(river).eq.2)then
-	Z(Ns(river),river)=(condd(3)-condd(2)*P(Ns(river)))/(condd(2)*R(Ns(river))+condd(1))
+	Z(Ns(river),river)=(condd(3)-condd(2)*P(Ns(river)))/&
+		(condd(2)*R(Ns(river))+condd(1))
 	do 130 i=Ns(river),2,-1
 	Q(i,river)=P(i)+R(i)*Z(i,river)
 	Z(i-1,river)=L(i)+M(i)*Z(i,river)