Changeset 4460

Timestamp:

2014-02-04T13:17:29+01:00 (10 years ago)

Author:

trackstand2

Message:

Addition of dbg SUM outputs and alteration of loop limits in dynldfiso and dynldf_bilapg

Location:

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

Files:

• dynldf_bilapg.F90 (modified) (7 diffs)
• dynldf_iso.F90 (modified) (1 diff)

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

@@ -121 +121 @@
       ENDIF
 
+#if defined ARPDBGSUM
+      WRITE (*,*) narea,': ARPDBG: start dyn_ldf_bilapg SUM(ua)= ',SUM(ua),' at step=',kt
+#endif
+
       ! Laplacian of (ub,vb) multiplied by ahm
       ! --------------------------------------  
       CALL ldfguv( ub, vb, zwk1, zwk2, 1 )      ! rotated harmonic operator applied to (ub,vb)
+#if defined ARPDBGSUM
+      WRITE (*,*) narea,': ARPDBG: dyn_ldf_bilapg SUM(zwk1,zwk2)= ',SUM(zwk1),SUM(zwk2),' at step=',kt
+#endif
       !                                         ! and multiply by ahmu, ahmv (output in (zwk1,zwk2) )
       CALL lbc_lnk( zwk1, 'U', -1. )   ;   CALL lbc_lnk( zwk2, 'V', -1. )     ! Lateral boundary conditions
@@ -151 +158 @@
       IF( wrk_not_released(3, 3,4) )   CALL ctl_stop('dyn_ldf_bilapg: failed to release workspace arrays')
       !
+#if defined ARPDBGSUM
+      WRITE (*,*) narea,': ARPDBG: end dyn_ldf_bilapg SUM(ua,va)= ',SUM(ua),SUM(va),' at step=',kt
+#endif
+
    END SUBROUTINE dyn_ldf_bilapg
 
@@ -235 +246 @@
          CALL ctl_stop('dyn:ldfguv: requested workspace arrays unavailable')   ;   RETURN
       END IF
-      CALL timing_start('ldfguv_1st')
+      !CALL timing_start('ldfguv_1st')
 
 #if defined key_z_first
@@ -605 +616 @@
       !                                                ! ===============
 #endif
-      CALL timing_stop('ldfguv_1st','section')
+      !CALL timing_stop('ldfguv_1st','section')
 
       !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,synchro,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
 
-      CALL timing_start('ldfguv_2nd')
+      !CALL timing_start('ldfguv_2nd')
       !                             ! ************ !   ! ===============
       DO jj = 2, jpjm1              !  Second step !   ! Horizontal slab
@@ -679 +690 @@
 #if defined key_z_first
          DO ji = 2, jpim1
-            DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1
+            DO jk = 2, mbkmax(ji,jj) ! -1 ! jpkm1
 #else
          DO jk = 2, jpkm1
@@ -712 +723 @@
 #if defined key_z_first
          DO ji = 2, jpim1
-            DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1
+            DO jk = 2, mbkmax(ji,jj) ! -1 ! jpkm1
 #else
          DO jk = 2, jpkm1
@@ -793 +804 @@
       END DO                                           !   End of slab
       !                                                ! ===============
-      CALL timing_stop('ldfguv_2nd','section')
+      !CALL timing_stop('ldfguv_2nd','section')
 
       IF( wrk_not_released(2, 1,2,3,4,5,6,7,8) )   CALL ctl_stop('dyn:ldfguv: failed to release workspace arrays')

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

@@ -1143 +1143 @@
 #endif
       !CALL timing_stop('dyn_ldf_iso_vslab','section')
+#if defined ARPDBGSUM
+      WRITE (*,*) narea,': ARPDBG: end dyn_ldf_iso SUM(ua)= ',SUM(ua),' at step=',kt
+#endif
 
       IF( wrk_not_released(2, 1,2,3,4,5,6,7,8) )   CALL ctl_stop('dyn_ldf_iso: failed to release workspace arrays')