Simulate Run-of-the-River (ROR) dams

- Introduced a conditional block that checks for ROR operation based on RULE = 0.
- For ROR dams, FLOWOUT is set equal to FLOWIN, constrained by minimum environmental flow (Qmin).
- Reservoir storage (VOL) is maintained constant.
- Power generation (ENERGYPRO) is calculated using turbine flow and effective head (REALHEAD).
- Reservoir water levels (HTK, HHO) are set to maximum operating head (HRESERMAX).

This update allows users to represent non-storage dams by setting operation strategy (RULE = 0) in the reservoir parameter file.

Author: HishamEldardiry

routing/model/reservoir.f

@@ -990,6 +990,9 @@ SUBROUTINE MAKE_CONVOLUTIONRS
                 READ(25,*) ENV_FLOW(J,W)  ! Read the percentage of environmental flow (% of maximum flow)
                 READ(25,*)
                 READ(25,*) RULE(J,W)
+		IF (RULE(J,W) .EQ. 0) THEN ! run-of-the-river dam
+                        READ(25,*) Dummy
+	        END IF
                 IF (RULE(J,W) .EQ. 1) THEN ! simplified rule curve
                     READ(25,*) HMAX(J,W), HMIN(J,W), OP1(J,1,W),OP1(J,2,W)
                     IF (HMAX(J,W) .LT. HMIN(J,W)) THEN
@@ -1337,6 +1340,27 @@ SUBROUTINE MAKE_CONVOLUTIONRS
                     ELSE
                         ENERGYPRO(J,I,W) = ENERGYPRO(J,I,W) *((HHO(J,I,W)+HHO(J,I+1,W))/2-(HRESERMAX(J,W)-REALHEAD(J,W)))/1000			! this part is for hydropower production estimation, ignore if work with irrigation reservoirs   
                     END IF 
+
+  		    ! Check for minimum environmental flow
+		    IF (FLOWOUT(J,I,W) .LE. Qmin(J,I,W)) THEN
+			FLOWOUT(J,I,W) = Qmin(J,I,W)
+		    END IF 
+
+      		    ! Set Hydrologic Budget (States & Fluxes) for ROR Dams (Run-of-the-River)
+                    IF (RULE(J,W) .EQ. 0) THEN
+                        print*,'RUN-OF-THE-RIVER (ROR) DAM........'
+                        FLOWOUT(J,I,W) = FLOWIN(J,I,W)
+                        IF (FLOWOUT(J,I,W) .LE. Qmin(J,I,W)) THEN
+                            FLOWOUT(J,I,W) = Qmin(J,I,W)
+                        END IF
+                        FLOWOUT_TURB(J,I,W) = FLOWOUT(J,I,W)
+                        VOL(J,I+1,W) = VRESER(J,W,I)
+                        ENERGYPRO(J,I,W) = 0.9 * 9.81 * FLOWOUT(J,I,W)
+                        ENERGYPRO(J,I,W) = ENERGYPRO(J,I,W) *(REALHEAD(J,W))/1000
+                        HTK(J,I,W) = REALHEAD(J,W)
+                        HHO(J,I,W)= HRESERMAX(J,W)
+                        HHO(J,I+1,W)= HRESERMAX(J,W)
+                    END IF
 c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx            
 c     STEP 6: Propagate water to the downstream reservoir, considering the time lag
 c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
@@ -1351,11 +1375,7 @@ SUBROUTINE MAKE_CONVOLUTIONRS
                         VRESER(J,W,I) = 0 
                         REALHEAD(J,W) = 0
                         QRESER(J,W) = 0
-                        FLOWOUT_TURB(J,I,W)=0.0  
-                    ELSE
-                        IF (FLOWOUT(J,I,W) .LE. Qmin(J,I,W)) THEN
-                            FLOWOUT(J,I,W) = Qmin(J,I,W)
-                        END IF   
+                        FLOWOUT_TURB(J,I,W)=0.0        
                     END IF
 
                     ! Reading UH_R files for reservoirs