[BUG] Stellarator constraint 65 infinity

[wrutten]

Describe the bug

Process encounters divide by zero error error. Temperature profiles are not computed (tprofile =0)

Steps to reproduce

Steps to reproduce the behaviour:

1. Using input file failing_IN.DAT.TXT
2. Add stellarator configuration file failing_stella_conf.json.TXT (remove .txt extension)
   ation)
3. Run process
4. See error
   'RuntimeWarning: divide by zero encountered in log
   zav_of_te = yi + c * (numpy.log(tprofile) - xi)'
   ID LEVEL MESSAGE
   14 3 CONSTRAINTS: NaN/infty error for constraint equation...
   Integer diagnostic values for this error:
   '1) 65
   Floating point diagnostic values for this error:
   '1) -Infinity

Expected behaviour

Screenshots / Tracebacks

This input file works on (commit 1d6905a68 origin/ARCHIVED_last_version_with_stellarator)

Environment (please complete the following inform

• OS: WSL: Ubuntu-20.04
• Version: v3.0.0-6-g01d10b0
• Branch: main

Additional context

This is an old stellarator input file, it worked on the archived stellarator version on GitLab.
(commit 1d6905a68 origin/ARCHIVED_last_version_with_stellarator)

Error occurred trying to get this input file working with new open source stellarator version

[timothy-nunn]

These divide by zeros don't appear to be causing this issue, however if you do ever need the profiles to be calculated you can add self.plasma_profile.run() in Stellarator.run which will calculate the profiles.

This issue appears to be with my VV stress calculations not being called, so defaulting to 0. I will try and add the calculation into the Stellarator, but i dont know how much sense it will make

[wrutten]

Discussed with Tim
It appears that the stellarator model does it's own VV stress calculation (using eq. 2.88 in J. Lion 2023 phd thesis). The resulting f_vv_actual is given in [NM/m3]. To use this in constraint 65 it should be converted to [Pa] using the vacuum thickness d_vv_in.

[je-cook]

This is now used in the refactor of stellarator in #4082 therefore it looks like its fixed