Changeset 3964 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC

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/OPA_SRC/SBC/sbcfwb.F90 (modified) (3 diffs)

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