declarations_khione.f

Go to the documentation of this file.

!                 **************************
                  MODULE declarations_khione
!                 **************************
!
!***********************************************************************
! KHIONE   V8P0                                              02/11/2016
!***********************************************************************
!
!brief    DECLARATION OF PRINCIPAL KHIONE VARIABLES
!
!history  F.SOUILLE (EDF)
!+        30/09/2019
!+        V8P0
!+        Added new variables
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief_def
!
!       NOTE: THIS MODULE IS ORGANIZED IN 6 PARTS
!
!       (1) VECTORS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!       (2) INTEGERS
!       (3) LOGICAL VALUES
!       (4) REALS
!       (5) STRINGS
!       (6) ALIASES (IF NECESSARY)
!
!
!=======================================================================
!
!  STRUCTURE USED FOR WRITTING OUT THE TIME SERIES FILE
!
!-----------------------------------------------------------------------
!
      TYPE str_cell_type
        INTEGER :: nval
        CHARACTER(LEN=24), ALLOCATABLE :: txt(:)
        CHARACTER(LEN=24), ALLOCATABLE :: unt(:)
        CHARACTER(LEN=24), ALLOCATABLE :: val(:)
      END TYPE str_cell_type
!
      TYPE str_line_type
        INTEGER :: nval
        CHARACTER(LEN=24) :: htxt
        CHARACTER(LEN=24) :: hunt
        CHARACTER(LEN=24) :: hval
        type(str_cell_type), ALLOCATABLE :: cells(:)
      END TYPE str_line_type
!
!=======================================================================
!
!       0) PARAMETERS
!
!-----------------------------------------------------------------------
!
!     MAXIMUM NUMBER OF VALUES OF KEY-WORDS OF ONE KIND
!       (INTEGER, LOGICAL, ETC.)
!
      INTEGER, PARAMETER :: maxkey = 300
!
!
!=======================================================================
!
!       1) OUTPUTS
!
!-----------------------------------------------------------------------
!
!     VARIABLES FOR GRAPHIC PRINTOUT
!
      CHARACTER(LEN=72) vardes,varimp
!
!     GRAPHIC PRINTOUT PERIOD KHIONE
!
      INTEGER leoprd
!
!     LISTING PRINTOUT PERIOD KHIONE
!
      INTEGER lisprd
!     MAXIMUM NUMBER OF OUTPUT VARIABLES
!
      INTEGER, PARAMETER :: maxvar = 50
!
!     TO KNOW IF A VARIABLE WILL BE EXITED ON FILE, ON LISTING
!
      LOGICAL sorleo(maxvar),sorimp(maxvar)
!
!     BLOCK OF VARIABLES FOR OUTPUT
!
      TYPE(bief_obj), TARGET :: varsor
!
!     NAMES AND UNITS OF VARIABLES
!       ( THESE TWO ARE SET IN LECDON_, THROUGH A CALL TO NOMVAR_ )
      CHARACTER(LEN=32) texte(maxvar),textpr(maxvar)
!
!     VARSOR FILTER FOR CONTINUED COMPUTATION
!
      INTEGER :: alire(maxvar)
!     BE ABSOLUTLY SURE THAT ALIRE IS SIZE MAXVAR
      parameter( alire = (/
!         1,2,3,4,5,6,7,8,9,0,1,2,3,4,5,6,7,8,9,0,1,2,3,4,5,6,7,8,9,0,
     &    0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
     &    0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 /) )
!
!     TITLE OF STUDY
!
      CHARACTER(LEN=72) titicecas
!
!     FOR TIME SERIES OUTPUT
!
      type(str_line_type) :: lines
!
!=======================================================================
!
!      2) FILES
!
!-----------------------------------------------------------------------
!
!
!     MAXIMUM RANK OF LOGICAL UNITS (OPENING FILES) AS DECLARED IN
!     SUBMIT STRINGS IN THE DICTIONARY
!
      INTEGER, PARAMETER :: maxlu_ice = 14
!
!     TELEMAC-2D FILES
!
      TYPE(bief_file) :: ice_files(maxlu_ice)
!
!     KHIONE RESULT FILE NUMBER
!
      INTEGER iceres
!
!     KHIONE BOUNDARY CONDIITON FILE NUMBER
!
      INTEGER icecli
!
!     KHIONE GEOMETRY FILE NUMBER
!
      INTEGER icegeo
!
!     KHIONE STEERING FILE NUMBER
!
      INTEGER icecas
!
!     KHIONE REFERENCE FILE NUMBER
!
      INTEGER iceref
!
!     KHIONE PREVIOUS ICE COVER FILE NUMBER (EULERIAN)
!
      INTEGER icecov
!
!     KHIONE PREVIOUS ICE PACKS FILE NUMBER (LAGRANGIAN)
!
      INTEGER iceblk
!
!     KHIONE CLOGGING ICE RESULTS FILE NUMBER
!
      INTEGER clgrfo
!
!=======================================================================
!
!      3) DICTIONARY
!
!-----------------------------------------------------------------------
!
!     COPY OF SUBMIT STRINGS IN THE DICTIONARY
!
      CHARACTER(LEN=PATH_LEN) submit(4,300)
