debugging runs with source term

erf-model · Jul 29, 2024 · 2ff39cf · 2ff39cf
1 parent 504f9a2
commit 2ff39cf
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 5 deletions.
diff --git a/model/src/w3iogomd.F90 b/model/src/w3iogomd.F90
@@ -2349,14 +2349,14 @@ SUBROUTINE W3OUTG ( A, FLPART, FLOUTG, FLOUTG2 )
     ! Receive theta array
     !call MPI_RECV(theta_values, n_elements, MPI_DOUBLE, 0, 4, MPI_COMM_WORLD, IERR_MPI)
     ! Print received values to a txt file
-    open(unit=6123, file='mpi_recv.txt', status='unknown', access='append', action="write")
+    open(unit=6123, file='ww3_mpi_recv.txt', status='unknown', access='append', action="write")
      DO JSEA=1, NSEAL
          CALL INIT_GET_ISEA(ISEA, JSEA)
          IX     = MAPSF(ISEA,1)
          IY     = MAPSF(ISEA,2)
          ! Need correct mapping of magnitude_values and theta_values
          COUNTER = IX + (IY-1) * NX 
-         WRITE(6123, *) "(", IX, IY, ")", ISEA, COUNTER, size(magnitude_values), magnitude_values(COUNTER), size(theta_values), theta_values(COUNTER), IERR_MPI
+         WRITE(6123, *) "(", IX, IY, ")", ISEA, JSEA, COUNTER, size(magnitude_values), magnitude_values(ISEA), size(theta_values), theta_values(ISEA), IERR_MPI
      END DO
     ! write(6123,*) 'Magnitude Values:', magnitude_values, 'Theta Values:', theta_values
     close(6123)

diff --git a/model/src/w3srcemd.F90 b/model/src/w3srcemd.F90
@@ -1099,7 +1099,6 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD,          &
 
 ! MY EDITS HERE
     ! open(unit=8123, file='sourceterm.txt', status='unknown', access='append', action="write")
-         print*, " This is the first ifdef for W3_ST6 in w3srce.md "
     ! WRITE(8123, *) "This is the first ifdef for W3_ST6 in w3srce.md"
     ! close(8123)
 
@@ -1854,7 +1853,6 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD,          &
 
     ! MY EDITS HERE
     COUNTER = 1
-    print*, " Source terms HERE "    
     IF (COUNTER.EQ.1) THEN
         OPEN(4121, file='output_SRC.txt', status='replace', action="write")
         WRITE(4121, *) U10DIR, "(x_vel, y_vel) = ", "(", COS(U10DIR), SIN(U10DIR), ") ", "Wind Vel = ", U10ABS

diff --git a/regtests/ww3_ts3/input/switch_ST6_PR1 b/regtests/ww3_ts3/input/switch_ST6_PR1
@@ -1 +1 @@
-NOGRB SHRD PR1 FLX4 LN1 ST6 NL1 BT0 DB0 TR0 BS0 IC0 IS0 REF0 WNT1 WNX1 RWND CRT1 CRX1 MGW MGP MGG O0 O1 O2 O3 O4 O5 O6 O7 O10 O11
+NOGRB DIST MPI PR1 FLX4 LN1 ST6 NL1 BT0 DB0 TR0 BS0 IC0 IS0 REF0 WNT1 WNX1 RWND CRT1 CRX1 MGW MGP MGG O0 O1 O2 O3 O4 O5 O6 O7 O10 O11
Original file line number	Diff line number	Diff line change
		@@ -1 +1 @@
		NOGRB SHRD PR1 FLX4 LN1 ST6 NL1 BT0 DB0 TR0 BS0 IC0 IS0 REF0 WNT1 WNX1 RWND CRT1 CRX1 MGW MGP MGG O0 O1 O2 O3 O4 O5 O6 O7 O10 O11
		NOGRB DIST MPI PR1 FLX4 LN1 ST6 NL1 BT0 DB0 TR0 BS0 IC0 IS0 REF0 WNT1 WNX1 RWND CRT1 CRX1 MGW MGP MGG O0 O1 O2 O3 O4 O5 O6 O7 O10 O11