Changeset 13179 for NEMO/branches/2020/ticket2487/src/ICE/iceistate.F90

Timestamp:

2020-06-29T20:04:05+02:00 (4 years ago)

Author:

smueller

Message:

Addition of initial sea-level compensation for non-zero sea-ice/snow mass (see ticket #2487)

File:

• NEMO/branches/2020/ticket2487/src/ICE/iceistate.F90 (modified) (2 diffs)

NEMO/branches/2020/ticket2487/src/ICE/iceistate.F90

@@ -92 +92 @@
       !!
       INTEGER  ::   ji, jj, jk, jl         ! dummy loop indices
-      REAL(wp) ::   ztmelts
+      REAL(wp) ::   ztmelts, z1_area, zsshadj
       INTEGER , DIMENSION(4)           ::   itest
       REAL(wp), DIMENSION(jpi,jpj)     ::   z2d
@@ -425 +425 @@
             END DO
          ENDIF
+      ELSE
+         z1_area = 1.0_wp / glob_sum( 'iceistate', e1e2t(:,:) )
+         zsshadj = glob_sum( 'iceistate', e1e2t(:,:) * snwice_mass(:,:) ) * r1_rau0 * z1_area
+         WRITE(ctmp1,'(A34,F10.6,A32)') 'iceistate:   mean ssh adjusted by ', -1.0_wp * zsshadj, ' m to compensate for the initial'
+         CALL ctl_warn( ctmp1,          '             ice+snow mass' )
+         sshn(:,:) = sshn(:,:) - zsshadj
+         sshb(:,:) = sshb(:,:) - zsshadj
       ENDIF