diff --git a/source/fortran/iteration_variables.f90 b/source/fortran/iteration_variables.f90 index 5ed451636c..2bbc3e4319 100755 --- a/source/fortran/iteration_variables.f90 +++ b/source/fortran/iteration_variables.f90 @@ -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 - !!
  • (174) DUMMY : Description + !!
  • (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 - !!
  • (175) DUMMY : Description + !!
  • (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 + !!
  • (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 + !!
  • (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 !------------------------------------------------------ @@ -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) @@ -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) diff --git a/source/fortran/numerics.f90 b/source/fortran/numerics.f90 index e534bc6bf5..3ad503c852 100755 --- a/source/fortran/numerics.f90 +++ b/source/fortran/numerics.f90 @@ -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