Skip to content

Commit

Permalink
change mole_dens subroutines to molar_conc
Browse files Browse the repository at this point in the history
  • Loading branch information
jcurtis2 committed Oct 23, 2024
1 parent 3463e7d commit 1c8f222
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 5 deletions.
8 changes: 4 additions & 4 deletions src/gas_state.F90
Original file line number Diff line number Diff line change
Expand Up @@ -177,7 +177,7 @@ end subroutine gas_state_ensure_nonnegative
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!> Convert (mol m^{-3}) to (ppb).
subroutine gas_state_mole_dens_to_ppb(gas_state, env_state)
subroutine gas_state_molar_conc_to_ppb(gas_state, env_state)

!> Gas state.
type(gas_state_t), intent(inout) :: gas_state
Expand All @@ -186,12 +186,12 @@ subroutine gas_state_mole_dens_to_ppb(gas_state, env_state)

call gas_state_scale(gas_state, 1d9 / env_state_air_molar_den(env_state))

end subroutine gas_state_mole_dens_to_ppb
end subroutine gas_state_molar_conc_to_ppb

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!> Convert (ppb) to (mol m^{-3}).
subroutine gas_state_ppb_to_mole_dens(gas_state, env_state)
subroutine gas_state_ppb_to_molar_conc(gas_state, env_state)

!> Gas state.
type(gas_state_t), intent(inout) :: gas_state
Expand All @@ -200,7 +200,7 @@ subroutine gas_state_ppb_to_mole_dens(gas_state, env_state)

call gas_state_scale(gas_state, env_state_air_molar_den(env_state) / 1e9)

end subroutine gas_state_ppb_to_mole_dens
end subroutine gas_state_ppb_to_molar_conc

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Expand Down
2 changes: 1 addition & 1 deletion src/scenario.F90
Original file line number Diff line number Diff line change
Expand Up @@ -228,7 +228,7 @@ subroutine scenario_update_gas_state(scenario, delta_t, env_state, &
call gas_state_interp_1d(scenario%gas_emission, &
scenario%gas_emission_time, scenario%gas_emission_rate_scale, &
env_state%elapsed_time, emissions, emission_rate_scale)
call gas_state_mole_dens_to_ppb(emissions, env_state)
call gas_state_molar_conc_to_ppb(emissions, env_state)
p = emission_rate_scale * delta_t / env_state%height
call gas_state_add_scaled(gas_state, emissions, p)
end if
Expand Down

0 comments on commit 1c8f222

Please sign in to comment.