lecdon_sisyphe.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE lecdon_sisyphe
!                   *************************
!
     &(motcar,file_desc,path,ncar,code,cas_file,dico_file)
!
!***********************************************************************
! SISYPHE   V7P2
!***********************************************************************
!
!brief    READS THE STEERING FILE BY CALL TO DAMOCLES.
!
!history  C.VILLARET + JMH (EDF-LNHE)
!+        02/05/2012
!+        V6P2
!+  File for liquid boundaries added
!+
!
!history  JWI
!+        31/05/2012
!+        V6P2
!+  added one increment to include wave orbital velocities
!+  (SORLEO(I+28+(NOMBLAY+4)*NSICLA+NOMBLAY).OR.
!
!history  PAT (LNHE)
!+        18/06/2012
!+        V6P2
!+   updated version with HRW's development
!
!history  Pablo Tassi PAT (EDF-LNHE)
!+        12/02/2013
!+        V6P3
!+ Preparing for the use of a higher NSICLM value
!+ (by Rebekka Kopmann)
!
!history  Pablo Tassi PAT (EDF-LNHE)
!+        12/02/2013
!+        V6P3
!+ Settling lag: determines choice between Rouse and Miles concentration profile
!+ (by Michiel Knaapen HRW)
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        18/05/2015
!+        V7P1
!+  Adding CHECK_MESH for the keyword 'CHECKING THE MESH'
!
!history  R. KOPMANN (BAW)
!+        13/07/2016
!+        V7P2
!+        Integrating liquid boundary file for QS
!
!history  S.E. BOURBAN (HRW)
!+        20/06/2016
!+        V7P2
!+   Now taking into account names of differentiators given by user.
!
!history  R.KOPMANN (BAW)
!+        07/12/2017
!+        V7P2
!+   VSPRES results file for CVSM handled by CAS-file.
!
!history  B.GLANDER (BAW)
!+        06/12/2018
!+        V7P2
!+  NEW VARIABLE: ZRL  REFERENCE LEVEL FOR NESTOR, CHANGE OF SORLEO, SORIMP
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FILE_DESC      |<--| STORES STRINGS 'SUBMIT' OF DICTIONARY
!| MOTCAR         |<--| VALUES OF KEY-WORDS OF TYPE CHARACTER
!| NCAR           |-->| NUMBER OF LETTERS IN STRING PATH
!| PATH           |-->| FULL PATH TO CODE DICTIONARY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      USE declarations_telemac
      USE declarations_sisyphe
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)               :: NCAR
      CHARACTER(LEN=24), INTENT(IN)     :: CODE
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: PATH
      CHARACTER(LEN=PATH_LEN), INTENT(INOUT)      :: MOTCAR(maxkeyword)
      CHARACTER(LEN=PATH_LEN), INTENT(INOUT) :: FILE_DESC(4,maxkeyword)
!     API
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: CAS_FILE
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: DICO_FILE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER            :: I,K,ERR
      INTEGER            :: MOTINT(maxkeyword)
      INTEGER            :: TROUVE(4,maxkeyword)
      INTEGER            :: ADRESS(4,maxkeyword)
      INTEGER            :: DIMENS(4,maxkeyword)
      DOUBLE PRECISION   :: SUMAVAI
      DOUBLE PRECISION   :: MOTREA(maxkeyword)
      LOGICAL            :: DOC,EFFPEN
      LOGICAL            :: MOTLOG(maxkeyword)
      CHARACTER(LEN=PATH_LEN) :: NOM_CAS
      CHARACTER(LEN=PATH_LEN) :: NOM_DIC
      CHARACTER(LEN=72)  :: MOTCLE(4,maxkeyword,2)

      CHARACTER(LEN=PATH_LEN) TEMPVAR
      INTEGER :: ID_DICO, ID_CAS
!
!-----------------------------------------------------------------------
!
      CHARACTER(LEN=2) CHAR2
!
!-----------------------------------------------------------------------
!
! INITIALISES THE VARIABLES FOR DAMOCLES CALL :
!
      DO k = 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(k)(1:1)=' '
!
        dimens(1,k) = 0
        dimens(2,k) = 0
        dimens(3,k) = 0
        dimens(4,k) = 0
      ENDDO
!
!     WRITES OUT INFO
      doc = .false.
!
!-----------------------------------------------------------------------
!     OPENS DICTIONNARY AND STEERING FILES
!-----------------------------------------------------------------------
!
      IF(ncar.GT.0) THEN
!
        nom_dic=path(1:ncar)//'SISDICO'
        nom_cas=path(1:ncar)//'SISCAS'
!
      ELSE
!
        nom_dic='SISDICO'
        nom_cas='SISCAS'
!
      ENDIF
      IF((cas_file(1:1).NE.' ').AND.(dico_file(1:1).NE.' ')) THEN
        WRITE(lu,*) 'FIXED DICO AND STEERING FILE PRESENT'
        nom_dic=dico_file
        nom_cas=cas_file
        WRITE(lu,*) 'NOM_DIC',nom_dic
        WRITE(lu,*) 'NOM_CAS',nom_cas
      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')
!
!-----------------------------------------------------------------------
!     CALLS DAMOCLES
!-----------------------------------------------------------------------
!
      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(sis_files,maxlu_sis,file_desc,maxkeyword)
