lecdon_artemis.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE lecdon_artemis
!                   *************************
!
     &(file_desc,path,ncar,
     & cas_file,dico_file)
!
!***********************************************************************
! ARTEMIS   V8P2
!***********************************************************************
!
!brief    READS THE STEERING FILE THROUGH A DAMOCLES CALL.
!
!history  J-M HERVOUET (LNH)
!+        17/08/1994
!+
!+   LINKED TO BIEF 5.0
!
!history  D. AELBRECHT (LNH)
!+        20/04/1999
!+        V6P0
!+
!
!history  C/PEYRARD (EDF)
!+        18/03/2014
!+        V6P1 - V7P0
!+   NEW KEY WORDS
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        18/05/2015
!+        V7P1
!+  Adding CHECK_MESH for the keyword 'CHECKING THE MESH'
!
!history Y AUDOUIN (LNHE)
!+       25/05/2015
!+       V7P0
!+       Modification to comply with the hermes module
!
!history  N.DURAND (HRW)
!+        November 2016
!+        V7P2
!+   Addition of new keywords in the steering file re: animation of the
!+   free surface. ANIMFS and ART_FILES(ARTAMP)%NAME read from file
!
!history  N.DURAND (HRW)
!+        August 2017
!+        V7P3
!+   1. Consolidation of nesting methods, meaning that CHAINTWC is now a
!+   choice of integers, and no longer a logical
!+   2. NDTWC (now NDIR) and NFTWC (now NF) are no longer read in the
!+   steering file but directly extracted from the spectrum file
!+   3. Addition of various checks before simulation can start
!+   4. Call to ARTEMIS_CONSTANTS added.
!+   5. Addition of new keywords in the steering file for nesting option 2:
!+   X_SFREF, Y_SFREF, ART_FILES(WACRES)%NAME
!
!history  N.DURAND (HRW)
!+        December 2018
!+        V8P0
!+   Addition of new keywords in the steering file for nesting option 2:
!+   ART_FILES(WACLQD)%NAME
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FILE_DESC      |<--| STORES STRINGS 'SUBMIT' OF DICTIONARY
!| NCAR           |-->| NUMBER OF LETTERS IN STRING PATH
!| PATH           |-->| FULL PATH TO CODE DICTIONARY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      USE declarations_telemac
      USE declarations_artemis
      USE interface_artemis, ex_lecdon_artemis => lecdon_artemis
!
      USE declarations_special
      IMPLICIT NONE
!
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      CHARACTER(LEN=PATH_LEN), INTENT(INOUT) :: FILE_DESC(4,maxkeyword)
      INTEGER, INTENT(IN)               :: NCAR
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: PATH
!     API
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: DICO_FILE
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: CAS_FILE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      LOGICAL LSTOP
      INTEGER I,KK
!
      CHARACTER(LEN=8) MNEMO(maxvar)
!
!-----------------------------------------------------------------------
!
! ARRAYS USED IN THE DAMOCLES CALL
!
      INTEGER ADRESS(4,maxkeyword),DIMENS(4,maxkeyword)
      DOUBLE PRECISION MOTREA(maxkeyword)
      INTEGER          MOTINT(maxkeyword)
      LOGICAL          MOTLOG(maxkeyword)
      CHARACTER(LEN=PATH_LEN) MOTCAR(maxkeyword)
      CHARACTER(LEN=72)  MOTCLE(4,maxkeyword,2)
      INTEGER          TROUVE(4,maxkeyword)
      LOGICAL DOC
      CHARACTER(LEN=PATH_LEN) :: NOM_CAS
      CHARACTER(LEN=PATH_LEN) :: NOM_DIC
      INTEGER :: ID_DICO, ID_CAS
! EXTRA DECLARATIONS FOR OUTPUT OF AMPLITUDE AND PHASE FILE
      CHARACTER(LEN=2)                  :: CLT      ! PERIOD COMPONENT NUMBER
!
! END OF DECLARATIONS FOR DAMOCLES CALL :
!
!-----------------------------------------------------------------------
!
! INITIALISES THE VARIABLES FOR DAMOCLES CALL :
!
      DO kk=1,maxkeyword
