Description: disable_implicit_type_coversion
Author: Youhei SASAKI <uwabami@gfd-dennou.org>

---
Origin: maintainer
Forwarded: not-needed
Reviewed-By: none
Last-Update: 2011-09-09

--- ispack-0.96.orig/stpack/src/stnm2l.f
+++ ispack-0.96/stpack/src/stnm2l.f
@@ -41,7 +41,7 @@
         N=L-1
         M=0
       ELSE
-        M=(2*MM+3-SQRT(1D0*((2*MM+3)*(2*MM+3)-4*(L+MM))))/2
+        M=INT((2*MM+3-SQRT(1D0*((2*MM+3)*(2*MM+3)-4*(L+MM))))*0.5D0)
         N=M+L-(M*(2*MM-M+3)-MM)
         IF(N.GT.MM) THEN
           N=N-(MM-M+1)
--- ispack-0.96.orig/sppack/src/spnini.f
+++ ispack-0.96/sppack/src/spnini.f
@@ -31,7 +31,7 @@
       RN(1,2)=1
 
       DO L=2,LM
-        N=SQRT(1D0*(L-1))
+        N=INT(SQRT(1D0*(L-1)))
         RN(L,1)=-N*(N+1)
         RN(L,2)=1D0/RN(L,1)
       END DO
--- ispack-0.96.orig/sppack/src/spswhi.f
+++ ispack-0.96/sppack/src/spswhi.f
@@ -32,7 +32,7 @@
       CL(1,4)=1
       CL(1,5)=0
       DO L=2,(MM+1)*(MM+1)
-        N=SQRT(1D0*(L-1))
+        N=INT(SQRT(1D0*(L-1)))
         DNUD=DNU*(N*(N+1)-2D0)**LEV
         C=-DNUD*((2-ALPHA)*(-N*(N+1))+2)/2
         S=-BARPHI*(-N*(N+1))
--- ispack-0.96.orig/sppack/src/spswli.f
+++ ispack-0.96/sppack/src/spswli.f
@@ -32,7 +32,7 @@
       CL(1,4)=1
       CL(1,5)=0
       DO L=2,(MM+1)*(MM+1)
-        N=SQRT(1D0*(L-1))
+        N=INT(SQRT(1D0*(L-1)))
         C=-DNU*((2-ALPHA)*(-N*(N+1))+2)/2
         S=-BARPHI*(-N*(N+1))
         FREQ2=S-C*C
--- ispack-0.96.orig/snpack/src/snnm2l.f
+++ ispack-0.96/snpack/src/snnm2l.f
@@ -31,7 +31,7 @@
 ************************************************************************
       SUBROUTINE SNL2NM(L,N,M)
 
-      N=SQRT(1D0*(L-1))
+      N=INT(SQRT(1D0*(L-1)))
       M=L-N*(N+1)-1
 
       END
--- ispack-0.96.orig/p3pack/src/p3espt.f
+++ ispack-0.96/p3pack/src/p3espt.f
@@ -31,7 +31,7 @@
       DO L=-LM,-1
         DO M=-MM,MM
           DO N=-NM,NM
