Skip to content

Calculate and add new fields to clm history #4288

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
slevis-lmwg wants to merge 1 commit into
base: b4b-dev
Choose a base branch
from
Draft

Calculate and add new fields to clm history #4288

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
19 changes: 18 additions & 1 deletion src/biogeophys/PhotosynthesisMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,6 +173,7 @@ module PhotosynthesisMod

real(r8), pointer, public :: rc14_canair_patch (:) ! patch C14O2/C12O2 in canopy air

real(r8), pointer, public :: psnsun_ln_patch (:) ! patch sunlit leaf photosynthesis at local noon (umol co2 /m**2/ s)
real(r8), pointer, public :: psnsun_patch (:) ! patch sunlit leaf photosynthesis (umol CO2/m**2/s)
real(r8), pointer, public :: psnsha_patch (:) ! patch shaded leaf photosynthesis (umol CO2/m**2/s)
real(r8), pointer, public :: c13_psnsun_patch (:) ! patch c13 sunlit leaf photosynthesis (umol 13CO2/m**2/s)
Expand Down Expand Up @@ -315,6 +316,7 @@ subroutine InitAllocate(this, bounds)
allocate(this%gb_mol_patch (begp:endp)) ; this%gb_mol_patch (:) = nan
allocate(this%rh_leaf_patch (begp:endp)) ; this%rh_leaf_patch (:) = nan
allocate(this%vpd_can_patch (begp:endp)) ; this%vpd_can_patch (:) = nan
allocate(this%psnsun_ln_patch (begp:endp)) ; this%psnsun_ln_patch (:) = nan
allocate(this%psnsun_patch (begp:endp)) ; this%psnsun_patch (:) = nan
allocate(this%psnsha_patch (begp:endp)) ; this%psnsha_patch (:) = nan
allocate(this%c13_psnsun_patch (begp:endp)) ; this%c13_psnsun_patch (:) = nan
Expand Down Expand Up @@ -419,6 +421,7 @@ subroutine Clean(this)
deallocate(this%gb_mol_patch )
deallocate(this%rh_leaf_patch )
deallocate(this%vpd_can_patch )
deallocate(this%psnsun_ln_patch )
deallocate(this%psnsun_patch )
deallocate(this%psnsha_patch )
deallocate(this%c13_psnsun_patch )
Expand Down Expand Up @@ -481,17 +484,19 @@ subroutine InitHistory(this, bounds)
!
! !USES:
use histFileMod , only: hist_addfld1d, hist_addfld2d
use clm_time_manager, only: is_near_local_noon
!
! !ARGUMENTS:
class(photosyns_type) :: this
type(bounds_type), intent(in) :: bounds
real(r8), pointer :: ptr_1d(:) ! pointer to 1d patch array
!
! !LOCAL VARIABLES:
integer :: begp, endp
integer :: begp, endp, begg
!---------------------------------------------------------------------

begp = bounds%begp; endp= bounds%endp
begg = bounds%begg

this%rh_leaf_patch(begp:endp) = spval
call hist_addfld1d (fname='RH_LEAF', units='fraction', &
Expand Down Expand Up @@ -548,6 +553,15 @@ subroutine InitHistory(this, bounds)
avgflag='A', long_name='sunlit leaf photosynthesis', &
ptr_patch=this%psnsun_patch)

this%psnsun_ln_patch(begp:endp) = spval
! write out local noon psn_sun (within +/- 1hr)
if ( is_near_local_noon( grc%londeg(begg), deltasec=3600 ) ) then
this%psnsun_ln_patch(begp:endp) = this%psnsun_patch(begp:endp)
end if
call hist_addfld1d (fname='PSNSUNLN', units='umolCO2/m^2/s', &
avgflag='A', long_name='sunlit leaf photosynthesis at local noon', &
ptr_patch=this%psnsun_ln_patch)

this%psnsha_patch(begp:endp) = spval
call hist_addfld1d (fname='PSNSHA', units='umolCO2/m^2/s', &
avgflag='A', long_name='shaded leaf photosynthesis', &
Expand Down Expand Up @@ -740,6 +754,7 @@ subroutine InitCold(this, bounds)
this%alphapsnsha_patch(p) = spval

if (lun%ifspecial(l)) then
this%psnsun_ln_patch(p) = 0._r8
this%psnsun_patch(p) = 0._r8
this%psnsha_patch(p) = 0._r8
if ( use_c13 ) then
Expand Down Expand Up @@ -1158,6 +1173,7 @@ subroutine TimeStepInit (this, bounds)
do p = bounds%begp, bounds%endp
l = patch%landunit(p)
if (.not. lun%lakpoi(l)) then
this%psnsun_ln_patch(p) = 0._r8
this%psnsun_patch(p) = 0._r8
this%psnsun_wc_patch(p) = 0._r8
this%psnsun_wj_patch(p) = 0._r8
Expand Down Expand Up @@ -1224,6 +1240,7 @@ subroutine NewPatchInit (this, p)
this%rc14_canair_patch(p) = spval
end if

this%psnsun_ln_patch(p) = 0._r8
this%psnsun_patch(p) = 0._r8
this%psnsha_patch(p) = 0._r8

Expand Down
Loading