Changeset 13811
- Timestamp:
- 2020-11-18T13:05:56+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/SI3-05_MP/src/ICE/icethd_pnd.F90
r13810 r13811 53 53 !! 54 54 !! ** Purpose : change melt pond fraction and thickness 55 !! 56 !! Note: Melt ponds affect currently affect radiative transfer 57 !! They carry no heat, and the melt water they carry is not 58 !! exchanged with the ocean 59 !! 60 !! This means freshwater is directly released after surface & bottom melt in ice_thd_dh 61 !! 62 !! A wfx_pnd has been coded for diagnostic purposes 63 !! Each time wfx_pnd is updated, wfx_sum / wfx_snw_sum must be updated 64 !! 65 !!! The current diagnostic lacks a contribution from drainage 55 66 !! 56 67 !!------------------------------------------------------------------- … … 226 237 !--- mass flux ---! 227 238 IF( zdv_mlt > 0._wp ) THEN 228 ! MV add comment on what that mass flux means 229 ! water removed from fw flux due to melt pond growth 239 ! MV add comment on what that mass flux means 230 240 zfac = zdv_mlt * rhow * r1_rdtice ! melt pond mass flux < 0 [kg.m-2.s-1] 231 241 wfx_pnd_1d(ji) = wfx_pnd_1d(ji) - zfac … … 277 287 278 288 ! Calculate the permeability of the ice (Assur 1958, see Flocco 2010) 279 ! zsbr = - ( t_i_1d(ji,jk) - rt0 ) / rTmlt280 289 DO jk = 1, nlay_i 281 290 zsbr = - 1.2_wp & … … 293 302 zdv_flush = MAX( zdv_flush, -v_ip_1d(ji) ) 294 303 v_ip_1d(ji) = v_ip_1d(ji) + zdv_flush 304 305 ! MV --- why pond drainage does not give back water into freshwater flux ? 295 306 296 307 !--- Set new pond area and depth ---! assuming linear relation between h_ip and a_ip_frac … … 407 418 408 419 REAL (wp), PARAMETER :: & 420 ! MV ouate de phoque!!! constants hard coded ???? 421 ! 917 = rhoice 422 ! 0.334 =Lfus!!!! 409 423 zr1_rlfus = 1._wp / 0.334e+6 / 917._wp , & ! (J/m^3) 410 424 zTp = -0.15_wp, & ! pond freezing temperature (C) … … 474 488 475 489 !js 23/04/19: rewired back to a fraction with a_i 476 zrfrac(:,:,:) = rn_pnd_fracmin + ( rn_pnd_fracmax - rn_pnd_fracmin ) * aicen(:,:,:) 490 !!! zrfrac(:,:,:) = rn_pnd_fracmin + ( rn_pnd_fracmax - rn_pnd_fracmin ) * aicen(:,:,:) 491 ! MV2020 492 zrfrac(:,:,:) = rn_apnd_min + ( rn_apnd_max - rn_apnd_min ) * at_i_1d(:) ! = ( 1 - r ) = fraction of melt water that is not flushed 477 493 zwfx_tmp(:,:) = 0._wp 478 479 ! MV ---> use expression from level-ice ponds, clarify what is needed480 494 481 495 !--- Add retained melt water to melt ponds … … 496 510 * ( dh_i_pnd(:,:,jl)*rhoic + dh_s_pnd(:,:,jl)*rhosn ) & 497 511 * z1_rhow * a_i(:,:,jl) 512 ! MV ---> use expression from level ice melt ponds (dv_mlt) 513 498 514 END DO 499 515 500 ! MV ---> use expression from level ice melt ponds (dv_mlt) 501 516 517 !!! MV ---> rewrite this accounting for level ice melt ponds 518 !!! 1) wfx_pnd_in is obsolete, ln_pnd_fw as well 519 !!! 2) 520 !!! wfx_sum 521 !!! 522 !!! 502 523 !js 03/05/19: we truncate negative values after calculating zvolp, in a 503 524 ! similar manner to the subroutine ice_thd_pnd_cesm. Variation dh_i_pnd and … … 531 552 wfx_snw_sum(:,:) = ( 1._wp - zwfx_tmp(:,:) ) * wfx_snw_sum(:,:) 532 553 ENDIF 554 555 533 556 534 557 !----------------------------------------------------------------- … … 568 591 END DO ! jj 569 592 593 594 595 !!! MV sounds like this should be replaced by lid 570 596 571 597 DO ij = 1, icells
Note: See TracChangeset
for help on using the changeset viewer.