!
!-----------------------------------------------------------------------
!
!     RETRIEVES FILES NUMBERS IN TELEMAC-2D FORTRAN PARAMETERS
!     AT THIS LEVEL LOGICAL UNITS ARE EQUAL TO THE FILE NUMBER
!
      DO i=1,maxlu_sis
        IF(sis_files(i)%TELNAME.EQ.'SISHYD') THEN
          sishyd=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISGEO') THEN
          sisgeo=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISCLI') THEN
          siscli=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISPRE') THEN
          sispre=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISRES') THEN
          sisres=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISREF') THEN
          sisref=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISCOU') THEN
          siscou=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISFON') THEN
          sisfon=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISSEC') THEN
          sissec=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISSEO') THEN
          sisseo=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISLIQ') THEN
          sisliq=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SISFLX') THEN
          sisflx=i
!       === NESTOR FILES
        ELSEIF(sis_files(i)%TELNAME.EQ.'SINACT') THEN
          sinact=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SINPOL') THEN
          sinpol=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SINREF') THEN
          sinref=i
        ELSEIF(sis_files(i)%TELNAME.EQ.'SINRST') THEN
          sinrst=i
!       ===
        ELSEIF(sis_files(i)%TELNAME.EQ.'VSPRES') THEN
          vspres=i
!
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
!
!     ALLOCATING ARRAYS DEPENDING ON MAXFRO
!
!-----------------------------------------------------------------------
!
!     MAXIMUM NUMBER OF BOUNDARIES
      maxfro   = motint( adress(1,58) )
!
      ALLOCATE(soldis(maxfro))
      ALLOCATE(okcgl(maxfro))
      ALLOCATE(okqgl(maxfro))
      ALLOCATE(cbor_classe(nsiclm*maxfro))
!
      DO k=1,maxfro
        okcgl(k)=.true.
        okqgl(k)=.true.
      ENDDO
!
!-----------------------------------------------------------------------
!
!     ASSIGNS THE STEERING FILE VALUES TO THE PARAMETER FORTRAN NAME
!
!-----------------------------------------------------------------------
!
!     OPTION OF MATRIX ASSEMBLY IS HARD-CODED
!
      optass = 1
      produc = 1
!
!     DISCRETISATIONS
!
      ielmt     = 11 ! SEDIMENTOLOGICAL VARIABLES
      ielmh_sis = 11 ! VARIABLES ASSOCIATED WITH WATER DEPTH
      ielmu_sis = 11 ! VARIABLES ASSOCIATED WITH VELOCITIES
!
!     FOR NOW PRINTOUTS START AT ZERO
!
      ptinig = 0
      ptinil = 0
!
!     NON-EQUILIBIRUM BEDLOAD
!
      loadmeth = 0
!
!     ICM           = MOTINT( ADRESS(1,  1) )
      icf           = motint( adress(1,  2) )
      npas          = motint( adress(1,  3) )
      nmaree        = motint( adress(1,  4) )
!     N1            = MOTINT( ADRESS(1,  5) )

      leopr         = motint( adress(1,  6) )
      lispr         = motint( adress(1,  7) )
      optban        = motint( adress(1, 11) )
      lvmac         = motint( adress(1, 12) )
      nsous         = motint( adress(1, 14) )
!
      mardat(1)     = motint( adress(1, 15) )
      mardat(2)     = motint( adress(1, 15) + 1 )
      mardat(3)     = motint( adress(1, 15) + 2 )
      martim(1)     = motint( adress(1, 16) )
      martim(2)     = motint( adress(1, 16) + 1 )
      martim(3)     = motint( adress(1, 16) + 2 )
!
      slvsed%SLV    = motint( adress(1, 17) )
      slvsed%KRYLOV = motint( adress(1, 18) )
      slvsed%PRECON = motint( adress(1, 19) )
      slvsed%NITMAX = motint( adress(1, 20) )
      choix         = motint( adress(1, 21) )
      dirflu        = motint( adress(1, 22) )
      npriv         = motint( adress(1, 23) )
      nadvar    = motint( adress(1,30) )
!     NUMBER OF DIRECTIONS FOR DIFFERENTIATED VARIABLES
      ad_numofdir  = motint( adress(1,59) )
!
!     NCSIZE        = MOTINT( ADRESS(1, 24) )
!     NUMBER OF PROCESSORS (ALREADY GIVEN IN INIT_FILES2;
!     MUST BE THE SAME, BUT WHEN USING COUPLED MODELS IT CAN
!     WRONGLY BE DIFFERENT)
      IF(ncsize.NE.motint(adress(1,24))) THEN
        WRITE(lu,*) 'DIFFERENT NUMBER OF PARALLEL PROCESSORS:'
        WRITE(lu,*) 'DECLARED BEFORE (CASE OF COUPLING ?):',ncsize
        WRITE(lu,*) 'SISYPHE :',motint(adress(1,24))
        WRITE(lu,*) 'VALUE ',ncsize,' IS KEPT'
      ENDIF
      resol         = motint( adress(1, 25) )
      slvtra%SLV    = motint( adress(1, 26) )
      slvtra%KRYLOV = motint( adress(1, 27) )
      slvtra%PRECON = motint( adress(1, 28) )
      slvtra%NITMAX = motint( adress(1, 29) )
      optdif        = motint( adress(1, 31) )
      optsup        = motint( adress(1, 32) )
      produc        = motint( adress(1, 33) )
      optass        = motint( adress(1, 34) )
      opdtra        = motint( adress(1, 35) )
      kfrot         = motint( adress(1, 37) )
      ncondis       = motint( adress(1, 38) )
      slopeff       = motint( adress(1, 39) )
      devia         = motint( adress(1, 40) )
      nomblay       = motint( adress(1,251) )
      nsicla        = motint( adress(1,252) )
      hidfac        = motint( adress(1,253) )
      icq           = motint( adress(1, 41) )
