Changeset 13642
- Timestamp:
- 2020-10-20T00:58:34+02:00 (4 years ago)
- Location:
- NEMO/releases/r4.0/r4.0-HEAD/src/ICE
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/ice1d.F90
r13284 r13642 145 145 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sst_1d 146 146 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: sss_1d 147 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: frq_m_1d 147 148 148 149 ! convergence check … … 225 226 ! 226 227 ii = ii + 1 227 ALLOCATE( sst_1d(jpij) , sss_1d(jpij) , STAT=ierr(ii) )228 ALLOCATE( sst_1d(jpij) , sss_1d(jpij) , frq_m_1d(jpij) , STAT=ierr(ii) ) 228 229 ! 229 230 ii = ii + 1 -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icethd.F90
r13640 r13642 442 442 CALL tab_2d_1d( npti, nptidx(1:npti), sst_1d(1:npti), sst_m ) 443 443 CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d(1:npti), sss_m ) 444 CALL tab_2d_1d( npti, nptidx(1:npti), frq_m_1d(1:npti), frq_m ) 444 445 ! 445 446 ! to update ice age -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icethd_dh.F90
r13589 r13642 139 139 ! 140 140 DO ji = 1, npti 141 zf_tt(ji) = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji) 141 zf_tt(ji) = qcn_ice_bot_1d(ji) + qsb_ice_bot_1d(ji) + fhld_1d(ji) + qtr_ice_bot_1d(ji) * frq_m_1d(ji) 142 142 zq_bot(ji) = MAX( 0._wp, zf_tt(ji) * rdt_ice ) 143 143 END DO -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/iceupdate.F90
r13589 r13642 138 138 ! solar flux transmitted thru the 1st level of the ocean (i.e. not used by sea-ice) 139 139 qsr(ji,jj) = ( 1._wp - at_i_b(ji,jj) ) * qsr_oce(ji,jj) * ( 1._wp - frq_m(ji,jj) ) & 140 ! + solar flux transmitted thru ice (also not used by sea-ice)141 & + SUM( a_i_b(ji,jj,:) * qtr_ice_bot(ji,jj,:) ) 140 ! + solar flux transmitted thru ice and the 1st ocean level (also not used by sea-ice) 141 & + SUM( a_i_b(ji,jj,:) * qtr_ice_bot(ji,jj,:) ) * ( 1._wp - frq_m(ji,jj) ) 142 142 ! 143 143 ELSE !-- cooling or no ice left
Note: See TracChangeset
for help on using the changeset viewer.