!
!       A FILENAME NOT GIVEN BY DAMOCLES WILL BE RECOGNIZED AS A WHITE SPACE
!       (IT MAY BE THAT NOT ALL COMPILERS WILL INITIALISE LIKE THAT)
!
        motcar(kk)(1:1)=' '
!
        dimens(1,kk) = 0
        dimens(2,kk) = 0
        dimens(3,kk) = 0
        dimens(4,kk) = 0
!
      ENDDO
!
!     WRITES OUT INFO
!
      doc = .false.
!
!-----------------------------------------------------------------------
!     OPENS DICTIONNARY AND STEERING FILES
!-----------------------------------------------------------------------
!
      IF(ncar.GT.0) THEN
!
        nom_dic=path(1:ncar)//'ARTDICO'
        nom_cas=path(1:ncar)//'ARTCAS'
!
      ELSE
!
        nom_dic='ARTDICO'
        nom_cas='ARTCAS'
!
      ENDIF
      IF((cas_file(1:1).NE.' ').AND.(dico_file(1:1).NE.' ')) THEN
        nom_dic=dico_file
        nom_cas=cas_file
      ENDIF
!
      CALL get_free_id(id_dico)
      OPEN(id_dico,file=nom_dic,form='FORMATTED',action='READ')
      CALL get_free_id(id_cas)
      OPEN(id_cas,file=nom_cas,form='FORMATTED',action='READ')
!
      CALL damocle( adress , dimens , maxkeyword, doc, lng    , lu ,
     &              motint , motrea , motlog , motcar  , motcle ,
     &              trouve , id_dico, id_cas , .false. , file_desc)
!
!-----------------------------------------------------------------------
!     CLOSES DICTIONNARY AND STEERING FILES
!-----------------------------------------------------------------------
!
      CLOSE(id_dico)
      CLOSE(id_cas)
!
!     DECODES 'SUBMIT' CHAINS
!
      CALL read_submit(art_files,maxlu_art,file_desc,maxkeyword)
!
!-----------------------------------------------------------------------
!
      DO i=1,maxlu_art
        IF(art_files(i)%TELNAME.EQ.'ARTGEO') THEN
          artgeo=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTCAS') THEN
          artcas=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTCLI') THEN
          artcli=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTFON') THEN
          artfon=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTRES') THEN
          artres=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTREF') THEN
          artref=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTBI1') THEN
          artbi1=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTBI2') THEN
          artbi2=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTFO1') THEN
          artfo1=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTFO2') THEN
          artfo2=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTRBI') THEN
          artrbi=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTRFO') THEN
          artrfo=i
        ELSEIF(art_files(i)%TELNAME.EQ.'ARTAMP') THEN
          artamp=i
        ELSEIF(art_files(i)%TELNAME.EQ.'WACSPE') THEN
          wacspe=i
        ELSEIF(art_files(i)%TELNAME.EQ.'WACRES') THEN
          wacres=i
        ELSEIF(art_files(i)%TELNAME.EQ.'WACLQD') THEN
          waclqd=i
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
!
!     SETTING CONSTANTS (PI AND RELATED)
!
      CALL artemis_constants
!
!-----------------------------------------------------------------------
!
! ASSIGNS THE STEERING FILE VALUES TO THE PARAMETER FORTRAN NAME :
!
!-----------------------------------------------------------------------
!
! INTEGER KEYWORDS :
!
      leoprd    = motint( adress(1, 1) )
      lisprd    = motint( adress(1, 2) )
      slvart%NITMAX    = motint( adress(1,3) )
      slvart%PRECON    = motint( adress(1,4) )
!     DISESP    = MOTINT( ADRESS(1, 5) )
!     STDGEO    = MOTINT( ADRESS(1, 6) )
!     STDRES    = MOTINT( ADRESS(1, 7) )
      slvart%SLV = motint( adress(1,8) )
      lisfon    = motint( adress(1, 9) )
      npale     = motint( adress(1,10) )
      ndale     = motint( adress(1,11) )
      ibreak    = motint( adress(1,12) )
      nitdis    = motint( adress(1,13) )
      regido    = motint( adress(1,14) )
      formfr    = motint( adress(1,15) )
      slvart%KRYLOV = motint( adress(1,16) )
      lvmac     = motint( adress(1,17) )
