- Timestamp:
- 2014-02-04T13:14:59+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r4427 r4438 225 225 IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt 226 226 227 #if defined ARPDBGSUM 227 228 WRITE(*,*)'ARPDBG, ssh_wzv: sum WWW of hdivn=',SUM(hdivn),' at step=',kt 228 229 WRITE(*,*)'ARPDBG, ssh_wzv: sum WWW of fse3t=',SUM(fse3t(:,:,:)),' at step=',kt 229 230 #endif 230 231 ! !------------------------------! 231 232 ! ! After Sea Surface Height ! … … 245 246 END DO 246 247 #endif 248 249 #if defined ARPDBGSUM 247 250 WRITE(*,*)'ARPDBG, ssh_wzv: sum XXX of zhdiv=',SUM(zhdiv),' at step=',kt 248 251 WRITE(*,*)'ARPDBG, ssh_wzv: sum XXX of ssha=',SUM(ssha),' at step=',kt 252 #endif 249 253 ! ! Sea surface elevation time stepping 250 254 ! In forward Euler time stepping case, the same formulation as in the leap-frog case can be used … … 256 260 ssha(:,:) = ( sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) ) ) * tmask(:,:,1) 257 261 #endif 262 263 #if defined ARPDBGSUM 258 264 WRITE(*,*)'ARPDBG, ssh_wzv: sum YYY of sshb=',SUM(sshb),' at step=',kt 259 265 WRITE(*,*)'ARPDBG, ssh_wzv: sum YYY of ssha=',SUM(ssha),' at step=',kt 266 #endif 267 260 268 #if defined key_agrif 261 269 CALL agrif_ssh( kt ) … … 424 432 ! 425 433 ELSE !** Leap-Frog time-stepping: Asselin filter + swap 434 435 #if defined ARPDBGSUM 426 436 WRITE(*,*) 'ARPDBG: ssh_nxt: SUM of sshb = ',SUM(sshb),' at step=',kt 427 437 WRITE(*,*) 'ARPDBG: ssh_nxt: SUM of sshn = ',SUM(sshn),' at step=',kt 438 #endif 439 428 440 zec = atfp * rdt / rau0 429 441 DO jj = 1, jpj
Note: See TracChangeset
for help on using the changeset viewer.