Ignore:
Timestamp:
2018-09-07T15:20:38+02:00 (3 years ago)
Author:
cbricaud
Message:

clean trcadv_crs.F90 and add ubs scheme for coarsenig

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv_crs.F90

    r6772 r10104  
    1818   USE trc             ! ocean passive tracers variables 
    1919   USE trcnam_trp      ! passive tracers transport namelist variables 
    20    USE traadv_cen2     ! 2nd order centered scheme (tra_adv_cen2   routine) 
    2120   USE traadv_tvd_crs  ! TVD      scheme           (tra_adv_tvd    routine) 
    22    USE traadv_muscl    ! MUSCL    scheme           (tra_adv_muscl  routine) 
    23    USE traadv_muscl2   ! MUSCL2   scheme           (tra_adv_muscl2 routine) 
    24    USE traadv_ubs      ! UBS      scheme           (tra_adv_ubs    routine) 
    25    USE traadv_qck      ! QUICKEST scheme           (tra_adv_qck    routine) 
    26    USE traadv_eiv      ! eddy induced velocity     (tra_adv_eiv    routine) 
     21   USE traadv_ubs_crs  ! TVD      scheme           (tra_adv_tvd    routine) 
    2722   USE ldftra_oce      ! lateral diffusion coefficient on tracers 
    2823   USE prtctl_trc      ! Print control 
     
    108103      zwn(:,:,jpk) = 0.e0                                 ! no transport trough the bottom 
    109104 
    110       IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif )   &  ! add the eiv transport (if necessary) 
    111          &              CALL tra_adv_eiv( kt, nittrc000, zun, zvn, zwn, 'TRC' ) 
    112105      ! 
    113106      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    114107      CASE ( 2 )   ;    CALL tra_adv_tvd_crs   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  TVD  
     108      CASE ( 5 )   ;    CALL tra_adv_ubs_crs   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  UBS  
    115109      ! 
    116110      CASE (-1 )                                      !==  esopa: test all possibility with control print  ==! 
Note: See TracChangeset for help on using the changeset viewer.