- Timestamp:
- 2015-10-01T14:48:08+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r5737 r5770 80 80 IF( nn_timing == 1 ) CALL timing_start('ssh_nxt') 81 81 ! 82 CALL wrk_alloc( jpi, jpj,zhdiv )82 CALL wrk_alloc( jpi,jpj, zhdiv ) 83 83 ! 84 84 IF( kt == nit000 ) THEN 85 !86 85 IF(lwp) WRITE(numout,*) 87 86 IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height' 88 87 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 89 !90 88 ENDIF 91 89 ! … … 159 157 !! Reference : Leclair, M., and G. Madec, 2009, Ocean Modelling. 160 158 !!---------------------------------------------------------------------- 161 ! 162 INTEGER, INTENT(in) :: kt ! time step 159 INTEGER, INTENT(in) :: kt ! time step 160 ! 161 INTEGER :: ji, jj, jk ! dummy loop indices 162 REAL(wp) :: z1_2dt ! local scalars 163 163 REAL(wp), POINTER, DIMENSION(:,: ) :: z2d 164 164 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3d, zhdiv 165 ! 166 INTEGER :: ji, jj, jk ! dummy loop indices 167 REAL(wp) :: z1_2dt ! local scalars 168 !!---------------------------------------------------------------------- 169 165 !!---------------------------------------------------------------------- 166 ! 170 167 IF( nn_timing == 1 ) CALL timing_start('wzv') 171 168 ! 172 169 IF( kt == nit000 ) THEN 173 !174 170 IF(lwp) WRITE(numout,*) 175 171 IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity ' … … 177 173 ! 178 174 wn(:,:,jpk) = 0._wp ! bottom boundary condition: w=0 (set once for all) 179 !180 175 ENDIF 181 176 ! !------------------------------! … … 216 211 217 212 #if defined key_bdy 218 IF (lk_bdy) THEN213 IF( lk_bdy ) THEN 219 214 DO jk = 1, jpkm1 220 215 wn(:,:,jk) = wn(:,:,jk) * bdytmask(:,:) … … 224 219 ! 225 220 IF( nn_timing == 1 ) CALL timing_stop('wzv') 226 227 221 ! 228 222 END SUBROUTINE wzv 223 229 224 230 225 SUBROUTINE ssh_swp( kt ) … … 265 260 sshb(:,:) = sshn(:,:) ! before <-- now 266 261 sshn(:,:) = ssha(:,:) ! now <-- after (before already = now) 262 ! 267 263 ELSE !** Leap-Frog time-stepping: Asselin filter + swap 268 264 sshb(:,:) = sshn(:,:) + atfp * ( sshb(:,:) - 2 * sshn(:,:) + ssha(:,:) ) ! before <-- now filtered
Note: See TracChangeset
for help on using the changeset viewer.