Skip to content

Workaround to fix issues with wrong namelists for nests in fv3atm #86

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
laurenchilutti merged 4 commits into from
Apr 9, 2021
Merged

Workaround to fix issues with wrong namelists for nests in fv3atm #86

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
306ff31
Add an option to use zero-gradient BC to reconstruct interface u/v fo…
XiaqiongZhou-NOAA Jan 25, 2021
4a27a04
Workaround to read the name of the CCPP suite from its current locati…
climbfuji Mar 26, 2021
97ccbe7
Merge branch 'dev/emc' of https://github.com/NOAA-GFDL/GFDL_atmos_cub…
climbfuji Mar 29, 2021
129a3ae
Merge branch 'dev/emc' of https://github.com/NOAA-GFDL/GFDL_atmos_cub…
climbfuji Mar 31, 2021
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
37 changes: 37 additions & 0 deletions driver/fvGFS/atmosphere.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -296,6 +296,24 @@ subroutine atmosphere_init (Time_init, Time, Time_step, Grid_box, area)
integer :: nthreads, ierr
integer :: nlunit = 9999
character (len = 64) :: fn_nml = 'input.nml'
! DH* 20210326
! This is a temporary workaround until the implementation
! of the generic interface to call GFDL or CCPP physics is
! completed. We need the name of the CCPP suite here in order
! to run the adiabatic init with fast physics turned on. All
! other vaiables are ignored (set to the same default value
! as in fv3atm's atmos_model.F90.
integer :: io
integer :: blocksize = 1
logical :: chksum_debug = .false.
logical :: dycore_only = .false.
logical :: debug = .false.
logical :: sync = .false.
integer, parameter :: maxhr = 4096
real, dimension(maxhr) :: fdiag = 0.
real :: fhmax=384.0, fhmaxhf=120.0, fhout=3.0, fhouthf=1.0,avg_max_length=3600.
namelist /atmos_model_nml/ blocksize, chksum_debug, dycore_only, debug, sync, fdiag, fhmax, fhmaxhf, fhout, fhouthf, ccpp_suite, avg_max_length
! *DH 20210326

!For regional
a_step = 0
Expand Down Expand Up @@ -450,6 +468,25 @@ subroutine atmosphere_init (Time_init, Time, Time_step, Grid_box, area)

! Do CCPP fast physics initialization before call to adiabatic_init (since this calls fv_dynamics)

! DH* 20210326
! First, read atmos_model_nml namelist section - this is a workaround to avoid
! unnecessary additional changes to the input namelists, in anticipation of the
! implementation of a generic interface for GFDL and CCPP fast physics soon
#ifdef INTERNAL_FILE_NML
read(input_nml_file, nml=atmos_model_nml, iostat=io)
ierr = check_nml_error(io, 'atmos_model_nml')
#else
unit = open_namelist_file ( )
ierr=1
do while (ierr /= 0)
read (unit, nml=atmos_model_nml, iostat=io, end=10)
ierr = check_nml_error(io,'atmos_model_nml')
enddo
10 call close_file (unit)
#endif
!write(0,'(a)') "It's me, and my physics suite is '" // trim(ccpp_suite) // "'"
! *DH 20210326

! For fast physics running over the entire domain, block
! and thread number are not used; set to safe values
cdata%blk_no = 1
Expand Down