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 2528 for trunk/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90 – NEMO

Ignore:
Timestamp:
2010-12-27T18:33:53+01:00 (13 years ago)
Author:
rblod
Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

    • Property svn:keywords set to Id
    r1271 r2528  
    1616 
    1717   !!---------------------------------------------------------------------- 
    18    !!   OPA 9.0 , LOCEAN-IPSL (2006) 
    19    !! $Id: agrif_opa_sponge.F90 1156 2008-06-26 16:06:45Z rblod $ 
    20    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
     18   !! NEMO/NST 3.3 , NEMO Consortium (2010) 
     19   !! $Id$ 
     20   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2121   !!---------------------------------------------------------------------- 
    2222 
     
    123123         DO jj = 1, jpjm1 
    124124            DO ji = 1, jpim1 
    125 #if defined key_zco 
    126                zabe1 = umask(ji,jj,jk) * spe1ur(ji,jj) 
    127                zabe2 = vmask(ji,jj,jk) * spe2vr(ji,jj) 
    128 #else 
    129125               zabe1 = umask(ji,jj,jk) * spe1ur(ji,jj) * fse3u(ji,jj,jk) 
    130126               zabe2 = vmask(ji,jj,jk) * spe2vr(ji,jj) * fse3v(ji,jj,jk) 
    131 #endif 
    132127               ztru(ji,jj,jk,jl) = zabe1 * ( trbdiff(ji+1,jj  ,jk,jl) - trbdiff(ji,jj,jk,jl) ) 
    133128               ztrv(ji,jj,jk,jl) = zabe2 * ( trbdiff(ji  ,jj+1,jk,jl) - trbdiff(ji,jj,jk,jl) ) 
     
    137132         DO jj = 2,jpjm1 
    138133            DO ji = 2,jpim1 
    139 #if defined key_zco 
    140                zbtr = spbtr2(ji,jj) 
    141 #else 
    142134               zbtr = spbtr2(ji,jj) / fse3t(ji,jj,jk) 
    143 #endif 
    144135               ! horizontal diffusive trends 
    145136               ztra = zbtr * (  ztru(ji,jj,jk,jl) - ztru(ji-1,jj,jk,jl)   & 
Note: See TracChangeset for help on using the changeset viewer.