homere_partel.f

Go to the documentation of this file.

!                   *********************
                    PROGRAM homere_partel
!                   *********************
!
!
!***********************************************************************
! PARALLEL   V6P2                                   21/08/2010
!***********************************************************************
!
!brief    PREPROCESSING STEP BEFORE A PARALLEL COMPUTATION
!
!history  Y. AUDOUIN (EDF)
!+        09/2012
!+        V7P0
!+        FIRST  VERSION
!
!history  V STOBIAC (EDF)
!+        09/2014
!+        V7P0
!+        READING OF THE MESH FORMAT
!
!history  Y. AUDOUIN (EDF)
!+        09/2012
!+        V7P0
!+        Reorganizing reading of parameters
!
!history  Y. AUDOUIN (EDF)
!+        21/05/2015
!+        V7P0
!+        Adding call to partel_resonly
!
!history C. COULET (ARTELIA)
!+       01/09/2016
!+       V7P2
!+       Modification to add the weir file management
!
!history  J. GRASSET (Daresbury Lab & EDF)
!+        01/05/2018
!+        Add code for managing concatenated mesh
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_partel
      USE declarations_special
      IMPLICIT NONE
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      CHARACTER(LEN=PATH_LEN) :: NAMEINP,NAMEGEO
      CHARACTER(LEN=PATH_LEN) :: NAMECLI
      INTEGER :: NPARTS
      INTEGER :: PMETHOD
      CHARACTER(LEN=PATH_LEN) :: NAMESEC
      CHARACTER(LEN=PATH_LEN) :: NAMEZFI
      CHARACTER(LEN=PATH_LEN) :: NAMESEU
      CHARACTER(LEN=8) :: GEOFORMAT,INPFORMAT
      CHARACTER(LEN=3) :: CONCATSTR
!
      LOGICAL :: RES_ONLY
      LOGICAL :: IS
      INTEGER I_LENCLI, I_LENINP, I_LENGEO
!
      CHARACTER(LEN=11) :: EXTENS
      EXTERNAL extens
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
!
      lng=2 ! JE NE PARLE PAS FRANCAIS, JE SUIS BARBARIEN
      lu=6  ! FORTRAN STANDARD OUPUT CHANNEL
      li=5  ! FORTRAN STANDARD INPUT CHANNEL
!
!----------------------------------------------------------------------
! INTRODUCE YOURSELF
!
      WRITE(lu,*) ' '
      WRITE(lu,*) '+-------------------------------------------------+'
      WRITE(lu,*) '  PARTEL/PARRES: TELEMAC METISOLOGIC PARTITIONER'
      WRITE(lu,*) '                                                   '
      WRITE(lu,*) '  REBEKKA KOPMANN & JACEK A. JANKOWSKI (BAW)'
      WRITE(lu,*) '                 JEAN-MICHEL HERVOUET (LNHE)'
      WRITE(lu,*) '                 CHRISTOPHE DENIS     (SINETICS) '
      WRITE(lu,*) '                 YOANN AUDOUIN        (LNHE) '
      WRITE(lu,*) '  PARTEL (C) COPYRIGHT 2000-2002 '
      WRITE(lu,*) '  BUNDESANSTALT FUER WASSERBAU, KARLSRUHE'
      WRITE(lu,*) ' '
      WRITE(lu,*) '  METIS 5.0.2 (C) COPYRIGHT 2012 '
      WRITE(lu,*) '  REGENTS OF THE UNIVERSITY OF MINNESOTA '
      WRITE(lu,*) ' '
      WRITE(lu,*) '  BIEF ',version,' (C) COPYRIGHT 2012 EDF'
      WRITE(lu,*) '+-------------------------------------------------+'
      WRITE(lu,*) ' '
      WRITE(lu,*) ' '
      WRITE(lu,*) '  MAXIMUM NUMBER OF PARTITIONS: ',maxnproc
      WRITE(lu,*) ' '
      WRITE(lu,*) '+--------------------------------------------------+'
      WRITE(lu,*) ' '
!
!----------------------------------------------------------------------
! NAMES OF THE INPUT FILES:
!
!
      WRITE(lu,*) '--INPUT FILE NAME <INPUT_NAME>: '
      READ(li,*) nameinp
