Call GWspatialParameters

Author: ccarouge

src/offline/cable_input.F90

@@ -2578,6 +2578,7 @@ SUBROUTINE load_parameters(met,air,ssnow,veg,climate,bgc,soil,canopy,rough,rad,
     USE POPmodule, ONLY: POP_INIT
     USE POPLUC_module, ONLY: POPLUC_INIT
     USE CABLE_LUC_EXPT, ONLY: LUC_EXPT_TYPE
+    USE cable_param_module, ONLY: GWspatialParameters
 
     IMPLICIT NONE
 
@@ -2661,6 +2662,11 @@ SUBROUTINE load_parameters(met,air,ssnow,veg,climate,bgc,soil,canopy,rough,rad,
          sum_flux,veg,mp)
     WRITE(logn,*) ' CABLE variables allocated with ', mp, ' patch(es).'
 
+! call below inserted by rk4417 - phase2
+    !CALL for gw_model false and true sets constants when false
+    CALL GWspatialParameters(logn,soil,ssnow)    ! MMY gw_model = True read var from gridinfo
+    ! MMY gw_model = False use default values
+
     IF (icycle > 0 .OR. CABLE_USER%CASA_DUMP_WRITE ) &
          CALL alloc_casavariable(casabiome,casapool,casaflux, &
          casamet,casabal,mp)

src/offline/cable_parameters.F90

@@ -67,7 +67,7 @@ MODULE cable_param_module
   PRIVATE
   PUBLIC get_default_params, write_default_params, derived_parameters,         &
        check_parameter_values, report_parameters, parID_type,                &
-       write_cnp_params, consistency_ice_veg_soil 
+       write_cnp_params, consistency_ice_veg_soil, GWspatialParameters 
   INTEGER :: patches_in_parfile=4 ! # patches in default global parameter
   ! file