Changeset 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

Timestamp:

2017-12-26T17:32:56+01:00 (6 years ago)

Author:

gm

Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90 (modified) (6 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

@@ -142 +142 @@
       !                                ! Parameter control
       IF( ln_dynldf_NONE ) THEN
-         IF(lwp) WRITE(numout,*) '   No viscous operator selected. ahmt and ahmf are not allocated'
+         IF(lwp) WRITE(numout,*) '   ==>>>   No viscous operator selected. ahmt and ahmf are not allocated'
          l_ldfdyn_time = .FALSE.
          RETURN
@@ -173 +173 @@
          !
          CASE(   0  )      !==  constant  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = constant '
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = constant '
             ahmt(:,:,:) = zah0 * tmask(:,:,:)
             ahmf(:,:,:) = zah0 * fmask(:,:,:)
             !
          CASE(  10  )      !==  fixed profile  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( depth )'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F( depth )'
             ahmt(:,:,1) = zah0 * tmask(:,:,1)                      ! constant surface value
             ahmf(:,:,1) = zah0 * fmask(:,:,1)
@@ -184 +184 @@
             !
          CASE ( -20 )      !== fixed horizontal shape read in file  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F(i,j) read in eddy_viscosity.nc file'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F(i,j) read in eddy_viscosity.nc file'
             CALL iom_open( 'eddy_viscosity_2D.nc', inum )
             CALL iom_get ( inum, jpdom_data, 'ahmt_2d', ahmt(:,:,1) )
@@ -198 +198 @@
             !
          CASE(  20  )      !== fixed horizontal shape  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)'
             IF( ln_dynldf_lap )   CALL ldf_c2d( 'DYN', 'LAP', zah0, ahmt, ahmf )    ! surface value proportional to scale factor
             IF( ln_dynldf_blp )   CALL ldf_c2d( 'DYN', 'BLP', zah0, ahmt, ahmf )    ! surface value proportional to scale factor^3
             !
          CASE( -30  )      !== fixed 3D shape read in file  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F(i,j,k) read in eddy_diffusivity_3D.nc file'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F(i,j,k) read in eddy_diffusivity_3D.nc file'
             CALL iom_open( 'eddy_viscosity_3D.nc', inum )
             CALL iom_get ( inum, jpdom_data, 'ahmt_3d', ahmt )
@@ -216 +216 @@
             !
          CASE(  30  )       !==  fixed 3D shape  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( latitude, longitude, depth )'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F( latitude, longitude, depth )'
             IF( ln_dynldf_lap )   CALL ldf_c2d( 'DYN', 'LAP', zah0, ahmt, ahmf )    ! surface value proportional to scale factor
             IF( ln_dynldf_blp )   CALL ldf_c2d( 'DYN', 'BLP', zah0, ahmt, ahmf )    ! surface value proportional to scale factor
@@ -223 +223 @@
             !
          CASE(  31  )       !==  time varying 3D field  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( latitude, longitude, depth , time )'
-            IF(lwp) WRITE(numout,*) '                                proportional to the velocity : |u|e/12 or |u|e^3/12'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F( latitude, longitude, depth , time )'
+            IF(lwp) WRITE(numout,*) '              proportional to the velocity : |u|e/12 or |u|e^3/12'
             !
             l_ldfdyn_time = .TRUE.     ! will be calculated by call to ldf_dyn routine in step.F90
             !
          CASE(  32  )       !==  time varying 3D field  ==!
-            IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( latitude, longitude, depth , time )'
-            IF(lwp) WRITE(numout,*) '             proportional to the local deformation rate and gridscale (Smagorinsky)'
-            IF(lwp) WRITE(numout,*) '                                                             : L^2|D| or L^4|D|/8'
+            IF(lwp) WRITE(numout,*) '   ==>>>   momentum mixing coef. = F( latitude, longitude, depth , time )'
+            IF(lwp) WRITE(numout,*) '              proportional to the local deformation rate and gridscale (Smagorinsky)'
+            IF(lwp) WRITE(numout,*) '                                                                : L^2|D| or L^4|D|/8'
             !
             l_ldfdyn_time = .TRUE.     ! will be calculated by call to ldf_dyn routine in step.F90