Hermes error when assembling free surface result

Hi,

V9.0 with MUMPS

+> res_SL_est_par

Traceback (most recent call last):

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/runcode.py", line 303, in <module>

    sys.exit(main())

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/runcode.py", line 287, in main

    run_study(cas_file, code_name, options)

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/execution/run_cas.py", line 170, in run_study

    run_local_cas(my_study, options)

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/execution/run_cas.py", line 78, in run_local_cas

    my_study.gather(options.sortie_file, options.nozip)

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/execution/study.py", line 744, in gather

    files = process_artnim(cas, [name removed]File "/opt/telemac-mascaret/v9.0.0/scripts/python3/execution/process.py", line 488, in process_artnim

    amp2wave(file_name_amp, file_name_wfs,

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/data_manip/computation/amp2wave.py", line 39, in amp2wave

    wave.write()

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/data_manip/extraction/telemac_file.py", line 2008, in write

    self.write_mesh()

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/data_manip/extraction/telemac_file.py", line 1935, in write_mesh

    self.set_header(self._title, self._nvar, self._varnames,

  File "/opt/telemac-mascaret/v9.0.0/scripts/python3/telapy/api/hermes.py", line 749, in set_header

    self.error = HermesFile._hermes.set_header(self._fformat, self.my_id,

_hermes.error: (shape(var_name, 0) == nvar) failed for 1st keyword nvar: set_header:nvar=1

This kind of error is quite obscure to me…

Python 3.8.10

Numpy 1.24.4

Full model + temp compute folder : [name removed]

cas_est_par.txt (3.6 KB)

Hello [name removed],

I have successfully run what you have provided for 3 configurations (you can read the installation configuration in $HOMETEL/configs/systel.edf.cfg):

  • 8 cores for S11.gfortran.dyn configuration (python 3.9.17 + numpy 1.21.6),

  • 32 cores for gaia.intel.dyn configuration (python 3.8.18 + numpy 1.24.4),

  • 32 cores for cronos.intel.dyn configuration (python 3.8.18 + numpy 1.24.4).

Do the 2 examples $HOMETEL/examples/artemis/creocean/art_creocean200-anim.cas and $HOMETEL/examples/artemis/bj78/art_bj78_animated.cas work for you or not?

Have you installed hermes?

Hope this helps,

[name removed]

Hello [name removed],

Thank you for testing. Even the gathering part works for you ? You have the file “ res_SL_est_par” ?

This is weird.

Yes Hermes is compiled and works with T2D results files.

I will try to run the 2 examples when I have access to our cluster.

Hello [name removed],

Yes, the computations end correctly with the gathering of the result file res_SL_est_par (+ “my work is done” in the listing file).

[name removed]

Both examples fail with same error as in first post…

In that case, there must be at least 1 difference in the way openTELEMAC is installed.

I suppose you have the same issue if running gretel.py --input-file gretel_AMP.par or gretel_ARTRES.par in the temporary directory?

Your error seems to appear when running /opt/telemac-mascaret/v9.0.0/scripts/python3/telapy/api/hermes.py

I wonder if you have installed the API for Hermes.

In your systel configuration file (see e.g. $HOMETEL/configs/systel.edf.cfg), for options, I use options: mpi hpc api (in particular api).

Hope this helps,

[name removed]

Gretel.py for ARTAMP and ARTRES runs fine.

I provide our systel.cfg (the one in production is slightly different but the essentiel is in the provided one) but I made it based on edf.cfg so…

“api” option is specified in config

I tired both f2py and f2py3 but I think one is a simlink of the other

Similar to https://opentelemac.org/index.php/kunena/scripts/14652-hermes-error-when-merging-2-files-with-kind-time

systel-2.cfg (4.67 KB)

[name removed], do you have any clue what could be the issue here ? Or any tips how to debug ?

Well…

I moved our cluster to lmod instead of relying on shaky environment variable management.

It works now using python 3.12, numpy 1.26.4