Changeset 9732

Timestamp:

2018-06-04T12:50:48+02:00 (7 years ago)

Author:

mathiot

Message:

replace the specific treatment of shallow ice shelves by a more general case

Location:

branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM

Files:

• CONFIG/SHARED/namelist_ref (modified) (1 diff)
• NEMO/OPA_SRC/DOM/domwri.F90 (modified) (1 diff)
• NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (5 diffs)

branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/CONFIG/SHARED/namelist_ref

@@ -121 +121 @@
 !-----------------------------------------------------------------------
    rn_isfdep_min    = 10.         ! minimum isf draft tickness (if lower, isf draft set to this value)
-   rn_isfhw_deep    = 1.e-3       ! minimum water column thickness to declare a column 'wet' in case of deep isf
-   rn_isfshallow    = 0.0         ! bathy between surface and rn_isfshallow are declared shallow
-   rn_isfhw_shallow = 1.e-3       ! minimum water column thickness to declare a column 'wet' in case of shallow isf
+   rn_glhw_min      = 1.e-3       ! minimum water column thickness to define the grounding line
+   rn_isfhw_min     = 10          ! minimum water column thickness in the cavity once the grounding line defined.
    ln_isfchannel    = .false.     ! remove channel (based on 2d mask build from isfdraft-bathy)
    ln_isfconnect    = .false.     ! force connection under the ice shelf (based on 2d mask build from isfdraft-bathy)

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

@@ -194 +194 @@
       zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp )
       CALL iom_rstput( 0, 0, inum4, 'isfdraft', zprt, ktype = jp_r4 )       !    ! nb of ocean T-points
+      zprt(:,:) = ssmask(:,:) * REAL( bathy(:,:) , wp )
+      CALL iom_rstput( 0, 0, inum4, 'bathy', zprt, ktype = jp_r4 )       !    ! nb of ocean T-points
+            
             
       IF( ln_sco ) THEN                                         ! s-coordinate

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

@@ -1269 +1269 @@
       !!                          1.0 = set iceshelf to the minimum depth allowed
       !!                          1.1 = ground ice shelf if water column less than X m
-      !!                          1.2 = ensure a minimum thickness for iceshelf cavity in shallow water
+      !!                          1.2 = ensure a minimum thickness for iceshelf cavity
       !!                          1.3 = remove channels and single point 'bay'
       !!                          1.4 = close single isolated point
@@ -1298 +1298 @@
       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_glhw_min   = 1.0e-3                             ! threshold on hw to define grounding line water
+      REAL(wp) :: rn_isfhw_min  = 1.0e-3                            ! threshold on hw to define isf draft into the cavity
       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
@@ -1304 +1305 @@
       !!---------------------------------------------------------------------
       NAMELIST/namzgr_isf/nn_kisfmax, rn_zisfmax, &
-              &           rn_isfdep_min, rn_isfhw_deep, rn_isfhw_shallow, rn_isfshallow, &
+              &           rn_isfdep_min, rn_glhw_min, rn_isfhw_min, &
               &           ln_isfcheminey, ln_isfconnect, ln_isfchannel, ln_isfsubgl, rn_isfsubgllon, rn_isfsubgllat
       !
@@ -1326 +1327 @@
          WRITE(numout,*) ' rn_zisfmax       = ',rn_zisfmax           !
          WRITE(numout,*) ' rn_isfdep_min    = ',rn_isfdep_min        !
-         WRITE(numout,*) ' rn_isfhw_deep    = ',rn_isfhw_deep        !
-         WRITE(numout,*) ' rn_isfhw_shallow = ',rn_isfhw_shallow     !
-         WRITE(numout,*) ' rn_isfshallow    = ',rn_isfshallow        !
+         WRITE(numout,*) ' rn_isfhw_min     = ',rn_isfhw_min         !
+         WRITE(numout,*) ' rn_glhw_min      = ',rn_glhw_min          !
          WRITE(numout,*) ' ln_isfcheminey   = ',ln_isfcheminey       !
          WRITE(numout,*) ' ln_isfconnect    = ',ln_isfconnect        !
@@ -1339 +1339 @@
       END WHERE
       !  
-      ! 1.1 ground ice shelf if water column less than X m
-      WHERE( bathy(:,:) - risfdep(:,:) < rn_isfhw_deep .AND. (bathy(:,:) > rn_isfshallow) ) 
+      ! 1.1 ground ice shelf if water column less than X m => set the grounding
+      ! line position
+      WHERE( bathy(:,:) - risfdep(:,:) < rn_glhw_min ) 
          risfdep(:,:)=0.0 ; misfdep(:,:) = 1
          bathy  (:,:)=0.0 ; mbathy (:,:) = 0
       END WHERE
       !
-      ! 1.2 ensure a minimum thickness for iceshelf cavity in shallow water
-      WHERE( (bathy(:,:) <= rn_isfshallow) .AND. (bathy(:,:) - risfdep(:,:) < rn_isfhw_shallow) ) 
-         risfdep(:,:)=0.0 ; misfdep(:,:) = 1
-         bathy  (:,:)=0.0 ; mbathy (:,:) = 0
-      END WHERE
+      ! 1.2 ensure a minimum thickness for iceshelf cavity => avoid to negative
+      ! e3t if ssh + sum(e3t*tmask) < 0
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            IF (bathy(ji,jj)-risfdep(ji,jj) < rn_isfhw_min) THEN
+               IF ((bathy(ji,jj) - rn_isfhw_min) > MAX(rn_isfdep_min,e3t_1d(1))) THEN 
+                  risfdep(ji,jj) = bathy(ji,jj) - rn_isfhw_min
+               ELSE
+                  risfdep(ji,jj)=0.0 ; misfdep(ji,jj) = 1
+                  bathy  (ji,jj)=0.0 ; mbathy (ji,jj) = 0
+               END IF
+            END IF
+         END DO
+      END DO
       !
       ! 1.3 Remove channels and single point 'bay'.