--- spml-g95.orig/src/aq_module.f90
+++ spml-g95/src/aq_module.f90
@@ -368,8 +368,6 @@
   save im, km, ra, alpha, beta, gamma, CF, CB, g_R, g_R_Weight, md, jmax
   save first_r2inv, first_Tau_D, first_Tau_N, first_Grid_N
 
-!$omp threadprivate(first_r2inv, first_Tau_D, first_Tau_N, first_Grid_N)
-
 contains
 
 ! --- 
@@ -794,8 +792,6 @@
 
     save R2MTX, kp
 
-!$omp threadprivate(R2MTX, kp)
-
     if ( size(aq_data,1) /= jmax ) then
        call MessageNotify('E','aq_r2Inv_aq', &
             '1st dim. of the spectral data should be same as dim. of MD.')
@@ -1012,8 +1008,6 @@
     integer :: k, j, kstr, kend
     save    :: alu, kp
 
-!$omp threadprivate(alu,kp)
-
     if ( size(aq_data,2)-1 < km ) then
        call MessageNotify('E','aq_BoundaryTau_D', &
             'The spectral dimension of input data too small.')
@@ -1117,8 +1111,6 @@
 
     save    :: alu, kp
 
-!$omp threadprivate(alu,kp)
-
     if ( size(aq_data,2)-1 < km ) then
        call MessageNotify('E','aq_BoundaryTau_N', &
             'The spectral dimension of input data too small.')
@@ -1269,8 +1261,6 @@
     integer :: i
     save    :: alu, kp
 
-!$omp threadprivate(alu,kp)
-
     if ( size(ag_data,2) < im ) then
        call MessageNotify('E','aq_BoundaryGrid_N', &
             'The dimension of input data too small.')
--- spml-g95.orig/src/at_module.f90
+++ spml-g95/src/at_module.f90
@@ -885,8 +885,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( size(at_data,2)-1 < km ) then
        call MessageNotify('E','at_BoundariesTau_DD', &
             'The Chebyshev dimension of input data too small.')
@@ -975,8 +973,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( size(at_data,2)-1 < km ) then
        call MessageNotify('E','at_BoundariesTau_DN', &
             'The Chebyshev dimension of input data too small.')
@@ -1065,8 +1061,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( size(at_data,2)-1 < km ) then
        call MessageNotify('E','at_BoundariesTau_ND', &
             'The Chebyshev dimension of input data too small.')
@@ -1156,8 +1150,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( size(at_data,2)-1 < km ) then
        call MessageNotify('E','at_BoundariesTau_NN', &
             'The Chebyshev dimension of input data too small.')
@@ -1316,8 +1308,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( im /= km ) then
        call MessageNotify('E','at_BoundariesGrid_DD', &
             'Chebyshev truncation and number of grid points should be same.')
@@ -1410,8 +1400,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( im /= km ) then
        call MessageNotify('E','at_BoundariesGrid_DN', &
             'Chebyshev truncation and number of grid points should be same.')
@@ -1505,8 +1493,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( im /= km ) then
        call MessageNotify('E','at_BoundariesGrid_ND', &
             'Chebyshev truncation and number of grid points should be same.')
@@ -1600,8 +1586,6 @@
     integer :: k
     save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
     if ( im /= km ) then
        call MessageNotify('E','at_BoundariesGrid_NN', &
             'Chebyshev truncation and number of grid points should be same.')
--- spml-g95.orig/src/eq_module.f90
+++ spml-g95/src/eq_module.f90
@@ -527,8 +527,6 @@
       integer :: k, l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(value)) then
          value1=0
       else
@@ -630,8 +628,6 @@
       save    :: alu, kp, first
       save    :: alub, kpb
 
-!$omp threadprivate(alu,kp,alub,kpb,first)
-
       if (.not. present(value)) then
          value1=0
       else
--- spml-g95.orig/src/et_galerkin_module.f90
+++ spml-g95/src/et_galerkin_module.f90
@@ -829,7 +829,6 @@
       logical :: first = .true.
       save LaplaInvMT, kpvot, first
 
-!$omp threadprivate(LaplaInvMT,kpvot,first)
 
       if ( .not. Set_RRFF ) &
            call MessageNotify('E','eq_Vor2Strm_et',&
--- spml-g95.orig/src/et_galerkin_module_fftj.f90
+++ spml-g95/src/et_galerkin_module_fftj.f90
@@ -827,7 +827,6 @@
       logical :: first = .true.
       save LaplaInvMT, kpvot, first
 
-!$omp threadprivate(LaplaInvMT,kpvot,first)
 
       if ( .not. Set_RRFF ) &
            call MessageNotify('E','eq_Vor2Strm_et',&
--- spml-g95.orig/src/et_module.f90
+++ spml-g95/src/et_module.f90
@@ -444,8 +444,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
@@ -528,8 +526,6 @@
       integer :: k,j
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
@@ -669,8 +665,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
@@ -790,8 +784,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
--- spml-g95.orig/src/et_module_fftj.f90
+++ spml-g95/src/et_module_fftj.f90
@@ -435,8 +435,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
@@ -519,8 +517,6 @@
       integer :: k,j
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
@@ -660,7 +656,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(values)) then
          value1=0 ; value2=0
@@ -781,8 +776,6 @@
       integer :: k,l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(values)) then
          value1=0 ; value2=0
       else