!     CONTROL SECTIONS
      ncp=dimens(1,42)
      ALLOCATE(ctrlsc(ncp),stat=err)
      CALL check_allocate(err, 'CTRLSC')
      IF(ncp.GE.1) THEN
        DO k=1,ncp
          ctrlsc(k) = motint( adress(1,42) + k-1 )
        ENDDO
      ENDIF
!     COORDINATES OF THE ORIGIN
      i_orig = motint( adress(1,43)   )
      j_orig = motint( adress(1,43)+1 )
      debug  = motint( adress(1,44)   )
!
      icr  =   motint(adress(1,46)   )
!
      iks  =   motint(adress(1,47)   )
!     CVSM MODEL
      pro_max_max =   max(nsicla*4+4,motint(adress(1,49) ))
      cvsmpperiod =   motint(adress(1,50)   )
      alt_model   =   motint(adress(1,52)   )
      vsmtype     =   motint(adress(1,53)   )
!
!     MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES
!
      maxadv = motint(adress(1,54))
!
!     SCHEME OPTION FOR ADVECTION
!
      optadv=1
      IF(resol.EQ.1) THEN
!       CHARACTERISTICS
!       OPTADV=OPTCHA (IN TELEMAC-2D, NOT IN SISYPHE YET)
      ELSEIF(resol.EQ.2) THEN
!       SUPG
        optadv=optsup
      ELSEIF(resol.EQ.5) THEN
!       PSI SCHEME
!       OPTADV_VI=OPTPSI (IN TELEMAC-2D, NOT IN SISYPHE YET)
      ENDIF
!     SCHEME OPTION FOR ADVECTION HAS PRIORITY WHEN PRESENT
      IF(trouve(1,55).EQ.2) THEN
        optadv = motint(adress(1,55))
      ENDIF
      nco_dist = motint( adress(1,56) )
      nsp_dist = motint( adress(1,57) )
!
! ############### !
! REAL KEYWORDS   !
! ############### !
!
      ! NON-EQUILIBRIRUM BEDLOAD
      ! ------------------------
      ls0         = 1.d0
!
      rc          = motrea( adress(2,  1) )
      xmve        = motrea( adress(2,  2) )
      xmvs        = motrea( adress(2,  3) )
      DO k=1,nsicla
        fdm(k)   = motrea( adress(2,  4) + k-1 )
      ENDDO
!     IF OLD NAME OF KEYWORD 28 HAS BEEN FOUND
      IF(trouve(2,28).EQ.2) THEN
        DO k=1,nsicla
          fdm(k) = motrea( adress(2,28) + k-1 )
        ENDDO
      ENDIF
      xkv         = motrea( adress(2,  5) )
!     SHIELDS NUMBERS
      DO k=1,dimens(2,6)
        ac(k)   = motrea( adress(2, 6) + k-1 )
      ENDDO
!     IF ALL SHIELDS NUMBERS ARE NOT GIVEN, THE LATEST
!     ONE IS TAKEN FOR THE FOLLOWING
!     FOR EXAMPLE IF ONLY ONE IS GIVEN, ALL WILL HAVE
!     THE SAME VALUE
      IF(dimens(2,6).LT.nsicla) THEN
        DO k=dimens(2,6)+1,nsicla
          ac(k) = motrea( adress(2, 6)+max(dimens(2,6),1)-1 )
        ENDDO
      ENDIF
      sfon        = motrea( adress(2,  7) )
      grav        = motrea( adress(2,  8) )
      zero        = motrea( adress(2,  9) )
      slvsed%ZERO = zero
      vce         = motrea( adress(2, 10) )
      hmin        = motrea( adress(2, 11) )
      delt        = motrea( adress(2, 12) )
      tprec       = motrea( adress(2, 13) )
      pmaree      = motrea( adress(2, 14) )
      teta        = motrea( adress(2, 15) )
      beta        = motrea( adress(2, 16) )
      slvsed%EPS  = motrea( adress(2, 17) )
      teta_susp   = motrea( adress(2, 18) )
      xkx         = motrea( adress(2, 19) )
      xky         = motrea( adress(2, 20) )
      slvtra%EPS  = motrea( adress(2, 21) )
!     SETTLING VELOCITIES (SAME TREATMENT THAN SHIELDS NUMBERS)
      DO k=1,dimens(2,22)
        xwc(k)   = motrea( adress(2, 22) + k-1 )
      ENDDO
      IF(dimens(2,22).LT.nsicla) THEN
        DO k=dimens(2,22)+1,nsicla
          xwc(k) = motrea( adress(2, 22)+dimens(2,22)-1 )
        ENDDO
      ENDIF
      crit_cfd    = motrea( adress(2, 23) )
      kspratio    = motrea( adress(2, 24) )
      phised      = motrea( adress(2, 25) )
      beta2       = motrea( adress(2, 26) )
      bijk        = motrea( adress(2, 27) )
