Changeset 3977

Timestamp:

2013-07-16T14:04:24+02:00 (11 years ago)

Author:

flavoni

Message:

add print of fwb correction value, and hminrhg for LIM2, in CNRS LIM3 branch

Location:

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO

Files:

• LIM_SRC_2/ice_2.F90 (modified) (1 diff)
• LIM_SRC_3/iceini.F90 (modified) (2 diffs)
• LIM_SRC_3/limsbc.F90 (modified) (3 diffs)
• OPA_SRC/SBC/sbcfwb.F90 (modified) (1 diff)

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90

@@ -51 +51 @@
    REAL(wp), PUBLIC ::   ahi0   = 350._wp      !: sea-ice hor. eddy diffusivity coeff. (m2/s)
    REAL(wp), PUBLIC ::   alphaevp = 1._wp      !: coefficient for the solution of EVP int. stresses
+   REAL(wp), PUBLIC ::   hminrhg = 0.05_wp     !: clem : ice thickness (in m) below which ice velocity is set to ocean velocity
 
    REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( ecc * ecc )

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

@@ -85 +85 @@
       CALL lim_itd_ini                 ! ice thickness distribution initialization
       !
-      CALL lim_sbc_init                ! ice surface boundary condition   
-
 
       !                                ! Initial sea-ice state
@@ -102 +100 @@
       ENDIF
       !
+      CALL lim_sbc_init                ! ice surface boundary condition   
+      !
       fr_i(:,:) = at_i(:,:)           ! initialisation of sea-ice fraction
       !

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

@@ -54 +54 @@
    REAL(wp)  ::   rzero  = 0._wp    
    REAL(wp)  ::   rone   = 1._wp
+   REAL(wp) ::   r1_rau0   ! = 1 / rau0
 
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau_oce, vtau_oce   ! air-ocean surface i- & j-stress     [N/m2]
@@ -409 +410 @@
       !!-------------------------------------------------------------------
       REAL(wp) :: zsum, zarea
+      r1_rau0 = 1._wp / rau0
       !
       IF(lwp) WRITE(numout,*)
@@ -436 +438 @@
          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
+         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 * r1_rau0
+         sshn(:,:) = sshn(:,:) - zsum 
+         sshb(:,:) = sshb(:,:) - zsum
          ENDIF
       !

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

@@ -142 +142 @@
             !                                         ! estimate from the previous year budget
          ENDIF
-         IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold
+         IF( MOD( kt, 5475 ) == 0 ) THEN
+            IF(lwp) WRITE(numout,*) 'correction of fresh water budget :', fwfold
+         ENDIF 
          ! 
          IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN      ! correct the freshwater fluxes