homere_gretel.f

Go to the documentation of this file.

!                   *********************
                    PROGRAM homere_gretel
!                   *********************
!
!
!***********************************************************************
! PARALLEL   V6P2                                   21/08/2010
!***********************************************************************
!
!brief    MERGES THE RESULTS OF A PARALLEL COMPUTATION
!+                TO WRITE A SINGLE FILE.
!
!history  Y. AUDOUIN (EDF)
!+        09/2012
!+        V7P0
!+        Reorganizing reading of parameters
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      USE declarations_partel
      IMPLICIT NONE
!
!-------------------------------------------------------------------------
!
      LOGICAL IS
      CHARACTER(LEN=PATH_LEN) :: GEO
      CHARACTER(LEN=PATH_LEN) :: BND
      CHARACTER(LEN=8)   :: GEOFORMAT,RESFORMAT
      CHARACTER(LEN=PATH_LEN) :: RES
      INTEGER            :: NPARTS, NPLAN, METHOD
!
!-------------------------------------------------------------------------
!
      li=5
      lu=6
      lng=2
!
!----------------------------------------------------------------------
! INTRODUCE YOURSELF
!
      WRITE(lu,*) ' '
      WRITE(lu,*) '+-------------------------------------------------+'
      WRITE(lu,*) '  GRETEL: TELEMAC MERGER'
      WRITE(lu,*) '                                                   '
      WRITE(lu,*) '  VERSION ',version
      WRITE(lu,*) '  HOLGER WEIL BEER (BAW)'
      WRITE(lu,*) '  JEAN-MICHEL HERVOUET (LNHE)'
      WRITE(lu,*) '  YOANN AUDOUIN        (LNHE) '
      WRITE(lu,*) '  GRETEL (C) COPYRIGHT 2003-2012 '
      WRITE(lu,*) '  BUNDESANSTALT FUER WASSERBAU, KARLSRUHE'
      WRITE(lu,*) ' '
      WRITE(lu,*) '+-------------------------------------------------+'
      WRITE(lu,*) ' '
      WRITE(lu,*) ' '
      WRITE(lu,*) '  MAXIMUM NUMBER OF PARTITIONS: ',maxnproc
      WRITE(lu,*) ' '
      WRITE(lu,*) '+--------------------------------------------------+'
      WRITE(lu,*) ' '
!
      WRITE (lu,*) '--GLOBAL GEOMETRY FILE: '
      READ(li,*) geo
!
      IF (geo.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(geo)
      ENDIF

      INQUIRE (file=geo,exist=is)
      IF (.NOT.is) THEN
        WRITE (lu,*)' FILE DOES NOT EXIST: ',trim(geo)
        CALL plante(1)
        stop
      ENDIF
!
      WRITE(lu,*)
     & '--GEOMETRY FILE FORMAT <FFORMAT> [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
!
      WRITE (lu,*) '--GLOBAL BOUNDARY FILE: '
      READ(li,*) bnd
!
      IF (bnd.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(bnd)
      ENDIF

      INQUIRE (file=bnd,exist=is)
      IF (.NOT.is) THEN
        WRITE (lu,*)' FILE DOES NOT EXIST: ',trim(bnd)
        CALL plante(1)
        stop
      ENDIF
!
      WRITE(lu,*) '--RESULT FILE: '
      READ(li,*) res
!
      IF (res.EQ.' ') THEN
        WRITE (lu,*) ' NO FILENAME '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) 'INPUT: ',trim(res)
      ENDIF
!
      WRITE(lu,*)
     & '--RESULT FILE FORMAT <FFORMAT> [MED,SERAFIN,SERAFIND]: '
      READ(li,*) resformat
      IF ( (resformat .NE. 'MED     ') .AND.
     &     (resformat(1:7) .NE. 'SERAFIN') ) THEN
        WRITE(lu,*)
     &  ' FILE FORMAT MUST BE "MED" OR "SERAFIN" OR "SERAFIND" '
        CALL plante(1)
        stop
      ELSE
        WRITE(lu,*) ' INPUT: ', resformat
      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,*) '--NUMBER OF PLANES: '
      READ (li,*) nplan
      WRITE (lu,*) ' INPUT: ',nplan

      WRITE (lu,*) '--METHOD (1: min memory, 2:max memory): '
      READ (li,*) method
      WRITE (lu,*) ' INPUT: ',method
!
      CALL gretel_autop(geo,geoformat,bnd,res,resformat,nparts,nplan,
     &                  method)
!
      stop 0
      END PROGRAM