!
!     INITIAL CONCENTRATIONS
!     ++++++++++++++++++++++
!
      DO k=1,nsiclm
!       DEFAULT VALUE
        cs0(k) = 0.d0
      ENDDO
      DO k=1,nsicla
        cs0(k)=motrea( adress(2,30) + k-1 )
      ENDDO
      DO k=1,10*maxfro
        cbor_classe(k)=0.d0
      ENDDO
      IF(dimens(2,31).GT.0) THEN
        DO k=1,dimens(2,31)
          cbor_classe(k)=motrea( adress(2,31) + k-1 )
        ENDDO
      ENDIF
!
!     COHESIVE SEDIMENT
!     +++++++++++++++++
!
      ncouch_tass = motint( adress(1,45)   )
!
!     DEFAULT VALUES, PREVIOUSLY IN THE DICTIONARY
!     EVEN IF THERE ARE NOT 10 LAYERS
      conc_vase(1)  =  50.d0
      conc_vase(2)  = 100.d0
      conc_vase(3)  = 150.d0
      conc_vase(4)  = 200.d0
      conc_vase(5)  = 250.d0
      conc_vase(6)  = 300.d0
      conc_vase(7)  = 350.d0
      conc_vase(8)  = 400.d0
      conc_vase(9)  = 450.d0
      conc_vase(10) = 500.d0
!
      IF(dimens(2,32).GT.0) THEN
        DO k=1,dimens(2,32)
          conc_vase(k)=motrea( adress(2,32) + k-1 )
        ENDDO
      ENDIF
!
!     OBSOLETE KEY WORD : CSF_VASE = MOTREA( ADRESS(2, 29) )
!
!      CSF_VASE = CONC_VASE(1)/XMVS
!
      IF(dimens(2,34).GT.0) THEN
        DO k=1,dimens(2,34)
          toce_vase(k)=motrea( adress(2,34) + k-1 )
        ENDDO
      ENDIF
!
!     KRONE AND PARTHENIADES EROSION AND DEPOSITION LAW
!
!     OBSOLETE KEY WORD : VITCE= MOTREA( ADRESS(2,35))
!
      vitce = sqrt(toce_vase(1)/xmve)
      vitcd= motrea( adress(2,36))
!     PARTHENIADES WITH CONVERSION TO M/S
      partheniades = motrea( adress(2,37))/xmvs
!
!     CONSOLIDATION MODEL
!
      tass = motlog(adress(3,23))
      itass  =   motint(adress(1,48)   )
!
!     MULTILAYER MODEL (WALTHER, 2008)
!     ITASS = 1
!
!     DEFAULT VALUES, PREVIOUSLY IN THE DICTIONARY
!     EVEN IF THERE ARE NOT 10LAYERS
      trans_mass(1)  = 5.d-5
      trans_mass(2)  = 4.5d-5
      trans_mass(3)  = 4.d-5
      trans_mass(4)  = 3.5d-5
      trans_mass(5)  = 3.d-5
      trans_mass(6)  = 2.5d-5
      trans_mass(7)  = 2.d-5
      trans_mass(8)  = 1.5d-5
      trans_mass(9)  = 1.d-5
      trans_mass(10) = 0.d0
!
      IF(dimens(2,33).GT.0) THEN
        DO k=1,dimens(2,33)
          trans_mass(k)=motrea( adress(2,33) + k-1 )
        ENDDO
      ENDIF
!
! V6P1
! THIEBOT MULTI LAYER MODEL
! ITASS=2
!
      conc_gel=motrea( adress(2,38))
      coef_n= motrea( adress(2,39))
      conc_max=motrea( adress(2,50))
!
!     PRESCRIBED SOLID DISCHARGES
!
      nsoldis=dimens(2,51)
      IF(nsoldis.GT.0) THEN
        DO k=1,nsoldis
          soldis(k)=motrea(adress(2,51)+k-1)
        ENDDO
      ENDIF
!
!     MINIMUM DEPTH FOR BEDLOAD
!
      hmin_bedload=motrea(adress(2,52))
!
!     HIDING EXPOSURE MULTI GRAIN MODEL
!
      DO k=1,nsicla
        hidi(k)  = motrea( adress(2,253) + k-1 )
        IF (trouve(2,255).EQ.1) THEN
          fd90(k)= fdm(k)
        ELSE
          fd90(k)= motrea( adress(2,255) + k-1 )
        ENDIF
        ava0(k)  = motrea( adress(2,258) + k-1 )
      ENDDO
      elay0       = motrea( adress(2,259) )
!
! UM: MPM-Factor
      mpm         = motrea( adress(2,260) )
! UM: ALPHA-Factor
      alpha       = motrea( adress(2,261) )
! UM: MOFAC-Factor
      mofac       = motrea( adress(2,262) )

      ! ################## !
      ! LOGICAL KEYWORDS !
      ! ################## !
! INDEX 99 IS ALREADY USED FOR KEYWORD 'LIST OF FILES'
! INDEX 54 IS ALREADY USED FOR KEYWORD 'DESCRIPTION OF LIBRARIES'
! INDEX 57 IS ALREADY USED FOR KEYWORD 'DEFAULT EXECUTABLE'
      ! SPHERICAL EQUATIONS HARD-CODED
      ! ----------------------------------
      spheri       = .false.


      ! COMPUTATION OF FALL VELOCITIES
      ! ------------------------------------------
      calwc = .false.
      ! IF TROUVE: VELOCITIES ARE USER-DEFINED
      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      IF (trouve(2, 22).EQ.2) calwc = .true.