!     KFROT     = MOTINT( ADRESS(1,18) )
      optass    = motint( adress(1,19) )
      produc    = motint( adress(1,20) )
      mardat(1) = motint( adress(1,21) )
      mardat(2) = motint( adress(1,21) + 1 )
      mardat(3) = motint( adress(1,21) + 2 )
      martim(1) = motint( adress(1,22) )
      martim(2) = motint( adress(1,22) + 1 )
      martim(3) = motint( adress(1,22) + 2 )
      npriv     = motint( adress(1,23) )
!     PLACEHOLDER FOR NCSIZE:
!     NCSIZE    = MOTINT( ADRESS(1,24) )
!     ORIGIN COORDINATES
      i_orig    = motint( adress(1,25)   )
      j_orig    = motint( adress(1,25)+1 )
!     RAPIDLY VARYING TOPOGRAPHY
      ipentco   = motint( adress(1,26) )
!     MAX ITERATION ON TETAP
      nittp   = motint( adress(1,27) )
!     DEBUGGER KEYWORD
      debug   = motint( adress(1,28) )
!     2 PLACEHOLDERS FOR NITFS AND NADVAR:
!     NITFS   = MOTINT( ADRESS(1,29) )
!     NADVAR  = MOTINT( ADRESS(1,30) )
      chaintwc= motint( adress(1,31) )

! FOR THE MOMENT WTITES TO FILE FROM THE START OF SIMULATION
      ptinig    = 0
      ptinil    = 0
!
! REAL KEYWORDS :
!
      per       = motrea( adress(2, 1) )
      tetah     = motrea( adress(2, 2) )
      grav      = motrea( adress(2, 3) )
      slvart%ZERO = motrea( adress(2, 4) )
      slvart%EPS = motrea( adress(2, 5) )
!     HMIN      = MOTREA( ADRESS(2, 6) )
      cotini    = motrea( adress(2, 7) )
      hautin    = motrea( adress(2, 8) )
      perdeb    = motrea( adress(2, 9) )
      perfin    = motrea( adress(2,10) )
      perpas    = motrea( adress(2,11) )
      perpic    = motrea( adress(2,12) )
      gamma     = motrea( adress(2,13) )
      tetmin    = motrea( adress(2,14) )
      tetmax    = motrea( adress(2,15) )
      expos     = motrea( adress(2,16) )
!     RELAX     = MOTREA( ADRESS(2,17) )
      epsdis    = motrea( adress(2,18) )
      reldis    = motrea( adress(2,19) )
      alfabj    = motrea( adress(2,20) )
      gammas    = motrea( adress(2,21) )
      kdally    = motrea( adress(2,22) )
      gdally    = motrea( adress(2,23) )
      visco     = motrea( adress(2,24) )
      diam90    = motrea( adress(2,25) )
      diam50    = motrea( adress(2,26) )
      mvsed     = motrea( adress(2,27) )
      mveau     = motrea( adress(2,28) )
      fwcoef    = motrea( adress(2,29) )
      ricoef    = motrea( adress(2,30) )
      ffon      = motrea( adress(2,31) )
      pmin      = motrea( adress(2,32) )
      pmax      = motrea( adress(2,33) )
      depref    = motrea( adress(2,34) )
      x_phref   = motrea( adress(2,35) )
      y_phref   = motrea( adress(2,35)+1 )
      epsdir    = motrea( adress(2,36) )
      epstp     = motrea( adress(2,37) )
      reltp     = motrea( adress(2,38) )
      tpstwc    = motrea( adress(2,39) )
!     PLACEHOLDERS FOR TINIFS AND DTFS:
!     TINIFS    = MOTREA( ADRESS(2,40) )
!     DTFS      = MOTREA( ADRESS(2,41) )
!     FOR THE F SPECTRUM REFERENCE POINT
      x_sfref   = motrea( adress(2,42) )
      y_sfref   = motrea( adress(2,42)+1 )
