Welcome, Guest
Username: Password: Remember me
  • Page:
  • 1
  • 2

TOPIC: Read non-erodible map from geometry file (noerod.f)

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10258

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Hello all,

I am trying to code noerod.f so as to read the depth of the sediment from a second variable in the geometry file. I am using v6p3r1 compiled with gfortran and want to run a tomawac-telemac2d-sisyphe coupled model in parallel.

I tried the noerod.f coding suggested on either of these threads:
Definition of non-erodable points
parallel runs with noerod.f

However, when the computation reaches the compilation of the noerod.f file, I get this error:
USE DECLARATIONS_TELEMAC2D
1
Fatal Error: Can't open module file 'declarations_telemac2d.mod' for reading (1): No such file or directory

Is it an installation issue? I have browsed the lib directory and the 'declarations_telemac2d.mod' exists. Is it because the structure of the installation has changed since the previous versions?

Any help would be appreciated.
Costas
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10266

  • cath
  • cath's Avatar
  • OFFLINE
  • Moderator
  • Posts: 52
  • Thank you received: 24
Hi

I assume that you are working in a coupled mode (telemac 2d and Siyphe).
In any case, it is not a good idea to have access to all the Telemac2d since there may be some compatibility issues (between telemac2d and Sisyphe, like for ex variables with same name and different types...).
It would be better to use only the variable you really need to access.
This restricted access to the Telemac-2d variables
USE DECLARATIONS_TELEMAC2D, ONLY: ...
Good luck Cath
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10268

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Dear Cath,

Yes, the aim is to eventually use it in coupled mode. The code I found use the following syntax:
USE BIEF
USE DECLARATIONS_TELEMAC
USE DECLARATIONS_TELEMAC2D, ONLY: T2D_FILES, T2DGEO
USE DECLARATIONS_SISYPHE , ONLY : IELMT, MESH
which employs the restrictions you suggest.

The problem however lies in the 'USE DECLARATIONS_TELEMAC2D'. It seems that the system cannot locate the relevant file.

Regards,
Costas
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10273

  • cath
  • cath's Avatar
  • OFFLINE
  • Moderator
  • Posts: 52
  • Thank you received: 24
If you are not in a couple mode, it is normal that the declarations of Telemac-2d are not recognised.
Why don't you use the Sisyphe geometry file instead.
Hope this helps, cath
The administrator has disabled public write access.
The following user(s) said Thank You: cyamin

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10279

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Hi Cath,

Thank you very much for the tip, I have no experience in programming and I am struggling to find my way around telemac routines.

I managed to sort this issue and the next one regarding the FFORMAT addition in FIND_IN_SEL for v6p3. Now the file compiles successfully but the program fails to run.
Debugger shows that it fails when the program calls the NOEROD subroutine. The error is:
LIT: READ ERROR
ONE INTENED TO READ A RECORD OF 1 VALUES OF TYPE: R8
ON LOGICAL UNIT: 1

This is the main routine for reading the non-erodible bed from the geometry file:
CALL BIEF_ALLVEC(1,NOER,'NOER ',IELMT,1,2,MESH)
ALLOCATE(W(NPOIN),STAT=ERR)
NAME='NOER '
CALL FIND_IN_SEL(NOER,NAME,SIS_FILES(SISGEO)%LU,'SERAFIND',W,OK)
CALL OV( 'X=Y-Z ',ZR,ZF,NOER%R,0.D0,NPOIN)
DEALLOCATE(W)
What am I doing wrong?

Many thanks in advance!
Costas
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10282

  • jmhervouet
  • jmhervouet's Avatar
Hello,

You are not far from success, what I see is that the name 'NOER ' must be given in 16 characters, so add spaces, and the format is most probably 'SERAFIN ' without D because nearly all geometry files are saved in single precision (hence the error message).

Regards,

Jean-Michel Hervouet
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10294

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Hello Jean-Michel,

Thank you very much for your help! I changed 'SERAFIND' to 'SERAFIN ' and it worked straight away. Results were as expected also.

In order to run it in coupled mode (Tomawac-Telemac2D-Sisyphe), apart from changing the declarations from SIS to T2D, I also had to specify the fortran file in the telemac2d cas file. Hence I assumed that all fortran sub-routines have to be assembled in this fortran file.

On the other hand, in the littoral example, there exist separate fortran files for Telemac2D and Tomawac. So one assumes that Tomawac has its own fortran file and Telemac2D and Sisyphe share the same (specified in Telemac2D). Is that correct?

Regards,
Costas
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10296

  • jmhervouet
  • jmhervouet's Avatar
Hello,

Yes, you may keep different fortran files and they are merged somewhere by the scripts. You may also merge them yourself in a single file for Telemac-2D (if I remember well the configuration "no file for telemac-2D and one file for Tomawac or Sisyphe" will not work, at least it did not work with Perl).

Regards,

Jean-Michel Hervouet
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 11 years 2 months ago #10297

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
It doesn't work. I had only a fortran file for sisyphe and I had to specify it in Telemac2D for the computation to start. It is much clearer now, thank you.

Costas
The administrator has disabled public write access.

Read non-erodible map from geometry file (noerod.f) 10 years 3 months ago #14028

  • cschwarz
  • cschwarz's Avatar
Dear all,

I have a problem with non-erodible beds related to the thread above.
First of all, I want to thank you very much for the nice explanations. I managed to changed my noerod.f-file according to it.

I am running a coupled telemac2d-sisyphe computation where I want to set part of the bed to be non-erodible. All seems to work fine, until the coupling when a Windows message appears(see WindosM.jpg) and the simulation stops(TelemacM.jpg).
I also added my Telemac2d steering file and the noerodBK.f file.

Do you have any suggestions?

Many thanks in advance,
Christian
Attachments:
The administrator has disabled public write access.
  • Page:
  • 1
  • 2
Moderators: Pablo, pavans

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.