- Timestamp:
- 2017-07-25T16:38:38+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r8342 r8369 631 631 IF( nn_icesal == 1 .OR. nn_icesal == 3 ) THEN 632 632 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 633 & - sm_i_1d(ji) 633 & - sm_i_1d(ji) * a_i_1d(ji) * dh_snowice(ji) * rhoic * r1_rdtice ! and get rn_icesal from the ocean 634 634 ENDIF 635 635 … … 665 665 666 666 ! --- ensure that a_i = 0 where ht_i = 0 --- 667 WHERE( ht_i_1d == 0._wp ) a_i_1d = 0._wp 668 667 DO ji = 1, nidx 668 IF( ht_i_1d(ji) == 0._wp ) a_i_1d(ji) = 0._wp 669 END DO 670 669 671 CALL wrk_dealloc( jpij, zqprec, zq_su, zq_bo, zf_tt, zq_rema, zsnw, zevap_rema ) 670 672 CALL wrk_dealloc( jpij, zdh_s_mel, zdh_s_pre, zdh_s_sub, zeh_i )
Note: See TracChangeset
for help on using the changeset viewer.