!
! LOGICAL KEYWORDS :
!
      listin    = motlog( adress(3, 1) )
      infogr    = motlog( adress(3, 2) )
      balaye    = motlog( adress(3, 3) )
      alemon    = motlog( adress(3, 4) )
      alemul    = motlog( adress(3, 5) )
      deferl    = motlog( adress(3, 6) )
      frotte    = motlog( adress(3, 7) )
      entfw     = motlog( adress(3, 8) )
      entreg    = motlog( adress(3, 9) )
      entrug    = motlog( adress(3, 10) )
      lishou    = motlog( adress(3, 11) )
      valid     = motlog( adress(3, 12) )
      courant   = motlog( adress(3, 13) )
      langauto  = motlog( adress(3, 14) )
      lphaseauto= motlog( adress(3, 15) )
!     SPHERICAL EQUATIONS, HARD-CODED
      spheri    = .false.
      check_mesh= motlog( adress(3, 16) )
      animfs    = motlog( adress(3, 17) )
      partel_concat = motlog(adress(3,18))
      IF(ncsize.LE.1) partel_concat=.false.!
! STRING KEYWORDS : SOME ARE USED BY THE LAUNCHING
!                   PROCEDURE
!
      titcas    = motcar( adress(4, 1) )(1:72)
      vardes    = motcar( adress(4, 2) )(1:72)
      CALL majus(vardes)
      varimp    = motcar( adress(4, 3) )(1:72)
      CALL majus(varimp)
!     FROM 4 TO 5: READ AND USED BY PRECOS
      art_files(artgeo)%NAME=motcar( adress(4,6) )
      IF(art_files(artgeo)%NAME(1:1).EQ.' ') THEN
        WRITE(lu,*) 'THE FOLLOWING KEYWORD IS MANDATORY:'
        WRITE(lu,*) 'GEOMETRY FILE (FICHIER DE GEOMETRIE)'
        CALL plante(1)
        stop
      ENDIF
!     PLACEHOLDERS:
!     NOMFOR    = MOTCAR( ADRESS(4, 7) )
!     NOMCAS    = MOTCAR( ADRESS(4, 8) )
!     ART_FILES(ARTCAS)%NAME=MOTCAR( ADRESS(4,8) )
      art_files(artcli)%NAME=motcar( adress(4,9) )
      IF(art_files(artcli)%NAME(1:1).EQ.' ') THEN
        WRITE(lu,*) 'THE FOLLOWING KEYWORD IS MANDATORY:'
        WRITE(lu,*) 'BOUNDARY CONDITIONS FILE '//
     &              '(FICHIER DES CONDITIONS AUX LIMITES)'
        CALL plante(1)
        stop
      ENDIF
!      WRITE(*,*) 'IN LECDON ',ART_FILES(ARTGEO)%NAME
      art_files(artres)%NAME=motcar( adress(4,10) )
      IF(art_files(artres)%NAME(1:1).EQ.' ') THEN
        WRITE(lu,*) 'THE FOLLOWING KEYWORD IS MANDATORY:'
        WRITE(lu,*) 'RESULTS FILE (FICHIER DES RESULTATS)'
        CALL plante(1)
        stop
      ENDIF
!     FROM 11 TO 14 : READ AND USED BY PRECOS
      art_files(artfon)%NAME=motcar( adress(4,15) )
      art_files(artbi1)%NAME=motcar( adress(4,16) )
      art_files(artbi2)%NAME=motcar( adress(4,17) )
      art_files(artfo1)%NAME=motcar( adress(4,18) )
      art_files(artfo2)%NAME=motcar( adress(4,19) )
      art_files(artrbi)%NAME=motcar( adress(4,20) )
      art_files(artrfo)%NAME=motcar( adress(4,21) )
!     FROM 22 TO 23 : READ AND USED BY PRECOS OR NOT USED
      cdtini    = motcar( adress(4,24) )(1:72)
      CALL majus(cdtini)
!     BINGEO    = MOTCAR( ADRESS(4,25) )(1:3)
!     CALL MAJUS(BINGEO)
!     BINRES    = MOTCAR( ADRESS(4,26) )(1:3)
!     CALL MAJUS(BINRES)
!     ACCOUNTNUMBER    = MOTCAR( ADRESS(4, 27) )
      art_files(artref)%NAME=motcar( adress(4,28) )
