Skip to content

Commit

Permalink
enaable rootbeta parametrization of froot (#453)
Browse files Browse the repository at this point in the history
# CABLE

## Description

enable parametrisation bases on rootbeta as default behaviour

Fixes #451 

## Type of change

Please delete options that are not relevant.

## Checklist

- [x] The new content is accessible and located in the appropriate
section.
- [x] I have checked that links are valid and point to the intended
content.
- [x] I have checked my code/text and corrected any misspellings

Please add a reviewer when ready for review.


<!-- readthedocs-preview cable start -->
----
📚 Documentation preview 📚:
https://cable--453.org.readthedocs.build/en/453/

<!-- readthedocs-preview cable end -->
  • Loading branch information
JhanSrbinovsky authored Dec 16, 2024
1 parent 2d5aedd commit 3879551
Show file tree
Hide file tree
Showing 2 changed files with 246 additions and 184 deletions.
45 changes: 29 additions & 16 deletions src/coupled/ESM1.5/CABLEfilesFromESM1.5/cable_um_init_subrs.F90
Original file line number Diff line number Diff line change
Expand Up @@ -328,28 +328,26 @@ END SUBROUTINE initialize_veg
!========================================================================

SUBROUTINE init_veg_pars_fr_vegin( soil_zse )
USE cable_params_mod, ONLY : vegin
USE cable_um_tech_mod, ONLY : veg, soil
USE cable_def_types_mod, ONLY : mp, ms

real, dimension(ms) :: soil_zse
USE cable_um_tech_mod, ONLY: veg, soil
USE cable_def_types_mod, ONLY: mp, ms

IMPLICIT NONE

REAL :: soil_zse(ms)

CALL init_veg_from_vegin(1, mp, veg, soil_zse)
CALL init_veg_from_vegin(1, mp, veg, soil_zse)

!froot fixed here for all vegetation types for ACCESS
!need more flexibility in next version to read in or parameterise
veg%froot(:,1) = 0.05
veg%froot(:,2) = 0.20
veg%froot(:,3) = 0.20
veg%froot(:,4) = 0.20
veg%froot(:,5) = 0.20
veg%froot(:,6) = 0.15

RETURN

END SUBROUTINE init_veg_pars_fr_vegin

SUBROUTINE init_veg_from_vegin(ifmp,fmp, veg, soil_zse )
use cable_def_types_mod, ONLY : veg_parameter_type, ms, mp
USE cable_params_mod, ONLY : vegin

USE cable_def_types_mod, ONLY: veg_parameter_type, ms, mp
USE cable_params_mod, ONLY: vegin
USE cable_common_module, ONLY: cable_user

integer :: ifmp, & ! start local mp, # landpoints (jhan:when is this not 1 )
fmp ! local mp, # landpoints
Expand Down Expand Up @@ -407,7 +405,20 @@ SUBROUTINE init_veg_from_vegin(ifmp,fmp, veg, soil_zse )
veg%clitt(h) = vegin%clitt(veg%iveg(h))
END DO ! over each veg patch in land point

!used in CM2 - overwritten here
IF (cable_user%access13roots) THEN

! prescribed values from vegin%froot so can be set to ESM1.5 values
! which were the same for all pfts (0.05, 0.2, 0.2, 0.2, 0.2, 0.15)
! or alternate prescribed root fractions could be used with pft dependence
DO is = 1, ms
DO h=1, mp
veg%froot(h,is) = vegin%froot(is,veg%iveg(h))
ENDDO
ENDDO

ELSE

! parametrized values used in CMIP6 ACCESS-CM2
! calculate vegin%froot from using rootbeta and soil depth
! (Jackson et al. 1996, Oceologica, 108:389-411)
totdepth = 0.0
Expand All @@ -420,7 +431,9 @@ SUBROUTINE init_veg_from_vegin(ifmp,fmp, veg, soil_zse )
veg%froot(:, is) = veg%froot(:, is)-veg%froot(:,is-1)
END DO

ENDIF

RETURN
END SUBROUTINE init_veg_from_vegin


Expand Down
Loading

0 comments on commit 3879551

Please sign in to comment.