SUBROUTINE get_idata (ng)
!
!git $Id$
!svn $Id: get_idata.F 1151 2023-02-09 03:08:53Z arango $
!================================================== Hernan G. Arango ===
! Copyright (c) 2002-2023 The ROMS/TOMS Group !
! Licensed under a MIT/X style license !
! See License_ROMS.md !
!=======================================================================
! !
! This routine reads input data that needs to be obtained only once. !
! !
! Currently, this routine is only executed in serial mode by the !
! main thread. !
! !
!=======================================================================
!
USE mod_param
USE mod_parallel
USE mod_grid
USE mod_iounits
USE mod_ncparam
USE mod_scalars
USE mod_sources
USE mod_stepping
USE mod_tides
!
USE nf_fread3d_mod, ONLY : nf_fread3d
USE nf_fread4d_mod, ONLY : nf_fread4d
USE strings_mod, ONLY : FoundError
!
implicit none
!
! Imported variable declarations.
!
integer, intent(in) :: ng
!
! Local variable declarations.
!
logical, save :: recordless = .TRUE.
logical, dimension(3) :: update = &
& (/ .FALSE., .FALSE., .FALSE. /)
!
integer :: LBi, UBi, LBj, UBj
integer :: itrc, is
real(r8) :: time_save = 0.0_r8
!
character (len=*), parameter :: MyFile = &
& "ROMS/Nonlinear/get_idata.F"
!
SourceFile=MyFile
!
! Lower and upper bounds for tiled arrays.
!
LBi=LBOUND(GRID(ng)%h,DIM=1)
UBi=UBOUND(GRID(ng)%h,DIM=1)
LBj=LBOUND(GRID(ng)%h,DIM=2)
UBj=UBOUND(GRID(ng)%h,DIM=2)
!
!-----------------------------------------------------------------------
! Turn on input data time wall clock.
!-----------------------------------------------------------------------
!
CALL wclock_on (ng, iNLM, 3, 118, MyFile)
!
!-----------------------------------------------------------------------
! Tide period, amplitude, phase, and currents.
!-----------------------------------------------------------------------
!
! Tidal Period.
!
IF (LprocessTides(ng)) THEN
IF (iic(ng).eq.0) THEN
CALL get_ngfld (ng, iNLM, idTper, TIDE(ng)%ncid, &
& 1, TIDE(ng), recordless, update(1), &
& 1, MTC, 1, 1, 1, NTC(ng), 1, &
& TIDES(ng) % Tperiod)
IF (FoundError(exit_flag, NoError, 138, MyFile)) RETURN
END IF
END IF
!
! Tidal elevation amplitude and phase. In order to read data as a
! function of tidal period, we need to reset the model time variables
! temporarily.
!
IF (LprocessTides(ng)) THEN
IF (iic(ng).eq.0) THEN
time_save=time(ng)
time(ng)=8640000.0_r8
tdays(ng)=time(ng)*sec2day
!
CALL get_2dfld (ng, iNLM, idTzam, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % SSH_Tamp)
IF (FoundError(exit_flag, NoError, 165, MyFile)) RETURN
!
CALL get_2dfld (ng, iNLM, idTzph, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % SSH_Tphase)
IF (FoundError(exit_flag, NoError, 177, MyFile)) RETURN
!
time(ng)=time_save
tdays(ng)=time(ng)*sec2day
END IF
END IF
!
! Tidal currents angle, phase, major and minor ellipse axis.
!
IF (LprocessTides(ng)) THEN
IF (iic(ng).eq.0) THEN
time_save=time(ng)
time(ng)=8640000.0_r8
tdays(ng)=time(ng)*sec2day
!
CALL get_2dfld (ng, iNLM, idTvan, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % UV_Tangle)
IF (FoundError(exit_flag, NoError, 205, MyFile)) RETURN
!
CALL get_2dfld (ng, iNLM, idTvph, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % UV_Tphase)
IF (FoundError(exit_flag, NoError, 217, MyFile)) RETURN
!
CALL get_2dfld (ng, iNLM, idTvma, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % UV_Tmajor)
IF (FoundError(exit_flag, NoError, 229, MyFile)) RETURN
!
CALL get_2dfld (ng, iNLM, idTvmi, TIDE(ng)%ncid, &
& 1, TIDE(ng), update(1), &
& LBi, UBi, LBj, UBj, MTC, NTC(ng), &
& GRID(ng) % rmask, &
& TIDES(ng) % UV_Tminor)
IF (FoundError(exit_flag, NoError, 241, MyFile)) RETURN
!
time(ng)=time_save
tdays(ng)=time(ng)*sec2day
END IF
END IF
!
!-----------------------------------------------------------------------
! Read in point Sources/Sinks position, direction, special flag, and
! mass transport nondimensional shape profile. Point sources are at
! U- and V-points.
!-----------------------------------------------------------------------
!
IF ((iic(ng).eq.0).and. &
& (LuvSrc(ng).or.LwSrc(ng).or.ANY(LtracerSrc(:,ng)))) THEN
CALL get_ngfld (ng, iNLM, idRxpo, SSF(ng)%ncid, &
& 1, SSF(ng), recordless, update(1), &
& 1, Nsrc(ng), 1, 1, 1, Nsrc(ng), 1, &
& SOURCES(ng) % Xsrc)
IF (FoundError(exit_flag, NoError, 754, MyFile)) RETURN
CALL get_ngfld (ng, iNLM, idRepo, SSF(ng)%ncid, &
& 1, SSF(ng), recordless, update(1), &
& 1, Nsrc(ng), 1, 1, 1, Nsrc(ng), 1, &
& SOURCES(ng) % Ysrc)
IF (FoundError(exit_flag, NoError, 763, MyFile)) RETURN
CALL get_ngfld (ng, iNLM, idRdir, SSF(ng)%ncid, &
& 1, SSF(ng), recordless, update(1), &
& 1, Nsrc(ng), 1, 1, 1, Nsrc(ng), 1, &
& SOURCES(ng) % Dsrc)
IF (FoundError(exit_flag, NoError, 772, MyFile)) RETURN
CALL get_ngfld (ng, iNLM, idRvsh, SSF(ng)%ncid, &
& 1, SSF(ng), recordless, update(1), &
& 1, Nsrc(ng), N(ng), 1, 1, Nsrc(ng), N(ng), &
& SOURCES(ng) % Qshape)
IF (FoundError(exit_flag, NoError, 782, MyFile)) RETURN
DO is=1,Nsrc(ng)
SOURCES(ng)%Isrc(is)= &
& MAX(1,MIN(NINT(SOURCES(ng)%Xsrc(is)),Lm(ng)+1))
SOURCES(ng)%Jsrc(is)= &
& MAX(1,MIN(NINT(SOURCES(ng)%Ysrc(is)),Mm(ng)+1))
END DO
END IF
!
!-----------------------------------------------------------------------
! Turn off input data time wall clock.
!-----------------------------------------------------------------------
!
CALL wclock_off (ng, iNLM, 3, 824, MyFile)
!
RETURN
END SUBROUTINE get_idata