!     GEOMETRY FILE FORMAT
      art_files(artgeo)%FMT = motcar( adress(4,29) )(1:8)
      CALL majus(art_files(artgeo)%FMT)
!     RESULTS FILE FORMAT
      art_files(artres)%FMT = motcar( adress(4,30) )(1:8)
      CALL majus(art_files(artres)%FMT)
!     REFERENCE FILE FORMAT
      art_files(artref)%FMT = motcar( adress(4,31) )(1:8)
      CALL majus(art_files(artref)%FMT)
!     BINARY FILE 1 FORMAT
      art_files(artbi1)%FMT = motcar( adress(4,32) )(1:8)
      CALL majus(art_files(artbi1)%FMT)
!     BINARY FILE 2 FORMAT
      art_files(artbi2)%FMT = motcar( adress(4,33) )(1:8)
      CALL majus(art_files(artbi2)%FMT)
      IF(animfs) THEN
!     PHASE AND AMPLITUDE FILENAME
        art_files(artamp)%NAME=motcar( adress(4,34) )
!     PHASE AND AMPLITUDE FILE FORMAT
        art_files(artamp)%FMT = motcar( adress(4,35) )(1:8)
        CALL majus(art_files(artamp)%FMT)
      ENDIF
      IF(chaintwc.GE.1) THEN
!     TOMAWAC OUTER SPECTRAL FILENAME
        art_files(wacspe)%NAME=motcar( adress(4,36) )
!     TOMAWAC OUTER SPECTRAL FILE FORMAT
        art_files(wacspe)%FMT = motcar( adress(4,37) )(1:8)
        CALL majus(art_files(wacspe)%FMT)
      ENDIF
!     PLACEHOLDERS:
!     NOMFST        = MOTCAR( ADRESS(4, 38) )
!     FREE SURFACE FILE FORMAT = MOTCAR( ADRESS(4, 39) )
!     NAME\_ADVAR   = MOTCAR( ADRESS(4, 40) )
      IF(chaintwc.EQ.2) THEN
!     TOMAWAC OUTER RESULT FILENAME
        art_files(wacres)%NAME=motcar( adress(4,41) )
!     TOMAWAC OUTER RESULT FILE FORMAT
        art_files(wacres)%FMT = motcar( adress(4,42) )(1:8)
        CALL majus(art_files(wacres)%FMT)
!     LIQUID BOUNDARY FILENAME, DERIVED FROM TOMAWAC OUTER RESULT FILE
        art_files(waclqd)%NAME=motcar( adress(4,43) )
!     LIQUID BOUNDARY FILE FORMAT
        art_files(waclqd)%FMT = motcar( adress(4,44) )(1:8)
        CALL majus(art_files(waclqd)%FMT)
      ENDIF
!     PLACEHOLDERS:
!     DICTIONARY    = MOTCAR( ADRESS(4, 4) )
!     PARTITIONING TOOL = MOTCAR( ADRESS(4, 5) )
!
      IF(listin) THEN
        WRITE(lu,102)
      ENDIF
102   FORMAT(1x,/,19x, '********************************************',/,
     &            19x, '*               LECDON:                    *',/,
     &            19x, '*        AFTER CALLING DAMOCLES            *',/,
     &            19x, '*        CHECKING OF DATA  READ            *',/,
     &            19x, '*         IN THE STEERING FILE             *',/,
     &            19x, '********************************************',/)
!
!-----------------------------------------------------------------------
!  NAME OF THE VARIABLES FOR THE RESULTS AND GEOMETRY FILES :
!-----------------------------------------------------------------------
!
      mnemo = '        '
      CALL nomvar_artemis(texte,textpr,mnemo)
!
! LOGICAL ARRAY FOR OUTPUT
!
      sorleo = .false.
      sorimp = .false.
      CALL sortie(vardes , mnemo , 100 , sorleo )
      CALL sortie(varimp , mnemo , 100 , sorimp )
!
!-----------------------------------------------------------------------
!
! IN THE CASE OF A PERIOD SWEEPING, THE FIRST PERIOD TO BE COMPUTED
! IS PERDEB
!
      IF (balaye) per = perdeb
