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 2107 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2010-09-21T14:30:20+02:00 (14 years ago)
Author:
mlelod
Message:

cosmetic changes again

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1837_MLF/NEMO/OPA_SRC/DYN/sshwzv.F90

    r2068 r2107  
    156156      ! In forward Euler time stepping case, the same formulation as in the leap-frog case can be used 
    157157      ! because emp_b field is initialized with the vlaues of emp field. Hence, 0.5 * ( emp + emp_b ) = emp 
    158       IF( neuler == 0 .AND. kt == nit000 ) THEN 
    159          z1_rau0 = 1.e0 / rau0 
    160          ssha(:,:) = (  sshb(:,:) - z2dt * ( z1_rau0 *   emp(:,:)                + zhdiv(:,:) )  ) * tmask(:,:,1) 
    161       ELSE 
    162          z1_rau0 = 0.5 / rau0 
    163          ssha(:,:) = (  sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) )  ) * tmask(:,:,1) 
    164       ENDIF 
     158      z1_rau0 = 0.5 / rau0 
     159      ssha(:,:) = (  sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) )  ) * tmask(:,:,1) 
    165160 
    166161#if defined key_obc 
Note: See TracChangeset for help on using the changeset viewer.