Welcome, Guest
Username: Password: Remember me

TOPIC: how to give two liquid boundary files

how to give two liquid boundary files 8 years 2 weeks ago #24214

  • zqzuoan
  • zqzuoan's Avatar
Thans very much.

It works now. I adjust the order of the frtran file, the subroutine TR in the front, the subroutine SL after and it works.

bestwishes
zqzuoan
The administrator has disabled public write access.

how to give two liquid boundary files 8 years 2 weeks ago #24271

  • jmhervouet
  • jmhervouet's Avatar
Hello,

One must be aware also that there is a system of priority in the different ways of giving the boundary conditions. It is keywords like PRESCRIBED VALUES OF TRACERS in the parameter file that trigger the reading of the LIQUID BOUNDARY FILES, otherwise the program reverts to reading the BOUNDARY CONDITIONS FILE where generally there are only zeros.

With best regards,

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

how to give two liquid boundary files 6 years 4 months ago #30809

  • yasmina
  • yasmina's Avatar
Dear zqzuoan,

Can you tell what we should mofidie in the subroutine SL.f, i was already create my boundary file, but i don't know what i should modifie in this subroutine and i must change something too in soubroutine of BORD.f???

Hope that you can help me.
Thank you.

BEST REGARDS.
YASMIN.
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30810

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi
To answer your question, we must know what you want to do exactly!
Why do you want 2 liquid boundary files!
Christophe
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30811

  • yasmina
  • yasmina's Avatar
Sorry i did't present my project.

So i had done a simulation of tide in a global model (Gulf Normand-Breton). And now i want to extract the result of tide the elevation for instance to impose them in a small model.

I have already publied a topic for this, but when i checked the other topics i found that we must modifie the subroutine bord.f , this topic is different and more recent and it's need a modification of SL.f , so actualy i try all solutions.

Actualy i work in the modification of the bord.f with following the advices of sebourban in the topic #749.

Thank you.
BEST REGARDS.
YASMINA.
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30812

  • yasmina
  • yasmina's Avatar
And for your question, no i need only one liquid boundary varying in time and space. This liquid boundary will represent my tide.
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30814

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
So you've got one method which works well for this kind of problems

Regards
Christophe
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30818

  • yasmina
  • yasmina's Avatar
The modification of the subroutine bord.f????

Okay i will trie this method.

Regards.
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30823

  • yasmina
  • yasmina's Avatar
Hello,
I have done what you adviced me, thank you for this.

So i had launched the simulation for 1000s just to verify that my subroutine works. The simulation works!!!! and the result seems correct juste i must launch it for more time.
I attached the output of the simulation:

