I think these variables should be changed to something more descriptive, particularly vgap2.
real(dp) :: vgap2
!! vertical gap between vacuum vessel and thermal shields (m)
real(dp) :: vgap
!! vertical gap between x-point and divertor (m) (if = 0, it is calculated)
maybe something like
vgap2: vgap_vv_thermalshield
vgap: vgap_xpoint_divertor
I think these variables should be changed to something more descriptive, particularly vgap2.
real(dp) :: vgap2
!! vertical gap between vacuum vessel and thermal shields (m)
real(dp) :: vgap
!! vertical gap between x-point and divertor (m) (if = 0, it is calculated)
maybe something like
vgap2: vgap_vv_thermalshield
vgap: vgap_xpoint_divertor