search

CABLE-LSM / CABLE-Trac-archive

Archive CABLE Trac contents as issues
Other
0 stars 0 forks source link

Added plant pool solution to analytical solution #242

Closed penguian closed 2 years ago

penguian commented 4 years ago

keyword_nogit owner:yxw599@nci.org.au resolution_fixed type_defect | by mgk576

I've added the solution for the steady-state plant pools to the analytical solution. Changes in spincasacnp.F90 and casa_cable.F90 in https://trac.nci.org.au/trac/cable/browser/branches/Users/mgk576/trunk_analytical

The main bits:

! Solve the C pools
casapool%cplant(npt,leaf) = (avgcnpp(npt) * avg_af(npt)) - &
                              (casapool%cplant(npt,leaf) * avg_lf(npt))

casapool%cplant(npt,wood) = (avgcnpp(npt) * avg_aw(npt)) - &
                              (casapool%cplant(npt,wood) * avg_lw(npt))

casapool%cplant(npt,froot) = (avgcnpp(npt) * avg_lr(npt)) - &
                              (casapool%cplant(npt,froot) * avg_lr(npt))

And

! Calculate average allocation fractions  for the plant pools
avg_af = avg_ar + casaflux%fracCalloc(:,leaf)
avg_aw = avg_ar + casaflux%fracCalloc(:,wood)
avg_ar = avg_aw + casaflux%fracCalloc(:,froot)

! Calculate average turnover fractions for the plant pools
avg_lf = avg_lr + (1.0 / casaflux%kplant(:,leaf))
avg_lw = avg_lr + (1.0 / casaflux%kplant(:,wood))
avg_lr = avg_lw + (1.0 / casaflux%kplant(:,froot))

Issue migrated from trac:242 at 2023-11-27 11:31:31 +1100

penguian commented 4 years ago

@mgk576@nci.org.au changed owner from bep599 to yxw599

penguian commented 4 years ago

@mgk576@nci.org.au changed component from MPI to biogeochem

penguian commented 2 years ago

@jxs599@nci.org.au changed status from new to closed

penguian commented 2 years ago

@jxs599@nci.org.au set resolution to fixed

penguian commented 2 years ago

@jxs599@nci.org.au set milestone to 1. Closed

penguian commented 2 years ago

@jxs599@nci.org.au commented

pushed @ 8689

penguian commented 1 year ago

@ccc561@nci.org.au set keywords to nogit