Skip to content

Add carbon molar mass physical constant #687

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

Merged
SeanBryan51 merged 1 commit into from
Feb 24, 2026
Merged

Add carbon molar mass physical constant #687

Show file tree
Hide file tree
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
2 changes: 1 addition & 1 deletion documentation/docs/user_guide/constants/physical_constants.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,4 +20,4 @@
| cswat | 4.218E3 | \( J \cdot kg^{-1} \cdot K^{-1} \) | Specific heat for water at \( 0^{\circ}C \) |
| density_liq | 1000.0 | \( kg \cdot m^{-3} \) | Density of liquid water |
| density_ice | 921.0 | \( kg \cdot m^{-3} \) | Density of ice |
| c_molar_mass | 1.201E-5 | \( \mu g / mol \) | Molar mass of carbon |
55 changes: 28 additions & 27 deletions src/offline/cable_output.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,7 @@ MODULE cable_output_module
USE netcdf
USE cable_common_module, ONLY: filename, calcsoilalbedo, CurYear,IS_LEAPYEAR, cable_user,&
gw_params
USE cable_phys_constants_mod, ONLY: c_molar_mass
IMPLICIT NONE
PRIVATE
PUBLIC open_output_file, write_output, close_output_file, create_restart, check_and_write, output_par_settings_type
Expand Down Expand Up @@ -1887,7 +1888,7 @@ SUBROUTINE write_output(dels, ktau, met, canopy, casaflux, casapool, casamet, ss
END IF
IF (output%NEE) THEN
! NEE: net ecosystem exchange [umol/m^2/s]
CALL generate_out_write_acc(ovid%NEE, 'NEE', out%NEE, REAL(canopy%fnee/1.201E-5, 4), ranges%NEE, patchout%NEE, out_settings)
CALL generate_out_write_acc(ovid%NEE, 'NEE', out%NEE, REAL(canopy%fnee/c_molar_mass, 4), ranges%NEE, patchout%NEE, out_settings)
END IF

!-----------------------WRITE SOIL STATE DATA-------------------------------
Expand Down Expand Up @@ -2087,25 +2088,25 @@ SUBROUTINE write_output(dels, ktau, met, canopy, casaflux, casapool, casamet, ss
! GPP: gross primary production C by veg [umol/m^2/s]
! added frday in the calculation of GPP (BP may08)

! temp_acc = REAL((-1.0 * canopy%fpn) / 1.201E-5, 4)
! temp_acc = REAL((-1.0 * canopy%fpn) / c_molar_mass, 4)
Comment thread
SeanBryan51 marked this conversation as resolved.
IF (output%GPP) THEN
CALL generate_out_write_acc(ovid%GPP, 'GPP', out%GPP, REAL((-1.0*canopy%fpn + canopy%frday) &
/1.201E-5, 4), ranges%GPP, patchout%GPP, out_settings)
/c_molar_mass, 4), ranges%GPP, patchout%GPP, out_settings)
END IF

! NPP: net primary production of C by veg [umol/m^2/s]
IF (output%NPP) THEN
! Add current timestep's value to total of temporary output variable:
!out%NPP = out%NPP + REAL((-1.0 * canopy%fpn - canopy%frp &
! - casaflux%clabloss/86400.0) / 1.201E-5, 4)
! - casaflux%clabloss/86400.0) / c_molar_mass, 4)
Comment thread
SeanBryan51 marked this conversation as resolved.
! vh ! expression below can be slightly different form that above in cases where
! leaf maintenance respiration is reduced in CASA
! (relative to its original value calculated in cable_canopy)
! in order to avoid negative carbon stores.
IF (output%casa) THEN
temp_acc = casaflux%cnpp/86400.0/1.201E-5
temp_acc = casaflux%cnpp/86400.0/c_molar_mass
ELSE
temp_acc = (-1.0*canopy%fpn - canopy%frp)/1.201E-5 ! &
temp_acc = (-1.0*canopy%fpn - canopy%frp)/c_molar_mass ! &
END IF
CALL generate_out_write_acc(ovid%NPP, 'NPP', out%NPP, temp_acc, ranges%NPP, patchout%NPP, out_settings)
END IF
Expand All @@ -2114,36 +2115,36 @@ SUBROUTINE write_output(dels, ktau, met, canopy, casaflux, casapool, casamet, ss
IF (output%AutoResp) THEN
! Add current timestep's value to total of temporary output variable:
!out%AutoResp = out%AutoResp + REAL((canopy%frp + canopy%frday + casaflux%clabloss/86400.0) &
! / 1.201E-5, 4)
! / c_molar_mass, 4)
! vh ! expression below can be slightly different form that above in cases where
! leaf maintenance respiration is reduced in CASA
! (relative to its original value calculated in cable_canopy)
! in order to avoid negative carbon stores.

