Changeset 14072 for NEMO/trunk/src/OCE/ZDF/zdfric.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/ZDF/zdfric.F90 (modified) (9 diffs)

NEMO/trunk/src/OCE/ZDF/zdfric.F90

@@ -12 +12 @@
    !!            3.3  !  2010-10  (C. Ethe, G. Madec) reorganisation of initialisation phase
    !!            3.3.1!  2011-09  (P. Oddo) Mixed layer depth parameterization
-   !!            4.0  !  2017-04  (G. Madec)  remove CPP ddm key & avm at t-point only 
+   !!            4.0  !  2017-04  (G. Madec)  remove CPP ddm key & avm at t-point only
    !!----------------------------------------------------------------------
 
@@ -28 +28 @@
    USE in_out_manager ! I/O manager
    USE iom            ! I/O manager library
-   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)  
+   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
 
 
@@ -43 +43 @@
    REAL(wp) ::   rn_alp      ! coefficient of the parameterization
    REAL(wp) ::   rn_ekmfc    ! Ekman Factor Coeff
-   REAL(wp) ::   rn_mldmin   ! minimum mixed layer (ML) depth    
+   REAL(wp) ::   rn_mldmin   ! minimum mixed layer (ML) depth
    REAL(wp) ::   rn_mldmax   ! maximum mixed layer depth
    REAL(wp) ::   rn_wtmix    ! Vertical eddy Diff. in the ML
@@ -61 +61 @@
       !!----------------------------------------------------------------------
       !!                 ***  ROUTINE zdf_ric_init  ***
-      !!                    
+      !!
       !! ** Purpose :   Initialization of the vertical eddy diffusivity and
       !!      viscosity coef. for the Richardson number dependent formulation.
@@ -109 +109 @@
       !!----------------------------------------------------------------------
       !!                 ***  ROUTINE zdfric  ***
-      !!                    
+      !!
       !! ** Purpose :   Compute the before eddy viscosity and diffusivity as
       !!                a function of the local richardson number.
       !!
-      !! ** Method  :   Local richardson number dependent formulation of the 
-      !!                vertical eddy viscosity and diffusivity coefficients. 
+      !! ** Method  :   Local richardson number dependent formulation of the
+      !!                vertical eddy viscosity and diffusivity coefficients.
       !!                The eddy coefficients are given by:
       !!                    avm = avm0 + avmb
@@ -122 +122 @@
       !!                    avm0= rn_avmri / (1 + rn_alp*Ri)**nn_ric
       !!                where ri is the before local Richardson number,
-      !!                rn_avmri is the maximum value reaches by avm and avt 
+      !!                rn_avmri is the maximum value reaches by avm and avt
       !!                and rn_alp, nn_ric are adjustable parameters.
       !!                Typical values : rn_alp=5. and nn_ric=2.
@@ -164 +164 @@
       END_3D
       !
-!!gm BUG <<<<====  This param can't work at low latitude 
+!!gm BUG <<<<====  This param can't work at low latitude
 !!gm               it provides there much to thick mixed layer ( summer 150m in GYRE configuration !!! )
       !
@@ -188 +188 @@
       !!---------------------------------------------------------------------
       !!                   ***  ROUTINE ric_rst  ***
-      !!                     
+      !!
       !! ** Purpose :   Read or write TKE file (en) in restart file
       !!
       !! ** Method  :   use of IOM library
-      !!                if the restart does not contain TKE, en is either 
-      !!                set to rn_emin or recomputed 
+      !!                if the restart does not contain TKE, en is either
+      !!                set to rn_emin or recomputed
       !!----------------------------------------------------------------------
       INTEGER         , INTENT(in) ::   kt     ! ocean time-step
@@ -202 +202 @@
       !!----------------------------------------------------------------------
       !
-      IF( TRIM(cdrw) == 'READ' ) THEN        ! Read/initialise 
+      IF( TRIM(cdrw) == 'READ' ) THEN        ! Read/initialise
          !                                   ! ---------------
          !           !* Read the restart file