Skip to content
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.

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

Possible out-dated LIBSTELL #6

Open
zhucaoxiang opened this issue May 13, 2020 · 1 comment
Open

Possible out-dated LIBSTELL #6

zhucaoxiang opened this issue May 13, 2020 · 1 comment
Assignees
Labels

Comments

@zhucaoxiang
Copy link
Collaborator

While compiling the code, I first used the full-size LIBSTELL that is up-to-date to the STELLOPT. I got the following error when running the test at examples/monoenergetic_geometryScheme5_ASCII.

 ****************************************************************************
 SFINCS: Stellarator Fokker-Plank Iterative Neoclassical Conservative Solver
 Version 3
 Using double precision.
 Parallel job (  16 processes) detected.
 Successfully read parameters from general namelist in input.namelist.
 Successfully read parameters from geometryParameters namelist in input.namelist
 .
 Successfully read parameters from speciesParameters namelist in input.namelist.
 Successfully read parameters from physicsParameters namelist in input.namelist.
 Successfully read parameters from resolutionParameters namelist in 
 input.namelist.
 Successfully read parameters from otherNumericalParameters namelist in 
 input.namelist.
 Successfully read parameters from preconditionerOptions namelist in 
 input.namelist.
 Successfully read parameters from export_f namelist in input.namelist.
 Setting Nxi_for_x_option=0, since RHSMode=3.
 Since RHSMode=3, ignoring the requested values of Zs, nHats, THats, nu_n, Er, a
 nd dPhiHatd*.
 mumps detected
 superlu_dist not detected
 ---- Physics parameters: ----
 Number of particle species =            1
 Delta (rho* at reference parameters)          =   4.569400000000000E-003
 alpha (e Phi / T at reference parameters)     =    1.00000000000000     
 nu_n (collisionality at reference parameters) =   8.330000000000001E-003
 Linear run
  Error No.           13  in READ_WOUT, iostat =           64
 Error opening wout file:
 ../../../../equilibria/wout_w7x_standardConfig.txt                             

Then I re-compiled the code the the native mini_libstell and everything is working correctly. My guess is that there are some updates in the read_wout_mod.f90 in STELLOPT and cause this incompatibility. I guess @lazersos might know more details.

@lazersos
Copy link
Collaborator

So I'm starting to dig into this and the issue appears to be that this text file is an oddball. I suspect it was generated by one of the old IPP routines and as such wasn't written by a supported version of VMEC, or may not have been written by VMEC at all. The issue seems to stem from a missing write of the lcurr variable before the names of the coils. I think the fix for now is to not use the text wout files.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

3 participants