SUBROUTINE sav_biomass (ng, Istr, Iend, &
& LBi, UBi, pmonth, &
& wtemp, PARz, DINwcr_loc, DINsed_loc, &
& DINwcr_sav_loc, DOwcr_loc, CO2wcr_loc, &
& LDeCwcr_loc, LDeNwcr_loc, &
& agb_loc, bgb_loc, pp_loc, agm_loc, &
& agar_loc, agbr_loc, sears_loc, agbg_loc, &
& bgag_loc, bgr_loc, bgm_loc)
!
!***********************************************************************
!****************************************** John C. Warner *************
!****************************************** Neil K. Ganju **************
!****************************************** Jeremy Testa ***************
!****************************************** Tarandeep S. Kalra *********
! !
! This routine computes equilibrium partial pressure of CO2 (pCO2) !
! in the surface seawater. !
! !
! On Input: !
! !
! Istr Starting tile index in the I-direction. !
! Iend Ending tile index in the I-direction. !
! LBi I-dimension lower bound. !
! UBi I-dimension upper bound. !
! t Water temperature (Celsius). !
! PArout PAR at depth (mu E m-2 s-1). !
! DINwcr_loc Dissolved Inorganic N in water col. (mu M) !
! DINsed_loc Dissolved Inorganic N in sediment col. (mu M) !
! DINwcr_sav_loc Change in dissolved Inorganic N in water col. !
! due to SAV model (mu M) !
! DOwcr_loc O2 interaction with bed !
! CO2wcr_loc CO2 interactions with bed !
! LdeCwcr_locLabile detrital carbon in fennel !
! LdeNwcr_locLabile detrital carbon in fennel !
! agb_loc Vector of above ground biomass (mmol N m-2) !
! bgb_loc Vector of below ground biomass (mmol N m-2) !
! !
! Estuarine SAV Model !
! Jeremy Testa, May 2015, Chesapeake Biological Laboratory !
! !
! Reference: !
! !
! Madden, C. J., Kemp, W. M. June 1996: Ecosystem Model of an !
! Estuarine Submersed Plant Community: Calibration and Simulation !
! of Eutrophication Responses: Estuarine Research Foundation !
! Vol. 19, No. 2B, p. 457-474 !
! !
!***********************************************************************
!
USE mod_kinds
USE mod_biology
!
implicit none
!
! Imported variable declarations.
!
integer, intent(in) :: ng, LBi, UBi
integer, intent(in) :: Istr, Iend
real(r8), intent(in) :: pmonth
# ifdef ASSUMED_SHAPE
real(r8), intent(in) :: wtemp(LBi:)
real(r8), intent(in) :: PARz(LBi:)
real(r8), intent(inout) :: DINwcr_loc(LBi:)
real(r8), intent(inout) :: DINsed_loc(LBi:)
real(r8), intent(inout) :: DINwcr_sav_loc(LBi:)
real(r8), intent(inout) :: DOwcr_loc(LBi:)
real(r8), intent(inout) :: CO2wcr_loc(LBi:)
real(r8), intent(inout) :: LdeCwcr_loc(LBi:)
real(r8), intent(inout) :: LdeNwcr_loc(LBi:)
real(r8), intent(inout) :: Agb_loc(LBi:)
real(r8), intent(inout) :: Bgb_loc(LBi:)
real(r8), intent(inout) :: pp_loc(LBi:)
real(r8), intent(inout) :: agm_loc(LBi:)
real(r8), intent(inout) :: agar_loc(LBi:)
real(r8), intent(inout) :: agbr_loc(LBi:)
real(r8), intent(inout) :: sears_loc(LBi:)
real(r8), intent(inout) :: agbg_loc(LBi:)
real(r8), intent(inout) :: bgag_loc(LBi:)
real(r8), intent(inout) :: bgr_loc(LBi:)
real(r8), intent(inout) :: bgm_loc(LBi:)
# else
real(r8), intent(in) :: wtemp(LBi:UBi)
real(r8), intent(in) :: PARz(LBi:UBi)
real(r8), intent(inout) :: DINwcr_loc(LBi:UBi)
real(r8), intent(inout) :: DINsed_loc(LBi:UBi)
real(r8), intent(inout) :: DINwcr_sav_loc(LBi:UBi)
real(r8), intent(inout) :: DOwcr_loc(LBi:UBi)
real(r8), intent(inout) :: CO2wcr_loc(LBi:UBi)
real(r8), intent(inout) :: LdeCwcr_loc(LBi:UBi)
real(r8), intent(inout) :: LdeNwcr_loc(LBi:UBi)
real(r8), intent(inout) :: Agb_loc(LBi:UBi)
real(r8), intent(inout) :: Bgb_loc(LBi:UBi)
real(r8), intent(inout) :: pp_loc(LBi:UBi)
real(r8), intent(inout) :: agm_loc(LBi:UBi)
real(r8), intent(inout) :: agar_loc(LBi:UBi)
real(r8), intent(inout) :: agbr_loc(LBi:UBi)
real(r8), intent(inout) :: sears_loc(LBi:UBi)
real(r8), intent(inout) :: agbg_loc(LBi:UBi)
real(r8), intent(inout) :: bgag_loc(LBi:UBi)
real(r8), intent(inout) :: bgr_loc(LBi:UBi)
real(r8), intent(inout) :: bgm_loc(LBi:UBi)
# endif
!
!
! Local variable declarations.
!
integer :: i, thresh, thresh2, nlocal
integer, parameter :: one=1
real(r8) :: ua, day_year
real(r8) :: llim, knt, nlim
real(r8) :: lmba
real(r8) :: ppbm
! real(r8) :: pp, agm, agar
! real(r8) :: agbr, sears, agbg, bgag
! real(r8) :: bgr, bgm
real(r8) :: shtdens, shtLen
real(r8) :: cff,cff1
real(r8) :: temp
real(r8) :: dtdays
real(r8), parameter :: eps=1.0e-10_r8
real(r8), parameter :: gr2mmol=1000.0_r8/14.007_r8
real(r8), parameter :: C2N_ratio=30.0_r8 ! Move this to the input file
real(r8), parameter :: gr2mmolC=1000.0_r8/12.011_r8
real(r8), parameter :: molNmolC=(14.007_r8/12.011_r8)*30.0_r8
real(r8), parameter :: pqrq=1.0_r8
!
! Initialize local variables and arrays
!
lmba = 1.0_r8
!
nlocal=INT((time(ng)-dstart*day2sec)/dt(ng))+1
! IF (nlocal==1) THEN
! DO i=Istr, Iend
! agb_loc(i)=11.9_r8
! bgb_loc(i)=71.3_r8
! DINsed_loc(i)=1000.0_r8
! END DO
! END IF
DO i=Istr, Iend
DINwcr_sav_loc(i)=0.0_r8
DOwcr_loc(i)=0.0_r8
CO2wcr_loc(i)=0.0_r8
LDeNwcr_loc(i)=0.0_r8
LDeCwcr_loc(i)=0.0_r8
END DO
!
!-----------------------------------------------------------------------
! pmonth contains the month information calculated in estuarybgc.h
!-----------------------------------------------------------------------
!
day_year = (pmonth - 52.0_r8)*365.0_r8
! Use temp to convert from grams carbon to mmol nitrogen if needed
temp=gr2mmol/C2N_ratio
dtdays=dt(ng)*sec2day
!
knt=knwc(ng)/knsed(ng)
!
DO i=Istr, Iend
IF (DINsed_loc(i).lt.1.0_r8) THEN
DINsed_loc(i)=1.0_r8
END IF
llim=PARz(i)/(PARz(i)+ki(ng))
nlim=(DINwcr_loc(i)+(knt*DINsed_loc(i)))/ &
& (knwc(ng)+(DINwcr_loc(i)+(knt*DINsed_loc(i))))
!
!-----------------------------------------------------------------------
! SAV Growth Rate
!-----------------------------------------------------------------------
!
IF ( GMODopt(ng).eq.1 ) THEN
ua=lmba*scl(ng)*(EXP(thta(ng)*(wtemp(i)-Topt(ng)))-1.0_r8)
ELSE
ua=scl2(ng)*thta2(ng)**(wtemp(i)-Topt(ng))
ENDIF
IF ( wtemp(i).gt.Tcrit(ng) ) THEN
thresh=1
ELSE
thresh=0
ENDIF
IF ( day_year.gt.273.0_r8 ) THEN
thresh2=1
ELSE
thresh2=0
ENDIF
!
!-----------------------------------------------------------------------
! Self-shading effect on growth and primary production rate
!-----------------------------------------------------------------------
!
lmba=1.0_r8-(agb_loc(i)/lmbamx(ng))**2
pp_loc(i)=agb_loc(i)*(ua*MIN(llim,nlim))
!
!-----------------------------------------------------------------------
! Above ground mortality and above ground active respiration
!-----------------------------------------------------------------------
!
agm_loc(i)=(kmag(ng)*agb_loc(i))**2
agar_loc(i)=pp_loc(i)*arsc(ng)*EXP(arc(ng)*wtemp(i))
!
!-----------------------------------------------------------------------
! Above ground basal respiration and seasonal root storage of carbon
!-----------------------------------------------------------------------
!
agbr_loc(i)=agb_loc(i)*(bsrc(ng)*EXP(rc(ng)*wtemp(i)))
sears_loc(i)=agb_loc(i)*RtSttl(ng)*thresh2
!
!-----------------------------------------------------------------------
! Translocation of above ground biomass to below ground
!-----------------------------------------------------------------------
!
agbg_loc(i)=pp_loc(i)*DOWNt(ng)
!
!-----------------------------------------------------------------------
! Translocation of below ground biomass to above ground
!-----------------------------------------------------------------------
!
bgag_loc(i)=bgb_loc(i)*trns(ng)*thresh
!
!-----------------------------------------------------------------------
! Below ground biomass respiration and Below ground biomass mortality
!-----------------------------------------------------------------------
!
bgr_loc(i)=bgb_loc(i)*bsrc(ng)*EXP(rc(ng)*wtemp(i))
bgm_loc(i)=bgb_loc(i)*(0.01_r8*EXP(km(ng)*wtemp(i)))
!
! If using carbon units, multiply this term by temp
DINsed_loc(i)=DINsed_loc(i)+(bgr_loc(i)+bgm_loc(i))*dtdays
!
!-----------------------------------------------------------------------
! Compute new AGB biomass (mmol N m-2)
!-----------------------------------------------------------------------
!
cff=pp_loc(i)*dtdays
!
!-----------------------------------------------------------------------
! If pp>(available nutrients), no growth
! assumes no new growth if growth exceeds available nutrients
!-----------------------------------------------------------------------
!
! IF(cff.gt.DINwcr_loc(i))THEN
! cff1=pp_loc(i)+bgag_loc(i)-agm_loc(i)-agar_loc(i) &
! & -agbr_loc(i)-sears_loc(i)-agbg_loc(i)
! ELSE
cff1=bgag_loc(i)-agm_loc(i)-agar_loc(i)-agbr_loc(i) &
& -sears_loc(i)-agbg_loc(i)
! ENDIF
!
!-----------------------------------------------------------------------
! Update AGB Biomass (mmol N m-2)
!-----------------------------------------------------------------------
!
agb_loc(i)=agb_loc(i)+cff1*dtdays
!
!-----------------------------------------------------------------------
! Updating Nitrogen in water column with N uptake by plant and N
! released from SAV respiration and mortality
! (temp--> converts gram Carbon units to mmol Nitrogen units)
! If using carbon units, multiply this term by temp
!-----------------------------------------------------------------------
!
IF(cff.gt.DINwcr_loc(i))THEN
DINwcr_sav_loc(i)=DINwcr_sav_loc(i)+(agar_loc(i)+ &
agbr_loc(i))*dtdays
ELSE
DINwcr_sav_loc(i)=DINwcr_sav_loc(i)+(agar_loc(i)+agbr_loc(i)- &
pp_loc(i))*dtdays
ENDIF
!
!-----------------------------------------------------------------------
! Compute new BGB Biomass (mmol N m-2)
!-----------------------------------------------------------------------
!
bgb_loc(i)=bgb_loc(i)+(sears_loc(i)+agbg_loc(i)-bgag_loc(i) &
& -bgm_loc(i)-bgr_loc(i))*dtdays
!
!-----------------------------------------------------------------------
! Compute Primary Production per unit biomass
!-----------------------------------------------------------------------
!
ppbm=(pp_loc(i)/(agb_loc(i)+eps))
!
!-----------------------------------------------------------------------
! Observed stems per square meter (average of 6 cores)
!-----------------------------------------------------------------------
!
shtdens=1100.0_r8
!
!-----------------------------------------------------------------------
! Empirical computation of shoot length from above ground biomass (cm)
!-----------------------------------------------------------------------
!
shtLen=2.27_r8*agb_loc(i)
!
!-----------------------------------------------------------------------
! O2 and CO2 interactions with bed
! If using carbon units, multiply DOwcr_loc and CO2wcr_loc by gr2mmolC
! (instead of molNmolC)
!-----------------------------------------------------------------------
!
DOwcr_loc(i)=DOwcr_loc(i)+(pp_loc(i)-agar_loc(i)-agbr_loc(i)) &
& *molNmolC*pqrq*dtdays
CO2wcr_loc(i)=CO2wcr_loc(i)+(agar_loc(i)+agbr_loc(i)-pp_loc(i)) &
& *molNmolC*pqrq*dtdays
!
!-----------------------------------------------------------------------
! Labile detrital carbon and nitrogen interactions with fennel
! If using carbon units, multiply LDeNwcr_loc by temp and
! LDeCwcr_loc by gr2mmolC (instead of molNmolC)
!-----------------------------------------------------------------------
!
LDeNwcr_loc(i)=LDeNwcr_loc(i)+(agm_loc(i))*dtdays
LDeCwcr_loc(i)=LDeCwcr_loc(i)+(agm_loc(i))*molNmolC*dtdays
END DO
!
RETURN
END SUBROUTINE SAV_BIOMASS_SUB
