- 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/traqsr.F90
r15027 r15028 244 244 ENDIF 245 245 ! 246 !!st IF( l_trdtra ) THEN ! qsr tracers trends saved for diagnostics247 !!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 ENDIF251 246 ! ! print mean trends (used for debugging) 252 247 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' qsr - Ta: ', mask1=tmask, clinfo3='tra-ta' ) … … 308 303 ENDIF 309 304 ! 310 DO_2D( 0, 0, 0, 0 ) ! pre-calculated expensive coefficient305 DO_2D( 0, 0, 0, 0 ) ! pre-calculated expensive coefficient 311 306 zlogc = LOG( MAX( 0.03_wp, MIN( sf_chl(1)%fnow(ji,jj,1) ,10._wp ) ) ) ! zlogc = log(zchl) with 0.03 <= Chl >= 10. 312 307 zc1 = 0.113328685307 + 0.803 * zlogc ! zc1 : log(zCze) = log (1.12 * zchl**0.803) … … 371 366 zeT(ji,jj) = zzeT ! total - - - 372 367 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 ! 378 369 END DO 379 370 !
Note: See TracChangeset
for help on using the changeset viewer.