Changeset 9513 for branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.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/dommsk.F90 (modified) (4 diffs)

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

@@ -25 +25 @@
    USE oce             ! ocean dynamics and tracers
    USE dom_oce         ! ocean space and time domain
+   USE domngb          ! find nearest wet point
+   USE domutl          ! fill closed 3d pool below isf
+   USE domzgr, ONLY : ln_isfsubgl, rn_isfsubgllon, rn_isfsubgllat ! import ln_isfsubgl to mask close sea below isf
+   !
    USE in_out_manager  ! I/O manager
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
@@ -133 +137 @@
       INTEGER  ::   iif, iil, ii0, ii1, ii   ! local integers
       INTEGER  ::   ijf, ijl, ij0, ij1       !   -       -
+      INTEGER  ::   jiseed, jjseed           !   -       -
       INTEGER  ::   ios
       INTEGER  ::   isrow                    ! index for ORCA1 starting row
@@ -187 +192 @@
       END DO  
       
-      ! (ISF) define barotropic mask and mask the ice shelf point
-      ssmask(:,:)=tmask(:,:,1) ! at this stage ice shelf is not masked
-      
       DO jk = 1, jpk
          DO jj = 1, jpj
@@ -199 +201 @@
          END DO  
       END DO  
-
+      !
+      IF ( ln_isfsubgl ) THEN
+         ! check closed wet pool
+         CALL dom_ngb(rn_isfsubgllon, rn_isfsubgllat, jiseed, jjseed, 'T', lwet=.TRUE.)
+         CALL fill_pool( jiseed, jjseed, 1, tmask, -1._wp )
+         ! at this point itab3d (:,1:ijmax,:) can have 3 different values :
+         !              0 where there where already 0
+         !              -1 where the ocean points are connected
+         !              1 where ocean points in tmask are not connected
+         IF (lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*)'dommsk : removal of subglacial lakes '
+            WRITE(numout,*)'~~~~~~~'
+            WRITE(numout,*)'Number of disconected points : ', COUNT(  (tmask(:,:,:) == 1) )
+            WRITE(numout,*)'lon/lat seed to detect main ocean is: ', rn_isfsubgllon, rn_isfsubgllat
+            WRITE(numout,*)'i/j     seed to detect main ocean is: ', jiseed, jjseed
+         END IF
+         DO jk = 1, jpk
+            WHERE (tmask(:,:,jk) > 0 .AND. misfdep(:,:) > 1) tmask(:,:,jk) = 0 ! remove only subglacial lake (ie similar to close sea only below an ice shelf 
+            WHERE (tmask(:,:,jk) < 0) tmask(:,:,jk) = 1 ! restore mask value
+         END DO
+      END IF
+
+      ! (ISF) define barotropic mask and mask the ice shelf point
+      DO jj = 1, jpj
+         DO ji = 1, jpi   ! vector loop
+            ssmask(ji,jj)  = MIN(1._wp,SUM(tmask(ji,jj,:)))
+         END DO
+      END DO
 !!gm  ????
 #if defined key_zdfkpp