Skip to content

Added casths_fraction and casths as iteration variables #3173

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 2 commits into from
May 8, 2024
Merged

Added casths_fraction and casths as iteration variables #3173

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
2102115
Added casths_fraction and casths as iteration variables
chris-ashe May 2, 2024
77ecc1d
comment update
chris-ashe May 8, 2024
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
47 changes: 25 additions & 22 deletions source/fortran/iteration_variables.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3804,52 +3804,55 @@ subroutine set_itv_170(ratio)
beta_div = ratio
end subroutine set_itv_170

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

subroutine init_itv_171
!! <LI> (171) DUMMY : Description
!! <LI> (171) casths_fraction: TF side case fraction of toroidal thickness
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(171) = 'DUMMY '
boundl(171) = 1.0d-99
boundu(171) = 1.0d99
lablxc(171) = 'casths_fraction '
boundl(171) = 0.01
boundu(171) = 0.99
end subroutine init_itv_171

real(kind(1.d0)) function itv_171()
use tfcoil_variables, only: casths_fraction
implicit none
itv_171 = DUMMY
itv_171 = casths_fraction
end function itv_171

subroutine set_itv_171(ratio)
use tfcoil_variables, only: casths_fraction
implicit none
real(kind(1.d0)) :: ratio
DUMMY = ratio
casths_fraction = ratio
end subroutine set_itv_171

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

subroutine init_itv_172
!! <LI> (172) DUMMY : Description
!! <LI> (172) casths: TF side case thickness [m]
use numerics, only: lablxc, boundl, boundu
implicit none
lablxc(172) = 'DUMMY '
boundl(172) = 1.0d-99
boundu(172) = 1.0d99
lablxc(172) = 'casths '
boundl(172) = 0.001
boundu(172) = 1.0
end subroutine init_itv_172

real(kind(1.d0)) function itv_172()
use tfcoil_variables, only: casths
implicit none
itv_172 = DUMMY
itv_172 = casths
end function itv_172

subroutine set_itv_172(ratio)
use tfcoil_variables, only: casths
implicit none
real(kind(1.d0)) :: ratio
DUMMY = ratio
casths = ratio
end subroutine set_itv_172


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

subroutine init_itv_173
Expand Down Expand Up @@ -4112,10 +4115,10 @@ subroutine loadxc
case (167); xcm(i) = itv_167()
case (168); xcm(i) = itv_168()
case (169); xcm(i) = itv_169()
! DUMMY Cases
case (170); xcm(i) = itv_170()
case (171); xcm(i) = itv_171()
case (172); xcm(i) = itv_172()
case (171); xcm(i) = itv_171()
case (172); xcm(i) = itv_172()
! DUMMY Cases
case (173); xcm(i) = itv_173()
case (174); xcm(i) = itv_174()
case (175); xcm(i) = itv_175()
Expand Down Expand Up @@ -4379,10 +4382,10 @@ subroutine convxc(xc,nn)
case (167); call set_itv_167(ratio)
case (168); call set_itv_168(ratio)
case (169); call set_itv_169(ratio)
! DUMMY Cases
case (170); call set_itv_170(ratio)
case (171); call set_itv_171(ratio)
case (172); call set_itv_172(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)
Expand Down
4 changes: 2 additions & 2 deletions source/fortran/numerics.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -372,8 +372,8 @@ module numerics
!! <LI> (168) fecrh_ignition: f-value for equation 91
!! <LI> (169) te0_ecrh_achievable: Max. achievable electron temperature at ignition point
!! <LI> (170) beta_div : field line angle wrt divertor target plate (degrees)
!! <LI> (171) EMPTY : Description
!! <LI> (172) EMPTY : Description
!! <LI> (171) casths_fraction : TF side case thickness as fraction of toridal case thickness
!! <LI> (172) casths : TF side case thickness [m]
!! <LI> (173) EMPTY : Description
!! <LI> (174) EMPTY : Description
!! <LI> (175) EMPTY : Description
Expand Down