- Timestamp:
- 2013-12-17T17:33:09+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r4328 r4338 110 110 END DO 111 111 ! ! 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 ! 114 115 z1_rau0 = 0.5_wp * r1_rau0 115 116 ssha(:,:) = ( sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) ) ) * tmask(:,:,1) … … 197 198 DO jk = 1, jpkm1 198 199 ! horizontal divergence of thickness diffusion transport ( velocity multiplied by e3t) 199 ! - ML - note: computation al lready 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) 200 201 DO jj = 2, jpjm1 201 202 DO ji = fs_2, fs_jpim1 ! vector opt.
Note: See TracChangeset
for help on using the changeset viewer.