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 10272 for NEMO/trunk/src/OCE/DYN/dynspg_ts.F90 – NEMO

Ignore:
Timestamp:
2018-11-02T12:15:39+01:00 (6 years ago)
Author:
deazer
Message:

Bug fix for wd case with implicit bottom friction.
(Currently only done for ln_dynadv_vec)
Trivial, just move lines so that the are done in all advection cases

see Ticket:

https://forge.ipsl.jussieu.fr/nemo/ticket/2149

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/DYN/dynspg_ts.F90

    r10256 r10272  
    11141114                            &   ) * ssvmask(ji,jj) 
    11151115  
    1116 !jth implicit bottom friction: 
    1117                   IF ( ll_wd ) THEN ! revert to explicit for bit comparison tests in non wad runs 
    1118                      ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rdtbt * zCdU_u(ji,jj) * hur_e(ji,jj)) 
    1119                      va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rdtbt * zCdU_v(ji,jj) * hvr_e(ji,jj)) 
    1120                   ENDIF 
    1121  
    11221116               END DO 
    11231117            END DO 
     
    11441138                            &               +    hv_n(ji,jj)  * zv_frc(ji,jj) ) & 
    11451139                            &   ) * zhvra 
     1140               END DO 
     1141            END DO 
     1142         ENDIF 
     1143!jth implicit bottom friction: 
     1144         IF ( ll_wd ) THEN ! revert to explicit for bit comparison tests in non wad runs 
     1145            DO jj = 2, jpjm1 
     1146               DO ji = fs_2, fs_jpim1   ! vector opt. 
     1147                     ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rdtbt * zCdU_u(ji,jj) * hur_e(ji,jj)) 
     1148                     va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rdtbt * zCdU_v(ji,jj) * hvr_e(ji,jj)) 
    11461149               END DO 
    11471150            END DO 
Note: See TracChangeset for help on using the changeset viewer.