Changeset 9076 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC

Timestamp:

2017-12-15T14:59:08+01:00 (6 years ago)

Author:

vancop

Message:

monocat to virtual_itd

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90 (modified) (3 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90

@@ -975 +975 @@
    
 
-   SUBROUTINE blk_ice_qcn( k_monocat, ptsu, ptb, phs, phi )
+   SUBROUTINE blk_ice_qcn( k_virtual_itd, ptsu, ptb, phs, phi )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE blk_ice_qcn  ***
@@ -990 +990 @@
       !!
       !!---------------------------------------------------------------------
-      INTEGER                   , INTENT(in   ) ::   k_monocat   ! single-category option
-      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   ptsu        ! sea ice / snow surface temperature
-      REAL(wp), DIMENSION(:,:)  , INTENT(in   ) ::   ptb         ! sea ice base temperature
-      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   phs         ! snow thickness
-      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   phi         ! sea ice thickness
+      INTEGER                   , INTENT(in   ) ::   k_virtual_itd   ! single-category option
+      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   ptsu            ! sea ice / snow surface temperature
+      REAL(wp), DIMENSION(:,:)  , INTENT(in   ) ::   ptb             ! sea ice base temperature
+      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   phs             ! snow thickness
+      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   phi             ! sea ice thickness
       !
       INTEGER , PARAMETER ::   nit = 10                  ! number of iterations
@@ -1013 +1013 @@
       !      I   Enhanced conduction factor  !
       ! -------------------------------------!
-      ! Emulates the enhancement of conduction by unresolved thin ice (k_monocat = 1/3)
+      ! Emulates the enhancement of conduction by unresolved thin ice (k_virtual_itd = 1/2)
       ! Fichefet and Morales Maqueda, JGR 1997
       !
       zgfac(:,:,:) = 1._wp
       
-      SELECT CASE ( k_monocat ) 
-      !
-      CASE ( 1 , 3 )
+      SELECT CASE ( k_virtual_itd )
+      !
+      CASE ( 1 , 2 )
          !
          zfac  = 1._wp /  ( rn_cnd_s + rcdic )