Changeset 4506 for branches/dev_r4028_CNRS_LIM3_MV2014/NEMOGCM/NEMO/LIM_SRC_3/limthd_ent.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/limthd_ent.F90 (modified) (3 diffs)

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

@@ -22 +22 @@
    USE domain         !
    USE phycst         ! physical constants
+   USE sbc_oce        ! Surface boundary condition: ocean fields
    USE ice            ! LIM variables
    USE par_ice        ! LIM parameters
@@ -86 +87 @@
          ztmelts        ,   &  ! ice melting point
          zqsnic         ,   &  ! enthalpy of snow ice layer
+         zsstK          ,   &  ! SST in Kelvin
          zhsnow         ,   &  ! temporary snow thickness variable
          zswitch        ,   &  ! dummy switch argument
@@ -540 +542 @@
       DO ji = kideb, kiut
          ! energy of the flooding seawater
-         zqsnic = rau0 * rcp * ( rtt - t_bo_b(ji) ) * dh_snowice(ji) * &
-            (rhoic - rhosn) / rhoic * REAL(snicswi(ji)) ! generally positive
+         ii = MOD( npb(ji) - 1, jpi ) + 1
+         ij =    ( npb(ji) - 1 ) / jpi + 1
+         zsstK  = sst_m(ii,ij) + rt0
+         zqsnic = ( rhosn - rhoic ) * dh_snowice(ji) * rcp * ( zsstK - rt0 ) * REAL ( snicswi(ji) ) ! MV HC 2014
+         
          ! Heat conservation diagnostic
          qt_i_in(ji,jl) = qt_i_in(ji,jl) + zqsnic 
-
-         qldif_1d(ji)   = qldif_1d(ji) + zqsnic * a_i_b(ji)
 
          ! enthalpy of the newly formed snow-ice layer