SUBROUTINE ad_get_data (ng) ! !git $Id$ !svn $Id: ad_get_data.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 in forcing, climatology and other data from ! ! NetCDF files. If there is more than one time-record, data is ! ! loaded into global two-time record arrays. The interpolation ! ! is carried elsewhere. ! ! ! ! Currently, this routine is only executed in serial mode by the ! ! main thread. ! ! ! !======================================================================= ! USE mod_param USE mod_boundary USE mod_coupling USE mod_clima USE mod_forces USE mod_grid USE mod_iounits USE mod_mixing USE mod_ncparam USE mod_ocean USE mod_scalars USE mod_sources USE mod_stepping ! USE strings_mod, ONLY : FoundError ! implicit none ! ! Imported variable declarations. ! integer, intent(in) :: ng ! ! Local variable declarations. ! logical, save :: recordless = .FALSE. logical, dimension(3) :: update = & & (/ .FALSE., .FALSE., .FALSE. /) ! integer :: ILB, IUB, JLB, JUB integer :: LBi, UBi, LBj, UBj integer :: i, ic, my_tile ! real(r8) :: scale ! character (len=*), parameter :: MyFile = & & "ROMS/Adjoint/ad_get_data.F" ! ! Lower and upper bounds for nontiled (global values) boundary arrays. ! my_tile=-1 ! for global values ILB=BOUNDS(ng)%LBi(my_tile) IUB=BOUNDS(ng)%UBi(my_tile) JLB=BOUNDS(ng)%LBj(my_tile) JUB=BOUNDS(ng)%UBj(my_tile) ! ! 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, iADM, 3, 92, MyFile) ! !======================================================================= ! Read in forcing data from FORCING NetCDF file. !======================================================================= ! !----------------------------------------------------------------------- ! Point Sources/Sinks time dependent data. !----------------------------------------------------------------------- ! ! Point Source/Sink vertically integrated mass transport. ! IF (LuvSrc(ng).or.LwSrc(ng)) THEN CALL get_ngfldr (ng, iADM, idRtra, SSF(ng)%ncid, & & 1, SSF(ng), recordless, update(1), & & 1, Nsrc(ng), 1, 2, 1, Nsrc(ng), 1, & & SOURCES(ng) % QbarG) IF (FoundError(exit_flag, NoError, 115, MyFile)) RETURN END IF ! ! Tracer Sources/Sinks. ! DO i=1,NT(ng) IF (LtracerSrc(i,ng)) THEN CALL get_ngfldr (ng, iADM, idRtrc(i), SSF(ng)%ncid, & & 1, SSF(ng), recordless, update(1), & & 1, Nsrc(ng), N(ng), 2, 1, Nsrc(ng), N(ng), & & SOURCES(ng) % TsrcG(:,:,:,i)) IF (FoundError(exit_flag, NoError, 131, MyFile)) RETURN END IF END DO ! !----------------------------------------------------------------------- ! Surface wind stress components from NLM forward file. !----------------------------------------------------------------------- ! CALL get_2dfldr (ng, iADM, idUsms, BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % umask, & & FORCES(ng) % sustrG) IF (FoundError(exit_flag, NoError, 218, MyFile)) RETURN ! CALL get_2dfldr (ng, iADM, idVsms, BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % vmask, & & FORCES(ng) % svstrG) IF (FoundError(exit_flag, NoError, 230, MyFile)) RETURN ! !----------------------------------------------------------------------- ! Surface air pressure from NLM forward file. !----------------------------------------------------------------------- ! CALL get_2dfldr (ng, iADM, idPair, BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & FORCES(ng) % PairG) IF (FoundError(exit_flag, NoError, 271, MyFile)) RETURN ! !----------------------------------------------------------------------- ! Surface solar shortwave radiation from NLM forward file. !----------------------------------------------------------------------- ! CALL get_2dfldr (ng, iADM, idSrad, BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & FORCES(ng) % srflxG) IF (FoundError(exit_flag, NoError, 439, MyFile)) RETURN ! !----------------------------------------------------------------------- ! Surface net heat flux from NLM forward file. !----------------------------------------------------------------------- ! CALL get_2dfldr (ng, iADM, idTsur(itemp), BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & FORCES(ng) % stfluxG(:,:,:,itemp)) IF (FoundError(exit_flag, NoError, 604, MyFile)) RETURN ! !----------------------------------------------------------------------- ! Surface net freshwater flux (E-P) from NLM forward file. !----------------------------------------------------------------------- ! CALL get_2dfldr (ng, iADM, idEmPf, BLK(ng)%ncid, & & 1, BLK(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & FORCES(ng) % stfluxG(:,:,:,isalt)) IF (FoundError(exit_flag, NoError, 705, MyFile)) RETURN ! !======================================================================= ! Read in open boundary conditions from BOUNDARY NetCDF file. !======================================================================= ! IF (LprocessOBC(ng)) THEN IF (ad_LBC(iwest,isFsur,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idZbry(iwest), & & BRYncid(idZbry(iwest),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 0, Mm(ng)+1, 1, & & BOUNDARY(ng) % zetaG_west) IF (FoundError(exit_flag, NoError, 810, MyFile)) RETURN END IF ! IF (ad_LBC(ieast,isFsur,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idZbry(ieast), & & BRYncid(idZbry(ieast),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 0, Mm(ng)+1, 1, & & BOUNDARY(ng) % zetaG_east) IF (FoundError(exit_flag, NoError, 823, MyFile)) RETURN END IF ! IF (ad_LBC(isouth,isFsur,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idZbry(isouth), & & BRYncid(idZbry(isouth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 0, Lm(ng)+1, 1, & & BOUNDARY(ng) % zetaG_south) IF (FoundError(exit_flag, NoError, 836, MyFile)) RETURN END IF ! IF (ad_LBC(inorth,isFsur,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idZbry(inorth), & & BRYncid(idZbry(inorth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 0, Lm(ng)+1, 1, & & BOUNDARY(ng) % zetaG_north) IF (FoundError(exit_flag, NoError, 849, MyFile)) RETURN END IF END IF ! IF (LprocessOBC(ng)) THEN IF (ad_LBC(iwest,isUbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU2bc(iwest), & & BRYncid(idU2bc(iwest),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 0, Mm(ng)+1, 1, & & BOUNDARY(ng) % ubarG_west) IF (FoundError(exit_flag, NoError, 867, MyFile)) RETURN END IF ! IF (ad_LBC(iwest,isVbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV2bc(iwest), & & BRYncid(idV2bc(iwest),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 1, Mm(ng)+1, 1, & & BOUNDARY(ng) % vbarG_west) IF (FoundError(exit_flag, NoError, 880, MyFile)) RETURN END IF ! IF (ad_LBC(ieast,isUbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU2bc(ieast), & & BRYncid(idU2bc(ieast),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 0, Mm(ng)+1, 1, & & BOUNDARY(ng) % ubarG_east) IF (FoundError(exit_flag, NoError, 893, MyFile)) RETURN END IF ! IF (ad_LBC(ieast,isVbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV2bc(ieast), & & BRYncid(idV2bc(ieast),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, 1, 2, 1, Mm(ng)+1, 1, & & BOUNDARY(ng) % vbarG_east) IF (FoundError(exit_flag, NoError, 906, MyFile)) RETURN END IF ! IF (ad_LBC(isouth,isUbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU2bc(isouth), & & BRYncid(idU2bc(isouth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 1, Lm(ng)+1, 1, & & BOUNDARY(ng) % ubarG_south) IF (FoundError(exit_flag, NoError, 919, MyFile)) RETURN END IF ! IF (ad_LBC(isouth,isVbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV2bc(isouth), & & BRYncid(idV2bc(isouth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 0, Lm(ng)+1, 1, & & BOUNDARY(ng) % vbarG_south) IF (FoundError(exit_flag, NoError, 932, MyFile)) RETURN END IF ! IF (ad_LBC(inorth,isUbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU2bc(inorth), & & BRYncid(idU2bc(inorth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 1, Lm(ng)+1, 1, & & BOUNDARY(ng) % ubarG_north) IF (FoundError(exit_flag, NoError, 945, MyFile)) RETURN END IF ! IF (ad_LBC(inorth,isVbar,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV2bc(inorth), & & BRYncid(idV2bc(inorth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, 1, 2, 0, Lm(ng)+1, 1, & & BOUNDARY(ng) % vbarG_north) IF (FoundError(exit_flag, NoError, 958, MyFile)) RETURN END IF END IF ! IF (LprocessOBC(ng)) THEN IF (ad_LBC(iwest,isUvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU3bc(iwest), & & BRYncid(idU3bc(iwest),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 0, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % uG_west) IF (FoundError(exit_flag, NoError, 977, MyFile)) RETURN END IF ! IF (ad_LBC(iwest,isVvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV3bc(iwest), & & BRYncid(idV3bc(iwest),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 1, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % vG_west) IF (FoundError(exit_flag, NoError, 990, MyFile)) RETURN END IF ! IF (ad_LBC(ieast,isUvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU3bc(ieast), & & BRYncid(idU3bc(ieast),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 0, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % uG_east) IF (FoundError(exit_flag, NoError, 1003, MyFile)) RETURN END IF ! IF (ad_LBC(ieast,isVvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV3bc(ieast), & & BRYncid(idV3bc(ieast),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 1, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % vG_east) IF (FoundError(exit_flag, NoError, 1016, MyFile)) RETURN END IF ! IF (ad_LBC(isouth,isUvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU3bc(isouth), & & BRYncid(idU3bc(isouth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 1, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % uG_south) IF (FoundError(exit_flag, NoError, 1029, MyFile)) RETURN END IF ! IF (ad_LBC(isouth,isVvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV3bc(isouth), & & BRYncid(idV3bc(isouth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 0, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % vG_south) IF (FoundError(exit_flag, NoError, 1042, MyFile)) RETURN END IF ! IF (ad_LBC(inorth,isUvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idU3bc(inorth), & & BRYncid(idU3bc(inorth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 1, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % uG_north) IF (FoundError(exit_flag, NoError, 1055, MyFile)) RETURN END IF ! IF (ad_LBC(inorth,isVvel,ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idV3bc(inorth), & & BRYncid(idV3bc(inorth),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 0, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % vG_north) IF (FoundError(exit_flag, NoError, 1068, MyFile)) RETURN END IF END IF ! IF (LprocessOBC(ng)) THEN DO i=1,NT(ng) IF (ad_LBC(iwest,isTvar(i),ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idTbry(iwest,i), & & BRYncid(idTbry(iwest,i),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 0, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % tG_west(:,:,:,i)) IF (FoundError(exit_flag, NoError, 1087, MyFile)) RETURN END IF END DO ! DO i=1,NT(ng) IF (ad_LBC(ieast,isTvar(i),ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idTbry(ieast,i), & & BRYncid(idTbry(ieast,i),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & JLB, JUB, N(ng), 2, 0, Mm(ng)+1, N(ng), & & BOUNDARY(ng) % tG_east(:,:,:,i)) IF (FoundError(exit_flag, NoError, 1102, MyFile)) RETURN END IF END DO ! DO i=1,NT(ng) IF (ad_LBC(isouth,isTvar(i),ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idTbry(isouth,i), & & BRYncid(idTbry(isouth,i),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 0, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % tG_south(:,:,:,i)) IF (FoundError(exit_flag, NoError, 1117, MyFile)) RETURN END IF END DO ! DO i=1,NT(ng) IF (ad_LBC(inorth,isTvar(i),ng)%acquire) THEN CALL get_ngfldr (ng, iADM, idTbry(inorth,i), & & BRYncid(idTbry(inorth,i),ng), & & nBCfiles(ng), BRY(1,ng), & & recordless, update(1), & & ILB, IUB, N(ng), 2, 0, Lm(ng)+1, N(ng), & & BOUNDARY(ng) % tG_north(:,:,:,i)) IF (FoundError(exit_flag, NoError, 1132, MyFile)) RETURN END IF END DO END IF ! !======================================================================= ! Read in data from Climatology NetCDF file. !======================================================================= ! ! Free-surface. ! IF (LsshCLM(ng)) THEN CALL get_2dfldr (ng, iADM, idSSHc, CLMncid(idSSHc,ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & CLIMA(ng) % sshG) IF (FoundError(exit_flag, NoError, 1158, MyFile)) RETURN END IF ! ! 2D momentum. ! IF (Lm2CLM(ng)) THEN CALL get_2dfldr (ng, iADM, idUbcl, CLMncid(idUbcl,ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % umask, & & CLIMA(ng) % ubarclmG) IF (FoundError(exit_flag, NoError, 1176, MyFile)) RETURN ! CALL get_2dfldr (ng, iADM, idVbcl, CLMncid(idVbcl,ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % vmask, & & CLIMA(ng) % vbarclmG) IF (FoundError(exit_flag, NoError, 1188, MyFile)) RETURN END IF ! ! 3D momentum. ! IF (Lm3CLM(ng)) THEN CALL get_3dfldr (ng, iADM, idUclm, CLMncid(idUclm,ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % umask, & & CLIMA(ng) % uclmG) IF (FoundError(exit_flag, NoError, 1207, MyFile)) RETURN ! CALL get_3dfldr (ng, iADM, idVclm, CLMncid(idVclm,ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % vmask, & & CLIMA(ng) % vclmG) IF (FoundError(exit_flag, NoError, 1219, MyFile)) RETURN END IF ! ! Tracers. ! ic=0 DO i=1,NT(ng) IF (LtracerCLM(i,ng)) THEN ic=ic+1 CALL get_3dfldr (ng, iADM, idTclm(i), & & CLMncid(idTclm(i),ng), & & nCLMfiles(ng), CLM(1,ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % rmask, & & CLIMA(ng) % tclmG(:,:,:,:,ic)) IF (FoundError(exit_flag, NoError, 1241, MyFile)) RETURN END IF END DO ! !----------------------------------------------------------------------- ! Read in forward state solution. !----------------------------------------------------------------------- ! ! Read in free-surface. ! CALL get_2dfldr (ng, iADM, idFsur, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % rmask, & & OCEAN(ng) % zetaG) IF (FoundError(exit_flag, NoError, 1265, MyFile)) RETURN ! ! Read 2D momentum. ! CALL get_2dfldr (ng, iADM, idUbar, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % umask, & & OCEAN(ng) % ubarG) IF (FoundError(exit_flag, NoError, 1279, MyFile)) RETURN CALL get_2dfldr (ng, iADM, idVbar, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % vmask, & & OCEAN(ng) % vbarG) IF (FoundError(exit_flag, NoError, 1291, MyFile)) RETURN ! ! Read in variables associated with time-averaged 2D momentum terms. ! CALL get_2dfldr (ng, iADM, idUfx1, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % umask, & & COUPLING(ng) % DU_avg1G) IF (FoundError(exit_flag, NoError, 1348, MyFile)) RETURN ! CALL get_2dfldr (ng, iADM, idUfx2, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % umask, & & COUPLING(ng) % DU_avg2G) IF (FoundError(exit_flag, NoError, 1360, MyFile)) RETURN ! CALL get_2dfldr (ng, iADM, idVfx1, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % vmask, & & COUPLING(ng) % DV_avg1G) IF (FoundError(exit_flag, NoError, 1372, MyFile)) RETURN ! CALL get_2dfldr (ng, iADM, idVfx2, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 2, 1, & & GRID(ng) % vmask, & & COUPLING(ng) % DV_avg2G) IF (FoundError(exit_flag, NoError, 1384, MyFile)) RETURN ! ! Read in 3D momentum. ! CALL get_3dfldr (ng, iADM, idUvel, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % umask, & & OCEAN(ng) % uG) IF (FoundError(exit_flag, NoError, 1398, MyFile)) RETURN ! CALL get_3dfldr (ng, iADM, idVvel, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % vmask, & & OCEAN(ng) % vG) IF (FoundError(exit_flag, NoError, 1410, MyFile)) RETURN ! ! Read in 3D tracers. ! DO i=1,NT(ng) CALL get_3dfldr (ng, iADM, idTvar(i), FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, & & GRID(ng) % rmask, & & OCEAN(ng) % tG(:,:,:,:,i)) IF (FoundError(exit_flag, NoError, 1478, MyFile)) RETURN END DO ! ! Read in vertical mixing variables. ! DO i=1,NAT scale=Fscale(idDiff(i),ng) ! save and rescale Fscale(idDiff(i),ng)=ad_Akt_fac(i,ng) CALL get_3dfldr (ng, iADM, idDiff(i), FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 0, N(ng), 2, 1, & & GRID(ng) % rmask, & & MIXING(ng) % AktG(:,:,:,:,i)) Fscale(idDiff(i),ng)=scale IF (FoundError(exit_flag, NoError, 1499, MyFile)) RETURN END DO ! scale=Fscale(idVvis,ng) ! save and rescale Fscale(idVvis,ng)=ad_Akv_fac(ng) CALL get_3dfldr (ng, iADM, idVvis, FWD(ng)%ncid, & & 1, FWD(ng), update(1), & & LBi, UBi, LBj, UBj, 0, N(ng), 2, 1, & & GRID(ng) % rmask, & & MIXING(ng) % AkvG) Fscale(idVvis,ng)=scale IF (FoundError(exit_flag, NoError, 1515, MyFile)) RETURN ! !----------------------------------------------------------------------- ! Turn off input data time wall clock. !----------------------------------------------------------------------- ! CALL wclock_off (ng, iADM, 3, 1783, MyFile) ! RETURN END SUBROUTINE ad_get_data