Changeset 14215 for NEMO/trunk/src/OCE/TRA/traqsr.F90
- Timestamp:
- 2020-12-18T14:49:22+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/TRA/traqsr.F90
r14189 r14215 151 151 ELSE ! No restart or Euler forward at 1st time step 152 152 z1_2 = 1._wp 153 DO_3D( is j, iej, isi, iei, 1, jpk )153 DO_3D( isi, iei, isj, iej, 1, jpk ) 154 154 qsr_hc_b(ji,jj,jk) = 0._wp 155 155 END_3D … … 157 157 ELSE !== Swap of qsr heat content ==! 158 158 z1_2 = 0.5_wp 159 DO_3D( is j, iej, isi, iei, 1, jpk )159 DO_3D( isi, iei, isj, iej, 1, jpk ) 160 160 qsr_hc_b(ji,jj,jk) = qsr_hc(ji,jj,jk) 161 161 END_3D … … 168 168 CASE( np_BIO ) !== bio-model fluxes ==! 169 169 ! 170 DO_3D( is j, iej, isi, iei, 1, nksr )170 DO_3D( isi, iei, isj, iej, 1, nksr ) 171 171 qsr_hc(ji,jj,jk) = r1_rho0_rcp * ( etot3(ji,jj,jk) - etot3(ji,jj,jk+1) ) 172 172 END_3D … … 190 190 ! most expensive calculations) 191 191 ! 192 DO_2D( is j, iej, isi, iei)192 DO_2D( isi, iei, isj, iej ) 193 193 ! zlogc = log(zchl) 194 194 zlogc = LOG ( MIN( 10. , MAX( 0.03, sf_chl(1)%fnow(ji,jj,1) ) ) ) … … 209 209 210 210 ! 211 DO_3D( is j, iej, isi, iei, 1, nksr + 1 )211 DO_3D( isi, iei, isj, iej, 1, nksr + 1 ) 212 212 ! zchl = ALOG( ze0(ji,jj) ) 213 213 zlogc = ze0(ji,jj) … … 239 239 ! 240 240 zcoef = ( 1. - rn_abs ) / 3._wp !* surface equi-partition in R-G-B 241 DO_2D( is j, iej, isi, iei)241 DO_2D( isi, iei, isj, iej ) 242 242 ze0(ji,jj) = rn_abs * qsr(ji,jj) 243 243 ze1(ji,jj) = zcoef * qsr(ji,jj) … … 250 250 ! 251 251 ! !* interior equi-partition in R-G-B depending on vertical profile of Chl 252 DO_3D( is j, iej, isi, iei, 2, nksr + 1 )252 DO_3D( isi, iei, isj, iej, 2, nksr + 1 ) 253 253 ze3t = e3t(ji,jj,jk-1,Kmm) 254 254 irgb = NINT( ztmp3d(ji,jj,jk) ) … … 264 264 END_3D 265 265 ! 266 DO_3D( is j, iej, isi, iei, 1, nksr ) !* now qsr induced heat content266 DO_3D( isi, iei, isj, iej, 1, nksr ) !* now qsr induced heat content 267 267 qsr_hc(ji,jj,jk) = r1_rho0_rcp * ( ztmp3d(ji,jj,jk) - ztmp3d(ji,jj,jk+1) ) 268 268 END_3D … … 274 274 zz0 = rn_abs * r1_rho0_rcp ! surface equi-partition in 2-bands 275 275 zz1 = ( 1. - rn_abs ) * r1_rho0_rcp 276 DO_3D( is j, iej, isi, iei, 1, nksr ) !* now qsr induced heat content276 DO_3D( isi, iei, isj, iej, 1, nksr ) !* now qsr induced heat content 277 277 zc0 = zz0 * EXP( -gdepw(ji,jj,jk ,Kmm)*xsi0r ) + zz1 * EXP( -gdepw(ji,jj,jk ,Kmm)*xsi1r ) 278 278 zc1 = zz0 * EXP( -gdepw(ji,jj,jk+1,Kmm)*xsi0r ) + zz1 * EXP( -gdepw(ji,jj,jk+1,Kmm)*xsi1r ) … … 292 292 ! 293 293 ! sea-ice: store the 1st ocean level attenuation coefficient 294 DO_2D( is j, iej, isi, iei)294 DO_2D( isi, iei, isj, iej ) 295 295 IF( qsr(ji,jj) /= 0._wp ) THEN ; fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rho0_rcp * qsr(ji,jj) ) 296 296 ELSE ; fraqsr_1lev(ji,jj) = 1._wp
Note: See TracChangeset
for help on using the changeset viewer.