Skip to content

Commit

Permalink
Bring in updates from SHiELD (#19)
Browse files Browse the repository at this point in the history 
* Copy code changes to atmos_model.F90.

* whitespace removal

---------

Co-authored-by: Linjiong Zhou <Linjiong.Zhou@noaa.gov>
  • Loading branch information
@laurenchilutti @linjiongzhou
laurenchilutti and linjiongzhou authored Feb 23, 2023
1 parent 6195f02 commit 7b8deb1
Showing 1 changed file with 60 additions and 13 deletions.
73 changes: 60 additions & 13 deletions SHiELD/atmos_model.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ module atmos_model_mod
use IPD_driver, only: IPD_initialize, IPD_setup_step, &
IPD_radiation_step, &
IPD_physics_step1, &
IPD_physics_step2
IPD_physics_step2, IPD_physics_end
#ifdef STOCHY
use stochastic_physics, only: init_stochastic_physics, &
run_stochastic_physics
Expand All @@ -97,6 +97,12 @@ module atmos_model_mod
use FV3GFS_io_mod, only: send_diag_manager_controlled_diagnostic_data
use fv_iau_mod, only: iau_external_data_type,getiauforcing,iau_initialize
use module_ocean, only: ocean_init
#if defined (USE_COSP)
use cosp2_test, only: cosp2_driver
#endif
#if defined (COSP_OFFLINE)
use cosp2_test, only: cosp2_offline
#endif
!-----------------------------------------------------------------------

implicit none
Expand Down Expand Up @@ -144,7 +150,7 @@ module atmos_model_mod
end type atmos_data_type
!</PUBLICTYPE >

integer :: fv3Clock, getClock, overrideClock, updClock, setupClock, radClock, physClock
integer :: fv3Clock, getClock, overrideClock, setupClock, radClock, physClock, diagClock, shieldClock

!-----------------------------------------------------------------------
integer :: blocksize = 1
Expand All @@ -168,6 +174,7 @@ module atmos_model_mod
!----------------
! IPD containers
!----------------
type(IPD_init_type) :: Init_parm
type(IPD_control_type) :: IPD_Control
type(IPD_data_type), allocatable :: IPD_Data(:) ! number of blocks
type(IPD_diag_type) :: IPD_Diag(250)
Expand Down Expand Up @@ -221,6 +228,8 @@ subroutine update_atmos_radiation_physics (Atmos)
integer :: nb, jdat(8)
integer :: nthrds

call mpp_clock_begin(shieldClock)

#ifdef OPENMP
nthrds = omp_get_max_threads()
#else
Expand Down Expand Up @@ -323,6 +332,8 @@ subroutine update_atmos_radiation_physics (Atmos)
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
! </SUBROUTINE>
Expand Down Expand Up @@ -360,7 +371,6 @@ subroutine atmos_model_init (Atmos, Time_init, Time, Time_step, iau_offset)
character(len=132) :: text
logical :: p_hydro, hydro, fexist
logical, save :: block_message = .true.
type(IPD_init_type) :: Init_parm
integer :: bdat(8), cdat(8)
integer :: ntracers
integer :: kdt_prev
Expand Down Expand Up @@ -557,17 +567,18 @@ subroutine atmos_model_init (Atmos, Time_init, Time, Time_step, iau_offset)
endif
#endif

setupClock = mpp_clock_id( 'GFS Step Setup ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
overrideClock = mpp_clock_id( 'GFS Override ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
radClock = mpp_clock_id( 'GFS Radiation ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
physClock = mpp_clock_id( 'GFS Physics ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
getClock = mpp_clock_id( 'Dynamics get state ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
updClock = mpp_clock_id( 'Dynamics update state ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
setupClock = mpp_clock_id( '---GFS Step Setup ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
overrideClock = mpp_clock_id( '---GFS Override ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
radClock = mpp_clock_id( '---GFS Radiation ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
physClock = mpp_clock_id( '---GFS Physics ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
diagClock = mpp_clock_id( '---GFS Diag ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
getClock = mpp_clock_id( '---GFS Get State ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
if (sync) then
fv3Clock = mpp_clock_id( 'FV3 Dycore ', flags=clock_flag_default+MPP_CLOCK_SYNC, grain=CLOCK_COMPONENT )
fv3Clock = mpp_clock_id( '--FV3 Dycore ', flags=clock_flag_default+MPP_CLOCK_SYNC, grain=CLOCK_COMPONENT )
else
fv3Clock = mpp_clock_id( 'FV3 Dycore ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
fv3Clock = mpp_clock_id( '--FV3 Dycore ', flags=clock_flag_default, grain=CLOCK_COMPONENT )
endif
shieldClock= mpp_clock_id( '--SHiELD Physics ', flags=clock_flag_default, grain=CLOCK_COMPONENT )

!-----------------------------------------------------------------------
end subroutine atmos_model_init
Expand Down Expand Up @@ -605,11 +616,12 @@ subroutine update_atmos_model_state (Atmos)

call set_atmosphere_pelist()
call mpp_clock_begin(fv3Clock)
call mpp_clock_begin(updClock)
call atmosphere_state_update (Atmos%Time, IPD_Data, IAU_Data, Atm_block)
call mpp_clock_end(updClock)
call mpp_clock_end(fv3Clock)

call mpp_clock_begin(shieldClock)
call mpp_clock_begin(diagClock)

if (chksum_debug) then
if (mpp_pe() == mpp_root_pe()) print *,'UPDATE STATE ', IPD_Control%kdt, IPD_Control%fhour
call FV3GFS_IPD_checksum(IPD_Control, IPD_Data, Atm_block)
Expand Down Expand Up @@ -650,6 +662,36 @@ subroutine update_atmos_model_state (Atmos)
time_intfull = time_intfull - Atmos%iau_offset*3600.
endif
endif
#if defined (USE_COSP)
!-----------------------------------------------------------------------
! The CFMIP Observation Simulator Package (COSP)
! Added by Linjiong Zhou
! May 2021
!-----------------------------------------------------------------------

if (IPD_Control%do_cosp) then

call cosp2_driver (IPD_Control, IPD_Data(:)%Grid, IPD_Data(:)%Statein, &
IPD_Data(:)%Stateout, IPD_Data(:)%Sfcprop, &
IPD_Data(:)%Radtend, IPD_Data(:)%Intdiag, Init_parm)

endif
#endif
#if defined (COSP_OFFLINE)
!-----------------------------------------------------------------------
! Output Variables for the Offline CFMIP Observation Simulator Package (COSP)
! Added by Linjiong Zhou
! Nov 2022
!-----------------------------------------------------------------------

if (IPD_Control%do_cosp) then

call cosp2_offline (IPD_Control, IPD_Data(:)%Statein, &
IPD_Data(:)%Stateout, IPD_Data(:)%Sfcprop, &
IPD_Data(:)%Radtend, IPD_Data(:)%Intdiag, Init_parm)

endif
#endif
call gfdl_diag_output(Atmos%Time, Atm_block, IPD_Data, IPD_Control%nx, IPD_Control%ny, fprint, &
IPD_Control%levs, 1, 1, 1.d0, time_int, time_intfull, &
IPD_Control%fhswr, IPD_Control%fhlwr, &
Expand All @@ -661,6 +703,9 @@ subroutine update_atmos_model_state (Atmos)
if (mod(isec,nint(3600*IPD_Control%fhzero)) == 0) diag_time = Atmos%Time
endif

call mpp_clock_end(diagClock)
call mpp_clock_end(shieldClock)

end subroutine update_atmos_model_state
! </SUBROUTINE>

Expand Down Expand Up @@ -692,6 +737,8 @@ subroutine atmos_model_end (Atmos)
!---local variables
integer :: idx

call IPD_physics_end (IPD_Control)

!-----------------------------------------------------------------------
!---- termination routine for atmospheric model ----

Expand Down

0 comments on commit 7b8deb1

Please sign in to comment.