added tglf_p_prime_elite_in and call to put_elite, made zgeev default

gafusion · Dec 12, 2024 · d78ebc4 · d78ebc4
1 parent ee9fb09
commit d78ebc4
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 2 deletions.
diff --git a/tglf/src/tglf_eigensolver.f90 b/tglf/src/tglf_eigensolver.f90
@@ -14,7 +14,7 @@ SUBROUTINE tglf_eigensolver
       IMPLICIT NONE
 !
       CHARACTER(1) :: rightvectors
-      LOGICAL :: use_zgeev = .false.
+      LOGICAL :: use_zgeev = .true.
       INTEGER :: is,js
       INTEGER ::  j1, j2, j, i
       INTEGER :: ifail,ib, jb, ia, ja, ia0, ja0

diff --git a/tglf/src/tglf_interface.f90 b/tglf/src/tglf_interface.f90
@@ -200,6 +200,7 @@ module tglf_interface
   integer  :: tglf_n_elite_in   = 700
   real     :: tglf_q_elite_in   = 2.0
   real     :: tglf_q_prime_elite_in = 16.0
+  real     :: tglf_p_prime_elite_in = 0.0
   real     :: tglf_R_elite_in(max_ELITE) 
   real     :: tglf_Z_elite_in(max_ELITE) 
   real     :: tglf_Bp_elite_in(max_ELITE) 

diff --git a/tglf/src/tglf_run.F90 b/tglf/src/tglf_run.F90
@@ -147,7 +147,16 @@ subroutine tglf_run()
           tglf_nfourier_in, &
           tglf_fourier_in)
 
-  else
+  elseif (tglf_geometry_flag_in == 3) then
+
+      call put_ELITE_geometry(tglf_n_elite_in,  &
+           tglf_q_elite_in,  &
+           tglf_q_prime_elite_in,  &
+           tglf_p_prime_elite_in,  &
+           tglf_R_elite_in,  &
+           tglf_Z_elite_in, &
+           tglf_Bp_elite_in  )
+  else   ! default to s-alpha
 
      call put_s_alpha_geometry(tglf_rmin_sa_in, &
           tglf_rmaj_sa_in, &