Changeset 3964 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2013-07-09T17:42:43+02:00 (11 years ago)

Author:

flavoni

Message:

correction of fresh water budget for LIM3, see ticket 1116

File:

• branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (4 diffs)

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -286 +286 @@
       !   mass of snow and ice per unit area          !
       !-----------------------------------------------!
-      IF( nn_ice_embd /= 0 ) THEN                               ! embedded sea-ice (mass required)
+     !SF IF( nn_ice_embd /= 0 ) THEN                               ! embedded sea-ice (mass required)
          snwice_mass_b(:,:) = snwice_mass(:,:)                  ! save mass from the previous ice time step
          !                                                      ! new mass per unit area
@@ -292 +292 @@
          !                                                      ! time evolution of snow+ice mass
          snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) * r1_rdtice
-      ENDIF
+      !SF ENDIF
 
       !-----------------------------------------------!
@@ -408 +408 @@
       !! ** input   : Namelist namicedia
       !!-------------------------------------------------------------------
+      REAL(wp) :: zsum, zarea
       !
       IF(lwp) WRITE(numout,*)
@@ -435 +436 @@
          snwice_mass  (:,:) = tms(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:)  )
          snwice_mass_b(:,:) = snwice_mass(:,:)
+         IF( .NOT.ln_rstart ) THEN           ! delete the initial ssh below sea-ice area
+            zarea     = glob_sum( e1e2t(:,:) )           ! interior global domain surface
+            zsum      = glob_sum( e1e2t(:,:) * ( snwice_mass(:,:) ) ) / zarea / rau0
+            sshn(:,:) = sshn(:,:) - zsum 
+            sshb(:,:) = sshb(:,:) - zsum
+         ENDIF
       !
       !