- Timestamp:
- 2017-10-18T19:14:32+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/LIM_SRC_3/icethd_dh.F90
r8586 r8637 148 148 ! Contribution to mass flux 149 149 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) + rhosn * h_s_1d(ji) * a_i_1d(ji) * r1_rdtice 150 dh_s_mlt(ji) = dh_s_mlt(ji) - h_s_1d(ji) 150 151 ! updates 151 h_s_1d(ji) = 0._wp152 h_s_1d(ji) = 0._wp 152 153 e_s_1d (ji,1) = 0._wp 153 154 t_s_1d (ji,1) = rt0 … … 211 212 ! snow melting only = water into the ocean (then without snow precip), >0 212 213 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) - rhosn * a_i_1d(ji) * zdeltah(ji,1) * r1_rdtice 214 dh_s_mlt(ji) = dh_s_mlt(ji) + zdeltah(ji,1) 213 215 ! updates available heat + precipitations after melting 214 216 zq_su (ji) = MAX( 0._wp , zq_su (ji) + zdeltah(ji,1) * zqprec(ji) ) … … 233 235 hfx_snw_1d(ji) = hfx_snw_1d(ji) - zdeltah(ji,jk) * a_i_1d(ji) * e_s_1d(ji,jk) * r1_rdtice 234 236 ! snow melting only = water into the ocean (then without snow precip) 235 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) - rhosn * a_i_1d(ji) * zdeltah(ji,jk) * r1_rdtice 237 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) - rhosn * a_i_1d(ji) * zdeltah(ji,jk) * r1_rdtice 238 dh_s_mlt(ji) = dh_s_mlt(ji) + zdeltah(ji,jk) 236 239 ! updates available heat + thickness 237 zq_su (ji) 240 zq_su (ji) = MAX( 0._wp , zq_su (ji) + zdeltah(ji,jk) * e_s_1d(ji,jk) ) 238 241 h_s_1d(ji) = MAX( 0._wp , h_s_1d(ji) + zdeltah(ji,jk) ) 239 242 END DO … … 571 574 zdeltah (ji,1) = MIN( 0._wp , MAX( zdeltah(ji,1) , - h_s_1d(ji) ) ) ! bound melting 572 575 dh_s_tot (ji) = dh_s_tot(ji) + zdeltah(ji,1) 573 h_s_1d (ji) = h_s_1d(ji) + zdeltah(ji,1)576 h_s_1d (ji) = h_s_1d(ji) + zdeltah(ji,1) 574 577 575 578 zq_rema(ji) = zq_rema(ji) + zdeltah(ji,1) * e_s_1d(ji,1) ! update available heat (J.m-2) … … 577 580 hfx_snw_1d(ji) = hfx_snw_1d(ji) - zdeltah(ji,1) * a_i_1d(ji) * e_s_1d(ji,1) * r1_rdtice ! W.m-2 (>0) 578 581 ! Contribution to mass flux 579 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) - rhosn * a_i_1d(ji) * zdeltah(ji,1) * r1_rdtice 582 wfx_snw_sum_1d(ji) = wfx_snw_sum_1d(ji) - rhosn * a_i_1d(ji) * zdeltah(ji,1) * r1_rdtice 583 dh_s_mlt(ji) = dh_s_mlt(ji) + zdeltah(ji,1) 580 584 ! 581 585 ! Remaining heat flux (W.m-2) is sent to the ocean heat budget … … 611 615 612 616 ! Case constant salinity in time: virtual salt flux to keep salinity constant 613 IF( nn_icesal == 1 .OR. nn_icesal == 3) THEN617 IF( nn_icesal /= 2 ) THEN 614 618 sfx_bri_1d(ji) = sfx_bri_1d(ji) - sss_1d (ji) * a_i_1d(ji) * zfmdt * r1_rdtice & ! put back sss_m into the ocean 615 619 & - s_i_1d(ji) * a_i_1d(ji) * dh_snowice(ji) * rhoic * r1_rdtice ! and get rn_icesal from the ocean
Note: See TracChangeset
for help on using the changeset viewer.