Skip to content

Commit

Permalink
ww3_nonlinear_cg: more debug
Browse files Browse the repository at this point in the history
  • Loading branch information
aronroland committed Oct 5, 2024
1 parent cacffb8 commit 2413701
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions model/src/w3profsmd_pdlib.F90
Original file line number Diff line number Diff line change
Expand Up @@ -5696,7 +5696,7 @@ SUBROUTINE PDLIB_JACOBI_GAUSS_SEIDEL_BLOCK(IMOD, FACX, FACY, DTG, VGX, VGY, LCAL
B_JAC = ZERO
ENDIF

IF (IAPROC == 1) WRITE(*,*) 'SUM B_JAC and ASPAR_DIAG_ALL', SUM(B_JAC), SUM(ASPAR_DIAG_ALL)
!IF (IAPROC == 1) WRITE(*,*) 'SUM B_JAC and ASPAR_DIAG_ALL', SUM(B_JAC), SUM(ASPAR_DIAG_ALL)
call print_memcheck(memunit, 'memcheck_____:'//' WW3_PROP SECTION 3')
!
! source terms
Expand Down Expand Up @@ -5838,7 +5838,7 @@ SUBROUTINE PDLIB_JACOBI_GAUSS_SEIDEL_BLOCK(IMOD, FACX, FACY, DTG, VGX, VGY, LCAL
#endif
END IF
END DO
IF (IAPROC == 1) WRITE(*,*) 'IMEM SUM ESUM', IP, SUM(ESUM)
!IF (IAPROC == 1) WRITE(*,*) 'IMEM SUM ESUM', IP, SUM(ESUM)
ENDIF ! IMEM

#ifdef W3_DEBUGSOLVERCOH
Expand Down Expand Up @@ -5926,7 +5926,7 @@ SUBROUTINE PDLIB_JACOBI_GAUSS_SEIDEL_BLOCK(IMOD, FACX, FACY, DTG, VGX, VGY, LCAL
PRE_VA(:, IP)=REAL(eSum)
#endif
eSum(1:NSPEC) = eSum(1:NSPEC) / ASPAR_DIAG(1:NSPEC)
IF (IAPROC == 1) WRITE(*,*) 'SUM ESUM AFTER DIV', IP, SUM(ESUM)
!IF (IAPROC == 1) WRITE(*,*) 'SUM ESUM AFTER DIV', IP, SUM(ESUM)
#ifdef W3_DEBUGFREQSHIFT
WRITE(740+IAPROC,*) 'JSEA=', JSEA, ' nbIter=', nbIter
DO ISP=1,NSPEC
Expand Down Expand Up @@ -6299,8 +6299,8 @@ SUBROUTINE PDLIB_JACOBI_GAUSS_SEIDEL_BLOCK(IMOD, FACX, FACY, DTG, VGX, VGY, LCAL
ENDIF ! FLSOU
END DO ! JSEA

IF (IAPROC == 1) WRITE(*,*) 'SUM VA', SUM(VA)
STOP
!IF (IAPROC == 1) WRITE(*,*) 'SUM VA', SUM(VA)
!STOP

#ifdef WEIGHTS
INQUIRE ( FILE='weights.ww3', EXIST = lexist )
Expand Down

0 comments on commit 2413701

Please sign in to comment.