Added units_mod

* Basic units module
* Implements some operator overloads
* Needs extension to more units

CMakeLists.txt

@@ -102,6 +102,7 @@ set(modules_src
 	"src/spline.f90"
 	"src/stats.f90"
 	"src/time.f90"
+	"src/units.f90"
 	"src/autoDiff/autoDiff.f90"
 	"src/autoDiff/autoDiffArray.f90"
 	"src/autoDiff/autoDiffExponential.f90"
@@ -170,7 +171,7 @@ function(addTest name)
 endfunction(addTest)
 
 addTest("testArray")
-# addTest("testAutoDiff")
+addTest("testAutoDiff")
 addTest("testConfig")
 addTest("testIterate")
 addTest("testKinds")
@@ -180,6 +181,7 @@ addTest("testStats")
 addTest("testTensor")
 addTest("testText")
 addTest("testTime")
+addTest("testUnits")
 addTest("testVtkIO")
 
 add_dependencies(testConfig symlinks)

src/test/testUnits.f90

@@ -0,0 +1,28 @@
+program testUnits_prg
+	!! Test program for units_mod
+	!! @todo
+	!! Needs serious improvements
+	use kinds_mod, only: wp
+	use units_mod
+	implicit none
+
+	call testUnits
+
+contains
+
+	subroutine testUnits
+		type(quantity_t)::dP
+		type(quantity_t)::rho
+		type(quantity_t)::V
+
+		rho = quantity_t(1.225_wp,'kg')/quantity_t(1.0_wp,'m')**3
+		dP  = quantity_t(7.25E-2_wp,'psi')
+		dP = dP%convert( quantity_t(1.0_wp,'Pa') )
+
+		V = sqrt(2.0_wp*dP/rho)
+
+		write(*,*) V%value,V%getChar()
+	end subroutine testUnits
+
+end program testUnits_prg 
+