Changeset 14072 for NEMO/trunk/src/OCE/LDF/ldfc1d_c2d.F90

Timestamp:

2020-12-04T08:48:38+01:00 (3 years ago)

Author:

laurent

Message:

Merging branch "2020/dev_r13648_ASINTER-04_laurent_bulk_ice", ticket #2369

File:

• NEMO/trunk/src/OCE/LDF/ldfc1d_c2d.F90 (modified) (7 diffs)

NEMO/trunk/src/OCE/LDF/ldfc1d_c2d.F90

@@ -2 +2 @@
    !!======================================================================
    !!                    ***  MODULE  ldfc1d_c2d  ***
-   !! Ocean physics:  profile and horizontal shape of lateral eddy coefficients 
+   !! Ocean physics:  profile and horizontal shape of lateral eddy coefficients
    !!=====================================================================
    !! History :  3.7  ! 2013-12  (G. Madec)  restructuration/simplification of aht/aeiv specification,
@@ -9 +9 @@
 
    !!----------------------------------------------------------------------
-   !!   ldf_c1d       : ah reduced by 1/4 on the vertical (tanh profile, inflection at 300m) 
+   !!   ldf_c1d       : ah reduced by 1/4 on the vertical (tanh profile, inflection at 300m)
    !!   ldf_c2d       : ah = F(e1,e2) (laplacian or = F(e1^3,e2^3) (bilaplacian)
    !!----------------------------------------------------------------------
@@ -29 +29 @@
    REAL(wp) ::   r1_4  = 0.25_wp          ! =1/4
    REAL(wp) ::   r1_12 = 1._wp / 12._wp   ! =1/12
- 
+
    !! * Substitutions
 #  include "do_loop_substitute.h90"
@@ -42 +42 @@
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE ldf_c1d  ***
-      !!              
+      !!
       !! ** Purpose :   1D eddy diffusivity/viscosity coefficients
       !!
       !! ** Method  :   1D eddy diffusivity coefficients F( depth )
-      !!                Reduction by zratio from surface to bottom 
-      !!                hyperbolic tangent profile with inflection point 
+      !!                Reduction by zratio from surface to bottom
+      !!                hyperbolic tangent profile with inflection point
       !!                at zh=500m and a width of zw=200m
       !!
@@ -95 +95 @@
          END_3D
          ! Lateral boundary conditions
-         CALL lbc_lnk_multi( 'ldfc1d_c2d', pah1, 'U', 1.0_wp , pah2, 'V', 1.0_wp )   
+         CALL lbc_lnk_multi( 'ldfc1d_c2d', pah1, 'U', 1.0_wp , pah2, 'V', 1.0_wp )
          !
       CASE DEFAULT                        ! error
@@ -107 +107 @@
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE ldf_c2d  ***
-      !!              
+      !!
       !! ** Purpose :   2D eddy diffusivity/viscosity coefficients
       !!
@@ -113 +113 @@
       !!       laplacian   operator :   ah proportional to the scale factor      [m2/s]
       !!       bilaplacian operator :   ah proportional to the (scale factor)^3  [m4/s]
-      !!       In both cases, pah0 is the maximum value reached by the coefficient 
+      !!       In both cases, pah0 is the maximum value reached by the coefficient
       !!       at the Equator in case of e1=ra*rad= ~111km, not over the whole domain.
       !!