! CV
      ! SHIELDS PARAMETER
      ! ------------------------------------------
      calac = .false.
      ! IF TROUVE
      ! ~~~~~~~~~~~~~
      IF (trouve(2, 6).EQ.2) calac = .true.


      bilma        = motlog( adress(3,  1) )
      perma        = motlog( adress(3,  2) )
      bandec       = motlog( adress(3,  3) )
      valid        = motlog( adress(3,  4) )
!     DTVAR        = MOTLOG( ADRESS(3,  5) )
      lumpi        = motlog( adress(3,  6) )
      susp         = motlog( adress(3,  7) )
      charr        = motlog( adress(3,  8) )
      houle        = motlog( adress(3, 10) )
      const_alayer = motlog( adress(3, 11) )
      lcondis      = motlog( adress(3, 12) )
      lgrafed      = motlog( adress(3, 13) )
!     USED TO CHECK SIS_FILES(SISPRE)%NAME
      debu         = motlog( adress(3, 14) )
      imp_inflow_c = motlog( adress(3, 15) )
      seccurrent   = motlog( adress(3, 16) )
      havesecfile  = motlog( adress(3, 59) )
      IF(code(1:9).EQ.'TELEMAC3D') seccurrent = .false.
      unit         = motlog( adress(3, 17) )
      vf           = motlog( adress(3,253) )
      corr_conv    = motlog( adress(3, 18) )
      DO k=1,nsicla
        sedco(k)   = .false.
      ENDDO
      IF(dimens(3,19).GT.0) THEN
        DO k=1,dimens(3,19)
          sedco(k) = motlog( adress(3,19) + k-1 )
        ENDDO
      ENDIF
      slide    = motlog( adress(3, 20) )
      dift     = motlog( adress(3, 21) )
      effpen   = motlog( adress(3, 22) )
      IF(.NOT.effpen) THEN
        slopeff=0
        devia=0
      ENDIF
!
      mixte=motlog(adress(3,24))
!     COUPLING WITH NESTOR
      nestor=motlog(adress(3,25))
!     V6P1
      kspred   =motlog(adress(3,26))
!
!     Settling lag: determines choice between Rouse and Miles concentration profile
!     SET_LAG = TRUE : Miles
!             = FALSE: Rouse
!
      set_lag  = motlog(adress(3,27) )
!     STATIONARY MODE: calculate sediment transport without updating the bed.
      stat_mode  = motlog(adress(3,28) )
!     Checking the mesh
      check_mesh = motlog(adress(3,29) )
!     NEW IMPLEMENTATION FOR CROSS-SECTION
      doflux = motlog(adress(3,61) )
      IF ( code(1:7) .NE. 'TELEMAC' ) THEN
!       SYMBOLIC LINEAR SOLVER FOR AD
        ad_symblinsolv  = motlog( adress(3,30) )
!       RESET TANGENTS ON ENTRY IN LINEAR SOLVER CG FOR AD
        ad_linsolv_resetderiv  = motlog( adress(3,31) )
!       ALLOW ADDITONAL INTERATIONS IN ITERATIVE LINEAR SOLVERS
!         SHOULD NOT BE USED IN PARALLEL MODE
        ad_linsolv_derivative_convergence  = motlog( adress(3,32) )
      ENDIF
      partel_concat = motlog(adress(3,62))
      IF(ncsize.LE.1) partel_concat=.false.