!
!
!=======================================================================
!
!       4) ALGORITHMIC DIFFERENTIATION
!
!-----------------------------------------------------------------------
!
!
!     NUMBER OF DIFFERENTIATING ARRAYS, AND THOSE WITH A GIVEN NAME
!
      INTEGER nadvar,n_names_advar
!
!     NAMES OF DIFFERENTING ARRAYS (GIVEN BY USER)
!
      CHARACTER(LEN=32) names_advar(maxvar)
!
!     BLOCK OF DIFFERENTIATED VECTORS
!
      TYPE(bief_obj), TARGET :: advar
!
!
!=======================================================================
!
!       5) SWITCHES AND OPTIONS
!
!-----------------------------------------------------------------------
!
!     NUMBER OF DIFFERENTIATING ARRAYS, AND THOSE WITH A GIVEN NAME
!
      INTEGER nicvar,n_names_icvar
!
!     DEBUGGER
      INTEGER debug
!
!     FRAZIL CONCENTRATION
      INTEGER ind_fra
!
!     WATER TEMPERATURE (DUPLICATED INDEX)
      INTEGER ind_t
!
!     SALINITY (DUPLICATED INDEX)
      INTEGER ind_s
!
!     FRAZIL PRECIPITATION (DUPLICATED INDEX)
      INTEGER ind_pre,ind_thi
!
!=======================================================================
!
!       1) HEAT FLUXES
!
!-----------------------------------------------------------------------
!
!     THERMAL BUDGET: NET SUM OF THERMAL FLUXES BETWEEN ATMO AND WATER
      TYPE(bief_obj), TARGET :: sumph
!
!     THERMAL BUDGET: NET SUM OF THERMAL FLUXES BETWEEN ATMO AND ICE
      TYPE(bief_obj), TARGET :: sumph_ice
!
!     THERMAL BUDGET: SOLAR RADIATION UNDER CLEAR SKY
      TYPE(bief_obj), TARGET :: phcl
!
!     THERMAL BUDGET: SOLAR RADIATION UNDER CLOUDY SKY
      TYPE(bief_obj), TARGET :: phri
!
!     THERMAL BUDGET: NET SOLAR RADIATION AFTER REFLECTION
      TYPE(bief_obj), TARGET :: phps
!
!     THERMAL BUDGET: EFFECTIVE BACK RADIATION
      TYPE(bief_obj), TARGET :: phib
!
!     THERMAL BUDGET: EVAPORATION HEAT FLUX
      TYPE(bief_obj), TARGET :: phie
!
!     THERMAL BUDGET: CONDUCTIVITY HEAT FLUX
      TYPE(bief_obj), TARGET :: phih
!
!     THERMAL BUDGET: PRECIPITATION
      TYPE(bief_obj), TARGET :: phip
!
!     THERMAL BUDGET: HEAT TRANSFER BETWEEN WATER AND ICE
      TYPE(bief_obj), TARGET :: phiw
!
!     WATER TEMPERATURE
!      TYPE(BIEF_OBJ), TARGET :: TW
!
!     HEAT TRANSFER COEFFICIENT BETWEEN SURFACE ICE AND WATER
!      TYPE(BIEF_OBJ), TARGET :: HWIN
!      DOUBLE PRECISION HWIN
!
!
!
!=======================================================================
!
!       2) PHYSICAL PROPERTIES AND CONSTANTS
!
!-----------------------------------------------------------------------
!
!     ICE EXPANDS UPON FREEZING, WHICH IS THE OPPOSITE OF MOST OTHER
!       SUBSTANCES. AT 0[oC] IT HAS A DENSITY OF 916.8[kg/m3] AS
!       COMPARED TO A DENSITY OF 999.8[kg/m3] FOR WATER AT THE SAME
!       TEMPERATURE. AS A RESULT, ICE FLOATS IN WATER.
!     NOTE: THE MAXIMUM WATER DENSITY IS 1000[kg/m3] AT 4[oC]
!
!     ICE DENSITY ( DEFAULT = 916.8 kg/m3 )
      DOUBLE PRECISION  rho_ice
!
!     AIR DENSITY ( DEFAULT = 1.225 kg/m3 )
      DOUBLE PRECISION  rho_air
!
!     WATER DENSITY
      DOUBLE PRECISION  ro0
!
!     SPECIFIC HEAT OF WATER (DEFAULT = 4.1855 kJ/kgK )
      DOUBLE PRECISION  cp_eau
!
!     SPECIFIC HEAT OF ICE (DEFAULT = 2.04 kJ/kgK )
      DOUBLE PRECISION  cp_ice
!
!     LATENT HEAT OF FREEZING (DEFAULT = 334. kJ/kgK )
      DOUBLE PRECISION  lh_ice
!
!     CONSTANT, COEFFICIENT FOR LINEARIZED HEAT FLUX EXCHANGES
!     .. BETWEEN WATER AND AIR
      DOUBLE PRECISION  cst_watair,lin_watair
!
!     .. BETWEEN ICE AND AIR
      DOUBLE PRECISION  cst_iceair,lin_iceair
!
!     COEFFICIENT FOR BACK RADIATION HEAT FLUX
      DOUBLE PRECISION coef_phib
!
!     COEFFICIENT FOR EVAPORATION HEAT FLUX
      DOUBLE PRECISION coef_phie
