Skip to content

Commit

Permalink
ww3_diffraction: minor edits ...
Browse files Browse the repository at this point in the history
  • Loading branch information
aronroland committed Feb 9, 2024
1 parent 305153c commit b837c0b
Showing 1 changed file with 3 additions and 4 deletions.
7 changes: 3 additions & 4 deletions model/src/w3profsmd_pdlib.F90
Original file line number Diff line number Diff line change
Expand Up @@ -7990,15 +7990,14 @@ SUBROUTINE PROP_REFRACTION_PR3(IP, ISEA, DTG, CAD, DoLimiter)
END IF
END DO
DO ISP=1,NSPEC
CAD(ISP)=DBLE(VCFLT(ISP))
CAD(ISP)=VCFLT(ISP)
END DO

IF (B_JGS_LDIFR) THEN
DO IK = 1, NK
DO ITH = 1, NTH
ISP = ITH + (IK-1)*NTH
CAD(ISP) = DIFRM(IP)*CAD(ISP)-CG(IK,IP)*(DIFRX(IP)*ESIN(ITH)-DIFRY(IP)*ECOS(ITH))
!CAD(ISP) = CAD(ISP)-CG(IK,IP)*(ESIN(ITH)-ECOS(ITH))
CAD(ISP) = DIFRM(IP)*CAD(ISP)-CG(IK,IP)*(DIFRX(IP)*ESIN(ISP)-DIFRY(IP)*ECOS(ISP))
END DO
END DO
END IF
Expand Down Expand Up @@ -8473,7 +8472,7 @@ SUBROUTINE DIFFRA_SIMPLE
ENDDO
END IF

IF (.FALSE.) THEN
IF (.false.) THEN
OPEN(555, FILE = 'ergdiffr.bin' , FORM = 'UNFORMATTED')
WRITE(555) 1.
WRITE(555) (SNGL(DIFRX(IP)), SNGL(DIFRY(IP)),SNGL(DIFRM(IP)), IP = 1, NP)
Expand Down

0 comments on commit b837c0b

Please sign in to comment.