Skip to content
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

Split update_atmos_radiation_physics in SHiELD/atmos_model.F90 #30

Merged
merged 2 commits into from
Nov 1, 2023
Merged

Split update_atmos_radiation_physics in SHiELD/atmos_model.F90 #30

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
8d126d3
Split update_atmos_radiation_physics
spencerkclark Oct 17, 2023
68832a8
Address comments from Rusty
spencerkclark Nov 1, 2023
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
38 changes: 35 additions & 3 deletions SHiELD/atmos_model.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -111,6 +111,9 @@ module atmos_model_mod
public update_atmos_radiation_physics
public update_atmos_model_state
public update_atmos_model_dynamics
public update_atmos_pre_radiation
public update_atmos_radiation
public update_atmos_physics
public atmos_model_init, atmos_model_end, atmos_data_type
public atmos_model_restart
public Atm_block, IPD_Control, IPD_Data
Expand Down Expand Up @@ -221,8 +224,22 @@ module atmos_model_mod
! compute/exchange fluxes with other component models. All fields in this
! variable type are allocated for the global grid (without halo regions).
! </INOUT>
subroutine update_atmos_radiation_physics(Atmos)
bensonr marked this conversation as resolved.
Show resolved Hide resolved
!--- This subroutine, which merely combines calls to finer-grained subroutines,
! is only required for backwards compatibility with
! FMScoupler/SHiELD/coupler_main.F90.
type (atmos_data_type), intent(in) :: Atmos

call update_atmos_pre_radiation(Atmos)

if (.not. dycore_only) then
call update_atmos_radiation(Atmos)
call update_atmos_physics(Atmos)
end if

end subroutine update_atmos_radiation_physics

subroutine update_atmos_radiation_physics (Atmos)
subroutine update_atmos_pre_radiation (Atmos)
!-----------------------------------------------------------------------
type (atmos_data_type), intent(in) :: Atmos
!--- local variables---
Expand Down Expand Up @@ -278,6 +295,14 @@ subroutine update_atmos_radiation_physics (Atmos)
#endif

call mpp_clock_end(setupClock)
endif
end subroutine update_atmos_pre_radiation

subroutine update_atmos_radiation (Atmos)
!-----------------------------------------------------------------------
type (atmos_data_type), intent(in) :: Atmos
!--- local variables---
integer :: nb, jdat(8), rc

if (mpp_pe() == mpp_root_pe() .and. debug) write(6,*) "radiation driver"
!--- execute the IPD atmospheric radiation subcomponent (RRTM)
Expand All @@ -296,6 +321,14 @@ subroutine update_atmos_radiation_physics (Atmos)
call FV3GFS_IPD_checksum(IPD_Control, IPD_Data, Atm_block)
endif

end subroutine update_atmos_radiation

subroutine update_atmos_physics (Atmos)
!-----------------------------------------------------------------------
type (atmos_data_type), intent(in) :: Atmos
!--- local variables---
integer :: nb, jdat(8), rc

if (mpp_pe() == mpp_root_pe() .and. debug) write(6,*) "physics driver"
!--- execute the IPD atmospheric physics step1 subcomponent (main physics driver)
call mpp_clock_begin(physClock)
Expand Down Expand Up @@ -331,12 +364,11 @@ subroutine update_atmos_radiation_physics (Atmos)
endif
call getiauforcing(IPD_Control,IAU_data)
if (mpp_pe() == mpp_root_pe() .and. debug) write(6,*) "end of radiation and physics step"
endif

call mpp_clock_end(shieldClock)

!-----------------------------------------------------------------------
end subroutine update_atmos_radiation_physics
end subroutine update_atmos_physics
! </SUBROUTINE>


Expand Down