!
!     COEFFICIENT FOR CONDUCTIVE HEAT FLUX
      DOUBLE PRECISION coef_phih
!
!     COEFFICIENT FOR PRECIPITATION HEAT FLUX
      DOUBLE PRECISION coef_phip
!
!     SOLAR CONSTANT (W/M^2)
      DOUBLE PRECISION sio
!
!     ALBEDO OF ICE
      DOUBLE PRECISION albe
!
!
!=======================================================================
!
!       ?) GEOGRAPHICAL CONSIDERATIONS
!
!     TODO: TO BE REPLACED BY A COMPUTATION FROM PHI0 (AND LAMBD0)
!
!     STANDARD LONGITUDE, IN DEG REES
      DOUBLE PRECISION alsm
!
!     LOCAL LONGITUDE, IN  DEGREES
      DOUBLE PRECISION allm
!
!     -1, FOR WEST LONGITUDE; + FOR EAST LONGITUDE
      DOUBLE PRECISION etadir

!
!
!=======================================================================
!
!       ?) FREEZE-UP RELATED PARAMETERS
!
!
!     NUMBER OF CLASSES FOR FRAZIL ICE
!
      INTEGER nc_fra
!
!     FRAZIL CRYSTAL RADIUS [ M ]
      DOUBLE PRECISION, ALLOCATABLE :: rk_frzl(:)
!
!     FRAZIL CRYSTAL THICKNESS [ M ]
      DOUBLE PRECISION, ALLOCATABLE :: ek_frzl(:)
!
!     FRAZIL CRYSTAL VOLUME [ M3 ]
      DOUBLE PRECISION, ALLOCATABLE :: vk_frzl(:)
!
!     NUSSELT NUMBER
      DOUBLE PRECISION nuss
!
!     MODEL FOR COMPUTATION OF THE NUSSELT NUMBER
      DOUBLE PRECISION nussi
!
!     TURBULENT NUSSELT NUMBER
      DOUBLE PRECISION, ALLOCATABLE :: nusst(:)
!
!     TURBULENT MODEL TO ESTIMATE K-EPS
      INTEGER iturb
!
!     MODEL FOR THE COMPUTATION OF THE BUOYANCY VELOCITY
      INTEGER ibuoy
!
!     FRAZIL PARTICLE BUOYANCY VELOCITY
      DOUBLE PRECISION, ALLOCATABLE :: vbk(:)
!
!     FRAZIL CRYSTAL DIAMETER THICKNESS RATIO [ - ]
      DOUBLE PRECISION de
!
!     WATER-FRAZIL THERMAL CONDUCTIVITY [ W/M/oC ]
      DOUBLE PRECISION tc_wt                        !  =>  REPLACES XKWP

!     THERMAL CONDUCTIVITY OF BLACK ICE ( W/M/oC )
      DOUBLE PRECISION tc_bi                        !  =>  REPLACES XKI

!     THERMAL CONDUCTIVITY OF SNOW ( W/M/oC )
      DOUBLE PRECISION tc_s                         !  =>  REPLACES XKS
!
!     SETTLING COEFFICIENT OF FRAZIL ON BARS
      DOUBLE PRECISION af
!
!     POROSITY OF SURFACE ICE
      DOUBLE PRECISION surf_ef
!
!     SECONDARY NUCLEATION MODEL
      INTEGER isnuc
!
!     SEEDDING MODEL
      INTEGER iseed
!
!     MINIMUM NUMBER OF FRAZIL PARTICLES PER UNIT VOLUME
!     (FOR SEEDING MODEL EQ 1)
      INTEGER minnk
!
!     FRAZIL ICE SEEDING RATE
!     (FOR SEEDING MODEL EQ 2)
      DOUBLE PRECISION seedr
!
!     SECONDARY NUCLEATION NMAX PARAMETER
!     FOR SVENSSON AND OMSTEDT COLLISION MODEL
      DOUBLE PRECISION snnmax
!
!     FLOCULATION - BREAKUP MODEL
      INTEGER ifloc
!
!     PARAMETER FOR NET EFFECT OF FLOCULATION - BREAKUP
!     RATIO OF PARTICLES THAT MOVE FROM CLASS I TO I+1
      DOUBLE PRECISION afloc
!
!     MEAN TURBULENT KINETIC ENERGY FOR FRAZIL THERMAL GROWTH
      TYPE(bief_obj), TARGET :: kgm
!
!     MEAN TURBULENT DISSIPATION RATE FOR FRAZIL THERMAL GROWTH
      TYPE(bief_obj), TARGET :: epsgm
!
!     TURBULENT INTENSITY FOR FRAZIL THERMAL GROWTH
      TYPE(bief_obj), TARGET :: alpgm
!
!     TURBULENT VISCOSITY FOR FRAZIL THERMAL GROWTH
      TYPE(bief_obj), TARGET :: nutgm
!
!     FRAZIL ICE SEEDING RATE
      TYPE(bief_obj), TARGET :: seed
!
!     PROBABILITY OF DEPOSITION OF FRAZIL PARTICLES REACHING THE
!       SURFACE LAYER - OPEN WATER
      TYPE(bief_obj), TARGET :: theta0
