Changeset 3402 for branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2012-05-25T18:43:49+02:00 (12 years ago)

Author:

acc

Message:

Branch: dev_r3385_NOCS04_HAMF; #665. Stage 2 of 2012 development: suppression of emps array and introduction of sfx (salt flux) array with associated code to setup the options for embedding the seaice into the ocean

File:

• branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (7 diffs)

branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -39 +39 @@
    USE prtctl           ! Print control
    USE cpl_oasis3, ONLY : lk_cpl
+   USE oce, ONLY : sshn, sshb, snwice_mass, snwice_mass_b, snwice_fmass 
 
    IMPLICIT NONE
@@ -92 +93 @@
       !!              - qns     : sea heat flux    : non solar (including heat content of the mass flux)
       !!              - emp     : freshwater budget: mass flux 
-      !!              - emps    : freshwater budget: salt flux due to Freezing/Melting
+      !!              - sfx     : freshwater budget: salt flux due to Freezing/Melting
       !!              - utau    : sea surface i-stress (ocean referential)
       !!              - vtau    : sea surface j-stress (ocean referential)
@@ -119 +120 @@
       CALL wrk_alloc( jpi, jpj, 1, zalb, zalbp )
 
-      zswitch = 1                                ! Default standard levitating sea-ice (salt exchanges only)
-!!gm ice embedment
-!     SELECT CASE( nn_ice_embd )                 ! levitating/embedded sea-ice option (not yet activated)
-!       CASE( 0    )   ;   zswitch = 1           ! standard levitating sea-ice : salt exchange only
-!       CASE( 1, 2 )   ;   zswitch = 0           ! other levitating sea-ice or embedded sea-ice : salt and volume fluxes
-!     END SELECT                                 !    
-!!gm end embedment
+      SELECT CASE( nn_ice_embd )                 ! levitating or embedded sea-ice option
+        CASE( 0    )   ;   zswitch = 1           ! (0) standard levitating sea-ice : salt exchange only
+        CASE( 1, 2 )   ;   zswitch = 0           ! (1) levitating sea-ice: salt and volume exchange but no pressure effect
+                                                 ! (2) embedded sea-ice : salt and volume fluxes and pressure
+      END SELECT                                 !    
+
       !------------------------------------------!
       !      heat flux at the ocean surface      !
@@ -147 +147 @@
 !!$            i1mfr   = 1.0 - AINT(  frld(ji,jj) )                   ! = 0. if ice-free ocean else 1. (after ice thermo)
 !!$
-!!$            IF( phicif(ji,jj) <= 0. ) THEN   ;   ifvt = zinda      ! = 1. if there was snow and ice before the ice thermo. which has been completely melted (possibly overmelted)
+!!$            IF( phicif(ji,jj) <= 0. ) THEN   ;   ifvt = zinda      ! = zinda if previous thermodynamic step overmelted the ice???
 !!$            ELSE                             ;   ifvt = 0.         ! 
 !!$            ENDIF
@@ -222 +222 @@
             ! salt flux only       : add concentration dilution term in salt flux  and no  F-M term in volume flux
             ! salt and mass fluxes : non concentartion dilution term in salt flux  and add F-M term in volume flux
-            emps(ji,jj) = zfsalt +                  zswitch  * zcd   ! salt flux (+ C/D if no ice/ocean mass exchange)
+            sfx (ji,jj) = zfsalt +                  zswitch  * zcd   ! salt flux (+ C/D if no ice/ocean mass exchange)
             emp (ji,jj) = zemp   + zemp_snw + ( 1.- zswitch) * zfmm  ! mass flux (- F/M mass flux if no ice/ocean mass exchange)
             !
          END DO
       END DO
+      !                                !------------------------------------------!
+      !                                !    mass of snow and ice per unit area    !
+      !                                !------------------------------------------!
+      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
+         snwice_mass  (:,:) = tms(:,:) * ( rhosn * hsnif(:,:) + rhoic * hicif(:,:)  ) * ( 1.0 - frld(:,:) )
+         !                                                      ! time evolution of snow+ice mass
+         snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / rdt_ice
+      ENDIF
 
       CALL iom_put( 'hflx_ice_cea', - fdtcn(:,:) )      
@@ -255 +265 @@
       IF(ln_ctl) THEN            ! control print
          CALL prt_ctl(tab2d_1=qsr   , clinfo1=' lim_sbc: qsr    : ', tab2d_2=qns   , clinfo2=' qns     : ')
-         CALL prt_ctl(tab2d_1=emp   , clinfo1=' lim_sbc: emp    : ', tab2d_2=emps  , clinfo2=' emps    : ')
+         CALL prt_ctl(tab2d_1=emp   , clinfo1=' lim_sbc: emp    : ', tab2d_2=sfx   , clinfo2=' sfx     : ')
          CALL prt_ctl(tab2d_1=utau  , clinfo1=' lim_sbc: utau   : ', mask1=umask,   &
             &         tab2d_2=vtau  , clinfo2=' vtau    : '        , mask2=vmask )
@@ -451 +461 @@
          END WHERE
       ENDIF
+      !                                      ! embedded sea ice
+      IF( nn_ice_embd /= 0 ) THEN            ! mass exchanges between ice and ocean (case 1 or 2) set the snow+ice mass
+         snwice_mass  (:,:) = tms(:,:) * ( rhosn * hsnif(:,:) + rhoic * hicif(:,:)  ) * ( 1.0 - frld(:,:) )
+         snwice_mass_b(:,:) = snwice_mass(:,:)
+      ELSE
+         snwice_mass  (:,:) = 0.e0           ! no mass exchanges
+         snwice_mass_b(:,:) = 0.e0           ! no mass exchanges
+      ENDIF
+      IF( nn_ice_embd == 2 .AND.          &  ! full embedment (case 2) & no restart : 
+         &   .NOT.ln_rstart ) THEN           ! deplete the initial ssh belew sea-ice area
+         sshn(:,:) = sshn(:,:) - snwice_mass(:,:) * r1_rau0
+         sshb(:,:) = sshb(:,:) - snwice_mass(:,:) * r1_rau0
+      ENDIF
       !
    END SUBROUTINE lim_sbc_init_2