Changeset 2107 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/DYN/sshwzv.F90
- Timestamp:
- 2010-09-21T14:30:20+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1837_MLF/NEMO/OPA_SRC/DYN/sshwzv.F90
r2068 r2107 156 156 ! In forward Euler time stepping case, the same formulation as in the leap-frog case can be used 157 157 ! 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) 165 160 166 161 #if defined key_obc
Note: See TracChangeset
for help on using the changeset viewer.