SUBROUTINE ana_m3clima (ng, tile, model) ! !! git $Id$ !! svn $Id: ana_m3clima.h 1151 2023-02-09 03:08:53Z arango $ !!====================================================================== !! Copyright (c) 2002-2023 The ROMS/TOMS Group ! !! Licensed under a MIT/X style license ! !! See License_ROMS.md ! !======================================================================= ! ! ! This routine sets analytical 3D momentum climatology fields. ! ! ! !======================================================================= ! USE mod_param USE mod_clima USE mod_ncparam ! ! Imported variable declarations. ! integer, intent(in) :: ng, tile, model ! ! Local variable declarations. ! character (len=*), parameter :: MyFile = & & __FILE__ ! #include "tile.h" ! CALL ana_m3clima_tile (ng, tile, model, & & LBi, UBi, LBj, UBj, & & IminS, ImaxS, JminS, JmaxS, & & CLIMA(ng) % uclm, & & CLIMA(ng) % vclm) ! ! Set analytical header file name used. ! #ifdef DISTRIBUTE IF (Lanafile) THEN #else IF (Lanafile.and.(tile.eq.0)) THEN #endif ANANAME(13)=MyFile END IF ! RETURN END SUBROUTINE ana_m3clima ! !*********************************************************************** SUBROUTINE ana_m3clima_tile (ng, tile, model, & & LBi, UBi, LBj, UBj, & & IminS, ImaxS, JminS, JmaxS, & & uclm, vclm) !*********************************************************************** ! USE mod_param USE mod_scalars ! USE exchange_3d_mod #ifdef DISTRIBUTE USE mp_exchange_mod, ONLY : mp_exchange3d #endif ! ! Imported variable declarations. ! integer, intent(in) :: ng, tile, model integer, intent(in) :: LBi, UBi, LBj, UBj integer, intent(in) :: IminS, ImaxS, JminS, JmaxS ! #ifdef ASSUMED_SHAPE real(r8), intent(out) :: uclm(LBi:,LBj:,:) real(r8), intent(out) :: vclm(LBi:,LBj:,:) #else real(r8), intent(out) :: uclm(LBi:UBi,LBj:UBj,N(ng)) real(r8), intent(out) :: vclm(LBi:UBi,LBj:UBj,N(ng)) #endif ! ! Local variable declarations. ! integer :: i, j, k #include "set_bounds.h" ! !----------------------------------------------------------------------- ! Set 3D momentum climatology. !----------------------------------------------------------------------- ! #if defined MY_APPLICATION DO k=1,N DO j=JstrT,JendT DO i=IstrP,IendT uclm(i,j,k)=??? END DO END DO DO j=JstrP,JendT DO i=IstrT,IendT vclm(i,j,k)=??? END DO END DO END DO #else ana_m3clima.h: No values provided for uclm and vclm. #endif ! ! Exchange boundary data. ! IF (EWperiodic(ng).or.NSperiodic(ng)) THEN CALL exchange_u3d_tile (ng, tile, & & LBi, UBi, LBj, UBj, 1, N(ng), & & uclm) CALL exchange_v3d_tile (ng, tile, & & LBi, UBi, LBj, UBj, 1, N(ng), & & vclm) END IF #ifdef DISTRIBUTE CALL mp_exchange3d (ng, tile, model, 2, & & LBi, UBi, LBj, UBj, 1, N(ng), & & NghostPoints, & & EWperiodic(ng), NSperiodic(ng), & & uclm, vclm) #endif ! RETURN END SUBROUTINE ana_m3clima_tile