- Timestamp:
- 2021-06-19T10:53:10+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traatf_qco.F90
r14618 r15028 203 203 DO jn = 1, kjpt 204 204 ! 205 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 205 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 206 !!st DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 206 207 ztn = pt(ji,jj,jk,jn,Kmm) 207 208 ztd = pt(ji,jj,jk,jn,Kaa) - 2._wp * ztn + pt(ji,jj,jk,jn,Kbb) ! time laplacian on tracers … … 262 263 IF( ( l_trdtra .AND. cdtype == 'TRA' ) .OR. ( l_trdtrc .AND. cdtype == 'TRC' ) ) THEN 263 264 ALLOCATE( ztrd_atf(jpi,jpj,jpk,kjpt) ) 264 ztrd_atf(:,:,:,:) = 0. 0_wp265 ztrd_atf(:,:,:,:) = 0._wp 265 266 ENDIF 266 267 zfact = 1._wp / p2dt … … 268 269 zfact2 = zfact1 * r1_rho0 269 270 DO jn = 1, kjpt 270 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 271 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 272 !!st DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 271 273 ze3t_b = e3t(ji,jj,jk,Kbb) 272 274 ze3t_n = e3t(ji,jj,jk,Kmm) … … 292 294 IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & 293 295 & ztc_f = ztc_f - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) ) 294 296 ! 295 297 ! 296 298 IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) &
Note: See TracChangeset
for help on using the changeset viewer.