Changeset 14252 for NEMO/trunk/src/ICE/icethd_pnd.F90
- Timestamp:
- 2020-12-24T09:07:31+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icethd_pnd.F90
r14072 r14252 511 511 !!------------------------------------------------------------------- 512 512 REAL(wp), PARAMETER :: & ! shared parameters for topographic melt ponds 513 zTd = 0.15_wp , & ! temperature difference for freeze-up (C)513 zTd = 273._wp , & ! temperature difference for freeze-up (K) 514 514 zvp_min = 1.e-4_wp ! minimum pond volume (m) 515 515 … … 530 530 zv_mlt ! total amount of meltwater produced 531 531 532 REAL(wp), DIMENSION(jpi,jpj) :: zvolp, & !! total melt pond water available before redistribution and drainage 533 zvolp_res !! remaining melt pond water available after drainage 532 REAL(wp), DIMENSION(jpi,jpj) :: zvolp_ini , & !! total melt pond water available before redistribution and drainage 533 zvolp , & !! total melt pond water volume 534 zvolp_res !! remaining melt pond water available after drainage 534 535 535 536 REAL(wp), DIMENSION(jpi,jpj,jpl) :: z1_a_i … … 624 625 END DO ! ji 625 626 627 zvolp_ini(:,:) = zvolp(:,:) 628 626 629 !-------------------------------------------------------------- 627 630 ! Redistribute and drain water from ponds … … 637 640 DO_2D( 1, 1, 1, 1 ) 638 641 639 IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. vt_ip(ji,jj) > zvp_min * at_i(ji,jj) ) THEN642 IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. zvolp_ini(ji,jj) > zvp_min * at_i(ji,jj) ) THEN 640 643 641 644 !-------------------------- … … 688 691 689 692 ! differential growth of base of surface floating ice layer 690 zdTice = MAX( - t_su(ji,jj,jl) - zTd, 0._wp ) ! > 0693 zdTice = MAX( - ( t_su(ji,jj,jl) - zTd ) , 0._wp ) ! > 0 691 694 zomega = rcnd_i * zdTice / zrhoi_L 692 695 zdHui = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_i(ji,jj,jl) )**2 ) &
Note: See TracChangeset
for help on using the changeset viewer.