!
      IF (nameinp.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(nameinp)
      ENDIF

      INQUIRE (file=nameinp,exist=is)
      IF (.NOT.is) THEN
        WRITE (lu,*)' FILE DOES NOT EXIST: ',trim(nameinp)
        CALL plante(1)
        stop
      ENDIF

      IF (code(1:1).EQ.' ') THEN
        code = nameinp(1:3)
      ENDIF
!
      WRITE(lu,*)
     & '--INPUT FILE FORMAT <INPFORMAT> [MED,SERAFIN,SERAFIND]: '
      READ(li,*) inpformat
      IF ( (inpformat .NE. 'MED     ') .AND.
     &     (inpformat(1:7) .NE. 'SERAFIN') ) THEN
        WRITE(lu,*)
     &  ' FILE FORMAT MUST BE "MED" OR "SERAFIN" OR "SERAFIND" '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) ' INPUT: ', inpformat
      ENDIF
!
      WRITE(lu, *) '--BOUNDARY CONDITIONS FILE NAME: '
      READ(li,*) namecli
      IF (namecli.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(namecli)
      ENDIF
!
      INQUIRE (file=namecli,exist=is)
      IF (.NOT.is) THEN
        WRITE (lu,*) 'FILE DOES NOT EXIST: ',trim(namecli)
        CALL plante(1)
        stop
      ENDIF
!
      WRITE(lu,fmt='(A,I6,A)')
     &  '--NUMBER OF PARTITIONS <NPARTS> [2 -',maxnproc,']: '
      READ(li,*) nparts
!
      IF ( (nparts > maxnproc) .OR. (nparts < 2) ) THEN
        WRITE(lu,fmt='(A,I6,A)')
     &  ' NUMBER OF PARTITIONS MUST BE IN [2 -',maxnproc,']'
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*)' INPUT: ', nparts
      ENDIF
!
      WRITE(lu,*) ' PARTITIONING METHOD <PMETHOD>
     &  [1 (METIS) OR 2 (SCOTCH)]: '
      READ(li,*) pmethod
      IF ( (pmethod > 2) .OR. (pmethod < 1) ) THEN
        WRITE(lu,
     &  '('' PARTITIONING METHOD MUST BE 1 OR 2'')')
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) '--INPUT: ', pmethod
      ENDIF

!
! #### THE SECTIONS FILE NAME
!
      namesec = ' '
      WRITE(lu,*)'--CONTROL SECTIONS FILE NAME (OR RETURN) : '
      READ(li,'(A)') namesec
!
      IF(namesec(1:1) .NE. ' ') THEN
        WRITE(lu,*) 'INPUT: ',trim(namesec)
        INQUIRE (file=namesec,exist=is)
        IF (.NOT.is) THEN
          WRITE (lu,*) ' FILE DOES NOT EXIST: ',trim(namesec)
          CALL plante(1)
          stop
        ENDIF
      ELSE
        WRITE(lu,*) ' NO SECTIONS '
      ENDIF
!
! #### THE ZONES FILE NAME
!
      namezfi = ' '
      WRITE(lu,*) '--CONTROL ZONES FILE NAME (OR RETURN) : '
      READ(li,'(A)') namezfi
!
      IF(namezfi(1:1) .NE. ' ') THEN
        WRITE(lu,*) 'INPUT: ',trim(namezfi)
        INQUIRE (file=namezfi,exist=is)
        IF (.NOT.is) THEN
          WRITE (lu,*) ' FILE DOES NOT EXIST: ', trim(namezfi)
          CALL plante(1)
          stop
        ENDIF
      ELSE
        WRITE(lu,*) ' NO ZONES '
      ENDIF
!
! #### THE WEIRS FILE NAME
!
      nameseu = ' '
      WRITE(lu,*) '--WEIR FILE NAME (OR RETURN) : '
      READ(li,'(A)') nameseu
!
      IF(nameseu(1:1) .NE. ' ') THEN
        WRITE(lu,*) 'INPUT: ',trim(nameseu)
        INQUIRE (file=nameseu,exist=is)
        IF (.NOT.is) THEN
          WRITE (lu,*) ' FILE DOES NOT EXIST: ', trim(nameseu)
          CALL plante(1)
          stop
        ENDIF
      ELSE
        WRITE(lu,*) ' NO WEIRS '
      ENDIF
!
!     Geometry file
!
      WRITE(lu,*) '--GEOMETRY FILE NAME <INPUT_NAME>: '
      READ(li,*) namegeo
