Changeset 13179 for NEMO/branches/2020/ticket2487/src/OCE/SBC/sbcice_cice.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/OCE/SBC/sbcice_cice.F90 (modified) (2 diffs)

NEMO/branches/2020/ticket2487/src/OCE/SBC/sbcice_cice.F90

@@ -155 +155 @@
       INTEGER, INTENT( in  ) ::   ksbc                ! surface forcing type
       REAL(wp), DIMENSION(jpi,jpj) :: ztmp1, ztmp2
-      REAL(wp) ::   zcoefu, zcoefv, zcoeff            ! local scalar
+      REAL(wp) ::   zcoefu, zcoefv, zcoeff            ! local scalars
+      REAL(wp) ::   z1_area, zsshadj                  !   "      "
       INTEGER  ::   ji, jj, jl, jk                    ! dummy loop indices
       !!---------------------------------------------------------------------
@@ -266 +267 @@
                END DO
             ENDIF
+         ELSE
+            z1_area = 1.0_wp / glob_sum( 'sbcice_cice', e1e2t(:,:) )
+            zsshadj = glob_sum( 'sbcice_cice', e1e2t(:,:) * snwice_mass(:,:) ) * r1_rau0 * z1_area
+            WRITE(ctmp1,'(A36,F10.6,A24)') 'sbcice_cice:   mean ssh adjusted by ', -1.0_wp * zsshadj, ' m to compensate for the'
+            CALL ctl_warn( ctmp1,          '               initial ice+snow mass' )
+            sshn(:,:) = sshn(:,:) - zsshadj
+            sshb(:,:) = sshb(:,:) - zsshadj
          ENDIF
       ENDIF