!
!     PROBABILITY OF DEPOSITION OF FRAZIL PARTICLES REACHING THE
!       SURFACE LAYER - ICE
      TYPE(bief_obj), TARGET :: theta1
!
!     COEFFICIENT QUANTIFYING THE RATE OF REENTRAINMENT OF SURFACE
!       ICE PER UNIT AREA
      TYPE(bief_obj), TARGET :: beta1
!
!     AVERAGE RISING VELOCITY OF FRAZIL ICE IN THE TURBULENT RIVER FLOW
      TYPE(bief_obj), TARGET :: vbb
!
!     ICE CONCENTRATION AS COMPUTED FROM SURFACE ICE PARTICLES
!     (ALSO USED TO DEFINE OPEN WATER VS. UNDER ICE COVER)
      TYPE(bief_obj), TARGET :: anfem
!
!     TOTAL ICE THICKNESS
      TYPE(bief_obj), TARGET :: thifem
!
!     SOLID ICE THICKNESS
      TYPE(bief_obj), TARGET :: thifems
!
!     FRAZIL ICE THICKNESS
      TYPE(bief_obj), TARGET :: thifemf
!
!     UNDERCOVER ICE THICKNESS
      TYPE(bief_obj), TARGET :: hun
!
!     FREEZING POINT OF WATER ( 0.oC BY DEFAULT )
      DOUBLE PRECISION cst_tmelt
      TYPE(bief_obj), TARGET :: tmelt

!     TOTAL NUMBER OF FRAZIL PARTICLE
      TYPE(bief_obj), TARGET :: ntot

!     TOTAL NUMBER OF FRAZIL PARTICLE
      TYPE(bief_obj), TARGET :: ctot
!
!     FRAZIL CONCENTRATION
      TYPE(bief_obj), TARGET :: frzl
!
!     FRAZIL NUMBER OF PARTICLE
      TYPE(bief_obj), TARGET :: nbp
!
!     WATER TEMPERATURE
      TYPE(bief_obj), TARGET :: temp
!
!     WATER SALINITY
      TYPE(bief_obj), TARGET :: sal
!
!=======================================================================
!
!       ?) ICE COVER RELATED PARAMETERS
!
!     ICE FRICTION
      TYPE(bief_obj), TARGET :: icestr

!     TURBULENT INTENSITY
      TYPE(bief_obj), TARGET :: vz
!
!     CRITICAL WATER TEMPERATURE FOR BORDER ICE
      TYPE(bief_obj), TARGET :: tcr
!
!     CRITICAL VELOCITY ABOVE WHICH STATIC BORDER ICE WILL NOT FORM
      DOUBLE PRECISION vcrbor
!
!     CRITICAL VELOCITY ABOVE WHICH DYNAMIC BORDER ICE WILL NOT FORM
      DOUBLE PRECISION vcrbom
!
!     INITIAL CONCENTRATION FOR BORDER / SKIM ICE FORMATION
      DOUBLE PRECISION anfem0      ! REPLACES ANMAXBORDER AND ANMAXSKIM
!
!     LAW OF FRICTION ON UNDER ICE COVER
      INTEGER ifrot
!
!     LAW FOR FRICTION COEFFICIENT
      INTEGER ifice
!
!     FRICTION COEFFICIENT OF THE UNDER ICE COVER
      DOUBLE PRECISION fice

!     MAX FRICTION COEFFICIENT OF THE UNDER ICE COVER
      DOUBLE PRECISION fice_max
!
!     EQUIVALENT SURFACE ICE THCIKNESS (LAW FOR FRICTION COEFFICIENT=1)
      DOUBLE PRECISION thie
!
!
!=======================================================================
!
!       ?) CLOGGING RELATED PARAMETERS
!
!     POROSITY OF ACCUMULATED ICE
      DOUBLE PRECISION clog_ef
!
!     ANGLE OF ACCUMULATED ICE
      DOUBLE PRECISION clog_theta
!
!     PHYSICAL CHARACTERISTICS OF THE INTAKE RACK
      DOUBLE PRECISION clog_vdist,clog_vdiam,clog_tdist,clog_tdiam
!
!     CLOGGED BOUNDARY NUMBERS
      INTEGER nfrclog,nseclog,ntotclog
      INTEGER, ALLOCATABLE :: numclog(:),seclog(:)
      INTEGER, ALLOCATABLE :: nseg_fs(:),liste_fs(:,:,:)
      DOUBLE PRECISION, ALLOCATABLE :: un1(:),un2(:)
!
!     GROWTH THICKNESS OF THE VERTICAL AND HORIZONTAL ACCRECION
      DOUBLE PRECISION, ALLOCATABLE :: clog_vwdth(:),clog_twdth(:)
!
!     PHYSICAL SIZE OF LENGTHS OF THE VERTICAL AND HORIZONTAL BARS
      DOUBLE PRECISION, ALLOCATABLE :: clog_vlgth(:),clog_tlgth(:)
!
!     TOTAL VOLUM OF ACCRECION
      DOUBLE PRECISION, ALLOCATABLE :: clog_volum(:)
!
!     CLOGGING SECTIONS
      TYPE(bief_obj), TARGET :: msksec
