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 4610 – NEMO

Changeset 4610


Ignore:
Timestamp:
2014-03-31T15:19:34+02:00 (10 years ago)
Author:
cetlod
Message:

bugfix : Add the MLE and ztilde parameterisation for passive tracers transport, see ticket #1290

Location:
trunk/NEMOGCM/NEMO/TOP_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90

    r3718 r4610  
    2424   USE traadv_qck      ! QUICKEST scheme           (tra_adv_qck    routine) 
    2525   USE traadv_eiv      ! eddy induced velocity     (tra_adv_eiv    routine) 
     26   USE traadv_mle      ! ML eddy induced velocity  (tra_adv_mle    routine) 
    2627   USE ldftra_oce      ! lateral diffusion coefficient on tracers 
    2728   USE prtctl_trc      ! Print control 
     
    100101         ! 
    101102      END DO 
    102       zwn(:,:,jpk) = 0.e0                                 ! no transport trough the bottom 
     103      ! 
     104      IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 
     105         zun(:,:,:) = zun(:,:,:) + un_td(:,:,:) 
     106         zvn(:,:,:) = zvn(:,:,:) + vn_td(:,:,:) 
     107      ENDIF 
     108      ! 
     109      zun(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom 
     110      zvn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom 
     111      zwn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom 
    103112 
    104113      IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif )   &  ! add the eiv transport (if necessary) 
    105114         &              CALL tra_adv_eiv( kt, nittrc000, zun, zvn, zwn, 'TRC' ) 
     115      ! 
     116      IF( ln_mle    )   CALL tra_adv_mle( kt, nittrc000, zun, zvn, zwn, 'TRC' )    ! add the mle transport (if necessary) 
    106117      ! 
    107118      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
  • trunk/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90

    r4306 r4610  
    5151   USE dom_oce  
    5252 
     53   USE domvvl, ONLY : un_td, vn_td          !: thickness diffusion transport 
     54   USE domvvl, ONLY : ln_vvl_ztilde         !: ztilde vertical coordinate 
     55   USE domvvl, ONLY : ln_vvl_layer          !: level  vertical coordinate 
    5356 
    5457   !* ocean fields: here now and after fields * 
     
    7174   USE oce , ONLY :   ssha    =>    ssha    !: sea surface height at t-point [m]    
    7275   USE oce , ONLY :   l_traldf_rot => l_traldf_rot  !: rotated laplacian operator for lateral diffusion 
    73 #if defined key_offline 
    74    USE oce , ONLY :   gtsu    =>    gtsu    !: t-, s- and rd horizontal gradient at u- and 
    75    USE oce , ONLY :   gtsv    =>    gtsv    !: 
    76    USE oce , ONLY :   gru     =>    gru     !: 
    77    USE oce , ONLY :   grv     =>    grv     !:  
    78 #endif 
    7976 
    8077   !* surface fluxes * 
Note: See TracChangeset for help on using the changeset viewer.