-            K=SQRT(1D0*(L*L+M*M+N*N))+0.5D0
+            K=INT(SQRT(1D0*(L*L+M*M+N*N))+0.5D0)
             IF(K.LE.KMAX) THEN
               ES(K)=ES(K)+((L*L+M*M)*Z(N,M,L,1)*Z(N,M,L,1)
      &          +(L*L+N*N)*Z(N,M,L,2)*Z(N,M,L,2)
@@ -43,7 +43,7 @@
       DO L=1,LM
         DO M=-MM,MM
           DO N=-NM,NM
-            K=SQRT(1D0*(L*L+M*M+N*N))+0.5D0
+            K=INT(SQRT(1D0*(L*L+M*M+N*N))+0.5D0)
             IF(K.LE.KMAX) THEN
               ES(K)=ES(K)+((L*L+M*M)*Z(N,M,L,1)*Z(N,M,L,1)
      &          +(L*L+N*N)*Z(N,M,L,2)*Z(N,M,L,2)
@@ -56,7 +56,7 @@
       L=0
       DO M=-MM,-1
         DO N=-NM,NM
-          K=SQRT(1D0*(M*M+N*N))+0.5D0
+          K=INT(SQRT(1D0*(M*M+N*N))+0.5D0)
           IF(K.LE.KMAX) THEN
             ES(K)=ES(K)+(M*M*Z(N,M,L,2)*Z(N,M,L,2)
      &        +(M*M+N*N)*Z(N,M,L,1)*Z(N,M,L,1))
@@ -66,7 +66,7 @@
       END DO
       DO M=1,MM
         DO N=-NM,NM
-          K=SQRT(1D0*(M*M+N*N))+0.5D0
+          K=INT(SQRT(1D0*(M*M+N*N))+0.5D0)
           IF(K.LE.KMAX) THEN
             ES(K)=ES(K)+(M*M*Z(N,M,L,2)*Z(N,M,L,2)
      &        +(M*M+N*N)*Z(N,M,L,1)*Z(N,M,L,1))
--- ispack-0.96.orig/sppack/src/spmini.f
+++ ispack-0.96/sppack/src/spmini.f
@@ -28,7 +28,7 @@
       LM=(MM+1)*(MM+1)
 
       DO L=1,LM
-        N=SQRT(1D0*(L-1))
+        N=INT(SQRT(1D0*(L-1)))
         M=L-N*(N+1)-1
         IRM(L,1)=-M+1+N*(N+1)
         IRM(L,2)=M
--- ispack-0.96.orig/sjpack/src/sjl2nm.f
+++ ispack-0.96/sjpack/src/sjl2nm.f
@@ -27,7 +27,7 @@
         M=0
         N=L-1
       ELSE
-        M=NN+1.5D0-SQRT(1D0*(NN+1)*(NN+1)-L+1)
+        M=INT(NN+1.5D0-SQRT(1D0*(NN+1)*(NN+1)-L+1))
         LD=1+NN+1+(M-1)*(2*NN+2-M)
         N=M+(L-LD)/2
         IF(L-LD.NE.2*(N-M)) THEN
--- ispack-0.96.orig/fepack/src/fepack.F
+++ ispack-0.96/fepack/src/fepack.F
@@ -184,13 +184,13 @@
       REAL*8 DD(ND)
 #ifdef LEN
       DO I=1,ND
-        RD=DD(I)
+        RD=REAL(DD(I))
         CALL FLCEN4(RD,IW)
         CALL FHUPUT(IU,4,IW)
       END DO
 #elif IEEE
       DO 10 I=1,ND
-        RD=DD(I)
+        RD=REAL(DD(I))
         CALL FHUPUT(IU,4,RD)
    10 CONTINUE
 #elif IBM
@@ -210,12 +210,12 @@
       DO I=1,ND
         CALL FHUGET(IU,4,IW)
         CALL FLCEN4(IW,RD)
-        DD(I)=RD
+        DD(I)=DBLE(RD)
       END DO
 #elif IEEE
       DO 10 I=1,ND
         CALL FHUGET(IU,4,RD)
-        DD(I)=RD
+        DD(I)=DBLE(RD)
    10 CONTINUE
 #elif IBM
       DO 10 I=1,ND
@@ -234,7 +234,7 @@
       INTEGER IW(2)
 
       DO I=1,ND
-        DD=RD(I)
+        DD=DBLE(RD(I))
         CALL FLCEN8(DD,IW)
         CALL FHUPUT(IU,8,IW)
       END DO
@@ -243,7 +243,7 @@
       REAL*8 DD
 
       DO 10 I=1,ND
-        DD=RD(I)
+        DD=DBLE(RD(I))
         CALL FHUPUT(IU,8,DD)
    10 CONTINUE
 #elif IBM
@@ -268,7 +268,7 @@
       DO 10 I=1,ND
         CALL FHUGET(IU,8,IW)
         CALL FLCEN8(IW,DD)
-        RD(I)=DD
+        RD(I)=REAL(DD)
    10 CONTINUE
 #elif IEEE
       REAL RD(ND)
@@ -276,7 +276,7 @@
 
       DO 10 I=1,ND
         CALL FHUGET(IU,8,DD)
-        RD(I)=DD
+        RD(I)=REAL(DD)
    10 CONTINUE
 #elif IBM
       REAL RD(ND)
--- ispack-0.96.orig/p3pack-mpi/src/p3empt.f
+++ ispack-0.96/p3pack-mpi/src/p3empt.f
@@ -49,7 +49,7 @@
         DO L=MAX(LS,1),LE
           DO M=-MM,MM
             DO N=-NM,NM
-              K=SQRT(1D0*(L*L+M*M+N*N))+0.5D0
+              K=INT(SQRT(1D0*(L*L+M*M+N*N))+0.5D0)
               IF(K.LE.KMAX) THEN
                 W(K)=W(K)+((L*L+M*M)*Z(N,M,1,L-LS)*Z(N,M,1,L-LS)
      &            +(L*L+N*N)*Z(N,M,2,L-LS)*Z(N,M,2,L-LS)
@@ -68,7 +68,7 @@
           L=0
           DO M=-MM,-1
             DO N=-NM,NM
-              K=SQRT(1D0*(M*M+N*N))+0.5D0
+              K=INT(SQRT(1D0*(M*M+N*N))+0.5D0)
               IF(K.LE.KMAX) THEN
                 W(K)=W(K)+(M*M*Z(N,M,2,0)*Z(N,M,2,0)
      &            +(M*M+N*N)*Z(N,M,1,0)*Z(N,M,1,0))
@@ -78,7 +78,7 @@
           END DO
           DO M=1,MM
             DO N=-NM,NM
-              K=SQRT(1D0*(M*M+N*N))+0.5D0
+              K=INT(SQRT(1D0*(M*M+N*N))+0.5D0)
               IF(K.LE.KMAX) THEN
                 W(K)=W(K)+(M*M*Z(N,M,2,0)*Z(N,M,2,0)
      &            +(M*M+N*N)*Z(N,M,1,0)*Z(N,M,1,0))
