Skip to content

modify full_ic_frq mod for gfortran compiler and start j loop from 1 … #128

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
cenlinhe merged 1 commit into from
Nov 14, 2023
Merged

modify full_ic_frq mod for gfortran compiler and start j loop from 1 … #128

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
13 changes: 11 additions & 2 deletions hrldas/HRLDAS_forcing/create_forcing.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,6 +91,7 @@ program create_forcing
integer :: ihour
integer :: numarg
integer :: unmatched
logical :: ic_flag
#ifndef _GFORTRAN_
! For some reason, gfortran does not like this declared external,
! but everyone else wants it declared external. So we use #ifndef
Expand Down Expand Up @@ -441,8 +442,16 @@ program create_forcing
ierr = nf90_close(ncid)
call error_handler(ierr, "Problem closing Netcdf file")

ic_flag = .false.
if ( full_ic_frq > -1 ) then ! full_ic_frq == -1 turns off extra processing for initial conditions.
if ( ( (full_ic_frq>0) .and. (mod(ihour, full_ic_frq)==0) ) .or. ( ihour==0 ) ) then
if ( ihour == 0) ic_flag = .true.
if ( full_ic_frq > 0 ) then
if (mod(ihour, full_ic_frq) == 0) then
ic_flag = .true.
endif
endif

if ( ic_flag .eqv. .true. ) then

call open_netcdf_for_output( trim(OutputDir)//"HRLDAS_setup_"//date(1:4)//date(6:7)//date(9:10)//date(12:13)//"_d"//hgrid, &
ncid, version, geo_em%idim, geo_em%jdim, geo_em, .true.)
Expand Down Expand Up @@ -2107,7 +2116,7 @@ subroutine read_era5_netcdf_unit(ncid, datastruct, ierr, current)
array = [string('Tair_f_inst'),string('z'),string('lsm')] ! Only z (for geopotential) and lsm (for land-sea mask) are needed
nullify(datastruct%data)
! Match the netcdf attribute with one of the Vtable entries
VARSEARCH : do j = 0, vtable_count
VARSEARCH : do j = 1, vtable_count
if (current%label == vtable(j)%name) then
VAR_NAME = array(j)%str
datastruct%field = vtable(j)%name
Expand Down