- Timestamp:
- 2020-02-12T13:06:25+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_dh.F90
r10888 r12369 131 131 DO ji = 1, npti 132 132 zq_top(ji) = MAX( 0._wp, qml_ice_1d(ji) * rdt_ice ) 133 134 IF ( to_print(ji) == 10 ) THEN 135 write(numout,*)'loading flux: zq_top(ji), qml_ice_1d(ji) = ',zq_top(ji), ' ', qml_ice_1d(ji) 136 ENDIF 133 137 END DO 134 138 ! … … 139 143 qml_ice_1d(ji) = zdum * MAX( 0._wp , SIGN( 1._wp, t_su_1d(ji) - rt0 ) ) 140 144 zq_top(ji) = MAX( 0._wp, qml_ice_1d(ji) * rdt_ice ) 145 146 IF ( to_print(ji) == 10 ) THEN 147 write(numout,*)'icethd_dh: qns_ice_1d(ji), qsr_ice_1d(ji), qtr_ice_top_1d(ji), qcn_ice_top_1d(ji) = ',qns_ice_1d(ji), ' ', qsr_ice_1d(ji), ' ', qtr_ice_top_1d(ji), ' ', qcn_ice_top_1d(ji) 148 write(numout,*)'icethd_dh: zdum, t_su_1d(ji), qml_ice_1d(ji), zq_top(ji) = ',zdum, ' ', t_su_1d(ji), qml_ice_1d(ji), zq_top(ji) 149 ENDIF 150 141 151 END DO 142 152 ! … … 146 156 zf_tt(ji) = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji) 147 157 zq_bot(ji) = MAX( 0._wp, zf_tt(ji) * rdt_ice ) 158 IF ( to_print(ji) == 1 ) THEN 159 write(numout,*)'icethd_dh: qcn_ice_bot_1d(ji), qsb_ice_bot_1d(ji) = ',qcn_ice_bot_1d(ji), qsb_ice_bot_1d(ji) 160 write(numout,*)'icethd_dh: fhld_1d(ji), zq_bot(ji) = ',fhld_1d(ji), ', ', zq_bot(ji) 161 ENDIF 148 162 END DO 149 163 … … 181 195 t_s_1d (ji,jk) = rt0 182 196 END IF 183 END DO 184 END DO 197 198 IF (to_print(ji) == 10) THEN 199 write(numout,*)'icethd_dh 8: t_s_1d(ji,jk), e_s_1d(ji,jk) = ',t_s_1d(ji,jk), ' ', e_s_1d(ji,jk) 200 END IF 201 END DO 202 END DO 203 204 DO ji = 1, npti 205 IF (to_print(ji) == 10) THEN 206 write(numout,*)'icethd_dh internal_melting: wfx_snw_sum_1d(ji) = ',wfx_snw_sum_1d(ji) 207 END IF 208 END DO 185 209 186 210 ! Snow precipitation … … 210 234 zq_top (ji) = MAX( 0._wp , zq_top (ji) + zdeltah(ji,1) * zqprec(ji) ) 211 235 zdh_s_pre(ji) = zdh_s_pre(ji) + zdeltah(ji,1) 236 237 IF ( to_print(ji) == 10 ) THEN 238 write(numout,*)'after precip: zq_top(ji), zdeltah(ji,1), zqprec(ji) = ',zq_top(ji), zdeltah(ji,1), zqprec(ji) 239 ENDIF 212 240 213 241 ! update thickness … … 220 248 ! 221 249 ENDIF 250 251 IF (to_print(ji) == 10) THEN 252 write(numout,*)'icethd_dh snow_precip: wfx_snw_sum_1d(ji) = ',wfx_snw_sum_1d(ji) 253 END IF 254 222 255 END DO 223 256 … … 253 286 ! 254 287 ENDIF 255 END DO 288 289 IF ( to_print(ji) == 10 ) THEN 290 write(numout,*)'after snow melt: zq_top(ji), zdeltah(ji,jk), e_s_1d(ji,jk), test = ',zq_top(ji), zdeltah(ji,jk), e_s_1d(ji,jk), ( zh_s(ji,jk) > 0._wp .AND. zq_top(ji) > 0._wp ) 291 ENDIF 292 293 IF (to_print(ji) == 10) THEN 294 write(numout,*)'icethd_dh loop: wfx_snw_sum_1d(ji), jk = ',wfx_snw_sum_1d(ji), ' ', jk 295 END IF 296 END DO 297 END DO 298 299 DO ji = 1, npti 300 IF (to_print(ji) == 10) THEN 301 write(numout,*)'icethd_dh snow_ground: wfx_snw_sum_1d(ji) = ',wfx_snw_sum_1d(ji) 302 END IF 256 303 END DO 257 304 … … 301 348 & ( ( zdh_s_pre(ji) ) * zqprec(ji) + & 302 349 & ( h_s_1d(ji) - zdh_s_pre(ji) ) * rhos * ( rcpi * ( rt0 - t_s_1d(ji,jk) ) + rLfus ) ) 350 351 IF (to_print(ji) == 10) THEN 352 write(numout,*)'icethd_dh 9: t_s_1d(ji,jk), e_s_1d(ji,jk) = ',t_s_1d(ji,jk), ' ', e_s_1d(ji,jk) 353 write(numout,*)'icethd_dh 9: zdh_s_pre(ji), zqprec(ji), h_s_1d(ji) = ',zdh_s_pre(ji), ' ', zqprec(ji), ' ', h_s_1d(ji) 354 END IF 303 355 END DO 304 356 END DO … … 362 414 363 415 END IF 416 417 IF ( to_print(ji) == 10 ) THEN 418 write(numout,*)'after top melt: zq_top(ji), test = ',zq_top(ji), ( t_i_1d(ji,jk) >= (ztmelts+rt0) ) 419 write(numout,*)'after top melt: wfx_sum_1d(ji), zdeltah(ji,jk), zfmdt, zdE = ', wfx_sum_1d(ji), zdeltah(ji,jk), zfmdt, zdE 420 ENDIF 364 421 365 422 ! Ice sublimation … … 510 567 zdeltah(ji,jk) = - zfmdt * r1_rhoi ! Gross thickness change 511 568 569 IF ( to_print(ji) == 1 ) THEN 570 write(numout,*)'icethd_dh: zfmdt, zq_bot(ji), zdE, zdeltah(ji,jk) = ',zfmdt,', ', zq_bot(ji), zdE, zdeltah(ji,jk) 571 ENDIF 572 512 573 zdeltah(ji,jk) = MIN( 0._wp , MAX( zdeltah(ji,jk), - zh_i(ji,jk) ) ) ! bound thickness change 513 574 … … 527 588 528 589 wfx_bom_1d(ji) = wfx_bom_1d(ji) - rhoi * a_i_1d(ji) * zdeltah(ji,jk) * r1_rdtice ! Mass flux 590 591 IF (to_print(ji) == 1) THEN 592 write(numout,*)'icethd_dh: wfx_bom_1d(ji), in meters, at_i_1d(ji) = ',wfx_bom_1d(ji),wfx_bom_1d(ji)*rdt_ice/(rhoi*at_i_1d(ji)), at_i_1d(ji) 593 ENDIF 529 594 530 595 ! update heat content (J.m-2) and layer thickness … … 536 601 END DO 537 602 END DO 603 604 538 605 539 606 ! Update temperature, energy … … 564 631 565 632 IF( ln_icectl .AND. zq_rema(ji) < 0. .AND. lwp ) WRITE(numout,*) 'ALERTE zq_rema <0 = ', zq_rema(ji) 633 634 IF (to_print(ji) == 10) THEN 635 write(numout,*)'icethd_dh snow_extra: wfx_snw_sum_1d(ji) = ',wfx_snw_sum_1d(ji) 636 END IF 637 566 638 END DO 567 639 … … 610 682 rswitch = 1._wp - MAX( 0._wp , SIGN( 1._wp , - h_i_1d(ji) ) ) 611 683 t_su_1d(ji) = rswitch * t_su_1d(ji) + ( 1._wp - rswitch ) * rt0 684 685 IF (to_print(ji) == 10) THEN 686 write(numout,*)'icethd_dh: t_su_1d(ji), t_i_1d(ji,1) = ',t_su_1d(ji), ' ', t_i_1d(ji,1) 687 END IF 612 688 END DO 613 689 … … 625 701 ! recalculate t_s_1d from e_s_1d 626 702 t_s_1d(ji,jk) = rt0 + rswitch * ( - e_s_1d(ji,jk) * r1_rhos * r1_rcpi + rLfus * r1_rcpi ) 703 704 IF (to_print(ji) == 10) THEN 705 write(numout,*)'icethd_dh 10: t_s_1d(ji,jk), e_s_1d(ji,jk) = ',t_s_1d(ji,jk), ' ', e_s_1d(ji,jk) 706 END IF 627 707 END DO 628 708 END DO
Note: See TracChangeset
for help on using the changeset viewer.
