Changeset 9513 for branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2018-04-27T14:30:02+02:00 (6 years ago)

Author:

mathiot

Message:

Add option to detect and remove subglacial lake (do not affect closed sea option)

File:

• branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (2 diffs)

branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -50 +50 @@
    PUBLIC   dom_zgr        ! called by dom_init.F90
 
+   ! remove subglacial lake under the ice sheet.
+   LOGICAL , PUBLIC  :: ln_isfsubgl
+   REAL(wp), PUBLIC  :: rn_isfsubgllon, rn_isfsubgllat
    !                              !!* Namelist namzgr_sco *
    LOGICAL  ::   ln_s_sh94         ! use hybrid s-sig Song and Haidvogel 1994 stretching function fssig1 (ln_sco=T)
@@ -1293 +1296 @@
       ! NAMELIST
       INTEGER  ::   ios
-      INTEGER  :: nn_kisfmax    = 999. !  
-      REAL(wp) :: rn_isfdep_min = 10.0_wp ! ice shelf minimal thickness 
+      INTEGER  :: nn_kisfmax    = 999.                              !  
+      REAL(wp) :: rn_isfdep_min = 10.0_wp                           ! ice shelf minimal thickness 
       REAL(wp) :: rn_isfhw_deep = 1.e-3 , rn_isfhw_shallow = 1.0e-3 ! threshold to define grounding line in deep/shallow water
-      REAL(wp) :: rn_isfshallow = 0.0_wp  ! threshold to define shallow ice shelf cavity
-      REAL(wp) :: rn_zisfmax    = 6000.0_wp ! maximun meter of ice we are allowed to dig to assure connectivity
-      LOGICAL  :: ln_isfcheminey= .FALSE. , ln_isfconnect= .FALSE. , ln_isfchannel= .FALSE. !  remove cheminey, assure connectivity
+      REAL(wp) :: rn_isfshallow = 0.0_wp                            ! threshold to define shallow ice shelf cavity
+      REAL(wp) :: rn_zisfmax    = 6000.0_wp                         ! maximun meter of ice we are allowed to dig to assure connectivity
+      LOGICAL  :: ln_isfcheminey= .FALSE. , ln_isfconnect= .FALSE. , ln_isfchannel= .FALSE. !  remove cheminey, assure connectivity, remove channel in water column (on z space not level space)
       !!---------------------------------------------------------------------
       NAMELIST/namzgr_isf/nn_kisfmax, rn_zisfmax, &
               &           rn_isfdep_min, rn_isfhw_deep, rn_isfhw_shallow, rn_isfshallow, &
-              &           ln_isfcheminey, ln_isfconnect, ln_isfchannel
+              &           ln_isfcheminey, ln_isfconnect, ln_isfchannel, ln_isfsubgl, rn_isfsubgllon, rn_isfsubgllat
       !
       IF( nn_timing == 1 )  CALL timing_start('zgr_isf')