!
! ################################### !
! CHARACTER STRING KEYWORDS           !
! ################################### !
!
      titca            = motcar( adress(4, 1) )(1:72)
      sortis           = motcar( adress(4, 2) )(1:72)
      varim            = motcar( adress(4, 3) )(1:72)
      sis_files(sisgeo)%NAME=motcar( adress(4,6) )
      IF(sis_files(sisgeo)%NAME(1:1).EQ.' ') THEN
        WRITE(lu,*) 'THE FOLLOWING KEYWORD IS MANDATORY:'
        WRITE(lu,*) 'GEOMETRY FILE (FICHIER DE GEOMETRIE)'
        CALL plante(1)
        stop
      ENDIF
      sis_files(siscli)%NAME=motcar( adress(4,9) )
      IF(sis_files(siscli)%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
      sis_files(sishyd)%NAME=motcar( adress(4,29) )
      sis_files(sispre)%NAME=motcar( adress(4,11) )
      sis_files(sisres)%NAME=motcar( adress(4,12) )
      sis_files(sisfon)%NAME=motcar( adress(4,16) )
      sis_files(sisres)%FMT = motcar( adress(4,31) )(1:8)
      IF(sis_files(sisres)%NAME(1:1).EQ.' ') THEN
        WRITE(lu,*) 'THE FOLLOWING KEYWORD IS MANDATORY:'
        WRITE(lu,*) 'RESULTS FILE (FICHIER DES RESULTATS)'
        CALL plante(1)
        stop
      ENDIF
      CALL majus(sis_files(sisres)%FMT)
!     RESULT FILE FORMAT FOR PREVIOUS SEDIMENTOLOGICAL
!     COMPUTATION...
      sis_files(sispre)%FMT = motcar( adress(4,34) )(1:8)
      CALL majus(sis_files(sispre)%FMT)
!     REFERENCE FILE FORMAT
      sis_files(sisref)%FMT = motcar( adress(4,33) )(1:8)
      CALL majus(sis_files(22)%FMT)
!     HYDRODYNAMIC FILE FORMAT
      sis_files(sishyd)%FMT = motcar( adress(4,32) )(1:8)
      CALL majus(sis_files(sishyd)%FMT)
!     WAVE FILE FORMAT (COUPLING WITH TOMAWAC)
      sis_files(siscou)%FMT = motcar( adress(4,35) )(1:8)
      CALL majus(sis_files(siscou)%FMT)
      bingeosis        = motcar( adress(4,18) )(1:3)
      binhydsis        = motcar( adress(4,19) )(1:3)
      binpresis        = motcar( adress(4,20) )(1:3)
      binressis        = motcar( adress(4,21) )(1:3)
      sis_files(sisref)%NAME=motcar( adress(4,22) )
      binrefsis        = motcar( adress(4,23) )(1:3)
!     NESTOR STEERING FILE
      sis_files(sinact)%NAME = motcar( adress(4,27) )  !  Nestor
      sis_files(sinpol)%NAME = motcar( adress(4,40) )  !  Nestor
      sis_files(sinref)%NAME = motcar( adress(4,41) )  !  Nestor
      sis_files(sinrst)%NAME = motcar( adress(4,64) )  !  Nestor
!     ******           = MOTCAR( ADRESS(4,28) )
!     WAVE FILE
      sis_files(siscou)%NAME=motcar( adress(4,30) )
!     SECTIONS
      sis_files(sissec)%NAME=motcar( adress(4,36) )
      sis_files(sisseo)%NAME=motcar( adress(4,37) )
!     FILE FOR LIQUID BOUNDARIES
      sis_files(sisliq)%NAME=motcar( adress(4,38) )
!     GEOMETRY FILE FORMAT
      sis_files(sisgeo)%FMT = motcar( adress(4,39) )(1:8)
      CALL majus(sis_files(sisgeo)%FMT)
!     NAMES OF PRIVATE VARIABLES
      n_names_priv = min(4,dimens(4,42))
      IF(n_names_priv.GT.0) THEN
        DO k=1,n_names_priv
          names_prive(k) = motcar(adress(4,42)+k-1)(1:32)
        ENDDO
      ENDIF
!
!     ADDITIONAL DIFFERENTIATED VARIABLES
!     NAMES OF THE DIFFERENTIATED VARIABLES
      n_names_advar = dimens(4,70)
      nadvar = max( nadvar,n_names_advar ) ! WARNING HERE ?
      IF(nadvar.GT.0) THEN
        DO i=1,nadvar
          WRITE(char2,'(I2)') i
          names_advar(i) =  'DERIVATIVE '//adjustl(char2)//'   '
     &                   // '??              '
        ENDDO
      ENDIF
      IF(n_names_advar.GT.0) THEN
        DO k=1,n_names_advar
          names_advar(k) = motcar(adress(4,70)+k-1)(1:32)
        ENDDO
      ENDIF
!
!     CVSM, But it's not Beautiful
      tempvar  =   motcar(adress(4,51)   )
      CALL lecdon_split_outputpoints(tempvar,cvsmoutput,cvsm_out_full)
!     FLUXLINEFILE
      sis_files(sisflx)%NAME=motcar( adress(4,69) )
!     C-VSM RESULTS FILE
      sis_files(vspres)%NAME=motcar( adress(4,53) )
      sis_files(vspres)%FMT =motcar( adress(4,55) )(1:8)
      CALL majus(sis_files(vspres)%FMT)
!
      WRITE(lu,102)
102   FORMAT(1x,/,19x, '********************************************',/,
     &            19x, '*               LECDON:                    *',/,
     &            19x, '*        AFTER CALLING DAMOCLES            *',/,
     &            19x, '*        CHECKING OF DATA  READ            *',/,
     &            19x, '*         IN THE STEERING FILE             *',/,
     &            19x, '********************************************',/)
!
!-----------------------------------------------------------------------
!
! LOGICALS FOR OUTPUT VARIABLES
!-----------------------------------------------------------------------
!
      nomblay=max(nomblay,ncouch_tass)
      ncouch_tass=nomblay
!
!     ADDITIONAL DIFFERENTIATED VARIABLES
      CALL nomvar_sisyphe(texte,textpr,mnemo,nsicla,unit,maxvar,
     &  npriv,nomblay,n_names_priv,names_prive,nadvar,names_advar)
!
      CALL sortie(sortis , mnemo , maxvar , sorleo )
      CALL sortie(varim  , mnemo , maxvar , sorimp )
!
      DO i = 1, 4
        IF ((npriv.LT.i).AND.
            ! JWI 31/05/2012 - added 1 to include wave orbital velocities
     &      (sorleo(i+29+(nomblay+4)*nsicla+nomblay).OR.
     &       sorimp(i+29+(nomblay+4)*nsicla+nomblay))) THEN
          npriv=max(npriv,i)
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
!
!     CANCELS OUTPUT OF VARIABLES WHICH ARE NOT BUILT IN THIS CASE
!
      IF(.NOT.susp) THEN
!V 7/09/2006 MIGHT WANT TO OUTPUT THE SUSPENDED COMPONENT IN BIJKER
!        SORIMP(24+4*NSICLA)=.FALSE.
!        SORIMP(25+4*NSICLA)=.FALSE.
!        SORIMP(26+4*NSICLA)=.FALSE.
      ENDIF
!
! JWI 31/05/2012 - added 1 to include wave orbital velocities
      IF(.NOT.charr) THEN
        sorleo(24+(nomblay+2)*nsicla)=.false.
        sorleo(25+(nomblay+2)*nsicla)=.false.
        sorleo(26+(nomblay+2)*nsicla)=.false.
        sorimp(24+(nomblay+2)*nsicla)=.false.
        sorimp(25+(nomblay+2)*nsicla)=.false.
        sorimp(26+(nomblay+2)*nsicla)=.false.
      ENDIF
! JWI END
!
!-----------------------------------------------------------------------
!
! CHECKS TETA VALUE
!
      IF( teta.LT.0.d0.OR.teta.GT.1.d0) THEN
          WRITE(lu,51)
51      FORMAT(/,1x,'BAD VALUE FOR TETA !                   ',/
     &          ,1x,'TETA MUST BE WITHIN 0 AND 1            ')
        CALL plante(1)
        stop
      ENDIF
!
!     INITIALISES MSK (MASKING VARIABLE)
!     FOR NOW MASKING IS ONLY DONE FOR ONE 'OPTION FOR THE TREATMENT
!     OF TIDAL FLATS'. SHOULD BE OFFERED AS AN OPTION FOR THE USER TO
!     CREATE ISLANDS IN THE FUTURE
      msk = .false.
      IF (.NOT.bandec) optban=0
      IF (optban.EQ.2) msk = .true.
!
!-----------------------------------------------------------------------
!
!     CHECKS WHETHER THERE IS A REFERENCE FILE
!
      IF (valid.AND.sis_files(sisref)%NAME.EQ.' ') THEN
          valid=.false.
        WRITE(lu,71)
        WRITE(lu,*)
71      FORMAT(/,1x,'VALIDATION IS NOT POSSIBLE :  ',/
     &          ,1x,'NO REFERENCE FILE  !                 ')
        CALL plante(1)
        stop
      ENDIF
!
!     CHECKS THE NUMBER OF NSICLA
!      IF(NSICLA.GT.10) THEN
      IF(nsicla.GT.nsiclm) THEN
        WRITE(lu,81) nsiclm
        WRITE(lu,*)
81      FORMAT(/,1x,'THE MAXIMUM NUMBER OF SEDIMENT CLASSES IS', i2)
        CALL plante(1)
        stop
      ENDIF
!     CHECKS THE SUM OF INITIAL AVAI
      sumavai = 0.d0
      DO i=1,nsicla
      sumavai = sumavai + ava0(i)
      ENDDO
      IF(abs(sumavai-1).GE.1.d-8) THEN
        WRITE(lu,91)
        WRITE(lu,*)
91      FORMAT(/,1x,'SUM OF SEDIMENT FRACTIONS IS NOT 1  ')
        CALL plante(1)
        stop
      ENDIF
!
!     WARNING FOR THE CHOICE OF RIGID BED METHOD
!
      IF(choix.GT.0.AND.choix.LT.4.AND.vf) THEN
        WRITE(lu,201)
        WRITE(lu,*)
201     FORMAT(/,1x,'FINITE VOLUMES CHOSEN: ',/
     &          ,1x,'METHOD 4 FOR RIGID BED WILL BE USED ')
        choix=4
      ENDIF
      IF (choix.EQ.4.AND..NOT.vf) THEN
        WRITE(lu,301)
        WRITE(lu,*)
301     FORMAT(/,1x,'FINITE ELEMENTS CHOSEN: ',/
     &          ,1x,'METHOD 4 FOR RIGID BED CAN NOT BE USED, METHOD 3 US
     &ED INSTEAD')
        choix=3
      ENDIF
!
!     CHECKS THAT THE BEDLOAD TRANSPORT FORMULATION AND THE HIDING
!     FACTOR FORMULATION CAN GO TOGETHER
!
      IF ((hidfac.EQ.3.AND.icf.NE.6).OR.
     &    (hidfac.EQ.1.AND.icf.NE.1).OR.
     &    (hidfac.EQ.2.AND.icf.NE.1)) THEN
        WRITE(lu,111)
        WRITE(lu,*)
111     FORMAT(/,1x,'CHOICE OF TRANSPORT FORMULA AND HIDING FACTOR FORMU
     &LATION NOT ALLOWED ')
        CALL plante(1)
        stop
      ENDIF
!
!     WITHOUT AND WITH COUPLING, SOME CORRECTIONS
!
      IF(code(1:7).EQ.'TELEMAC'.AND.
     &   sis_files(sishyd)%NAME(1:1).NE.' ') THEN
        sis_files(sishyd)%NAME(1:1)=' '
        WRITE(lu,113)
113     FORMAT(/,1x,'COUPLING: HYDRODYNAMIC FILE IGNORED')
      ENDIF
!
!     COMPUTATION CONTINUED
!
      IF(debu) THEN
        IF(sis_files(sispre)%NAME(1:1).EQ.' ') THEN
          WRITE(lu,313)
313       FORMAT(/,1x,'COMPUTATION CONTINUED:',/,
     &             1x,'PREVIOUS SEDIMENTOLOGICAL FILE MISSING')
          CALL plante(1)
          stop
        ENDIF
      ELSE
        IF(sis_files(sispre)%NAME(1:1).NE.' ') THEN
          sis_files(sispre)%NAME(1:1)=' '
          WRITE(lu,213)
213       FORMAT(/,1x,'NO COMPUTATION CONTINUED:',/,
     &             1x,'PREVIOUS SEDIMENTOLOGICAL FILE IGNORED')
        ENDIF
      ENDIF
!
! METHODS NOT CODED UP FOR SUSPENSION
! -------------------------------------------
!
      IF(susp) THEN
        IF(resol.NE.adv_car   .AND.resol.NE.adv_sup   .AND.
     &     resol.NE.adv_lpo   .AND.resol.NE.adv_nsc   .AND.
     &     resol.NE.adv_psi   .AND.resol.NE.adv_lpo_tf.AND.
     &     resol.NE.adv_nsc_tf.AND.resol.NE.adv_psi_tf      ) THEN
          WRITE(lu,303) resol
303       FORMAT(1x,'RESOLVING METHOD NOT IMPLEMENTED : ',1i6)
          CALL plante(1)
          stop
        ENDIF
      ENDIF
!
      IF(.NOT.houle) sis_files(siscou)%NAME(1:1)=' '
      IF(houle) THEN
        IF(icf.NE.4.AND.icf.NE.5.AND.icf.NE.8.AND.icf.NE.9) THEN
          WRITE(lu,1304) icf
1304      FORMAT(' TRANSPORT FORMULA',1i3,1x,
     &       'DOES NOT TAKE WAVES INTO ACCOUNT,',/,1x,
     &       'TRY 4, 5, 8 OR 9')
          CALL plante(1)
          stop
        ENDIF
      ENDIF
!
! BEDLOAD AND SUSPENDED TRANSPORT COUPLING
! ---------------------------------
!
      IF((icf==30.OR.icf==3.OR.icf==9).AND.susp.AND.charr) THEN
        WRITE(lu,1302) icf
        CALL plante(1)
        stop
      ENDIF
1302  FORMAT('FOR THE FORMULA',1i3,/,1x,
     &       'THE SUSPENSION TERM IS CALCULATED TWICE,'
     &      ,' WITH TOTAL LOAD FORMULA AND SUSPENSION ')
!
! REFERENCE CONCENTRATION
!
! MODIFICATION CV 31/12      IF(ICQ.EQ.2.AND.(PERCOU.NE.1.OR..NOT.CHARR)) THEN
!
      IF(icq.EQ.2.AND.(percou.GT.1.OR..NOT.charr)) THEN
        WRITE(lu,1402) icq
1402  FORMAT('FOR THE BIJKER REFERENCE CONCENTRATION',1i3,/,1x,
     &       'BEDLOAD MUST BE COMPUTED, CHOOSE:',/,1x,
     &       'BEDLOAD = YES')
        CALL plante(1)
        stop
      ENDIF
!
!     CHECKS CONSISTENCY OF BEDLOAD LAWS
!
!     SOULSBY SLOPE EFFECT : REQUIRES A THRESHOLD FORMULA
!
      IF(slopeff.EQ.2) THEN
        IF(icf.NE.1) THEN
          WRITE(lu,1404) icf
1404      FORMAT('BED-LOAD TRANSPORT FORMULA, HERE ICF=',1i3,/,1x,
     &           'MUST HAVE A THRESHOLD',/,1x,
     &           'IF FORMULA FOR SLOPE EFFECT=2 (SOULSBY)')
        ENDIF
      ENDIF
!
! V6P0 : COHERENCE IF CONSOLIDATION MODEL IS USED
! VITCE AND CSF_VASE STEM FROM THE FIRST LAYER OF THE MULTI-LAYER MODEL
! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
! Mots cles supprimes vitesse critique d'erosion et concentration du lit
! CV : si premiere couche est  vide cela n'est pas correct

!
      IF(mixte) THEN
!
!       FILLS VOIDS WITH MUD:
! CV: verifier que la concentration en cohesif est non nulle
!
        csf_sable= 1.d0
!V: verrouiller les options
        nsicla=2
        sedco(1)=.false.
        sedco(2)=.true.
        charr=.false.
        susp=.true.
!V
      ELSE
        csf_sable= (1.d0-xkv)
      ENDIF
!
      IF((.NOT.mixte).AND.sedco(1)) THEN
        charr=.false.
        !SUSP=.TRUE. (In general, but not necessary)
      ENDIF
!
      IF(nomblay.GT.nlaymax) THEN
        WRITE (lu,*) 'NUMBER OF BED LOAD MODEL LAYERS LARGER THAN '
        WRITE (lu,*) 'THE MAXIMUM PROGRAMMED VALUE OF ', nlaymax
        CALL plante(1)
        stop
      ENDIF
!      IF(NOMBLAY.LT.2) THEN
!        WRITE (LU,*) 'BEWARE: NUMBER OF BED LOAD MODEL LAYERS'
!        WRITE (LU,*) '======= LOWER THAN THE DEFAULT VALUE OF 2'
!      ENDIF
!
!----------------------------------------------------------------
!
!  V6P1: FOR THE BED FRICTION PREDICTOR USE LAW OF FRICTION 5 (NIKURADSE)
!
      IF(kspred) kfrot=5
!
!----------------------------------------------------------------
!
      RETURN
      END