!
!-----------------------------------------------------------------------
!
! IF NOTHING IS TO BE WRITTEN TO THE LISTING FILE, THERE SHOULD BE
! NO INFORMATION WRITTEN ABOUT THE SOLVEUR
!
      IF (.NOT.listin) infogr = .false.
!
!-----------------------------------------------------------------------
!
! FOR RANDOM SEA COMPUTATIONS, INHIBITS THE GRAPHICAL AND LISTING
! OUTPUT OF THE PHASES, SPEEDS, FREE SURFACE ELEVATION, AND
! POTENTIAL BECAUSE THEY DO NOT MEAN ANYTHING.
!
! BUT WRITES OUT AN AVERAGE WAVE NUMBER, AN AVERAGE PHASE CELERITY
! AND AN AVERAGE GROUP VELOCITY, COMPUTED FROM THE MEAN WAVE
! PERIOD T01. ALSO WRITES OUT AN AVERAGE INCIDENCE
!
! MOREOVER, CHECKS THAT THE NUMBER OF DISCRETISED PERIODS AND
! DIRECTIONS IS AT LEAST 1
!
      IF(alemon .OR. alemul) THEN
!
!       2 : PHASE
!
        sorleo( 2) = .false.
        sorimp( 2) = .false.
!
!       3 AND 4 : U0 AND V0
!
        sorleo( 3) = .false.
        sorimp( 3) = .false.
        sorleo( 4) = .false.
        sorimp( 4) = .false.
!
!       5 : FREE SURFACE
!
        sorleo( 5) = .false.
        sorimp( 5) = .false.
!
!       11 AND 12 : REAL AND IMAGINARY PARTS OF THE POTENTIAL
!
        sorleo(11) = .false.
        sorimp(11) = .false.
        sorleo(12) = .false.
        sorimp(12) = .false.
!
        IF (npale.LE.0) npale = 1
        IF (ndale.LE.0) ndale = 1
!
      ELSE
!
!       17, 18 AND 19 : T01, T02 AND TM
!
        sorleo(17) = .false.
        sorimp(17) = .false.
        sorleo(18) = .false.
        sorimp(18) = .false.
        sorleo(19) = .false.
        sorimp(19) = .false.
      ENDIF
!
! NO PERIOD SWEEPING FOR RANDOM SEAS
!
      IF (alemon .OR. alemul) balaye = .false.
!
! NDALE=1 FOR MONO-DIRECTIONAL SEAS
!
      IF (alemon .AND. .NOT.alemul) ndale = 1
!
! NPALE=1 FOR REGULAR SEAS
!
      IF (.NOT.alemon .AND. .NOT.alemul .AND. .NOT.balaye) npale = 1
!
! IF AMPLITUDE AND PHASE FILE IS REQUIRED, NAME OF VARIABLES IN TEXTANIM
! AND LOGICAL ARRAY FOR OUTPUT IN SORNIM
!
      IF (animfs) THEN
        DO i=1,ndale
          WRITE(clt,205) i
205       FORMAT(i2.2)
          textanim(2*i-1) = 'WAVE HEIGHT D'//clt//'              (M)'
          textanim(2*i)   = 'WAVE PHASE D'//clt//'             (RAD)'
!
          sornim(2*i-1) = .true.
          sornim(2*i) = .true.
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
!  DOES NOT WRITE OUT 'PRIVE' VARIABLES IF THEY'VE NOT BEEN ALLOCATED
!
      lstop = .false.
      DO i=1,4
        IF ((sorleo(12+i).OR.sorimp(12+i)).AND.(npriv.LT.i)) THEN
          WRITE(lu,17) i,npriv
          lstop = .true.
        ENDIF
      ENDDO
      IF (lstop) THEN
        CALL plante(1)
        stop
      ENDIF
 17   FORMAT(1x,'PRIVATE ARRAY ',1i1,' CANNOT BE USED '
     &      ,1x,'BECAUSE IT WAS NOT ALLOCATED.',/
     &      ,1x,'CHECK ''NPRIV''  (AT THE TIME BEING ',1i1,' ).',/)
!
!-----------------------------------------------------------------------
!
!  WRITES OUT THE TITLE
!
      IF(listin) THEN
        WRITE(lu,3001) titcas
