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 15189 for NEMO/branches/2021 – NEMO

Changeset 15189 for NEMO/branches/2021


Ignore:
Timestamp:
2021-08-13T14:51:28+02:00 (3 years ago)
Author:
techene
Message:

#2715 remove advection & time filtering and add domain lateral boundary management (previously done in trc_atf)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/TOP/TRP/trctrp.F90

    r14086 r15189  
    8282         IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_trc       ! tracers sponge 
    8383#endif 
    84                                 CALL trc_adv    ( kt, Kbb, Kmm, tr, Krhs )      ! horizontal & vertical advection  
     84#if ! defined key_RK3 
     85         !                                                 ! MLF only: add the advection trend to the RHS 
     86                                CALL trc_adv    ( kt, Kbb, Kmm, tr, Krhs )      ! horizontal & vertical advection 
     87#endif 
    8588                                CALL trc_ldf    ( kt, Kbb, Kmm,       tr, Krhs )  ! lateral mixing 
    8689                                CALL trc_zdf    ( kt, Kbb, Kmm, Krhs, tr, Kaa  )  ! vert. mixing & after tracer   ==> after 
    87                                 CALL trc_atf    ( kt, Kbb, Kmm, Kaa , tr )        ! time filtering of "now" tracer fields     
     90#if defined key_RK3 
     91         !                                                 ! RK3: only manage lateral boundary 
     92# if defined key_agrif 
     93                                CALL Agrif_trc  ( kt )                            ! AGRIF zoom boundaries 
     94# endif 
     95         !                                                                        ! Update after tracer on domain lateral boundaries 
     96                                CALL lbc_lnk( 'stprk3_stg', tr(:,:,:,:,Kaa), 'T', 1._wp )    
     97         ! 
     98         IF( ln_bdy )           CALL trc_bdy    ( kt, Kbb, Kmm, Kaa ) 
     99#else 
     100         !                                                 ! MLF: apply Asselin time filter and manage lateral boundary 
     101                                CALL trc_atf    ( kt, Kbb, Kmm, Kaa , tr )        ! time filtering of "now" tracer fields 
     102#endif 
    88103         ! 
    89104         ! Subsequent calls use the filtered values: Kmm and Kaa  
Note: See TracChangeset for help on using the changeset viewer.