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 4338 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2013-12-17T17:33:09+01:00 (10 years ago)
Author:
acc
Message:

Branch dev_MERGE_2013. #1209. Fix to compute after scale factors before wn computation in time-splitting case. dom_vvl_sf_nxt is now called twice but correctly computes the baroclinic contribution only once. Still need to sort out the asselin filtering of the separate components

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r4328 r4338  
    110110      END DO 
    111111      !                                                ! Sea surface elevation time stepping 
    112       ! In forward Euler time stepping case, the same formulation as in the leap-frog case can be used 
    113       ! because emp_b field is initialized with the vlaues of emp field. Hence, 0.5 * ( emp + emp_b ) = emp 
     112      ! In time-split case we need a first guess of the ssh after (using the baroclinic timestep) in order to 
     113      ! compute the vertical velocity which can be used to compute the non-linear terms of the momentum equations. 
     114      !  
    114115      z1_rau0 = 0.5_wp * r1_rau0 
    115116      ssha(:,:) = (  sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) )  ) * tmask(:,:,1) 
     
    197198         DO jk = 1, jpkm1 
    198199            ! horizontal divergence of thickness diffusion transport ( velocity multiplied by e3t) 
    199             ! - ML - note: computation allready done in dom_vvl_sf_nxt. Could be optimized (not critical and clearer this way) 
     200            ! - ML - note: computation already done in dom_vvl_sf_nxt. Could be optimized (not critical and clearer this way) 
    200201            DO jj = 2, jpjm1 
    201202               DO ji = fs_2, fs_jpim1   ! vector opt. 
Note: See TracChangeset for help on using the changeset viewer.