diff --git a/model/src/w3profsmd_pdlib.F90 b/model/src/w3profsmd_pdlib.F90 index 3d4636ef6..e9866a363 100644 --- a/model/src/w3profsmd_pdlib.F90 +++ b/model/src/w3profsmd_pdlib.F90 @@ -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 @@ -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)