- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/TOP/TRP/trcsbc.F90
r13237 r13899 18 18 USE oce_trc ! ocean dynamics and active tracers variables 19 19 USE trc ! ocean passive tracers variables 20 USE prtctl _trc! Print control for debbuging20 USE prtctl ! Print control for debbuging 21 21 USE iom 22 22 USE trd_oce … … 88 88 zfact = 0.5_wp 89 89 DO jn = 1, jptra 90 CALL iom_get( numrtr, jpdom_auto glo, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) ) ! before tracer content sbc90 CALL iom_get( numrtr, jpdom_auto, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) ) ! before tracer content sbc 91 91 END DO 92 92 ELSE ! No restart or restart not found: Euler forward time stepping … … 121 121 ! 122 122 DO jn = 1, jptra 123 DO_2D _01_00123 DO_2D( 0, 1, 0, 0 ) 124 124 sbc_trc(ji,jj,jn) = zsfx(ji,jj) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 125 125 END_2D … … 129 129 ! 130 130 DO jn = 1, jptra 131 DO_2D _01_00131 DO_2D( 0, 1, 0, 0 ) 132 132 sbc_trc(ji,jj,jn) = ( zsfx(ji,jj) + fmmflx(ji,jj) ) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 133 133 END_2D … … 137 137 ! 138 138 DO jn = 1, jptra 139 DO_2D _01_00139 DO_2D( 0, 1, 0, 0 ) 140 140 zse3t = 1. / e3t(ji,jj,1,Kmm) 141 141 ! tracer flux at the ice/ocean interface (tracer/m2/s) … … 161 161 IF( l_trdtrc ) ztrtrd(:,:,:) = ptr(:,:,:,jn,Krhs) ! save trends 162 162 ! 163 DO_2D _01_00163 DO_2D( 0, 1, 0, 0 ) 164 164 zse3t = zfact / e3t(ji,jj,1,Kmm) 165 165 ptr(ji,jj,1,jn,Krhs) = ptr(ji,jj,1,jn,Krhs) + ( sbc_trc_b(ji,jj,jn) + sbc_trc(ji,jj,jn) ) * zse3t … … 187 187 ! 188 188 IF( sn_cfctl%l_prttrc ) THEN 189 WRITE(charout, FMT="('sbc ')") ; CALL prt_ctl_ trc_info(charout)190 CALL prt_ctl _trc( tab4d=ptr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' )189 WRITE(charout, FMT="('sbc ')") ; CALL prt_ctl_info( charout, cdcomp = 'top' ) 190 CALL prt_ctl( tab4d_1=ptr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm, clinfo3='trd' ) 191 191 ENDIF 192 192 IF( l_trdtrc ) DEALLOCATE( ztrtrd )
Note: See TracChangeset
for help on using the changeset viewer.