--- a/fem/src/modules/HelmholtzBEM.src	2013-01-01 14:20:36.418566986 +0800
+++ b/fem/src/modules/HelmholtzBEM.src	2013-01-01 14:19:30.962572469 +0800
@@ -78,7 +78,7 @@
 
      LOGICAL :: AllocationsDone = .FALSE., GotIt
  
+     COMPLEX(KIND=dp), POINTER, CONTIGUOUS :: Potential(:),ForceVector(:), &
-     COMPLEX(KIND=dp), POINTER CONTIG :: Potential(:),ForceVector(:), &
                Diagonal(:)
      INTEGER, POINTER :: PotentialPerm(:), BoundaryPerm(:)
 
@@ -680,7 +680,7 @@
 !------------------------------------------------------------------------------
        INTEGER ::  N
  
+       COMPLEX(KIND=dp), CONTIGUOUS ::  A(:,:),x(:),b(:)
-       COMPLEX(KIND=dp) CONTIG ::  A(:,:),x(:),b(:)
 !------------------------------------------------------------------------------
 
        SELECT CASE( ListGetString( Solver % Values, 'Linear System Solver' ) )
@@ -709,7 +709,7 @@
 !------------------------------------------------------------------------------
        TYPE(Solver_t) :: SolverParam
        INTEGER :: N
+       COMPLEX(KIND=dp), DIMENSION(:), CONTIGUOUS :: x,b
-       COMPLEX(KIND=dp), DIMENSION(:) CONTIG :: x,b
 !------------------------------------------------------------------------------
        REAL(KIND=dp) :: dpar(50)
 
--- a/fem/src/modules/CoupledNSUpdate.src	2012-12-31 23:55:38.000000000 +0800
+++ b/fem/src/modules/CoupledNSUpdate.src	2013-01-01 14:32:44.266506017 +0800
@@ -91,7 +91,7 @@
   INTEGER, ALLOCATABLE :: Indexes(:)
   TYPE(Variable_t), POINTER :: FlowSol
 
-  REAL(KIND=dp), POINTER CONTIG :: Mx(:),Mb(:),Mr(:)
+  REAL(KIND=dp), POINTER, CONTIGUOUS :: Mx(:),Mb(:),Mr(:)
   TYPE(ValueList_t), POINTER :: SolverParams
 
   SAVE STIFF, LOAD, FORCE, rho, mu, Velocity, AllocationsDone, &
--- a/fem/src/modules/PoissonBEM.src	2012-12-31 23:55:38.000000000 +0800
+++ b/fem/src/modules/PoissonBEM.src	2013-01-01 14:33:47.036500759 +0800
@@ -658,7 +658,7 @@
 !------------------------------------------------------------------------------
        TYPE(Solver_t) :: SolverParam
        INTEGER :: N
-       REAL(KIND=dp), DIMENSION(:) CONTIG :: x,b
+       REAL(KIND=dp), DIMENSION(:), CONTIGUOUS :: x,b
 !------------------------------------------------------------------------------
        REAL(KIND=dp) :: dpar(50)