diff --git a/SourceCode/PotentWavForce.f90 b/SourceCode/PotentWavForce.f90
index ff94e52..78e54ba 100644
--- a/SourceCode/PotentWavForce.f90
+++ b/SourceCode/PotentWavForce.f90
@@ -113,7 +113,7 @@ SUBROUTINE EFORCE(WK,W1,TP,BETA,AMP,EXFC)
        MOD=ABS(EXFC(MD))/(RHO*G*AMP) !SQRT(REL**2+MDMG**2)
        NREL=-IMG
        NIMG=-REL
-       PHS(MD)=ATAN2D(NIMG,NREL)
+       PHS(MD)=ATAN2(NIMG,NREL)*180.0D0/PI
        WRITE(20+MD,1010) WK,W1,REAL(EXFC(MD)),IMAG(EXFC(MD))
        
        IF (ABS(TP+1.D0).GT.1.E-6.AND.ABS(TP).GT.1.E-6) THEN
@@ -221,4 +221,4 @@ SUBROUTINE RFORCE(WK,W1,TP,AMAS,BDMP)
       END SUBROUTINE RFORCE
 !-------------------------------------------------------------------------------
 END MODULE PotentWavForce
-!*******************************************************************************
\ No newline at end of file
+!*******************************************************************************
diff --git a/SourceCode/PressureElevation.f90 b/SourceCode/PressureElevation.f90
index 152e43d..98e62fe 100644
--- a/SourceCode/PressureElevation.f90
+++ b/SourceCode/PressureElevation.f90
@@ -374,7 +374,7 @@ SUBROUTINE OutputPressureElevation_Diffraction(NFILE)
        REL=REAL(NVCP)
        IMG=IMAG(NVCP)
        MOD=ABS(NVCP)
-       PHS=ATAN2D(IMG,REL)
+       PHS=ATAN2(IMG,REL)*180.0D0/PI
        WRITE(NFILE,1020) OUFR,BETA*180.0D0/PI,IPT,MOD,PHS,REL,IMG
        
       ENDDO
diff --git a/SourceCode/PrintOutput.f90 b/SourceCode/PrintOutput.f90
index e7bafd3..3c6c11f 100644
--- a/SourceCode/PrintOutput.f90
+++ b/SourceCode/PrintOutput.f90
@@ -177,7 +177,7 @@ SUBROUTINE PrintBody_CmplxVal(NFILE,W1,NBETA,RAOType,CVAB)
         MDL(II)=CDABS(CVAB(II))
         NREL(II)=-IMG(II)
         NIMG(II)=-REL(II)
-        PHS(II)=ATAN2D(NIMG(II),NREL(II))
+        PHS(II)=ATAN2(NIMG(II),NREL(II))*180.0D0/PI
         IF (PHS(II).LT.0.D0) PHS(II)=PHS(II)+360.D0
       ENDDO
        
@@ -208,7 +208,7 @@ SUBROUTINE PrintBody_Exfc(NFILE,W1,NBETA,RAOType,CVAB)
         REL(II)=DREAL(CVAB(II))
         IMG(II)=DIMAG(CVAB(II))
         MDL(II)=CDABS(CVAB(II))
-        PHS(II)=ATAN2D(REL(II),IMG(II))
+        PHS(II)=ATAN2(REL(II),IMG(II))*180.0D0/PI
       ENDDO
        
       WRITE(FMT,*) '(F8.4,',NBETA,'(ES14.6),',NBETA,'(F12.4))'
diff --git a/SourceCode/SolveMotion.f90 b/SourceCode/SolveMotion.f90
index 2ae692d..d0fb233 100644
--- a/SourceCode/SolveMotion.f90
+++ b/SourceCode/SolveMotion.f90
@@ -101,7 +101,7 @@ SUBROUTINE SolveMotion(W1,TP,OUFR,BETA,AMP,AMAS,BDMP,&
        MOD=SQRT(REL**2+IMG**2) !ABS(EXFC(IMD))/NFAC
        NREL= REL
        NIMG=-IMG
-       PHS=ATAN2D(NIMG,NREL)
+       PHS=ATAN2(NIMG,NREL)*180.0D0/PI
        
        IF (ABS(TP+1.D0).GT.1.E-6.AND.ABS(TP).GT.1.E-6) THEN
         WRITE(63,1030)  OUFR,BETA*180.0D0/PI,MD,MOD,PHS,NREL,NIMG