Changeset 4438 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN

Timestamp:

2014-02-04T13:14:59+01:00 (10 years ago)

Author:

trackstand2

Message:

Protect debug SUM output in dynnxt and sshwzv with ARPDBGSUM

Location:

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN

Files:

• dynnxt.F90 (modified) (3 diffs)
• sshwzv.F90 (modified) (4 diffs)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

@@ -168 +168 @@
       CALL lbc_lnk( va, 'V', -1. ) 
       !
+#if defined ARPDBGSUM
       WRITE (*,*) 'ARPDBG: dyn_nxt: after lbc update sum(ua)=',SUM(ua),' at step=',kt
+#endif
+
 # if defined key_obc
       !                                !* OBC open boundaries
@@ -225 +228 @@
 #endif
 
+#if defined ARPDBGSUM
       WRITE(*,*) 'ARPDBG: dyn_nxt b4 time filter SUM(ua)=',SUM(ua),'at step=',kt
+#endif
 
       ! Time filter and swap of dynamics arrays
@@ -404 +409 @@
       ENDIF
 
+#if defined ARPDBGSUM
       WRITE(*,*) 'ARPDBG: dyn_nxt end SUM(un)=',SUM(un),'at step=',kt
+#endif
 
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=un, clinfo1=' nxt  - Un: ', mask1=umask,   &

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

@@ -225 +225 @@
       IF( neuler == 0 .AND. kt == nit000 )   z2dt = rdt
 
+#if defined ARPDBGSUM
       WRITE(*,*)'ARPDBG, ssh_wzv: sum WWW of hdivn=',SUM(hdivn),' at step=',kt
       WRITE(*,*)'ARPDBG, ssh_wzv: sum WWW of fse3t=',SUM(fse3t(:,:,:)),' at step=',kt
-
+#endif
       !                                           !------------------------------!
       !                                           !   After Sea Surface Height   !
@@ -245 +246 @@
       END DO
 #endif
+
+#if defined ARPDBGSUM
       WRITE(*,*)'ARPDBG, ssh_wzv: sum XXX of zhdiv=',SUM(zhdiv),' at step=',kt
       WRITE(*,*)'ARPDBG, ssh_wzv: sum XXX of ssha=',SUM(ssha),' at step=',kt
+#endif
       !                                                ! Sea surface elevation time stepping
       ! 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)
 #endif
+
+#if defined ARPDBGSUM
       WRITE(*,*)'ARPDBG, ssh_wzv: sum YYY of sshb=',SUM(sshb),' at step=',kt
       WRITE(*,*)'ARPDBG, ssh_wzv: sum YYY of ssha=',SUM(ssha),' at step=',kt
+#endif
+
 #if defined key_agrif
       CALL agrif_ssh( kt )
@@ -424 +432 @@
             !
          ELSE                                         !** Leap-Frog time-stepping: Asselin filter + swap
+
+#if defined ARPDBGSUM
             WRITE(*,*) 'ARPDBG: ssh_nxt: SUM of sshb = ',SUM(sshb),' at step=',kt
             WRITE(*,*) 'ARPDBG: ssh_nxt: SUM of sshn = ',SUM(sshn),' at step=',kt
+#endif
+
             zec = atfp * rdt / rau0
             DO jj = 1, jpj