New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14789 for NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/OCE/LDF/ldfc1d_c2d.F90 – NEMO

Ignore:
Timestamp:
2021-05-05T13:18:04+02:00 (3 years ago)
Author:
mcastril
Message:

[2021/HPC-11_mcastril_HPDAonline_DiagGPU] Update externals

Location:
NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
         5^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8^/vendors/PPR@HEAD            ext/PPR 
        89 
        910# SETTE 
        10 ^/utils/CI/sette@13559        sette 
         11^/utils/CI/sette@14244        sette 
  • NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/OCE/LDF/ldfc1d_c2d.F90

    r13497 r14789  
    22   !!====================================================================== 
    33   !!                    ***  MODULE  ldfc1d_c2d  *** 
    4    !! Ocean physics:  profile and horizontal shape of lateral eddy coefficients  
     4   !! Ocean physics:  profile and horizontal shape of lateral eddy coefficients 
    55   !!===================================================================== 
    66   !! History :  3.7  ! 2013-12  (G. Madec)  restructuration/simplification of aht/aeiv specification, 
     
    99 
    1010   !!---------------------------------------------------------------------- 
    11    !!   ldf_c1d       : ah reduced by 1/4 on the vertical (tanh profile, inflection at 300m)  
     11   !!   ldf_c1d       : ah reduced by 1/4 on the vertical (tanh profile, inflection at 300m) 
    1212   !!   ldf_c2d       : ah = F(e1,e2) (laplacian or = F(e1^3,e2^3) (bilaplacian) 
    1313   !!---------------------------------------------------------------------- 
     
    2929   REAL(wp) ::   r1_4  = 0.25_wp          ! =1/4 
    3030   REAL(wp) ::   r1_12 = 1._wp / 12._wp   ! =1/12 
    31   
     31 
    3232   !! * Substitutions 
    3333#  include "do_loop_substitute.h90" 
     
    4242      !!---------------------------------------------------------------------- 
    4343      !!                  ***  ROUTINE ldf_c1d  *** 
    44       !!               
     44      !! 
    4545      !! ** Purpose :   1D eddy diffusivity/viscosity coefficients 
    4646      !! 
    4747      !! ** Method  :   1D eddy diffusivity coefficients F( depth ) 
    48       !!                Reduction by zratio from surface to bottom  
    49       !!                hyperbolic tangent profile with inflection point  
     48      !!                Reduction by zratio from surface to bottom 
     49      !!                hyperbolic tangent profile with inflection point 
    5050      !!                at zh=500m and a width of zw=200m 
    5151      !! 
     
    9595         END_3D 
    9696         ! Lateral boundary conditions 
    97          CALL lbc_lnk_multi( 'ldfc1d_c2d', pah1, 'U', 1.0_wp , pah2, 'V', 1.0_wp )    
     97         CALL lbc_lnk( 'ldfc1d_c2d', pah1, 'U', 1.0_wp , pah2, 'V', 1.0_wp ) 
    9898         ! 
    9999      CASE DEFAULT                        ! error 
     
    107107      !!---------------------------------------------------------------------- 
    108108      !!                  ***  ROUTINE ldf_c2d  *** 
    109       !!               
     109      !! 
    110110      !! ** Purpose :   2D eddy diffusivity/viscosity coefficients 
    111111      !! 
     
    113113      !!       laplacian   operator :   ah proportional to the scale factor      [m2/s] 
    114114      !!       bilaplacian operator :   ah proportional to the (scale factor)^3  [m4/s] 
    115       !!       In both cases, pah0 is the maximum value reached by the coefficient  
     115      !!       In both cases, pah0 is the maximum value reached by the coefficient 
    116116      !!       at the Equator in case of e1=ra*rad= ~111km, not over the whole domain. 
    117117      !! 
     
    140140         END_2D 
    141141      CASE( 'TRA' )                       ! U- and V-points 
    142          DO_2D( 1, 1, 1, 1 ) 
     142         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    143143            pah1(ji,jj,1) = pUfac * MAX( e1u(ji,jj), e2u(ji,jj) )**knn 
    144144            pah2(ji,jj,1) = pUfac * MAX( e1v(ji,jj), e2v(ji,jj) )**knn 
Note: See TracChangeset for help on using the changeset viewer.