Changeset 11423 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90

Timestamp:

2019-08-08T16:02:49+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo : add UKESM ice sheet coupling method (ticket #2142)

File:

• NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90 (modified) (1 diff)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfload.F90

@@ -90 +90 @@
       !
       !                                !- Surface value + ice shelf gradient
-      risfload = 0._wp                       ! compute pressure due to ice shelf load 
+      risfload(:,:) = 0._wp                       ! compute pressure due to ice shelf load 
       DO jj = 1, jpj                         ! (used to compute hpgi/j for all the level from 1 to miku/v)
          DO ji = 1, jpi                      ! divided by 2 later
             ikt = mikt(ji,jj)
             !
-            ! top layer of the ice shelf
-            risfload(ji,jj) = risfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1)
-            !
-            ! core layers of the ice shelf
-            DO jk = 2, ikt-1
-               risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk)
-            END DO
-            !
-            ! deepest part of the ice shelf (between deepest T point and ice/ocean interface
-            risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) &
-               &                                              * ( risfdep(ji,jj) - gdept_n(ji,jj,ikt-1) )
-            !
+            IF ( ikt > 1 ) THEN
+               !
+               ! top layer of the ice shelf
+               risfload(ji,jj) = risfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w_n(ji,jj,1)
+               !
+               ! core layers of the ice shelf
+               DO jk = 2, ikt-1
+                  risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w_n(ji,jj,jk)
+               END DO
+               !
+               ! deepest part of the ice shelf (between deepest T point and ice/ocean interface
+               risfload(ji,jj) = risfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) &
+                  &                                              * ( risfdep(ji,jj) - gdept_n(ji,jj,ikt-1) )
+               !
+            END IF
          END DO
       END DO
-      !
-      ! mask it 
-      risfload(:,:) = risfload(:,:) * mskisf_cav(:,:)  ! need to be saved for diaar5
       !
    END SUBROUTINE isf_load_isomip