!
!=======================================================================
!
!       ?) ICE CHARACTERISATION
!
!     REPLACEMENT OF:
!       ICE ISLAND INDEX:   TYPE(BIEF_OBJ), TARGET :: JAMFEM
!       ICE TYPE INDEX:     TYPE(BIEF_OBJ), TARGET :: ICEREGION
!     ICETYPE IS:
!       1 - OPEN WATER
!       2 - STATIC BORDER ICE, FULLY WITHIN THE BORDER ICE COVER
!       3 - STATIC BORDER ICE THAT HAS GROWN A SUFFICIENTLY THICK COVER
!     ICELOC IS:
!       1 - OPEN WATER
!       2 - BOUNDARY NODE
!       3 - EDGE NODE
      TYPE(bief_obj), TARGET :: icetype
      TYPE(bief_obj), TARGET :: icetypep
      TYPE(bief_obj), TARGET :: iceloc
!
!     INTEGERS WORKING ARRAYS
      TYPE(bief_obj), TARGET :: it1,it2
!
!-----------------------------------------------------------------------
!
!     IMAGE ELEMENT
      INTEGER, ALLOCATABLE :: imgelm(:,:)

      !PARCEL'S NUMBER OF BOOM IMAGE
      INTEGER, ALLOCATABLE ::naphabm(:,:)

!     TRASH RACK CROSS SECTION NODES
      INTEGER, ALLOCATABLE :: ncgno(:),nfzno(:)
!
      INTEGER ntwbc
!
!     NUMBER OF INTAKES
      INTEGER nitk

!     SWTICH FOR FRAZIL CLOGGING
      INTEGER nclog

      INTEGER, PARAMETER :: maxfro = 50

!     MAXIMUM TOTAL NEIGHBORING PARCEL NUMBER
      INTEGER maxkrn

!      NUNBER OF ICE BOOM
      INTEGER nbm

!      MAXIMUM NUNBER OF ICE BOOM
      INTEGER maxbm

!     NUMBER OF ICE BOOM SPANS
      INTEGER nbmsp

!     ICE BOOM TYPES
      INTEGER bmtype

!      NUMBER OF SLOT TO ISSUE ICE
      INTEGER itbdf

!      NUMBER OF BREAKUP REACHES
      INTEGER ncov

!      NUMBER OF SOLID BOUNDARY
      INTEGER nmlbn

!      NUMBER OF OPEN BOUNDARY
      INTEGER nmobn

!      NUMBER OF UNDERCOVER BOUNDARY
      INTEGER nuqsec

!     TYEP OF UNDERCOVER BOUDNARY CONDITION
      INTEGER iuqbc

!     COUNTER FOR UNDERCOVER BOUNDARY CONDITION
      INTEGER iundt

      INTEGER itmax
      !     TO KNOW IF A VARIABLE WILL BE EXITED ON FILE, ON LISTING
!      NUMBER OF ICE BOUNDARY
      INTEGER msbd

!      ICE CONCENTRATION BE SPECIFIED
      INTEGER iqbjb

!      NUMBER OF SOLID BOUNDARY
      INTEGER nlb

!     NUMBER OF OPEN BOUNDARY
      INTEGER nob

!     NUMBER OF TIMESTEP FOR UNDERCOVER
      INTEGER iund

      INTEGER ntbn,nhbn,nqbn,n0bn,nhtn,nfbn,nqb1,nq0,nq01,
     &  nel,nsbn

!     NUMBER OF PARCEL IN THE DOMAIN
      INTEGER nparin

!     NUMBER OF PARCEL MOVING OUT OF THE DOMAIN
      INTEGER nparout

!     DELETE THE BAD PARCELS
      INTEGER nbadparcels

!     TOTAL NUMBER OF PARCELS
      INTEGER ntotalparcels

!     FILE NUMBER FOR UNDERCOVER BOUNDARY FILE
      INTEGER iudc

!     FILE NUMBER FOR ICE BOOM FILE
      INTEGER inbm

!     FILE NUMBER FOR ICE BOOM FORCE FILE
      INTEGER ibmf

!     FILE NUMBER FOR BREAK UP FILE
      INTEGER ibrk

!     FILE NUMBER FOR ICE CONDITION ON FEM NODE FILE
      INTEGER ihdw0

!     FILE NUMBER FOR FRAZIL DEPOSITION AND ENTRAINMENT
      INTEGER ifrz

!     FILE NUMBER FOR TRASH TRACK AND CROSS SECTIONS
      INTEGER icsg

!     ICE FLUX TIME FRAME
      INTEGER kflxi

      INTEGER kprt
      INTEGER kdmph
      INTEGER kdmpi
      INTEGER kdmphi

      INTEGER kwnd
      INTEGER kelv,kflx,
     &  khdw,khdi,khot,kstr,kthr,irnd
!          ! UPDATE BOUNDARY AND OUTPUT
     & jobflx,nsec,ndisc,nelev,istrt,
     & ihdw,ihot0,ihot,istr0,istr,ilnk,
     & ipro,iplt,iclg,iqbd

!     NUMBER OF ICE PARCELS
      INTEGER knp

!     STEPS TO CHECK BREAKUP
      INTEGER nbkdtm

!     NUMBER OF REACHES FOR BREAKUP
      INTEGER nrch

!     ATMOSPHERE-WATER EXCHANGE MODEL
!
      INTEGER atmoexch
