- Timestamp:
- 2017-12-04T10:41:40+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_dh.F90
r8882 r8885 270 270 !------------------------------ 271 271 ! qla_ice is always >=0 (upwards), heat goes to the atmosphere, therefore snow sublimates 272 ! clemcomment: not counted in mass/heat exchange in iceupdate.F90 since this is an exchange with atm. (not ocean)272 ! comment: not counted in mass/heat exchange in iceupdate.F90 since this is an exchange with atm. (not ocean) 273 273 zdeltah(1:npti,:) = 0._wp 274 274 DO ji = 1, npti … … 332 332 hfx_res_1d(ji) = hfx_res_1d(ji) + zfmdt * a_i_1d(ji) * zEi * r1_rdtice 333 333 334 ! Contribution to salt flux ( clem:using s_i_1d and not sz_i_1d(jk) is ok)334 ! Contribution to salt flux (using s_i_1d and not sz_i_1d(jk) is ok) 335 335 sfx_res_1d(ji) = sfx_res_1d(ji) - rhoic * a_i_1d(ji) * zdeltah(ji,jk) * s_i_1d(ji) * r1_rdtice 336 336 … … 358 358 zQm = zfmdt * zEw ! Energy of the melt water sent to the ocean [J/m2, <0] 359 359 360 ! Contribution to salt flux >0 ( clem:using s_i_1d and not sz_i_1d(jk) is ok)360 ! Contribution to salt flux >0 (using s_i_1d and not sz_i_1d(jk) is ok) 361 361 sfx_sum_1d(ji) = sfx_sum_1d(ji) - rhoic * a_i_1d(ji) * zdeltah(ji,jk) * s_i_1d(ji) * r1_rdtice 362 362 … … 377 377 zdeltah(ji,jk) = zdeltah(ji,jk) + zdum 378 378 dh_i_sub(ji) = dh_i_sub(ji) + zdum 379 ! Salt flux > 0 (clem 2016: flux is sent to the ocean for simplicity but salt should remain in the ice except if all ice is melted.379 ! Salt flux > 0 (clem: flux is sent to the ocean for simplicity but salt should remain in the ice except if all ice is melted. 380 380 ! It must be corrected at some point) 381 381 sfx_sub_1d(ji) = sfx_sub_1d(ji) - rhoic * a_i_1d(ji) * zdum * s_i_1d(ji) * r1_rdtice … … 526 526 hfx_res_1d(ji) = hfx_res_1d(ji) + zfmdt * a_i_1d(ji) * zEi * r1_rdtice 527 527 528 ! Contribution to salt flux ( clem:using s_i_1d and not sz_i_1d(jk) is ok)528 ! Contribution to salt flux (using s_i_1d and not sz_i_1d(jk) is ok) 529 529 sfx_res_1d(ji) = sfx_res_1d(ji) - rhoic * a_i_1d(ji) * zdeltah(ji,jk) * s_i_1d(ji) * r1_rdtice 530 530 … … 559 559 hfx_thd_1d(ji) = hfx_thd_1d(ji) + zfmdt * a_i_1d(ji) * zEw * r1_rdtice 560 560 561 ! Contribution to salt flux ( clem:using s_i_1d and not sz_i_1d(jk) is ok)561 ! Contribution to salt flux (using s_i_1d and not sz_i_1d(jk) is ok) 562 562 sfx_bom_1d(ji) = sfx_bom_1d(ji) - rhoic * a_i_1d(ji) * zdeltah(ji,jk) * s_i_1d(ji) * r1_rdtice 563 563
Note: See TracChangeset
for help on using the changeset viewer.