IF (output%casa) THEN
temp_acc = canopy%frday/1.201E-5 + &
temp_acc = canopy%frday/c_molar_mass + &
(casaflux%crmplant(:, 2)/86400.0 + casaflux%crmplant(:, 3)/86400.0 + &
casaflux%crgplant/86400.0 + casaflux%clabloss/86400.)/1.201E-5
casaflux%crgplant/86400.0 + casaflux%clabloss/86400.)/c_molar_mass
ELSE
temp_acc = (canopy%frp + canopy%frday)/1.201E-5
temp_acc = (canopy%frp + canopy%frday)/c_molar_mass
END IF
CALL generate_out_write_acc(ovid%AutoResp, 'AutoResp', out%AutoResp, temp_acc, ranges%AutoResp, patchout%AutoResp, out_settings)
IF (output%casa) THEN
! rootresp alt: REAL(0.3*casaflux%crmplant(:,2)/86400.0/ 1.201E-5, 4)
! rootresp alt: REAL(0.3*casaflux%crmplant(:,2)/86400.0/ c_molar_mass, 4)
Comment thread
SeanBryan51 marked this conversation as resolved.
CALL generate_out_write_acc(ovid%RootResp, 'RootResp', out%RootResp, &
REAL(casaflux%crmplant(:, 3)/86400.0/1.201E-5, 4), ranges%AutoResp, patchout%AutoResp, out_settings)
REAL(casaflux%crmplant(:, 3)/86400.0/c_molar_mass, 4), ranges%AutoResp, patchout%AutoResp, out_settings)
CALL generate_out_write_acc(ovid%StemResp, 'StemResp', out%StemResp, &
REAL(casaflux%crmplant(:, 2)/86400.0/1.201E-5, 4), ranges%AutoResp, patchout%AutoResp, out_settings)
REAL(casaflux%crmplant(:, 2)/86400.0/c_molar_mass, 4), ranges%AutoResp, patchout%AutoResp, out_settings)
END IF
END IF

IF (output%LeafResp) THEN
! LeafResp: Leaf respiration [umol/m^2/s]
CALL generate_out_write_acc(ovid%LeafResp, 'LeafResp', out%LeafResp, REAL(canopy%frday/1.201E-5, 4), ranges%LeafResp, patchout%LeafResp, out_settings)
CALL generate_out_write_acc(ovid%LeafResp, 'LeafResp', out%LeafResp, REAL(canopy%frday/c_molar_mass, 4), ranges%LeafResp, patchout%LeafResp, out_settings)
END IF
IF (output%HeteroResp) THEN
! HeteroResp: heterotrophic respiration [umol/m^2/s]
CALL generate_out_write_acc(ovid%HeteroResp, 'HeteroResp', out%HeteroResp, REAL(canopy%frs/1.201E-5, 4), ranges%HeteroResp, patchout%HeteroResp, out_settings)
CALL generate_out_write_acc(ovid%HeteroResp, 'HeteroResp', out%HeteroResp, REAL(canopy%frs/c_molar_mass, 4), ranges%HeteroResp, patchout%HeteroResp, out_settings)
END IF

! output patch area
Expand All @@ -2163,37 +2164,37 @@ SUBROUTINE write_output(dels, ktau, met, canopy, casaflux, casapool, casamet, ss
IF (cable_user%POPLUC) THEN
temp_acc = -(casaflux%Crsoil - casaflux%cnpp &
- casapool%dClabiledt)/86400.0 &
/1.201E-5 !- &
/c_molar_mass !- &
!REAL((casaflux%FluxCtohwp + casaflux%FluxCtoclear )/86400.0 &
!/ 1.201E-5, 4)
!/ c_molar_mass, 4)
ELSE
temp_acc = -(casaflux%Crsoil - casaflux%cnpp &
- casapool%dClabiledt)/86400.0 &
/1.201E-5
/c_molar_mass
END IF
CALL generate_out_write_acc(ovid%NBP, 'NBP', out%NBP, temp_acc, ranges%NEE, patchout%NBP, out_settings)
END IF

