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 5878 for branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90 – NEMO

Ignore:
Timestamp:
2015-11-11T19:11:01+01:00 (8 years ago)
Author:
deazer
Message:

Adds simple ramp of diffusion coefficient in BDY rim using TANH weighting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90

    r5876 r5878  
    3232   USE wrk_nemo        ! Memory Allocation 
    3333   USE timing          ! Timing 
     34#if defined key_bdy  
     35   USE bdy_oce  
     36#endif  
    3437 
    3538   IMPLICIT NONE 
     
    110113      REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d 
    111114      REAL(wp), POINTER, DIMENSION(:,:,:) ::  zdkt, zdk1t, zdit, zdjt, ztfw  
     115      REAL(wp), DIMENSION(jpi,jpj)     ::   zfactor  ! multiplier for diffusion 
    112116      !!---------------------------------------------------------------------- 
    113117      ! 
     
    124128      ENDIF 
    125129      ! 
     130#if defined key_bdy 
     131      zfactor(:,:) = sponge_factor(:,:) 
     132#else 
     133      zfactor(:,:) = 1.0 
     134#endif 
    126135      !                                                          ! =========== 
    127136      DO jn = 1, kjpt                                            ! tracer loop 
     
    200209            DO jj = 1 , jpjm1 
    201210               DO ji = 1, fs_jpim1   ! vector opt. 
    202                   zabe1 = ( fsahtu(ji,jj,jk) + pahtb0 ) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
    203                   zabe2 = ( fsahtv(ji,jj,jk) + pahtb0 ) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
     211                  zabe1 = zfactor(ji,jj) * ( fsahtu(ji,jj,jk) + pahtb0 ) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
     212                  zabe2 = zfactor(ji,jj) * ( fsahtv(ji,jj,jk) + pahtb0 ) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
    204213                  ! 
    205214                  zmsku = 1. / MAX(  tmask(ji+1,jj,jk  ) + tmask(ji,jj,jk+1)   & 
Note: See TracChangeset for help on using the changeset viewer.