Computes spectral geopotential from spectral temperature T and spectral topography phis, as in GFDL Climate Group GCM. The result is updated in the state%phi variable.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(ModelState_t), | intent(inout) | :: | state | |||
| integer, | intent(in) | :: | time_level |
subroutine set_geopotential(state, time_level)
use physical_constants, only : rgas
use model_state, only : ModelState_t
type(ModelState_t), intent(inout) :: state
integer, intent(in) :: time_level
state%phi = get_geopotential(state%t(:, :, :, time_level), &
state%phis, state%xgeop1, state%xgeop2, &
state%mod_geometry%hsg, state%mod_geometry%fsg)
end subroutine set_geopotential