!
!     ! INDEX FOR DISCHARGE USED FOR FRAZIL ACCRETION
!      INTEGER, PARAMETER :: KCGQ = 1
!
!
!
!
!-----------------------------------------------------------------------
!
!       3) LOGICAL
!
!     WITH ICEDYNAMICS CALCULATION
      LOGICAL icedynamics
!
!     IF YES, MASS-BALANCE
      LOGICAL icebilmas
!
!     IF YES, VALIDATION
      LOGICAL icevalid
!
!     STEADY HYDRODYNAMICS
      LOGICAL  perma
!
!     IF TEMPERATURE IS VARIABLE
      LOGICAL yatemp

!      IF YES, TECPLOT OUTPUT
      LOGICAL ioutput

!     IF YES, ICE INPUT
      LOGICAL iciput
!
!     HEAT EXCHANGES
      LOGICAL thermal_budget
!
!     ICE COVER IMPACT ON HYDRO
      LOGICAL icover_impact
!
!     CLOGGING
      LOGICAL clogging
!
!     BORDER_ICE
      LOGICAL bd_ice
!
!     SALINITY
      LOGICAL salinity
!
!     FRZAIL PRECIPITATION
      LOGICAL prec
!
!     SCHEME OPTION FOR THERMAL GROWTH
      INTEGER itgm
!
!     ENERGY BALANCE VERSION
      INTEGER inrjb
!
!-----------------------------------------------------------------------
!
!       4) REALS
!
!-----------------------------------------------------------------------

      DOUBLE PRECISION roi
!  ICE LATENT HEAT
!      DOUBLE PRECISION XLATEN    !       => MOVED TO KHIONE UNDER LH_ICE

      !  CLOUDINESS
      DOUBLE PRECISION cc

!     HEIGHT THE WIND DATA IS MEASUR ED(M)
      DOUBLE PRECISION z1

!     ELEVATION MEASURED FROM SEA  LEVEL(M)
      DOUBLE PRECISION zh

!--------------------------------------------------------------
!  GROUP FOR SIMULATING WATER TEMPERATURE AND THERMAL GROWTH AND DECAY OF ICE

!     COEFFICIENT FOR LINEARIZED HEAT FLUX CALCULATION ON WATER-AIR SURFACE
      DOUBLE PRECISION hwa

!     COEFFICIENT FOR LINEARIZED HEAT FLUX CALCULATION ON ICE-AIR SURFACE
      DOUBLE PRECISION hia

!     COEFFICIENT FOR LINEARIZED HEAT FLUX CALCULATION ON ICE-AIR SURFACE
      DOUBLE PRECISION alp

!     COEFFICIENT FOR LINEARIZED HEAT FLUX CALCULATION ON WATER-AIR SURFACE
      DOUBLE PRECISION alpw
!
!     CONSTANT FOR CALCULATING HEAT TRANSFER COEFFICIENT BETWEEN RIVER WATER AND ICE FOR TURBULENT FLOW
      DOUBLE PRECISION cwi1
!
!     CONSTANT FOR CALCULATING HEAT TRANSFER COEFFICIENT BETWEEN RIVER WATER AND ICE FOR SUPERCOOLED TURBULENT FLOW
      DOUBLE PRECISION ciw1
!
!     NUSSELT NUMBER FOR CALCULATING HEAT TRANSFER COEFFICIENT BETWEEN RIVER WATER AND ICE FOR LAMINAR FLOW
      DOUBLE PRECISION ata

!     HEAT TRANSFER COEFFICIENT FROM WATER TO ICE
      DOUBLE PRECISION cwi
!
!     KINEMATIC VISCOSITY OF WATER
      DOUBLE PRECISION xnu

!------------------------------------------------------------------------
!   GROUP FOR SIMULATING BORDER ICE FORMATION
!
!     CRITICAL WATER SURFACE TEMPERATURE FOR BORDER ICE FORMATION
      DOUBLE PRECISION tc
!
!     CHANNEL WIDTH FOR THE COMPUTATION OF SURFACE WATER TEMPERATURE
      DOUBLE PRECISION bch
!
!     MAXIMUM CONCENTRATION FOR BORDER ICE FORMATION
!      DOUBLE PRECISION ANMAXBORDER
!
!----------------------------------------------------------------------------
! GROUP FOR SKIM ICE SIMULATION
!
!     MAXIMUM CONCENTRATION FOR SKIM ICE RUN
!      DOUBLE PRECISION ANMAXSKIM                 => REPLACED BY ANFEM0
!
!--------------------------------------------------------------------------
! GROUP FOR MASS EXCHANGE BETWEEN SUSPENDED LAYER, SURFACE LAYER AND ANCHOR ICE

!     POROSITY OF SURFACE ICE
      DOUBLE PRECISION ef

!     POROSITY OF FRAZIL ICE
      DOUBLE PRECISION ef0

!     VERTICAL BAR LENGTH
      DOUBLE PRECISION hbar
!
!-------------------------------------------------------------------------------
!  GROUP FOR ICE COVER FORMATION
!
!     INITIAL AND CURRENT WIDTH FOR VERTICAL BARS
      DOUBLE PRECISION wfv0,wfv1
