Changeset 719 for trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
- Timestamp:
- 2007-10-16T16:59:56+02:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
- Property svn:keywords changed from Id to Author Date Id Revision
r710 r719 14 14 !!---------------------------------------------------------------------- 15 15 !! OPA 9.0 , LOCEAN-IPSL (2005) 16 !! $ Id$16 !! $Header$ 17 17 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 18 18 !!---------------------------------------------------------------------- … … 20 20 USE oce ! ocean dynamics and tracers 21 21 USE dom_oce ! ocean space and time domain 22 USE sbc_oce ! surface boundary condition: ocean22 USE phycst ! physical constants 23 23 USE zdf_oce ! ocean vertical physics 24 USE phycst ! physical constants25 24 USE in_out_manager ! I/O manager 25 USE taumod ! surface ocean stress 26 26 27 27 IMPLICIT NONE … … 36 36 !!---------------------------------------------------------------------- 37 37 !! OPA 9.0 , LOCEAN-IPSL (2005) 38 !! $ Id$38 !! $Header$ 39 39 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 40 40 !!---------------------------------------------------------------------- … … 141 141 !!! change les resultats (derniers digit, pas significativement + rapide 1* de moins) 142 142 !!! ua(ji,jj,1) = ub(ji,jj,1) & 143 !!! + p2dt * ( ua(ji,jj,1) + utau(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) )143 !!! + p2dt * ( ua(ji,jj,1) + taux(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) ) 144 144 z2dtf = p2dt / ( fse3u(ji,jj,1)*rau0 ) 145 145 ua(ji,jj,1) = ub(ji,jj,1) & 146 + p2dt * ua(ji,jj,1) + z2dtf * utau(ji,jj)146 + p2dt * ua(ji,jj,1) + z2dtf * taux(ji,jj) 147 147 END DO 148 148 END DO … … 236 236 !!! change les resultats (derniers digit, pas significativement + rapide 1* de moins) 237 237 !!! va(ji,jj,1) = vb(ji,jj,1) & 238 !!! + p2dt * ( va(ji,jj,1) + vtau(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) )238 !!! + p2dt * ( va(ji,jj,1) + tauy(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) ) 239 239 z2dtf = p2dt / ( fse3v(ji,jj,1)*rau0 ) 240 240 va(ji,jj,1) = vb(ji,jj,1) & 241 + p2dt * va(ji,jj,1) + z2dtf * vtau(ji,jj)241 + p2dt * va(ji,jj,1) + z2dtf * tauy(ji,jj) 242 242 END DO 243 243 END DO
Note: See TracChangeset
for help on using the changeset viewer.