refactor: remove BaseTmethod parameter and always use weekday/weekend base temperatures

Removes the BaseTmethod parameter which controlled whether to use a single
BaseT_HC value or separate weekday/weekend base temperatures for heating/
cooling degree day (HDD/CDD) calculations in the anthropogenic heat flux
module.

Changes:
- Remove BaseTMethod parameter from suews_ctrl_const.f95
- Remove BaseTMethod and BaseT_HC from update_DailyState_Day subroutine
- Always use weekday/weekend-specific BaseT_Heating and BaseT_Cooling arrays
- Simplify HDD/CDD calculation logic by removing SELECT CASE branching

This change aligns with the data model which already provides weekday/weekend
base temperature profiles, and follows the removal of BaseT_HC from input
files in SUEWS v2023a.

Note: BaseT and BaseTe remain as separate parameters for vegetation LAI
phenology calculations (growing degree days and senescence degree days).

Addresses #144

Author: sunt05

src/suews/src/suews_ctrl_const.f95

@@ -1130,7 +1130,6 @@ MODULE data_in
 
    ! ---- Model options set in RunControl --------------------------------------------------------
    INTEGER :: EmissionsMethod, & ! anthropogenic emissions method
-              BaseTMethod, & ! base temperature method for HDD/CDD calculations used in QF module
               CBLuse, & !CBL slab model used (1) or not used (0)
               MultipleMetFiles, & !Indicates whether a single met file is used for all grids (0) or one for each grid (1)
               MultipleInitFiles, & !Indicates whether a single initial conditions file is used for all grids (0) or one for each grid (1)

src/suews/src/suews_phys_dailystate.f95

@@ -82,8 +82,6 @@ SUBROUTINE SUEWS_cal_DailyState( &
       TYPE(SUEWS_STATE), INTENT(INOUT) :: modState
 
       ! INTEGER :: WaterUseMethod
-      INTEGER, PARAMETER :: BaseTMethod = 2 ! base t method [-]
-      REAL(KIND(1D0)), PARAMETER :: BaseT_HC = 18.2 !base temperature for heating degree dayb [degC] ! to be fully removed TODO
 
       ! TYPE(IRRIGATION_PRM), INTENT(IN) :: irrPrm
       ! INTEGER :: Ie_start !Starting time of water use (DOY)
@@ -503,12 +501,10 @@ SUBROUTINE SUEWS_cal_DailyState( &
                ! regular update at all timesteps of a day
                IF (execute_subroutines) THEN
                   CALL update_DailyState_Day( &
-                     BaseTMethod, &
                      DayofWeek_id, &
                      avkdn, & !input
                      Tair, &
                      Precip, &
-                     BaseT_HC, &
                      BaseT_Heating, BaseT_Cooling, &
                      nsh_real, &
                      Tmin_id, Tmax_id, lenDay_id, & !inout
@@ -746,26 +742,22 @@ END SUBROUTINE SUEWS_cal_DailyState
    ! END SUBROUTINE update_DailyState_End
 
    SUBROUTINE update_DailyState_Day( &
-      BaseTMethod, &
       DayofWeek_id, &
       avkdn, & !input
       Tair, &
       Precip, &
-      BaseT_HC, &
       BaseT_Heating, BaseT_Cooling, &
       nsh_real, &
       Tmin_id, Tmax_id, lenDay_id, & !inout
       HDD_id) !inout
       ! use time, only: id, id_prev_t
       IMPLICIT NONE
 
-      INTEGER, INTENT(IN) :: BaseTMethod
       INTEGER, DIMENSION(3), INTENT(in) :: DayofWeek_id
 
       REAL(KIND(1D0)), INTENT(IN) :: avkdn
       REAL(KIND(1D0)), INTENT(IN) :: Tair ! Ambient air temperature [degC], this can be from either forcing or diagnostic
       REAL(KIND(1D0)), INTENT(IN) :: Precip
-      REAL(KIND(1D0)), INTENT(IN) :: BaseT_HC
       REAL(KIND(1D0)), DIMENSION(2), INTENT(IN) :: BaseT_Heating
       REAL(KIND(1D0)), DIMENSION(2), INTENT(IN) :: BaseT_Cooling
       REAL(KIND(1D0)), INTENT(IN) :: nsh_real
@@ -793,18 +785,9 @@ SUBROUTINE update_DailyState_Day( &
       iu = 1 !Set to 1=weekday
       IF (DayofWeek_id(1) == 1 .OR. DayofWeek_id(1) == 7) iu = 2 !Set to 2=weekend
 
-      SELECT CASE (BaseTMethod)
-      CASE (1)
-         BaseT_Heating_use = BaseT_HC
-         BaseT_Cooling_use = BaseT_HC
-      CASE (2)
-         BaseT_Heating_use = BaseT_Heating(iu)
-         BaseT_Cooling_use = BaseT_Cooling(iu)
-
-      CASE default
-         CALL ErrorHint(75, "RunControl.nml", -999, -999, -999)
-
-      END SELECT
+      ! Use weekday/weekend-specific base temperatures for heating/cooling degree day calculations
+      BaseT_Heating_use = BaseT_Heating(iu)
+      BaseT_Cooling_use = BaseT_Cooling(iu)
 
       ! Daily min and max temp (these get updated through the day) ---------------------
       Tmin_id = MIN(Tair, Tmin_id) !Daily min T in column 3