!
!     INITIAL AND CURRENT WIDTH FOR TRANSVERSE BARS
      DOUBLE PRECISION wft0,wft1
!
!     INITIAL AND CURRENT WIDTH FOR VERTICAL BARS
      DOUBLE PRECISION wfv0dw,wfv1dw
!
!     INITIAL AND CURRENT WIDTH FOR TRANSVERSE BARS
      DOUBLE PRECISION wft0dw,wft1dw
!
!
      DOUBLE PRECISION time,corio,pi,sqpi
      DOUBLE PRECISION  rho,rhoi,raddig

!     GEOMETRY MAX & MIN
      DOUBLE PRECISION xmin,ymin,xmax,ymax,xgmin,ygmin

!     INITIAL ICE PARCEL CONCENTRATION
      DOUBLE PRECISION an0

!     INITIAL FRAZIL ICE THICKNESS
      DOUBLE PRECISION thi0f

!     INITIAL SIZE OF ICE PARCEL
      DOUBLE PRECISION hpi0
!
!     A MINIMUM VALUE
      DOUBLE PRECISION vsm
!
!     NUMBER OF ICE PARCELS
!     DOUBLE PRECISION KNP
!
!     ICE PARCEL MASS
      DOUBLE PRECISION pm1
!
      DOUBLE PRECISION angle   ! MAY NOT NEED;   WIND ANGLE
!
      DOUBLE PRECISION refel   ! MAY NOT NEED; REFERENCE LEVEL FOR HYDRO
!
      DOUBLE PRECISION optim  ! MAY NOT NEED; OPTIMAL TIME STEP FOR HYDRO
!
!     ICE MASS DENSITY AND ICE MASS AT TN
      DOUBLE PRECISION umi0,pmi0,pmi0f

!     FRICTION COEFFICIENT BETWEEN ICE-RIVERBANK; ICE-RIVERBED
      DOUBLE PRECISION fric1,fric2
!
!     AIR-ICE DRAG COEFFICIENT
      DOUBLE PRECISION wdca

!     ICE INTERNAL FRICTION ANGLE
      DOUBLE PRECISION phi
!
!     EMPIRICAL CONSTANT
      DOUBLE PRECISION pj
!
!     MAX WATER-ICE DRAG COEFFICIENT
      DOUBLE PRECISION cwmax

!     BOLTZMANN CONSTANT (WM-2K-4)
      DOUBLE PRECISION sgma

!     SPECIFIC HEAT OF WATER
      DOUBLE PRECISION cp

!     TIME FOR ICE DYNAMICS
      DOUBLE PRECISION timice0
      DOUBLE PRECISION tmn
      DOUBLE PRECISION endtime

!      UNDERCOVER TIME STEP
      DOUBLE PRECISION undt

!     STARTING TIME OF ICE DYNAMICS
      DOUBLE PRECISION ticebg

!     BOOM CURTAIN DEPTH
      DOUBLE PRECISION dwall

      ! WEATHER TIME
      DOUBLE PRECISION tima0,tima1,tair0,tair1


!     VARIABLES FOR BOUNDARY CONDITIONS AND NODAL INFORMATION FOR
!     HYDRODYNAMICS MAY NOT BE NEEDED
      DOUBLE PRECISION sinbn,cosbn,hb,
     & qb,stemp,ctemp,disch,dischu,
     & wndtime,prttime,dmptimh,dmptimi,
     & fac,qbt,qbti,qin,xbgn,ybgn,xend,yend

!     ICE BOOM COORDINATE
      DOUBLE PRECISION, ALLOCATABLE :: xbm(:,:),ybm(:,:)

!     ICE BOOM SPAN COORDINATE
      DOUBLE PRECISION, ALLOCATABLE :: xbmsp(:,:),ybmsp(:,:)

!     FORCE PER UNIT LENGTH IN SPAN J
      DOUBLE PRECISION, ALLOCATABLE :: asignp(:)

!     ICE BOOM FLOAT OR SUBMERGE
      DOUBLE PRECISION, ALLOCATABLE :: frnum1(:)

!     ICE THICKNESS AT THE ICE BOOM J
      DOUBLE PRECISION, ALLOCATABLE :: thisp(:)

!     CRITICAL SUBMERGENCE LINE LOAD OF BOOM
      DOUBLE PRECISION, ALLOCATABLE :: bmlineld(:)

!     CRITICAL ICE THICKNESS ABOVE SUBMERGED BOOM FOR RE-EMERGENCE
      DOUBLE PRECISION, ALLOCATABLE :: thiemerg(:)

!     BOTTOM ELEVATION OF THE GATE
      DOUBLE PRECISION, ALLOCATABLE :: gate(:)

!     BUOYANT FORCE OF THE ICE BOOM WHEN IT IS SUBMERGED
      DOUBLE PRECISION, ALLOCATABLE :: bmbuoy(:)

!     FRICTION COEFFICIENT BETWEEN THE ICE BOOM AND THE BOTTOM OF THE ICE
      DOUBLE PRECISION, ALLOCATABLE :: bmfric(:)

!     DIRECTION OF THE SOLID BOUNDARY
      DOUBLE PRECISION, ALLOCATABLE :: papha1(:),papha2(:)

