diff --git a/test/Ave_SS.py b/test/Ave_SS.py
index 3ea8d1c7..b43e2edd 100644
--- a/test/Ave_SS.py
+++ b/test/Ave_SS.py
@@ -14,30 +14,30 @@ def main():
     #[e] read modpara
 
     #
-    Ref_ave_Temp     = np.zeros([max_eigen],dtype=float)
-    Ref_err_Temp     = np.zeros([max_eigen],dtype=float)
-    Ref_ave_Ene      = np.zeros([max_eigen],dtype=float)
-    Ref_err_Ene      = np.zeros([max_eigen],dtype=float)
+    Ref_ave_Temp     = np.zeros([max_eigen],dtype=np.float64)
+    Ref_err_Temp     = np.zeros([max_eigen],dtype=np.float64)
+    Ref_ave_Ene      = np.zeros([max_eigen],dtype=np.float64)
+    Ref_err_Ene      = np.zeros([max_eigen],dtype=np.float64)
     # 
-    ave_Temp     = np.zeros([max_eigen],dtype=float)
-    err_Temp     = np.zeros([max_eigen],dtype=float)
+    ave_Temp     = np.zeros([max_eigen],dtype=np.float64)
+    err_Temp     = np.zeros([max_eigen],dtype=np.float64)
     # 
-    InvTemp      = np.zeros([max_set,max_eigen],dtype=float)
-    ave_InvTemp  = np.zeros([max_eigen],dtype=float)
-    err_InvTemp  = np.zeros([max_eigen],dtype=float)
-    Ene          = np.zeros([max_set,max_eigen],dtype=float)
-    ave_Ene      = np.zeros([max_eigen],dtype=float)
-    err_Ene      = np.zeros([max_eigen],dtype=float)
-    Ene2         = np.zeros([max_set,max_eigen],dtype=float)
-    ave_Ene2     = np.zeros([max_eigen],dtype=float)
-    err_Ene2     = np.zeros([max_eigen],dtype=float)
+    InvTemp      = np.zeros([max_set,max_eigen],dtype=np.float64)
+    ave_InvTemp  = np.zeros([max_eigen],dtype=np.float64)
+    err_InvTemp  = np.zeros([max_eigen],dtype=np.float64)
+    Ene          = np.zeros([max_set,max_eigen],dtype=np.float64)
+    ave_Ene      = np.zeros([max_eigen],dtype=np.float64)
+    err_Ene      = np.zeros([max_eigen],dtype=np.float64)
+    Ene2         = np.zeros([max_set,max_eigen],dtype=np.float64)
+    ave_Ene2     = np.zeros([max_eigen],dtype=np.float64)
+    err_Ene2     = np.zeros([max_eigen],dtype=np.float64)
     #
-    Spc          = np.zeros([max_set,max_eigen],dtype=float)
-    ave_Spc      = np.zeros([max_eigen],dtype=float)
-    err_Spc      = np.zeros([max_eigen],dtype=float)
-    Ent          = np.zeros([max_set,max_eigen-1],dtype=float)
-    ave_Ent      = np.zeros([max_eigen-1],dtype=float)
-    err_Ent      = np.zeros([max_eigen-1],dtype=float)
+    Spc          = np.zeros([max_set,max_eigen],dtype=np.float64)
+    ave_Spc      = np.zeros([max_eigen],dtype=np.float64)
+    err_Spc      = np.zeros([max_eigen],dtype=np.float64)
+    Ent          = np.zeros([max_set,max_eigen-1],dtype=np.float64)
+    ave_Ent      = np.zeros([max_eigen-1],dtype=np.float64)
+    err_Ent      = np.zeros([max_eigen-1],dtype=np.float64)
     #
 
     for cnt_set in range(0,max_set):