!------------------------WRITE REMAINING CASA DATA----------------------------------
CALL generate_out_write_acc(ovid%dCdt, 'dCdt', out%dCdt, &
REAL((casapool%ctot - casapool%ctot_0)/86400.0/1.201E-5, 4), ranges%NEE, patchout%dCdt, out_settings)
REAL((casapool%ctot - casapool%ctot_0)/86400.0/c_molar_mass, 4), ranges%NEE, patchout%dCdt, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnover, 'PlantTurnover', out%PlantTurnover, &
REAL((SUM(casaflux%Cplant_turnover, 2))/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnover, out_settings)
REAL((SUM(casaflux%Cplant_turnover, 2))/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnover, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverLeaf, 'PlantTurnoverLeaf', out%PlantTurnoverLeaf, &
REAL((casaflux%Cplant_turnover(:, 1))/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverLeaf, out_settings)
REAL((casaflux%Cplant_turnover(:, 1))/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverLeaf, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverFineRoot, 'PlantTurnoverFineRoot', out%PlantTurnoverFineRoot, &
REAL((casaflux%Cplant_turnover(:, 3))/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverFineRoot, out_settings)
REAL((casaflux%Cplant_turnover(:, 3))/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverFineRoot, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverWood, 'PlantTurnoverWood', out%PlantTurnoverWood, &
REAL((casaflux%Cplant_turnover(:, 2))/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverWood, out_settings)
REAL((casaflux%Cplant_turnover(:, 2))/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverWood, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverWoodDist, 'PlantTurnoverWoodDist', out%PlantTurnoverWoodDist, &
REAL(casaflux%Cplant_turnover_disturbance/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverWoodDist, out_settings)
REAL(casaflux%Cplant_turnover_disturbance/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverWoodDist, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverWoodCrowding, 'PlantTurnoverWoodCrowding', out%PlantTurnoverWoodCrowding, &
REAL(casaflux%Cplant_turnover_crowding/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverWoodCrowding, out_settings)
REAL(casaflux%Cplant_turnover_crowding/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverWoodCrowding, out_settings)
CALL generate_out_write_acc(ovid%PlantTurnoverWoodResourceLim, 'PlantTurnoverWoodResourceLim', out%PlantTurnoverWoodResourceLim, &
REAL((casaflux%Cplant_turnover_resource_limitation)/86400.0/1.201E-5, 4), ranges%NEE, patchout%PlantTurnoverWoodResourceLim, out_settings)
REAL((casaflux%Cplant_turnover_resource_limitation)/86400.0/c_molar_mass, 4), ranges%NEE, patchout%PlantTurnoverWoodResourceLim, out_settings)
IF (cable_user%POPLUC) THEN
CALL generate_out_write_acc(ovid%LandUseFlux, 'LandUseFlux', out%LandUseFlux, &
REAL((casaflux%FluxCtohwp + casaflux%FluxCtoclear)/86400.0/1.201E-5, 4), ranges%NEE, patchout%LandUseFlux, out_settings)
REAL((casaflux%FluxCtohwp + casaflux%FluxCtoclear)/86400.0/c_molar_mass, 4), ranges%NEE, patchout%LandUseFlux, out_settings)
END IF

! plant carbon [kg C m-2]
Expand Down
1 change: 1 addition & 0 deletions src/params/cable_phys_constants_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,6 +42,7 @@ MODULE cable_phys_constants_mod
REAL, PARAMETER :: cswat = 4.218e3 ! specific heat for water at 0°C (J/kg/K)
REAL, PARAMETER :: density_liq = 1000.0 ! density of liquid water
REAL, PARAMETER :: density_ice = 921.0 ! density of ice
REAL, PARAMETER :: c_molar_mass = 1.201e-5 ! molar mass of carbon (ug/mol)

! Teten coefficients
REAL, PARAMETER :: tetena = 6.106 ! Magnus Tetans (Murray 1967)
Expand Down