!     DISTANCE BETWEEN ICE PARCEL AND SOLID BOUNDARY
      DOUBLE PRECISION, ALLOCATABLE :: pardis(:)

!     DIRECTION OF THE ICE BOOM
      DOUBLE PRECISION, ALLOCATABLE :: paphabm1(:),paphabm2(:)

!     DISTANCE BETWEEN ICE APARCEL AND ICE BOOM
      DOUBLE PRECISION, ALLOCATABLE :: pardisbm(:)

!!     DIRECTION OF THE ICE BOOM (SIN,COS)
!      DOUBLE PRECISION, ALLOCATABLE :: ALPHA11(:),ALPHA12(:)

!!     DISTANCE BETWEEN ICE PARCEL AND ICE BOOM
!      DOUBLE PRECISION, ALLOCATABLE :: PARDS1(:)

!     TIME FOR ICE COVER BREAKUP
      DOUBLE PRECISION, ALLOCATABLE :: brktim(:)

!     CRITICAL VALUE OF WATER SURFACE ELEVATION AT REFERENCE NODE
      DOUBLE PRECISION, ALLOCATABLE :: brkwse(:)

!     CRITICAL VALUE OF WATER SURFACE ELEVATION INCREASE AT REFERENCE NODE
      DOUBLE PRECISION, ALLOCATABLE :: brkdwse(:)

!     WATER SURFACE ELEVATON RECORD FOR BREAKUP
      DOUBLE PRECISION, ALLOCATABLE :: wserd(:,:)

!     INITIAL ICE COVER THICKNESS
      DOUBLE PRECISION, ALLOCATABLE :: covti(:)

!      FINITE ELEMENT ORIENTATION AND HEIGHT
      DOUBLE PRECISION, ALLOCATABLE :: alafae(:),elmht(:)

!      FINITE ELEMENT MAX EDGE IN X,Y
      DOUBLE PRECISION, ALLOCATABLE :: xmaxh(:),ymaxh(:)

!      FINITE ELEMENT INFORMAITON
      DOUBLE PRECISION, ALLOCATABLE :: a1(:),a2(:),a3(:)
      DOUBLE PRECISION, ALLOCATABLE :: b1(:),b2(:),b3(:)
      DOUBLE PRECISION, ALLOCATABLE :: c1(:),c2(:),c3(:)
      DOUBLE PRECISION, ALLOCATABLE :: clmp(:)
!
!-----------------------------------------------------------------------
!      LINK FILE
      CHARACTER(LEN=50) lnkfile

!      ICE BOUNDARY FILE
      CHARACTER(LEN=50) iqbfile

!      ICE BOOM FILE AND BOOM FORCE FILE
      CHARACTER(LEN=50) bmfile,bffile

!      BREAKUP FILE
      CHARACTER(LEN=50) bkfile

!      UNDERCOVER TRANSPORT FILE
      CHARACTER(LEN=50) udfile

!     FRAZIL COEFFICIENTS FILE
      CHARACTER(LEN=50) fzfile

!     TRASH RACK CROSS SECTION FILE
      CHARACTER(LEN=50) trfile

!     FRAZIL CLOGGING OUTPUT FILE
      CHARACTER(LEN=50) clgfile

!     EXISTING ICE CONDITION FILE
      CHARACTER(LEN=50) ht0file,st0file,hd0file

      CHARACTER(LEN=50) casename,geofile,timfile,prtfile,
     &  ptmfile,flxfile,elvfile,wndfile,dscfile,
     &  ho0file,errfile,state,thrfile,delfile,
     &  qnifile,qnwfile,dbmfile,
     &  lelfile,nbmfile,profile,tbcfile,linkefile,swifile,mcrfile,
     &  weafile,icevolumefile,netfile,covfile,parfile,
     &  clgdwfile,dscloc,elvloc,dsciloc

      INTEGER, PARAMETER :: max_ice_file = 1000

!      TYPE(BIEF_FILE) :: PLTFILE(MAX_ICE_FILE)

      CHARACTER(LEN=50) pltfile(max_ice_file)

      CHARACTER(LEN=50) hdwfile(max_ice_file)

      CHARACTER(LEN=50) hotfile(max_ice_file)

      CHARACTER(LEN=50) strfile(max_ice_file)
!
!=======================================================================
!
!       5) BLOCKS
!
!-----------------------------------------------------------------------
!
!     BLOCK OF ICE VARIABLES
!
      TYPE(bief_obj), TARGET :: icvar
!
!
!     BLOCK OF WORKING ARRAYS
!
      TYPE(bief_obj), TARGET :: tb
!
!=======================================================================
!
!       X) POINTERS
!
!-----------------------------------------------------------------------
!
!     DECLARATION OF POINTERS FOR ALIASES TO WORKING VECTORS IN TB
!     TARGETS ARE DEFINED IN POINT_KHIONE.
!
      TYPE(bief_obj),POINTER :: t1,t2,t3,t4,t5,t6
!
!=======================================================================
!
!       X) WORKING MATRIX AND ARRAY
!
!-----------------------------------------------------------------------
!
      TYPE(bief_obj),TARGET :: bm1,bm2
      TYPE(bief_obj),TARGET :: cv1,cv2,cv3,cv4
!
!=======================================================================
!
      SAVE

      END MODULE declarations_khione