Merge commit 'aa9c0a0e32928a2e40c7476e99c53d008f10d3c4'

sce/Makefile

@@ -28,6 +28,7 @@ RANLIB = ranlib
 OBJ = o
 FFLAGS = $(OPT)  $(DEBUG)  -o 
 FFLAGS2 =  $(DEBUG) -O3 -o 
+VAR_TYPE = 1 # 0 single, 1 double
 
 
 

sce/osnf/Makefile

@@ -22,7 +22,7 @@ RANLIB = ranlib
 OBJ = o
 FFLAGS = $(OPT)  $(DEBUG)  -o 
 FFLAGS2 =  $(DEBUG) -O3 -o 
-
+VAR_TYPE = 1 # 0 single, 1 double
 
 
 main.exe	:  osnf_lib.a  main.$(OBJ) 
@@ -41,7 +41,7 @@ osnf_lib.a	:  numerics.$(OBJ) zeroin.$(OBJ) sfmin.$(OBJ) \
             svode.$(OBJ) slinpk.$(OBJ) vode.$(OBJ) dlinpk.$(OBJ) vode_integrate.$(OBJ) \
             erfinv.$(OBJ) tridiagonal.$(OBJ) hygfx.$(OBJ) random.$(OBJ)
 numerics_type.$(OBJ)	: numerics_type.f90 
-	$(FOR) numerics_type.f90 $(FFLAGS)numerics_type.$(OBJ)
+	$(FOR) numerics_type.f90 -cpp -DVAR_TYPE=$(VAR_TYPE) $(FFLAGS)numerics_type.$(OBJ)
 numerics.$(OBJ)	: numerics.f90 numerics_type.$(OBJ)
 	$(FOR) numerics.f90 $(FFLAGS)numerics.$(OBJ)
 zeroin.$(OBJ)	: zeroin.f

sce/osnf/numerics_type.f90

@@ -6,9 +6,15 @@ module numerics_type
     integer, parameter :: i1b = selected_int_kind(2)
     integer, parameter :: sp = kind(1.0)
     integer, parameter :: dp = kind(1.0d0)
+#if VAR_TYPE == 1
     integer, parameter :: wp = kind(1.0d0)
+	integer, parameter :: wpc = kind((1.0d0,1.0d0))
+#endif
+#if VAR_TYPE == 0
+    integer, parameter :: wp = kind(1.0)
+	integer, parameter :: wpc = kind((1.0,1.0))
+#endif
 	integer, parameter :: spc = kind((1.0,1.0))
 	integer, parameter :: dpc = kind((1.0d0,1.0d0))
-	integer, parameter :: wpc = kind((1.0d0,1.0d0))
     real(wp), parameter :: pi=3.141592653589793238462643383279502884197_wp
 end module numerics_type