Changeset 3977 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3

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/LIM_SRC_3

Files:

• iceini.F90 (modified) (2 diffs)
• limsbc.F90 (modified) (3 diffs)

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
       !