search

msotov / SPECIES

Automatic code for the computation of stellar parameters using echelle spectra
MIT License
6 stars 0 forks source link

SPECIES Can't Find Temporary Files #1

Open spencerhurt opened 3 years ago

spencerhurt commented 3 years ago

I've tried running SPECIES on several different spectra, including the example solar spectrum. After calculating the line widths, the code always throws an error with the following output in the log:

2021-04-15 18:51:11 INFO run_iteration: sun01_harps 2021-04-15 18:51:11 INFO run_iteration: parent process: 87956 2021-04-15 18:51:11 INFO run_iteration: process id: 90785 2021-04-15 18:51:11 INFO run_iteration: Created EW file. 2021-04-15 18:51:11 INFO run_iteration: Alias for sun01_harps is vevwkg 2021-04-15 18:51:11 INFO run_iteration: Attempting to find colors for sun01_harps 2021-04-15 18:51:11 INFO run_iteration: Using colors only from Vizier. 2021-04-15 18:51:11 INFO run_iteration: Photometry information is: 2021-04-15 18:51:11 INFO run_iteration: name : sun01 2021-04-15 18:51:11 INFO run_iteration: Luminosity class is dwarf, from photometry. 2021-04-15 18:51:11 INFO run_iteration: Initial metallicity is 0.000000 2021-04-15 18:51:11 WARNING run_iteration: No valid temperature from photometry was computed. Check the photometric information of your star 2021-04-15 18:51:11 INFO write_log: hold_m=False, hold_t=False, hold_p=False, hold_v=False 2021-04-15 18:51:11 INFO write_log: Boundaries are: metallicity = (-3.00, 1.00), temperature = (3500.00, 9000.00), gravity = (0.50, 4.90), velocity = (0.00, 5.00) 2021-04-15 18:51:11 INFO write_log: Initial values are: feh=0.00, T=5500, logg=4.36, vt=1.23 2021-04-15 18:51:12 ERROR : sun01_harps: ERROR IN CODE, STOPPING THE COMPUTATION 2021-04-15 18:51:12 ERROR : line 322: [Errno 2] No such file or directory: './output/vevwkg_out.test' 2021-04-15 18:51:12 ERROR : <class 'FileNotFoundError'> 2021-04-15 18:51:12 ERROR : [Errno 2] No such file or directory: './output/vevwkg_out.test' 2021-04-15 18:51:12 ERROR : <traceback object at 0x7fbfa4c0e6e0>

While running the code, nothing appears in the output directory. I've tried pinning down why the ./output/{alias}_out.test file isn't generated but haven't found anything and was hoping you could help. Thanks in advance!

msotov commented 3 years ago

Hi Spencer,

If no output file is being generated then it could be a problem with the MOOG installation. Please type 'MOOGSILENT' in your console, if it executes something then it means the problem is somewhere else, and if it doesn't then it means MOOG is not correctly installed. I'll do my best to help with both cases.

spencerhurt commented 3 years ago

Typing MOOGSILENT into my console doesn't return anything. I've tried going through the steps outlined in the wiki to build MOOG but after running make -f Makefile.maclapsilent, the following error is returned:

OpacHydrogen.f:163:50:

  163 |          call linter (wfflog,fflog(1,itheta),nnnn,wavelog,fftlog,1)
      |                                                  1
Error: Rank mismatch in argument 'xnew' at (1) (rank-1 and scalar)
OpacHydrogen.f:169:47:

  169 |      .              maxwave = map1(wbf,bf,nnnn,wave,hminbf,1)
      |                                               1
Error: Rank mismatch in argument 'xnew' at (1) (rank-1 and scalar)
make: *** [OpacHydrogen.o] Error 1
afgupta commented 3 years ago

Spencer,

I encountered the same issue when installing SPECIES & MOOG a few weeks ago. It seems the syntax in the OpacHydrogen.f file is no longer supported with some fortran compilers. If you're using gfortran, you should be able to get it to work by adding the -fallow-argument-mismatch flag to the make -f Makefile.maclapsilent call.

See the related issue & solution here: https://github.com/scipy/scipy/issues/11611