Skip to content

Commit

Permalink
fix indentation in w3wave
Browse files Browse the repository at this point in the history
  • Loading branch information
JessicaMeixner-NOAA committed Dec 30, 2024
1 parent e82df78 commit 9e7c923
Showing 1 changed file with 63 additions and 63 deletions.
126 changes: 63 additions & 63 deletions model/src/w3wavemd.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1492,79 +1492,79 @@ SUBROUTINE W3WAVE ( IMOD, ODAT, TEND, STAMP, NO_OUT &

IF ((IOBP_LOC(JSEA).eq.1..or.IOBP_LOC(JSEA).eq. 3).and.IOBDP_LOC(JSEA).eq.1.and.IOBPA_LOC(JSEA).eq.0) THEN

IX = MAPSF(ISEA,1)
IY = MAPSF(ISEA,2)
DELA=1.
DELX=1.
DELY=1.
IX = MAPSF(ISEA,1)
IY = MAPSF(ISEA,2)
DELA=1.
DELX=1.
DELY=1.

#ifdef W3_REF1
IF (GTYPE.EQ.RLGTYPE) THEN
DELX=SX*CLATS(ISEA)/FACX
DELY=SY/FACX
DELA=DELX*DELY
END IF
IF (GTYPE.EQ.CLGTYPE) THEN
! Maybe what follows works also for RLGTYPE ... to be verified
DELX=HPFAC(IY,IX)/ FACX
DELY=HQFAC(IY,IX)/ FACX
DELA=DELX*DELY
END IF
REFLEC=REFLC(:,ISEA)
REFLEC(4)=BERG(ISEA)*REFLEC(4)
REFLED=REFLD(:,ISEA)
IF (GTYPE.EQ.RLGTYPE) THEN
DELX=SX*CLATS(ISEA)/FACX
DELY=SY/FACX
DELA=DELX*DELY
END IF
IF (GTYPE.EQ.CLGTYPE) THEN
! Maybe what follows works also for RLGTYPE ... to be verified
DELX=HPFAC(IY,IX)/ FACX
DELY=HQFAC(IY,IX)/ FACX
DELA=DELX*DELY
END IF
REFLEC=REFLC(:,ISEA)
REFLEC(4)=BERG(ISEA)*REFLEC(4)
REFLED=REFLD(:,ISEA)
#endif

#ifdef W3_BT4
D50=SED_D50(ISEA)
PSIC=SED_PSIC(ISEA)
D50=SED_D50(ISEA)
PSIC=SED_PSIC(ISEA)
#endif
!
!
#ifdef W3_DEBUGSRC
IF (IX .eq. DEBUG_NODE) THEN
WRITE(740+IAPROC,*) 'NODE_SRCE_IMP_PRE : IX=', IX, ' JSEA=', JSEA
END IF
WRITE(740+IAPROC,*) 'IT/IX/IY/IMOD=', IT, IX, IY, IMOD
WRITE(740+IAPROC,*) 'ISEA/JSEA=', ISEA, JSEA
WRITE(740+IAPROC,*) 'Before sum(VA)=', sum(VA(:,JSEA))
FLUSH(740+IAPROC)
#endif
CALL W3SRCE(srce_imp_pre, IT, ISEA, JSEA, IX, IY, IMOD, &
VAold(:,JSEA), VA(:,JSEA), &
VSioDummy, VDioDummy, SHAVETOT(JSEA), &
ALPHA(1:NK,JSEA), WN(1:NK,ISEA), &
CG(1:NK,ISEA), CLATS(ISEA), DW(ISEA), U10(ISEA), &
U10D(ISEA), &
IF (IX .eq. DEBUG_NODE) THEN
WRITE(740+IAPROC,*) 'NODE_SRCE_IMP_PRE : IX=', IX, ' JSEA=', JSEA
END IF
WRITE(740+IAPROC,*) 'IT/IX/IY/IMOD=', IT, IX, IY, IMOD
WRITE(740+IAPROC,*) 'ISEA/JSEA=', ISEA, JSEA
WRITE(740+IAPROC,*) 'Before sum(VA)=', sum(VA(:,JSEA))
FLUSH(740+IAPROC)
#endif
CALL W3SRCE(srce_imp_pre, IT, ISEA, JSEA, IX, IY, IMOD, &
VAold(:,JSEA), VA(:,JSEA), &
VSioDummy, VDioDummy, SHAVETOT(JSEA), &
ALPHA(1:NK,JSEA), WN(1:NK,ISEA), &
CG(1:NK,ISEA), CLATS(ISEA), DW(ISEA), U10(ISEA), &
U10D(ISEA), &
#ifdef W3_FLX5
TAUA(ISEA), TAUADIR(ISEA), &
#endif
AS(ISEA), UST(ISEA), &
USTDIR(ISEA), CX(ISEA), CY(ISEA), &
ICE(ISEA), ICEH(ISEA), ICEF(ISEA), &
ICEDMAX(ISEA), &
REFLEC, REFLED, DELX, DELY, DELA, &
TRNX(IY,IX), TRNY(IY,IX), BERG(ISEA), &
FPIS(ISEA), DTDYN(JSEA), &
FCUT(JSEA), DTGpre, TAUWX(JSEA), TAUWY(JSEA), &
TAUOX(JSEA), TAUOY(JSEA), TAUWIX(JSEA), &
TAUWIY(JSEA), TAUWNX(JSEA), &
TAUWNY(JSEA), PHIAW(JSEA), CHARN(JSEA), &
TWS(JSEA), PHIOC(JSEA), TMP1, D50, PSIC, TMP2, &
PHIBBL(JSEA), TMP3, TMP4, PHICE(JSEA), &
TAUOCX(JSEA), TAUOCY(JSEA), WNMEAN(JSEA), &
RHOAIR(ISEA), ASF(ISEA))
IF (.not. LSLOC) THEN
VSTOT(:,JSEA) = VSioDummy
VDTOT(:,JSEA) = VDioDummy
ENDIF
TAUA(ISEA), TAUADIR(ISEA), &
#endif
AS(ISEA), UST(ISEA), &
USTDIR(ISEA), CX(ISEA), CY(ISEA), &
ICE(ISEA), ICEH(ISEA), ICEF(ISEA), &
ICEDMAX(ISEA), &
REFLEC, REFLED, DELX, DELY, DELA, &
TRNX(IY,IX), TRNY(IY,IX), BERG(ISEA), &
FPIS(ISEA), DTDYN(JSEA), &
FCUT(JSEA), DTGpre, TAUWX(JSEA), TAUWY(JSEA), &
TAUOX(JSEA), TAUOY(JSEA), TAUWIX(JSEA), &
TAUWIY(JSEA), TAUWNX(JSEA), &
TAUWNY(JSEA), PHIAW(JSEA), CHARN(JSEA), &
TWS(JSEA), PHIOC(JSEA), TMP1, D50, PSIC, TMP2, &
PHIBBL(JSEA), TMP3, TMP4, PHICE(JSEA), &
TAUOCX(JSEA), TAUOCY(JSEA), WNMEAN(JSEA), &
RHOAIR(ISEA), ASF(ISEA))
IF (.not. LSLOC) THEN
VSTOT(:,JSEA) = VSioDummy
VDTOT(:,JSEA) = VDioDummy
ENDIF
#ifdef W3_DEBUGSRC
WRITE(740+IAPROC,*) 'After sum(VA)=', sum(VA(:,JSEA))
WRITE(740+IAPROC,*) ' sum(VSTOT)=', sum(VSTOT(:,JSEA))
WRITE(740+IAPROC,*) ' sum(VDTOT)=', sum(VDTOT(:,JSEA))
WRITE(740+IAPROC,*) ' SHAVETOT=', SHAVETOT(JSEA)
FLUSH(740+IAPROC)
WRITE(740+IAPROC,*) 'After sum(VA)=', sum(VA(:,JSEA))
WRITE(740+IAPROC,*) ' sum(VSTOT)=', sum(VSTOT(:,JSEA))
WRITE(740+IAPROC,*) ' sum(VDTOT)=', sum(VDTOT(:,JSEA))
WRITE(740+IAPROC,*) ' SHAVETOT=', SHAVETOT(JSEA)
FLUSH(740+IAPROC)
#endif
ENDIF
ENDIF
END DO ! JSEA
END IF ! PDLIB
#endif
Expand Down

0 comments on commit 9e7c923

Please sign in to comment.