- Timestamp:
- 2019-06-06T16:29:54+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GO6_mixing/src/ICE/icethd_zdf_bl99.F90
r10888 r11083 206 206 ! 207 207 l_T_converged(:) = .FALSE. 208 ! !============================!209 208 ! Convergence calculated until all sub-domain grid points have converged 210 209 ! Calculations keep going for all grid points until sub-domain convergence (vectorisation optimisation) 211 210 ! but values are not taken into account (results independant of MPI partitioning) 212 211 ! 212 ! !============================! 213 213 DO WHILE ( ( .NOT. ALL (l_T_converged(1:npti)) ) .AND. iconv < iconv_max ) ! Iterative procedure begins ! 214 ! !============================!214 ! !============================! 215 215 iconv = iconv + 1 216 216 ! … … 742 742 zdti_max = MAX ( zdti_max , MAXVAL( ABS( t_s_1d(ji,1:nlay_s) - ztsb(ji,1:nlay_s) ) ) ) 743 743 ! t_i 744 DO jk = 0, nlay_i744 DO jk = 1, nlay_i 745 745 ztmelts = -rTmlt * sz_i_1d(ji,jk) + rt0 746 746 t_i_1d(ji,jk) = MAX( MIN( t_i_1d(ji,jk), ztmelts ), rt0 - 100._wp ) … … 856 856 t_i_1d (1:npti,:) = ztiold (1:npti,:) 857 857 qcn_ice_1d(1:npti) = qcn_ice_top_1d(1:npti) 858 859 !!clem 860 ! remettre t_su_1d, qns_ice_1d et dqns_ice_1d comme avant puisqu'on devrait faire comme si on avant conduction = input 861 !clem 858 862 ENDIF 859 863 !
Note: See TracChangeset
for help on using the changeset viewer.