--- spml-g95.orig/src/wq_module.f90
+++ spml-g95/src/wq_module.f90
@@ -2010,7 +2010,6 @@
       integer :: n, l, lend
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid=.TRUE.
@@ -2158,7 +2157,6 @@
       integer :: k
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid=.TRUE.
@@ -2276,7 +2274,6 @@
       save    :: alu, kp, first
       save    :: alub, kpb
 
-!$omp threadprivate(alu,kp,alub,kpb,first)
 
       if (.not. present(cond)) then
          rigid=.TRUE. 
@@ -2469,7 +2466,6 @@
       save    :: alu, kp, first
       save    :: alub, kpb
 
-!$omp threadprivate(alu,kp,alub,kpb,first)
 
       if (.not. present(cond)) then
          rigid=.TRUE. 
@@ -2601,7 +2597,6 @@
       integer :: n, l, lend
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2738,7 +2733,6 @@
       integer  :: l, n, nn(2), lend
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2853,7 +2847,6 @@
       integer  :: k, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
--- spml-g95.orig/src/wt_galerkin_module.f90
+++ spml-g95/src/wt_galerkin_module.f90
@@ -303,7 +303,6 @@
       logical :: first = .true.
       save LaplaInvMT, kpvot, first
 
-!$omp threadprivate(LaplaInvMT,kpvot,first)
 
       if ( .not. Set_RRFF ) &
            call MessageNotify('E','wq_LaplaPol2Pol_wt',&
--- spml-g95.orig/src/wt_galerkin_module_sjpack.f90
+++ spml-g95/src/wt_galerkin_module_sjpack.f90
@@ -352,7 +352,6 @@
       logical :: first = .true.
       save LaplaInvMT, kpvot, first
 
-!$omp threadprivate(LaplaInvMT,kpvot,first)
 
       if ( .not. Set_RRFF ) &
            call MessageNotify('E','wq_LaplaPol2Pol_wt',&
--- spml-g95.orig/src/wt_module.f90
+++ spml-g95/src/wt_module.f90
@@ -2037,8 +2037,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
       else
@@ -2184,8 +2182,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
       else
@@ -2338,8 +2334,6 @@
       integer :: k,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
       else
@@ -2480,8 +2474,6 @@
       integer :: l,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
       else
@@ -2610,8 +2602,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(new)) then
          new_matrix=.false.
       else 
@@ -2698,8 +2688,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(new)) then
          new_matrix=.false.
       else 
@@ -2787,8 +2775,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(new)) then
          new_matrix=.false.
       else 
@@ -2879,8 +2865,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(new)) then
          new_matrix=.false.
       else 
--- spml-g95.orig/src/wt_module_sjpack.f90
+++ spml-g95/src/wt_module_sjpack.f90
@@ -2025,7 +2025,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2172,7 +2171,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2326,7 +2324,6 @@
       integer :: k,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2468,7 +2465,6 @@
       integer :: l,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2598,7 +2594,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2686,7 +2681,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2775,7 +2769,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2867,7 +2860,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
--- spml-g95.orig/src/wt_zonal_galerkin_module.f90
+++ spml-g95/src/wt_zonal_galerkin_module.f90
@@ -353,8 +353,6 @@
       logical :: first = .true.
       save LaplaInvMT, kpvot, first
 
-!$omp threadprivate(LaplaInvMT,kpvot,first)
-
       if ( .not. Set_RRFF ) &
            call MessageNotify('E','wq_LaplaPol2Pol_wt',&
                     'at_aq_galerkin_RRFF_module not initialized.')
--- spml-g95.orig/src/wt_zonal_module.f90
+++ spml-g95/src/wt_zonal_module.f90
@@ -2018,8 +2018,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
-
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
       else
@@ -2165,7 +2163,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2319,7 +2316,6 @@
       integer :: k,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2461,7 +2457,6 @@
       integer :: l,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2591,7 +2586,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2679,7 +2673,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2768,7 +2761,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2860,7 +2852,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
--- spml-g95.orig/src/wt_zonal_module_sjpack.f90
+++ spml-g95/src/wt_zonal_module_sjpack.f90
@@ -1997,7 +1997,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2144,7 +2143,6 @@
       integer :: l
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2298,7 +2296,6 @@
       integer :: k,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2440,7 +2437,6 @@
       integer :: l,n
       save    :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(cond)) then
          rigid1=.TRUE. ; rigid2=.TRUE.
@@ -2570,7 +2566,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2658,7 +2653,6 @@
       integer  :: l, n
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2747,7 +2741,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -2839,7 +2832,6 @@
       integer  :: l, n, nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
--- spml-g95.orig/src/wtu_module.f90
+++ spml-g95/src/wtu_module.f90
@@ -680,7 +680,6 @@
       integer  :: l
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -805,7 +804,6 @@
       integer  :: l
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -932,7 +930,6 @@
       integer  :: n,l,nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
@@ -1062,7 +1059,6 @@
       integer  :: n,l,nn(2)
       save     :: alu, kp, first
 
-!$omp threadprivate(alu,kp,first)
 
       if (.not. present(new)) then
          new_matrix=.false.
