Skip to content

updated fv_regional_bc.F90 to read levsp using fms2_io primitives #152

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 1 commit into from
Oct 6, 2021
Merged

updated fv_regional_bc.F90 to read levsp using fms2_io primitives #152

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
34 changes: 12 additions & 22 deletions model/fv_regional_bc.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -40,7 +40,8 @@ module fv_regional_mod
register_global_attribute, read_data, &
register_field, FmsNetcdfFile_t, &
FmsNetcdfDomainFile_t, write_data, &
get_global_attribute, global_att_exists
get_global_attribute, global_att_exists,&
get_dimension_size
use tracer_manager_mod,only: get_tracer_index,get_tracer_names
use field_manager_mod, only: MODEL_ATMOS
use time_manager_mod, only: get_time &
Expand Down Expand Up @@ -1352,40 +1353,31 @@ subroutine start_regional_restart(Atm &
!*** Local variables
!---------------------
!
type(FmsNetcdfFile_t) :: Gfs_ctl
type(FmsNetcdfFile_t) :: Grid_input
integer, allocatable, dimension(:) :: pes !< Array of the pes in the current pelist
integer :: ierr, ios
real, allocatable :: wk2(:,:)
!
logical :: filtered_terrain = .true.
logical :: gfs_dwinds = .true.
integer :: levp = 64
logical :: checker_tr = .false.
integer :: nt_checker = 0
namelist /external_ic_nml/ filtered_terrain, levp, gfs_dwinds &
,checker_tr, nt_checker
! variables for reading the dimension from the gfs_ctrl
integer ncid, levsp
integer :: levp, levsp
!-----------------------------------------------------------------------
!***********************************************************************
!-----------------------------------------------------------------------
!
!-----------------------------------------------------------------------
!*** Read the number of model layers in the external forecast (=levp).
!-----------------------------------------------------------------------
!
read (input_nml_file,external_ic_nml,iostat=ios)
ierr = check_nml_error(ios,'external_ic_nml')
if(ierr/=0)then
write(0,11011)ierr
11011 format(' start_regional_restart failed to read external_ic_nml ierr=',i3)
allocate(pes(mpp_npes()))
call mpp_get_current_pelist(pes)
if( open_file(Gfs_ctl, 'INPUT/gfs_ctrl.nc', "read", pelist=pes) ) then
!--- read in the number of levsp
call get_dimension_size(Gfs_ctl, 'levsp', levsp)
call close_file(Gfs_ctl)
else
call mpp_error(FATAL,'==> Error in fv_regional::start_regional_restart file INPUT/gfs_ctl.nc does not exist')
endif

!--- read in ak and bk from the control file using fms_io read_data ---
call open_ncfile( 'INPUT/gfs_ctrl.nc', ncid ) ! open the file
call get_ncdim1( ncid, 'levsp', levsp )
call close_ncfile( ncid )

levp = levsp-1
!
!-----------------------------------------------------------------------
Expand All @@ -1409,8 +1401,6 @@ subroutine start_regional_restart(Atm &
allocate (wk2(levp+1,2))
allocate (ak_in(levp+1)) !<-- Save the input vertical structure for
allocate (bk_in(levp+1)) ! remapping BC updates during the forecast.
allocate(pes(mpp_npes()))
call mpp_get_current_pelist(pes)
if (Atm%flagstruct%hrrrv3_ic) then
if (open_file(Grid_input, 'INPUT/hrrr_ctrl.nc', "read", pelist=pes)) then
call read_data(Grid_input,'vcoord',wk2)
Expand Down