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

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

@@ -184 +184 @@
                &                            + fdtcn(ji,jj)                        &   ! turbulent ice-ocean heat
                &                            + fsbbq(ji,jj) * ( 1.0 - zinda )  )   &   ! residual heat from previous step
+               ! MV HC 2014 check that this is good.
+               ! We should remove the heat content of precip that has fallen on sea ice
                & - pfrld(ji,jj)**betas * sprecip(ji,jj) * lfus                    )   ! latent heat of sprecip melting
+               ! MV HC 2014 partie heat content manque
             !
             ! Positive heat budget is used for bottom ablation
@@ -282 +285 @@
             CALL tab_2d_1d( nbpb, qldif_1d   (1:nbpb), qldif           , jpi, jpj, npb(1:nbpb) )
             CALL tab_2d_1d( nbpb, rdm_ice_1d (1:nbpb), rdm_ice         , jpi, jpj, npb(1:nbpb) )
+            CALL tab_2d_1d( nbpb, rdq_ice_1d (1:nbpb), rdq_ice         , jpi, jpj, npb(1:nbpb) )
             CALL tab_2d_1d( nbpb, rdm_snw_1d (1:nbpb), rdm_snw         , jpi, jpj, npb(1:nbpb) )
+            CALL tab_2d_1d( nbpb, rdq_snw_1d (1:nbpb), rdq_snw         , jpi, jpj, npb(1:nbpb) )
             CALL tab_2d_1d( nbpb, dmgwi_1d   (1:nbpb), dmgwi           , jpi, jpj, npb(1:nbpb) )
             CALL tab_2d_1d( nbpb, qlbbq_1d   (1:nbpb), zqlbsbq         , jpi, jpj, npb(1:nbpb) )
@@ -349 +354 @@
                CALL tab_1d_2d( nbpb, qfvbq         , npb, qfvbq_1d  (1:nbpb)   , jpi, jpj )
                CALL tab_1d_2d( nbpb, rdm_ice       , npb, rdm_ice_1d(1:nbpb)   , jpi, jpj )
+               CALL tab_1d_2d( nbpb, rdq_ice       , npb, rdq_ice_1d(1:nbpb)   , jpi, jpj )
                CALL tab_1d_2d( nbpb, rdm_snw       , npb, rdm_snw_1d(1:nbpb)   , jpi, jpj )
+               CALL tab_1d_2d( nbpb, rdq_snw       , npb, rdq_snw_1d(1:nbpb)   , jpi, jpj )
                CALL tab_1d_2d( nbpb, dmgwi         , npb, dmgwi_1d  (1:nbpb)   , jpi, jpj )
                CALL tab_1d_2d( nbpb, rdvosif       , npb, dvsbq_1d  (1:nbpb)   , jpi, jpj )