New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15028 for NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traatf_qco.F90 – NEMO

Ignore:
Timestamp:
2021-06-19T10:53:10+02:00 (3 years ago)
Author:
techene
Message:

#2605 : additionnal cleaning sette ~OK (as stated previously) for MLF and RK3

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traatf_qco.F90

    r14618 r15028  
    203203      DO jn = 1, kjpt 
    204204         ! 
    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 ) 
    206207            ztn = pt(ji,jj,jk,jn,Kmm) 
    207208            ztd = pt(ji,jj,jk,jn,Kaa) - 2._wp * ztn + pt(ji,jj,jk,jn,Kbb)  ! time laplacian on tracers 
     
    262263      IF( ( l_trdtra .AND. cdtype == 'TRA' ) .OR. ( l_trdtrc .AND. cdtype == 'TRC' ) )   THEN 
    263264         ALLOCATE( ztrd_atf(jpi,jpj,jpk,kjpt) ) 
    264          ztrd_atf(:,:,:,:) = 0.0_wp 
     265         ztrd_atf(:,:,:,:) = 0._wp 
    265266      ENDIF 
    266267      zfact = 1._wp / p2dt 
     
    268269      zfact2 = zfact1 * r1_rho0 
    269270      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 ) 
    271273            ze3t_b = e3t(ji,jj,jk,Kbb) 
    272274            ze3t_n = e3t(ji,jj,jk,Kmm) 
     
    292294            IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr )                            & 
    293295               &     ztc_f  = ztc_f  - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) ) 
    294                ! 
     296            ! 
    295297            ! 
    296298            IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) )                                          & 
Note: See TracChangeset for help on using the changeset viewer.