Changeset 3964

Timestamp:

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

Author:

flavoni

Message:

correction of fresh water budget for LIM3, see ticket 1116

Location:

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO

Files:

• LIM_SRC_3/limsbc.F90 (modified) (4 diffs)
• OPA_SRC/SBC/sbcfwb.F90 (modified) (3 diffs)
• OPA_SRC/oce.F90 (modified) (1 diff)

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
       !
       !

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

@@ -36 +36 @@
    REAL(wp) ::   fwfold    ! fwfold to be suppressed
    REAL(wp) ::   area      ! global mean ocean surface (interior domain)
+   REAL(wp) ::   r1_rau0   ! = 1 / rau0 
 
    !! * Substitutions
@@ -130 +131 @@
          ikty = 365 * 86400 / rdttra(1)    !!bug  use of 365 days leap year or 360d year !!!!!!!
          IF( MOD( kt, ikty ) == 0 ) THEN
+            !
+            r1_rau0 = 1._wp / rau0
+            !
             a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow
                                                       ! sum over the global domain
@@ -138 +142 @@
             !                                         ! estimate from the previous year budget
          ENDIF
+         IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold
          ! 
          IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN      ! correct the freshwater fluxes

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/oce.F90

@@ -85 +85 @@
          &     spgu  (jpi,jpj)   , spgv(jpi,jpj)   ,                       &
          &     gtsu(jpi,jpj,jpts), gtsv(jpi,jpj,jpts),                     &
-         &     gru(jpi,jpj)      , grv(jpi,jpj)                      , STAT=ierr(2) )
+         &     gru(jpi,jpj)      , grv(jpi,jpj),                           &
+         &     snwice_mass(jpi,jpj), snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj)    , STAT=ierr(2) )
          !
       oce_alloc = MAXVAL( ierr )