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/traqsr.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/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.