Skip to content

Added triang and kappa as iteration variables #3563

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
timothy-nunn merged 1 commit into from
Mar 6, 2025
Merged

Added triang and kappa as iteration variables #3563

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
88 changes: 67 additions & 21 deletions source/fortran/iteration_variables.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3870,51 +3870,93 @@ subroutine set_itv_173(ratio)
f_tritium = ratio
end subroutine set_itv_173

!---------------------------------
! DUMMY variables below here
!---------------------------------

subroutine init_itv_174
!! <LI> (174) DUMMY : Description
!! <LI> (174) triang : Plasma triangularity
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(174) = 'DUMMY '
boundl(174) = 1.0d-99
boundu(174) = 1.0d99
lablxc(174) = 'triang '
boundl(174) = 0.00
boundu(174) = 1.00
end subroutine init_itv_174

real(kind(1.d0)) function itv_174()
implicit none
itv_174 = DUMMY
use physics_variables, only: triang
itv_174 = triang
end function itv_174

subroutine set_itv_174(ratio)
implicit none
use physics_variables, only: triang
real(kind(1.d0)) :: ratio
DUMMY = ratio
triang = ratio
end subroutine set_itv_174

!---------------------------------

subroutine init_itv_175
!! <LI> (175) DUMMY : Description
!! <LI> (175) kappa : Plasma elongation
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(175) = 'DUMMY '
boundl(175) = 1.0d-99
boundu(175) = 1.0d99
lablxc(175) = 'kappa '
boundl(175) = 0.00
boundu(175) = 10.00
end subroutine init_itv_175

real(kind(1.d0)) function itv_175()
implicit none
itv_175 = DUMMY
use physics_variables, only: kappa
itv_175 = kappa
end function itv_175

subroutine set_itv_175(ratio)
use physics_variables, only: kappa
real(kind(1.d0)) :: ratio
kappa = ratio
end subroutine set_itv_175

!---------------------------------
! DUMMY variables below here
!---------------------------------

subroutine init_itv_176
!! <LI> (174) DUMMY : Description
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(176) = 'DUMMY '
boundl(176) = 1.0d-99
boundu(176) = 1.0d99
end subroutine init_itv_176

real(kind(1.d0)) function itv_176()
implicit none
itv_176 = DUMMY
end function itv_176

subroutine set_itv_176(ratio)
implicit none
real(kind(1.d0)) :: ratio
DUMMY = ratio
end subroutine set_itv_175
end subroutine set_itv_176

!---------------------------------

subroutine init_itv_177
!! <LI> (175) DUMMY : Description
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(177) = 'DUMMY '
boundl(177) = 1.0d-99
boundu(177) = 1.0d99
end subroutine init_itv_177

real(kind(1.d0)) function itv_177()
implicit none
itv_177 = DUMMY
end function itv_177

subroutine set_itv_177(ratio)
implicit none
real(kind(1.d0)) :: ratio
DUMMY = ratio
end subroutine set_itv_177

!------------------------------------------------------

Expand Down Expand Up @@ -4115,9 +4157,11 @@ subroutine loadxc
case (171); xcm(i) = itv_171()
case (172); xcm(i) = itv_172()
case (173); xcm(i) = itv_173()
! DUMMY Cases
case (174); xcm(i) = itv_174()
case (175); xcm(i) = itv_175()
! DUMMY Cases
case (176); xcm(i) = itv_176()
case (177); xcm(i) = itv_177()

case default
idiags(1) = i ; idiags(2) = ixc(i)
Expand Down Expand Up @@ -4377,10 +4421,12 @@ subroutine convxc(xc,nn)
case (170); call set_itv_170(ratio)
case (171); call set_itv_171(ratio)
case (172); call set_itv_172(ratio)
! DUMMY Cases
case (173); call set_itv_173(ratio)
case (174); call set_itv_174(ratio)
case (175); call set_itv_175(ratio)
! DUMMY Cases
case (176); call set_itv_176(ratio)
case (177); call set_itv_177(ratio)

case default
call report_error(57)
Expand Down
2 changes: 1 addition & 1 deletion source/fortran/numerics.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@ module numerics

public

integer, parameter :: ipnvars = 175
integer, parameter :: ipnvars = 177
!! ipnvars FIX : total number of variables available for iteration

integer, parameter :: ipeqns = 92
Expand Down