!
      IF (namegeo.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(namegeo)
      ENDIF

      INQUIRE (file=namegeo,exist=is)
      IF (.NOT.is) THEN
        WRITE (lu,*)' FILE DOES NOT EXIST: ',trim(namegeo)
        CALL plante(1)
        stop
      ENDIF
!
      WRITE(lu,*)
     & '--GEOMETRY FILE FORMAT <GEOFORMAT> [MED,SERAFIN,SERAFIND]: '
      READ(li,*) geoformat
      IF ( (geoformat .NE. 'MED     ') .AND.
     &     (geoformat(1:7) .NE. 'SERAFIN') ) THEN
        WRITE(lu,*)
     &  ' FILE FORMAT MUST BE "MED" OR "SERAFIN" OR "SERAFIND" '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) ' INPUT: ', geoformat
      ENDIF
!
!     Concatenate files
!
      WRITE(lu,*) '--CONCATENATE FILES <YES-NO>: '
      READ(li,*) concatstr
!
      IF (concatstr.NE.'YES'.AND.concatstr.NE.'NO') THEN
        WRITE (lu,*) 'DONT KNOW IF I SHOUD CONCATENATE:'
     &               //trim(concatstr)
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'CONCATENATE: ',trim(concatstr)
        IF(concatstr=='YES')THEN
          partel_concat=.true.
        ELSE
          partel_concat=.false.
        ENDIF
      ENDIF
!
      ! Check if the geometry has already been partitionned
      IF(partel_concat)THEN
        INQUIRE(file=trim(namegeo)//'-CONCAT', exist=res_only)
      ELSE
        INQUIRE(file=trim(namegeo)//extens(nparts-1,0), exist=res_only)
      ENDIF
      IF(namegeo.EQ.nameinp) res_only = .false.

!
! FIND THE INPUT FILE CORE NAME LENGTH
!
      i_leninp = len(trim(nameinp))
!
      IF (i_leninp > path_len) THEN
        WRITE(lu,*) ' '
        WRITE(lu,*) 'ATTENTION:'
        WRITE(lu,*) 'THE NAME OF THE INPUT FILE:'
        WRITE(lu,*) nameinp
        WRITE(lu,*) 'IS LONGER THAN ',path_len,' CHARACTERS'
        WRITE(lu,*) 'WHICH IS THE LONGEST APPLICABLE NAME FOR TELEMAC '
        WRITE(lu,*) 'INPUT AND OUTPUT FILES. STOPPED. '
        CALL plante(1)
        stop
      ENDIF
!
! CORE NAME LENGTH
!
      i_lencli = len(trim(namecli))
!
      IF (i_lencli > path_len) THEN
        WRITE(lu,*) ' '
        WRITE(lu,*) 'ATTENTION:'
        WRITE(lu,*) 'THE NAME OF THE BOUNDARY CONDITIONS FILE:'
        WRITE(lu,*) namecli
        WRITE(lu,*) 'IS LONGER THAN ',path_len,' CHARACTERS'
        WRITE(lu,*) 'WHICH IS THE LONGEST APPLICABLE NAME FOR TELEMAC '
        WRITE(lu,*) 'INPUT AND OUTPUT FILES. STOPPED. '
        CALL plante(1)
        stop
      ENDIF
!
! FIND THE GEOMETRY FILE CORE NAME LENGTH
!
      i_lengeo = len(trim(namegeo))
!
      IF (i_lengeo > path_len) THEN
        WRITE(lu,*) ' '
        WRITE(lu,*) 'ATTENTION:'
        WRITE(lu,*) 'THE NAME OF THE GEOMETRY FILE:'
        WRITE(lu,*) namegeo
        WRITE(lu,*) 'IS LONGER THAN ',path_len,' CHARACTERS'
        WRITE(lu,*) 'WHICH IS THE LONGEST APPLICABLE NAME FOR TELEMAC '
        WRITE(lu,*) 'INPUT AND OUTPUT FILES. STOPPED. '
        CALL plante(1)
        stop
      ENDIF
!
      IF(res_only) THEN
        ! Reading geometry from xxxGEO
        CALL parres(namegeo, nameinp, nparts, geoformat, inpformat)
      ELSE
        CALL partel(nameinp, namecli, nparts, pmethod, inpformat,
     &              namesec, namezfi, nameseu)
      ENDIF
!
      stop 0
      END PROGRAM