MASTER PROCESSOR NUMBER            0  OF THE GROUP OF            4
 EXECUTABLE FILE: C:\Stage_yasmina\Modélisation_hydrodynamique\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\A.EXE
 BARRIER PASSED
 LISTING OF TELEMAC2D------------------------------------------------------------------------------

                TTTTT  EEEEE  L      EEEEE  M   M  AAAAA  CCCCC
                  T    E      L      E      MM MM  A   A  C
                  T    EEE    L      EEE    M M M  AAAAA  C
                  T    E      L      E      M   M  A   A  C
                  T    EEEEE  LLLLL  EEEEE  M   M  A   A  CCCCC

                        2D    VERSION V7P3   FORTRAN 2003

                        ~^~^~^~^~^~^~^~^~^~^~^^~^~^~^~^~^~
                          ~                            ~
                               \   '    o      '
                               /\ o       \  o
                             >=)'>    '   /\ '
                               \/   \   >=)'>        ~
                               /    /\    \/
                        ~         >=)'>   /     .
                                    \/                   )
                                    /                   (
                                          ~          )   )
                          }     ~              (    (   (
                         {                      )    )   )
                          }  }         .       (    (   (
                         {  {               /^^^^^^^^^^^^
                        ^^^^^^^^^\         /
                                  ^^^^^^^^^

 WARNING IN DICTIONARY:
 FOR KEYWORD: VARIABLES FOR GRAPHIC PRINTOUTS
 THE NUMBER OF DEFAULT VALUES            1  IS DIFFERENT FROM THE DECLARED SIZE            4

 END OF FILE FOR DAMOCLES

                   ********************************************
                   *               LECDON:                    *
                   *        AFTER CALLING DAMOCLES            *
                   *        CHECKING OF DATA  READ            *
                   *         IN THE STEERING FILE             *
                   ********************************************

 EXITING LECDON. NAME OF THE STUDY:

 OPENING FILES FOR TELEMAC2D
  - FILE    C:\Stage_yasmina\Modélisation_hydrodynamique\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\T2DGEO00003-00000 OF TYPE BIG_ENDIAN
  - FILE    C:\Stage_yasmina\Modélisation_hydrodynamique\medium_box\Telemac(large_to_small)\init2_baie.cas_2018-07-05-17h46min50s\T2DGEO00003-00000 IS SINGLE PRECISION

                          *****************************
                          *    MEMORY ORGANIZATION    *
                          *****************************

 READ_MESH_INFO: TITLE= newSelafin
            NUMBER OF ELEMENTS:   102997
            NUMBER OF POINTS:    53446

            FORMAT NOT INDICATED IN TITLE

 MXPTEL (BIEF) : MAXIMUM NUMBER OF ELEMENTS AROUND A POINT:   9
                 MAXIMUM NUMBER OF POINTS AROUND A POINT:  10
 (GLOBAL MESH)
 SEGBOR (BIEF) : NUMBER OF BOUNDARY SEGMENTS =   3911
 INCLUDING THOSE DUE TO DOMAIN DECOMPOSITION
 CORRXY (BIEF):NO MODIFICATION OF COORDINATES

 MESH: MESH   ALLOCATED

                     *************************************
                     *    END OF MEMORY ORGANIZATION:    *
                     *************************************

 INITIALISING TELEMAC2D FOR
 INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)
 STRCHE (BIEF): NO MODIFICATION OF FRICTION

 NUMBER OF LIQUID BOUNDARIES:           1

 CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM

================================================================================
 ITERATION        0    TIME:   0.0000 S
 TELEMAC2D INITIALISED
 USING STREAMLINE VERSION 7.3 FOR CHARACTERISTICS
 PROSOU : IN CARTESIAN COORDINATES, THE CORIOLIS
          PARAMETER IS READ IN THE STEERING FILE
          IT IS THE KEY WORD 'CORIOLIS
          COEFFICIENT', IT IS UNIFORM IN SPACE
 GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS:     100 RELATIVE PRECISION:   0.1404821E-03
 GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS:     100 RELATIVE PRECISION:   0.1539796E-03
 GRACJG (BIEF) : EXCEEDING MAXIMUM ITERATIONS:     100 RELATIVE PRECISION:   0.1154087E-03

================================================================================
 ITERATION       50    TIME:  8 MN  20.0000 S   (      500.0000 S)
--------------------------------------------------------------------------------
                          ADVECTION STEP
--------------------------------------------------------------------------------
                    DIFFUSION-PROPAGATION STEP
 GRACJG (BIEF) :       34 ITERATIONS, RELATIVE PRECISION:   0.9383331E-04
--------------------------------------------------------------------------------
                       BALANCE OF WATER VOLUME
     VOLUME IN THE DOMAIN :   0.9516157E+08 M3
     FLUX BOUNDARY    1:    -64721.91     M3/S  ( >0 : ENTERING  <0 : EXITING )
     RELATIVE ERROR IN VOLUME AT T =        500.0     S :    0.2752341E-06

================================================================================
 ITERATION      100    TIME: 16 MN  40.0000 S   (     1000.0000 S)
--------------------------------------------------------------------------------
                          ADVECTION STEP
--------------------------------------------------------------------------------
                    DIFFUSION-PROPAGATION STEP
 GRACJG (BIEF) :       35 ITERATIONS, RELATIVE PRECISION:   0.9594290E-04
--------------------------------------------------------------------------------
                       BALANCE OF WATER VOLUME
     VOLUME IN THE DOMAIN :   0.7236964E+08 M3
     FLUX BOUNDARY    1:    -31629.34     M3/S  ( >0 : ENTERING  <0 : EXITING )
     RELATIVE ERROR IN VOLUME AT T =        1000.     S :    0.1036767E-05
--------------------------------------------------------------------------------
                   FINAL BALANCE OF WATER VOLUME

     RELATIVE ERROR CUMULATED ON VOLUME:   -0.2310381E-04

     INITIAL VOLUME              :    0.1713690E+09 M3
     FINAL VOLUME                :    0.7236964E+08 M3
     VOLUME THAT ENTERED THE DOMAIN:   -0.9900332E+08 M3  ( IF <0 EXIT )
     TOTAL VOLUME LOST             :    -3959.276     M3

 END OF TIME LOOP

 EXITING MPI

 CORRECT END OF RUN
STOP 0

STOP 0
 ELAPSE TIME :
STOP 0
                             30  SECONDS
STOP 0



The problem after the modification of the time step from 100 to 2000, the simulation don't work. In the begenning it's indicate me that the liquid boundary file has a tabulation replace it with space, wich i was done. But now it's indicate that there is some problems in the subroutine SL.


ITERATION        0    TIME:   0.0000 S
 TELEMAC2D INITIALISED
 THE LIQUID BOUNDARIES FILE CONTAINSapplication called MPI_Abort(MPI_COMM_WORLD,
 2) - process 0

           6  LINES WITH:
 ETA(34331ETA(34288ETA(34240ETA(34197ETA(34147ETA(34094ETA(34043ETA(33989ETA(339
33ETA(33876ETA(33821ETA(33764ETA(33707ETA(33651ETA(33592ETA(33533ETA(33473ETA(33
414ETA(33353ETA(33291ETA(33230ETA(33168ETA(33104ETA(33041ETA(32974ETA(32907ETA(3
2839ETA(32773ETA(32707ETA(32639ETA(32573ETA(32506ETA(32439ETA(32372ETA(32306ETA(
32240ETA(32173ETA(32107ETA(32044ETA(31979ETA(31914ETA(31051ETA(31785ETA(31720ETA
(31054ETA(31588ETA(31521ETA(31459ETA(31398ETA(31338ETA(31278ETA(31220ETA(31162ET
A(31104ETA(31048ETA(30992ETA(30936ETA(30883ETA(30830ETA(30779ETA(30727ETA(30677E
TA(30629ETA(30582ETA(30535ETA(30490ETA(30447ETA(30403ETA(30361ETA(30321ETA(30281
ETA(30243ETA(30205ETA(30170ETA(30137ETA(30103ETA(30074ETA(30043ETA(30014ETA(2998
6ETA(29960ETA(29938ETA(29913ETA(20801ETA(29872ETA(29852ETA(29834ETA(29819ETA(298
44ETA(29826ETA(29809ETA(29802ETA(29806

 SL: MORE PRESCRIBED ELEVATIONS ARE REQUIRED
      IN THE PARAMETER FILE
      AT LEAST      1 MUST BE GIVEN

 PLANTE: PROGRAM STOPPED AFTER AN ERROR
 RETURNING EXIT CODE:            2
job aborted:
rank: node: exit code[: error message]

Now i don't understand what it can be the solution.
I attached my liquid boundary file:


File Attachment:

File Name: liq.txt
File Size: 12 KB


And this is my fortran file:

File Attachment:

File Name: bord.f
File Size: 30 KB


Hope that you can help me in this error.

Best regards.
YASMINA.
The administrator has disabled public write access.

how to give two liquid boundary files 6 years 4 months ago #30825

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
This is not possible without any other change.
The error message you obtain:
SL: MORE PRESCRIBED ELEVATIONS ARE REQUIRED
IN THE PARAMETER FILE
AT LEAST 1 MUST BE GIVEN
is commented in your fortran file so it can't appear except if the fortran is not used!

Regards
Christophe
The administrator has disabled public write access.
Moderators: pham

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