Changeset 4506 for branches/dev_r4028_CNRS_LIM3_MV2014/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2014-02-21T15:20:39+01:00 (10 years ago)

Author:

vancop

Message:

[Heat conservation in LIM3, part HC1 (LIM_SRC_3_HC17)]

File:

• branches/dev_r4028_CNRS_LIM3_MV2014/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (4 diffs)

branches/dev_r4028_CNRS_LIM3_MV2014/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -106 +106 @@
       INTEGER  ::   iflt, ial , iadv , ifral, ifrdv   !   -      -
       REAL(wp) ::   zinda, zemp, zemp_snow, zfmm      ! local scalars
-      REAL(wp) ::   zemp_snw                          !   -      -
+      REAL(wp) ::   zemp_snw, zqmass,   zcd           !   -      -  
+      REAL(wp) ::   zswitch                           !   -      -  
       REAL(wp) ::   zfcm1 , zfcm2                     !   -      -
       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb, zalbp     ! 2D/3D workspace
@@ -113 +114 @@
       
       IF( lk_cpl )   CALL wrk_alloc( jpi, jpj, jpl, zalb, zalbp )
+
+      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                                 !                                                                      
 
       !------------------------------------------!
@@ -172 +179 @@
                &    + ifrdv   * (       qfvbq(ji,jj) +             qdtcn(ji,jj) ) * r1_rdtice   &
                &    + fhmec(ji,jj)                                                              & ! snow melt when ridging
-               &    + fheat_mec(ji,jj)                                                          & ! ridge formation
                &    + fheat_res(ji,jj)                                                            ! residual heat flux
             ! qcmif   Energy needed to bring the ocean surface layer until its freezing (ok)
@@ -184 +190 @@
             fsbbq(ji,jj) = ( 1._wp - ( ifvt + iflt ) ) * fscmbq(ji,jj)     
 
-            ! used to compute the oceanic heat flux at the next time step
+            ! heat flux associated with ice-ocean mass exchange
+            zqmass = ( rdq_snw(ji,jj)                                         & 
+                 & + rdq_ice(ji,jj) * ( 1.- zswitch) ) * r1_rdtice   !  heat flux due to snow & ice heat content 
             qsr(ji,jj) = zfcm1                                       ! solar heat flux 
-            qns(ji,jj) = zfcm2 - fdtcn(ji,jj)                        ! non solar heat flux
-            !                           ! fdtcn : turbulent oceanic heat flux
+            qns(ji,jj) = zfcm2 - fdtcn(ji,jj) + zqmass               ! non solar heat flux   
+                               ! fdtcn : turbulent oceanic heat flux
          END DO
       END DO