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 3952 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90 – NEMO

Ignore:
Timestamp:
2013-07-02T15:20:13+02:00 (11 years ago)
Author:
flavoni
Message:

ok for LIM2, still reproducibility error on LIM3 case

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

    r3938 r3952  
    176176                     zerp = zsrp * ( 1. - 2.*rnfmsk(ji,jj) )   &      ! No damping in vicinity of river mouths 
    177177                        &        * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj,1) )    
     178#if defined key_lim3 
    178179                     emps(ji,jj) = emps(ji,jj) + zerp 
     180#else 
     181                     emps(ji,jj) = emps(ji,jj) + zerp / MAX( sss_m(ji,jj) , 1.e-20 )      ! lim2 case emps is a mass flux 
     182#endif 
    179183                     erp (ji,jj) = zerp / MAX( sss_m(ji,jj), 1.e-20 ) ! converted into an equivalent volume flux (diagnostic only) 
    180184                  END DO 
     
    193197                     IF( ln_sssr_bnd )   zerp = SIGN( 1., zerp ) * MIN( zerp_bnd, ABS(zerp) ) 
    194198                     emp(ji,jj) = emp (ji,jj) + zerp 
     199#if defined key_lim2 
     200                     emps(ji,jj) = emps(ji,jj) + zerp      ! lim2 case: emps is a mass flux 
     201#endif 
    195202                     erp(ji,jj) = zerp 
    196203                  END DO 
Note: See TracChangeset for help on using the changeset viewer.