Changeset 11099 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/TOP/TRP/trcatf.F90
- Timestamp:
- 2019-06-11T15:59:58+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/TOP/TRP/trcatf.F90
r11057 r11099 30 30 USE trd_oce 31 31 USE trdtra 32 USE tra nxt32 USE traatf 33 33 USE bdy_oce , ONLY: ln_bdy 34 34 USE trcbdy ! BDY open boundaries … … 143 143 ENDIF 144 144 ! ! Leap-Frog + Asselin filter time stepping 145 IF( .NOT.( (neuler == 0 .AND. kt == nittrc000) .OR. ln_top_euler ) ) THEN ! Only time filter if not an Euler timestep 145 IF( (neuler == 0 .AND. kt == nittrc000) .OR. ln_top_euler ) THEN ! Euler time-stepping 146 ! 147 IF (l_trdtrc .AND. .NOT. ln_linssh ) THEN ! Zero Asselin filter contribution must be explicitly written out since for vvl 148 ! ! Asselin filter is output by tra_nxt_vvl that is not called on this time step 149 ztrdt(:,:,:,:) = 0._wp 150 DO jn = 1, jptra 151 CALL trd_tra( kt, Kmm, Kaa, 'TRC', jn, jptra_atf, ztrdt(:,:,:,jn) ) 152 ENDDO 153 END IF 154 ! 155 ELSE 146 156 IF( .NOT. l_offline ) THEN ! Leap-Frog + Asselin filter time stepping 147 IF( ln_linssh ) THEN ; CALL tra_ nxt_fix( kt, Kbb, Kmm, Kaa, nittrc000, 'TRC', ptr, jptra ) ! linear ssh148 ELSE ; CALL tra_ nxt_vvl( kt, Kbb, Kmm, Kaa, nittrc000, rdttrc, 'TRC', ptr, sbc_trc, sbc_trc_b, jptra ) ! non-linear ssh157 IF( ln_linssh ) THEN ; CALL tra_atf_fix( kt, Kbb, Kmm, Kaa, nittrc000, 'TRC', ptr, jptra ) ! linear ssh 158 ELSE ; CALL tra_atf_vvl( kt, Kbb, Kmm, Kaa, nittrc000, rdttrc, 'TRC', ptr, sbc_trc, sbc_trc_b, jptra ) ! non-linear ssh 149 159 ENDIF 150 160 ELSE 151 CALL trc_atf_off( kt, Kbb, Kmm, Kaa ) ! offline161 CALL trc_atf_off( kt, Kbb, Kmm, Kaa, ptr ) ! offline 152 162 ENDIF 153 163 ! … … 260 270 !! Default option Empty module 261 271 !!---------------------------------------------------------------------- 272 USE par_oce 273 USE par_trc 262 274 CONTAINS 263 275 SUBROUTINE trc_atf( kt, Kbb, Kmm, Kaa, ptr ) 264 INTEGER , INTENT(in) :: kt265 INTEGER , INTENT(in) :: Kbb, Kmm, Kaa! time level indices266 REAL(wp), DIMENSION(jpi,jpj,jpk,jptra,jpt), INTENT(inout) :: ptr ! passive tracers276 INTEGER , INTENT(in) :: kt 277 INTEGER, INTENT(in ) :: Kbb, Kmm, Kaa ! time level indices 278 REAL(wp), DIMENSION(jpi,jpj,jpk,jptra,jpt), INTENT(inout) :: ptr ! passive tracers and RHS of tracer equation 267 279 WRITE(*,*) 'trc_atf: You should not have seen this print! error?', kt 268 280 END SUBROUTINE trc_atf
Note: See TracChangeset
for help on using the changeset viewer.