Changeset 14789 for NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/OCE/LDF/ldfc1d_c2d.F90
- Timestamp:
- 2021-05-05T13:18:04+02:00 (3 years ago)
- 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 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev _r12970_AGRIF_CMEMSext/AGRIF5 ^/vendors/AGRIF/dev@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 ^/vendors/PPR@HEAD ext/PPR 8 9 9 10 # SETTE 10 ^/utils/CI/sette@1 3559sette11 ^/utils/CI/sette@14244 sette
-
- Property svn:externals
-
NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/OCE/LDF/ldfc1d_c2d.F90
r13497 r14789 2 2 !!====================================================================== 3 3 !! *** 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 5 5 !!===================================================================== 6 6 !! History : 3.7 ! 2013-12 (G. Madec) restructuration/simplification of aht/aeiv specification, … … 9 9 10 10 !!---------------------------------------------------------------------- 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) 12 12 !! ldf_c2d : ah = F(e1,e2) (laplacian or = F(e1^3,e2^3) (bilaplacian) 13 13 !!---------------------------------------------------------------------- … … 29 29 REAL(wp) :: r1_4 = 0.25_wp ! =1/4 30 30 REAL(wp) :: r1_12 = 1._wp / 12._wp ! =1/12 31 31 32 32 !! * Substitutions 33 33 # include "do_loop_substitute.h90" … … 42 42 !!---------------------------------------------------------------------- 43 43 !! *** ROUTINE ldf_c1d *** 44 !! 44 !! 45 45 !! ** Purpose : 1D eddy diffusivity/viscosity coefficients 46 46 !! 47 47 !! ** 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 50 50 !! at zh=500m and a width of zw=200m 51 51 !! … … 95 95 END_3D 96 96 ! 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 ) 98 98 ! 99 99 CASE DEFAULT ! error … … 107 107 !!---------------------------------------------------------------------- 108 108 !! *** ROUTINE ldf_c2d *** 109 !! 109 !! 110 110 !! ** Purpose : 2D eddy diffusivity/viscosity coefficients 111 111 !! … … 113 113 !! laplacian operator : ah proportional to the scale factor [m2/s] 114 114 !! 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 116 116 !! at the Equator in case of e1=ra*rad= ~111km, not over the whole domain. 117 117 !! … … 140 140 END_2D 141 141 CASE( 'TRA' ) ! U- and V-points 142 DO_2D( 1, 1, 1, 1)142 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 143 143 pah1(ji,jj,1) = pUfac * MAX( e1u(ji,jj), e2u(ji,jj) )**knn 144 144 pah2(ji,jj,1) = pUfac * MAX( e1v(ji,jj), e2v(ji,jj) )**knn
Note: See TracChangeset
for help on using the changeset viewer.