3001    FORMAT(/1x,'NAME OF THE STUDY :',1x,a72,/)
      ENDIF
!
!-----------------------------------------------------------------------
!
! NO TIDAL FLATS ==> MSK = .FALSE.
!
      msk = .false.
!
!-----------------------------------------------------------------------
!
! COMPULSORY CHOICE OF KEYWORD FOR THE DIRECT SOLVEUR
!
      IF( (slvart%SLV.EQ.8).AND.(optass.NE.3) ) THEN
        WRITE(lu,3003)
!
3003    FORMAT(1x,'WITH DIRECT SYSTEM SOLVER, EDGE-BASED STORAGE',/,1x,
     &             'IS MANDATORY',///)
        CALL plante(1)
        stop
      ENDIF
!
! USE OF DIRECT SOLVEUR IS NOT POSSIBLE WITH PARALLELISM
!
      IF(ncsize.GT.1) THEN
        IF(slvart%SLV.EQ.8) THEN
          WRITE(lu,3005)
3005      FORMAT(1x,'WITH PARALLELISM,',/,1x,
     &            'NO DIRECT SYSTEM SOLVER 8 USE 9',///)
          CALL plante(1)
          stop
        ENDIF
      ELSE
        IF(slvart%SLV.EQ.9) THEN
          WRITE(lu,3006)
3006      FORMAT(1x,'WITH NO PARALLELISM,',/,1x,
     &            'NO DIRECT SYSTEM SOLVER 9 USE 8',///)
          CALL plante(1)
          stop
        ENDIF
      ENDIF
!
!-----------------------------------------------------------------------
!
!  TESTS PMIN%MAX and TETMIN%MAX
!
!     CHECKS IF TETMIN<TETMAX
      IF(tetmax.LE.tetmin) THEN
        WRITE(lu,*) 'LECDON_ARTEMIS:'
        WRITE(lu,*) 'MAXIMUM ANGLE OF PROPAGATION < MINIMUM ANGLE'
        WRITE(lu,*) '=> CHANGE MAXIMUM OR MINIMUM ANGLE OF PROPAGATION'
        WRITE(lu,*) '   IN THE USER PARAMETER FILE.                   '
        CALL plante(1)
        stop
      ENDIF
!     CHECKS IF PMIN<PMAX
      IF(pmax.LE.pmin) THEN
        WRITE(lu,*) 'LECDON_ARTEMIS:'
        WRITE(lu,*) 'MAXIMUM SPECTRAL PERIOD < MINIMUM SPECTRAL PERIOD'
        WRITE(lu,*) '=> CHANGE MAXIMUM OR MINIMUM SPECTRAL PERIOD     '
        WRITE(lu,*) '   IN THE USER PARAMETER FILE.                   '
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
!  TESTS SUPPLY OF REQUIRED INPUT FOR NESTING
!
      IF(chaintwc.GE.1) THEN
        IF (art_files(wacspe)%NAME .EQ. '') THEN
          WRITE(lu,3007)
3007      FORMAT(/1x,'PLEASE SUPPLY TOMAWAC OUTER SPECTRAL FILE',/)
          stop
        ENDIF
      ENDIF
      IF(chaintwc.EQ.2) THEN
        IF (art_files(wacres)%NAME .EQ. '') THEN
          WRITE(lu,3009)
3009      FORMAT(/1x,'PLEASE SUPPLY TOMAWAC OUTER RESULT FILE',/)
          stop
        ENDIF
        IF (art_files(waclqd)%NAME .EQ. '') THEN
          WRITE(lu,3011)
3011      FORMAT(/1x,'PLEASE SUPPLY TOMAWAC LIQUID BOUNDARY FILE',/)
          stop
        ENDIF
        IF (x_sfref .LT. -99999.) THEN
          WRITE(lu,3013)
3013      FORMAT(/1x,'PLEASE SUPPLY X FOR THE REFERENCE F SPECTRUM',/)
          stop
        ENDIF
        IF (y_sfref .LT. -99999.) THEN
          WRITE(lu,3015)
3015      FORMAT(/1x,'PLEASE SUPPLY Y FOR THE REFERENCE F SPECTRUM',/)
          stop
        ENDIF
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END