Changeset 15028


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

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

Location:
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA
Files:
2 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) )                                          & 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90

    r15027 r15028  
    244244      ENDIF 
    245245      ! 
    246 !!st      IF( l_trdtra ) THEN     ! qsr tracers trends saved for diagnostics 
    247 !!st         ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) - ztrdt(:,:,:) 
    248 !!st         CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_tem, jptra_qsr, ztrdt ) 
    249 !!st         DEALLOCATE( ztrdt ) 
    250 !!st      ENDIF 
    251246      !                       ! print mean trends (used for debugging) 
    252247      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
     
    308303      ENDIF 
    309304      ! 
    310       DO_2D( 0, 0, 0, 0 )                          ! pre-calculated expensive coefficient 
     305       DO_2D( 0, 0, 0, 0 )                          ! pre-calculated expensive coefficient 
    311306         zlogc = LOG(  MAX( 0.03_wp, MIN( sf_chl(1)%fnow(ji,jj,1) ,10._wp ) )  ) ! zlogc = log(zchl)   with 0.03 <= Chl >= 10.  
    312307         zc1   = 0.113328685307 + 0.803 * zlogc                               ! zc1 : log(zCze)  = log (1.12  * zchl**0.803) 
     
    371366            zeT(ji,jj) = zzeT                                   ! total          -        -      - 
    372367         END_2D 
    373 !!stbug         IF( jk == 1 ) THEN                        !- sea-ice : store the 1st level attenuation coefficient 
    374 !!stbug            WHERE( qsr(A2D(0)) /= 0._wp )   ;   fraqsr_1lev(A2D(0)) = 1._wp - zeT(A2D(0)) / qsr(A2D(0)) 
    375 !!stbug            ELSEWHERE                       ;   fraqsr_1lev(A2D(0)) = 1._wp 
    376 !!stbug            END WHERE 
    377 !!stbug         ENDIF 
     368         ! 
    378369      END DO 
    379370      ! 
Note: See TracChangeset for help on using the changeset viewer.