The TELEMAC-MASCARET system  trunk
deall_gaia.f
Go to the documentation of this file.
1 ! **************************
2  SUBROUTINE deall_gaia
3 ! **************************
4 !
5 !
6 !***********************************************************************
7 ! GAIA
8 !***********************************************************************
9 !
11 !
12 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
14 !
15  USE bief
18 !
19  IMPLICIT NONE
20 
21 !
22 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
23 !
24  INTEGER :: I,K
25 !
26 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
27 !
28 !
29  CALL outbief(mesh)
30 ! Deallocate the mesh structure
31  CALL dealmesh(mesh)
32 
33  NULLIFY(ikle)
34  NULLIFY(x)
35  NULLIFY(y)
36  NULLIFY(nelem)
37  NULLIFY(nelmax)
38  NULLIFY(nptfr)
39  NULLIFY(nptfrx)
40  NULLIFY(typelm)
41  NULLIFY(npoin)
42  NULLIFY(npmax)
43  NULLIFY(mxptvs)
44  NULLIFY(lv)
45 !
46 
47 
48  CALL bief_deallobj(varsor)
49 ! Deallocating all the blocks first
50  CALL bief_deallobj(s )
51  CALL bief_deallobj(e )
52  CALL bief_deallobj(z )
53  CALL bief_deallobj(esomt )
54  CALL bief_deallobj(cumbe )
55  CALL bief_deallobj(emax )
59  CALL bief_deallobj(q )
60  CALL bief_deallobj(qu )
61  CALL bief_deallobj(qv )
62  CALL bief_deallobj(u2d )
63  CALL bief_deallobj(v2d )
64  CALL bief_deallobj(thetac)
65  CALL bief_deallobj(qs )
66  CALL bief_deallobj(qsx )
67  CALL bief_deallobj(qsy )
68  CALL bief_deallobj(qs_c )
69  CALL bief_deallobj(qsxc )
70  CALL bief_deallobj(qsyc )
71  CALL bief_deallobj(hiding)
72  CALL bief_deallobj(zf )
73  CALL bief_deallobj(zr )
74  CALL bief_deallobj(radsec)
75  CALL bief_deallobj(zref )
76  CALL bief_deallobj(coefpn)
77  CALL bief_deallobj(cf )
78  CALL bief_deallobj(tob )
79  CALL bief_deallobj(tobw )
82  CALL bief_deallobj(taup )
83  CALL bief_deallobj(mu )
84  CALL bief_deallobj(ksp )
85  CALL bief_deallobj(ks )
86  CALL bief_deallobj(ksr )
87  CALL bief_deallobj(thetaw)
88  CALL bief_deallobj(fw )
89  CALL bief_deallobj(uw )
90  CALL bief_deallobj(hw )
91  CALL bief_deallobj(tw )
92  CALL bief_deallobj(acladm)
93  CALL bief_deallobj(unladm)
94  CALL bief_deallobj(hcpl )
95  CALL bief_deallobj(ecpl )
96  CALL bief_deallobj(elay )
97  CALL bief_deallobj(estrat)
98  CALL bief_deallobj(kx )
99  CALL bief_deallobj(ky )
100  CALL bief_deallobj(kz )
103  CALL bief_deallobj(unorm )
104  CALL bief_deallobj(maskb )
105  CALL bief_deallobj(mask )
106  CALL bief_deallobj(afbor )
107  CALL bief_deallobj(bfbor )
108  CALL bief_deallobj(flbor )
109  CALL bief_deallobj(q2bor )
110 ! BOUNDARY FLUX FOR CALL TO CVDFTR
111  CALL bief_deallobj(flbor_gai )
112  CALL bief_deallobj(flbortra )
113  CALL bief_deallobj(cstaeq)
114 ! MAK ADDITION
115  CALL bief_deallobj(csratio)
116  CALL bief_deallobj(hn )
117  CALL bief_deallobj(hprop )
118  CALL bief_deallobj(volu2d)
119  CALL bief_deallobj(v2dpar)
120  CALL bief_deallobj(unsv2d)
121  CALL bief_deallobj(mpm_aray)
123 
124  ! MSK
125  CALL bief_deallobj(maskel)
126  CALL bief_deallobj(msktmp)
127  CALL bief_deallobj(maskpt)
128 ! FOR MIXED SEDIMENTS
129  CALL bief_deallobj(ms_sable )
130  CALL bief_deallobj(ms_vase )
131  CALL bief_deallobj( liebor)
132  CALL bief_deallobj( liqbor)
133  CALL bief_deallobj( limtec)
134  CALL bief_deallobj( numliq)
135  CALL bief_deallobj( clt )
136  CALL bief_deallobj( clu )
137  CALL bief_deallobj( clv )
138  CALL bief_deallobj( limdif)
139  CALL bief_deallobj( lihbor)
141  CALL bief_deallobj( limpro)
142  CALL bief_deallobj( it1 )
143  CALL bief_deallobj( it2 )
144  CALL bief_deallobj( it3 )
145  CALL bief_deallobj( it4 )
146 ! NUMBER OF LAYERS
147  CALL bief_deallobj( nlayer)
148  CALL bief_deallobj(breach)
149  CALL bief_deallobj(ifamas)
150 
151  DEALLOCATE(avail)
152  DEALLOCATE(es)
153  DEALLOCATE(m2t)
154  DEALLOCATE(mpa2t)
155 
156  DEALLOCATE(ivide)
157 
158  DEALLOCATE(sanfra)
159 
160  CALL bief_deallobj(masktr)
161  CALL bief_deallobj(ebor )
162  CALL bief_deallobj(qbor )
163  DO k=1,nomblay
164  ALLOCATE(laythi%ADR(k)%P%R(1))
165  ENDDO
166  CALL bief_deallobj(laythi)
167  CALL bief_deallobj(layconc)
169  CALL bief_deallobj(tocemud)
170  CALL bief_deallobj(parthe)
171 !
172  CALL bief_deallobj(qscl )
173  CALL bief_deallobj(qscl_c)
174  CALL bief_deallobj(qsclxc)
175  CALL bief_deallobj(qsclyc)
176  CALL bief_deallobj(mudb)
177  CALL bief_deallobj(f_mudb)
178 !
179  CALL bief_deallobj(qscl_s)
180  CALL bief_deallobj(zfcl_c)
181  CALL bief_deallobj(fludp )
182  CALL bief_deallobj(fludpt)
183  CALL bief_deallobj(fluer )
184  CALL bief_deallobj(fluert)
185 !
186  CALL bief_deallobj(zfcl_ms)
187  CALL bief_deallobj(evcl_mb)
188  CALL bief_deallobj(calfa_cl)
189  CALL bief_deallobj(salfa_cl)
190  CALL bief_deallobj(flbcla)
191  ! Parallel to what was done in point_gaia
192  CALL bief_deallobj(ratios)
193  CALL bief_deallobj(mass_s)
194  CALL bief_deallobj(ratiom)
195  CALL bief_deallobj(mass_m)
196 
197  CALL bief_deallobj(am1_s)
198  CALL bief_deallobj(am2_s)
199  CALL bief_deallobj(mbor )
200 
201  CALL bief_deallobj( w1 )
202  CALL bief_deallobj( te1)
203  CALL bief_deallobj( te2)
204  CALL bief_deallobj( te3)
205 
206  CALL bief_deallobj(varcl)
207  CALL bief_deallobj(prive)
208  CALL bief_deallobj(tb )
209 
210  IF(nadvar.GT.0) THEN
211  CALL bief_deallobj(advar )
212  ENDIF
213 ! to check: all variables of new model deallocated?
214 !
215  DEALLOCATE(ratio_mud_sand)
216  DEALLOCATE(mass_sand)
217  DEALLOCATE(mass_mud)
218  DEALLOCATE(mass_mud_tot)
219  DEALLOCATE(mass_sand_tot)
220  DEALLOCATE(mass_mix_tot)
221  DEALLOCATE(ratio_sand)
222  DEALLOCATE(ratio_mud)
223  DEALLOCATE(ratio_debimp)
224  DEALLOCATE(toce_sand)
225  DEALLOCATE(toce_mud)
226  DEALLOCATE(conc_mud)
227  DEALLOCATE(trans_mass)
228  DEALLOCATE(partheniades)
229  DEALLOCATE(qer_mud)
230  DEALLOCATE(qer_sand)
231  DEALLOCATE(time)
232  DEALLOCATE(toce_mix)
233  DEALLOCATE(cae_ilay)
234  DEALLOCATE(qe_moy)
235  DEALLOCATE(fluer_pur_mud)
236  DEALLOCATE(fluer_pur_sand)
237  DEALLOCATE(fluer_mix)
238  DEALLOCATE(conc_mud_activ_layer)
239  DEALLOCATE(mass_mud_activ_layer)
240  DEALLOCATE(mass_sand_active_layer)
241  DEALLOCATE(mass_sand_masked)
242  DEALLOCATE(evcl_m_tot_sand)
243  DEALLOCATE(ratio_evol_tot_sand)
244  DEALLOCATE(flux_mass_mud)
245  DEALLOCATE(flux_mass_mud_tot)
246  DEALLOCATE(flux_mass_sand)
247  DEALLOCATE(num_transf)
248  DEALLOCATE(flux_neg_mud_activ_layer)
249  DEALLOCATE(flux_neg_sand_activ_layer)
250  DEALLOCATE(flux_pos_mud_activ_layer)
251  DEALLOCATE(flux_pos_sand_activ_layer)
252  IF(bilma) THEN
253  DEALLOCATE(bedload_b_flux)
254  DEALLOCATE(sumbedload_b)
255  ENDIF
256  IF(bilma.OR.vsmtype==1) THEN
257  DEALLOCATE(mcumucla)
258  ENDIF
259 !
260  ! Resetting variable
261  init_fluxpr = .true.
262  IF(deja_rfc) THEN
263  deja_rfc = .false.
264  DEALLOCATE(infic_rfc)
265  DEALLOCATE(time_rfc)
266  ENDIF
267  IF(deja_flusec) THEN
268  deja_flusec = .false.
269  old_method_flusec = .false.
270  DEALLOCATE(nseg)
271  DEALLOCATE(liste)
272  DEALLOCATE(volnegs)
273  DEALLOCATE(volposs)
274  DEALLOCATE(volnegc)
275  DEALLOCATE(volposc)
276  DEALLOCATE(flx)
277  DEALLOCATE(volneg)
278  DEALLOCATE(volpos)
279  DEALLOCATE(flxs)
280  DEALLOCATE(flxc)
281  ENDIF
282  IF(deja_flusec2) THEN
283  old_method_flusec = .false.
284  deja_flusec2 = .false.
285  DO i = 1,numberoflines_flusec2
286  DEALLOCATE(fluxlinedata_flusec2(i)%SECTIONIDS)
287  DEALLOCATE(fluxlinedata_flusec2(i)%DIRECTION)
288  ENDDO
289  DEALLOCATE(fluxlinedata_flusec2)
290  DEALLOCATE(flux_flusec2)
291  DEALLOCATE(volflux_flusec2)
292  ENDIF
293 
294  ! from lecdon
295  DEALLOCATE(soldis)
296  DEALLOCATE(okcgl)
297  DEALLOCATE(okqgl)
298  DEALLOCATE(ctrlsc)
299  ! from lecdon_gaia
300  IF (bilma) THEN
301  DEALLOCATE(sumbedload_b_flux)
302  ENDIF
303 
304  ! from lecdon_telemac2d_gaia and lecdon_telemac3d_gaia
305  IF (ALLOCATED(schadvsed)) THEN
306  DEALLOCATE(schadvsed)
307  DEALLOCATE(optadv_sed)
308  DEALLOCATE(slvsed)
309  DEALLOCATE(sed0)
310  DEALLOCATE(presed)
311  DEALLOCATE(sedsce)
312  ENDIF
313  IF (ALLOCATED(verprosed)) THEN
314  DEALLOCATE(verprosed)
315  DEALLOCATE(dnusedh)
316  DEALLOCATE(dnusedv)
317  ENDIF
318 
319 
320  DEALLOCATE(type_sed)
321  DEALLOCATE(num_imud_icla)
322  DEALLOCATE(num_icla_imud)
323  DEALLOCATE(num_isand_icla)
324  DEALLOCATE(num_icla_isand)
325  DEALLOCATE(num_isusp_icla)
326 !
327 !-----------------------------------------------------------------------
328 !
329 !=======================================================================
330 !
331 ! WRITES OUT TO LISTING :
332 !
333 ! IF(LISTIN) THEN
334 ! IF(LNG.EQ.1) WRITE(LU,22)
335 ! IF(LNG.EQ.2) WRITE(LU,23)
336 ! ENDIF
337 !22 FORMAT(1X,///,21X,'****************************************',/,
338 ! &21X, '* FIN DE LA DEALLOCATION DE LA MEMOIRE : *',/,
339 ! &21X, '****************************************',/)
340 !23 FORMAT(1X,///,21X,'*************************************',/,
341 ! &21X, '* END OF MEMORY ORGANIZATION: *',/,
342 ! &21X, '*************************************',/)
343 !
344 !-----------------------------------------------------------------------
345 !
346  RETURN
347  END SUBROUTINE deall_gaia
double precision, dimension(:,:,:), allocatable cae_ilay
integer, pointer nptfr
Number of boundary points.
type(bief_obj), target esomt
Cumulated bed evolution over time [kg/m2].
type(bief_obj), target hcpl
Water depth saved for constant flow discharge.
type(bief_obj), target kx
Coefficients of the dispersion tensor (dim. npoin)
type(bief_obj), target uw
Orbital wave velocity.
type(bief_obj), target qs
Total solid discharge (bedload+suspension)
double precision, dimension(:), allocatable ratio_evol_tot_sand
Ratio between EVCL_M_TOT_SAND and the mass in the active layer.
type(bief_obj), target s
Void structure.
logical bilma
Mass balance.
type(bief_obj), target ms_vase
type(bief_obj), target zfcl_c
Bed evolution per class (due to bedload)
double precision, dimension(:,:), allocatable, target trans_mass
Mass transfer for consolidation between layers, for each bed layer, for each point. It is variable in space for the active layer (layer 1) only. It is a choice related to the way the bed model is done: the active layer is treated along with the others, there are no specialised arrays for the active layer.
type(bief_obj), target liqbor
Type of boundary conditions on sand transport rate.
type(bief_obj), target breach
Array saying whether the non-erodable bottom has been reached (vf)
type(bief_obj), target ecpl
Evolution saved for constant flow discharge.
type(bief_obj), target e
Evolution of the bed mass at each point for each time step [kg/m2].
double precision, dimension(:), allocatable mass_sand_active_layer
Mass of sand in active layer (array defined for temporary work in some subroutines) ...
type(bief_obj), target laythi
Layer thicknesses (for output variables)
type(bief_obj), target fludpt
Deposition flux for implicitation.
type(bief_obj), target maskel
Mask.
type(bief_obj), target acladm
Mean diameter of active-layer.
type(bief_obj), target ky
type(bief_obj), pointer ikle
Connectivity table.
type(bief_obj), target zfcl_ms
Evolution for each class due to sloping bed effects.
double precision, dimension(:,:), allocatable fluer_mix
double precision, dimension(:,:), allocatable ratio_mud_sand
Ratio of mud to sand, for ilayer,ipoin.
double precision, dimension(:), allocatable flxc
type(bief_obj), target limtec
Type of boundary conditions.
double precision, dimension(:), allocatable evcl_m_tot_sand
Evolution of total (sum from each class) mass of sand for each point.
type(bief_obj), target u2d
Components of depth-averaged velocity.
type(bief_obj), target uconv_gai
Components of velocity vectors.
type(bief_obj), target tob
Bed shear stress [n/m2].
double precision, dimension(:), allocatable flxs
double precision, dimension(:,:), allocatable, target toce_mud
Critical erosion shear stress of the mud, for each bed layer, for each point. It is variable in space...
subroutine outbief(MESH)
Definition: outbief.f:6
type(bief_obj), target q
Flow rate.
double precision, dimension(:,:,:), allocatable, target mass_mud
Surface mass of mud (kg/m2), for imud,ilayer,ipoin.
type(bief_obj), target calfa_cl
Cosinus of the angle between mean flow q and transport qs.
type(bief_obj), target nlayer
Number of layers for each point.
integer, pointer npmax
Maximum number of 2d points in the mesh.
double precision, dimension(:,:), allocatable infic_rfc
type(bief_obj), target v2d
type(bief_obj), target kz
type(bief_obj), target evcl_mb
Mass evolution for class (due to bedload)
integer, target nomblay
Number of bed load model layers = NUMSTRAT+1 to take the active layer into account.
type(bief_obj), target it3
double precision, dimension(:), pointer x
2d coordinates of the mesh
integer, dimension(:), allocatable verprosed
subroutine dealmesh(MESH)
Definition: dealmesh.f:6
type(bief_obj), target zr
Non erodable (rigid) bottom elevation.
type(bief_obj), target evol_mb
Total mass evolution due to bedload [kg/m2].
type(bief_obj), target te3
type(bief_obj), target volu2d
Integral of bases.
double precision, dimension(:), allocatable qer_sand
double precision, dimension(:), allocatable flx
integer, pointer lv
Vector length (for vectorisation) from mesh object.
type(bief_obj), target qscl_s
Suspended load transport rate for a sediment class.
type(bief_obj), target ratiom
Sediment composition : mud (for output variables)
double precision, dimension(:,:), allocatable, target es
Layer thicknesses as double precision.
type(bief_obj), target fluert
Erosion flux for implicitation.
double precision, dimension(:,:), allocatable mass_sand_tot
Surface total mass of sand (kg/m2), for ilayer,ipoin.
double precision, dimension(:,:), allocatable flux_neg_sand_activ_layer
Negative sand transfer to apply in case we have an active layer.
type(bief_obj), target varcl
Block of clandestine variables.
type(bief_obj), target mbor
Boundary matrix.
type(bief_obj), target qscl
Total transport rate for a sediment class : bedload+suspended load.
subroutine bief_deallobj(OBJ)
Definition: bief_deallobj.f:7
double precision, dimension(:), allocatable volneg
type(bief_obj), target afbor
Flux condition nu df/dn=afbor * f + bfbor.
type(bief_obj), target clu
double precision, dimension(:,:), allocatable fluer_pur_mud
double precision, dimension(:), allocatable ratio_debimp
Ratio of sand in the prescribed solid discharge.
type(bief_obj), target qsxc
Solid discharge (bedload), along x and y.
double precision, dimension(:), allocatable dnusedh
Coefficient for horizontal diffusion of suspended sediments in 3D.
type(bief_obj), target am1_s
Matrices.
type(bief_obj), target qbor
Imposed solid transport at the boundary In m3/s, for every class.
integer, pointer typelm
Type of element.
double precision, dimension(:), allocatable volnegc
double precision, dimension(:,:,:), allocatable flux_mass_mud
Flux of mud transfert inbetween layers.
type(bief_obj), target flbor_gai
Flux at the boundaries.
double precision, dimension(:,:), allocatable sumbedload_b
Cumulated bedload on boundary for every class (kg): variable for mass balance.
type(bief_obj), target ms_sable
For mixed sediments.
type(bief_obj), target tobcw_max
Maximum of total current + wave shear stress.
type(bief_obj), target qsx
Solid discharge, along x and y.
type(bief_obj), target ks
Total bed roughness.
double precision, dimension(:), pointer y
type(bief_obj), target evol_ms
Total mass evolution due to suspension [kg/m2].
double precision, dimension(:,:,:), allocatable, target mass_sand
Surface mass of sand (kg/m2), for isand,ilayer,ipoin.
type(bief_obj), target it4
type(bief_obj), target numliq
Liquid boundary numbering.
type(bief_obj), target thetaw
Wave direction (deg wrt ox axis) !!!!!some say oy axis!!!!!
integer, dimension(:), allocatable num_isusp_icla
Tables to switch from suspended sediment number to class number.
type(bief_obj), target flulim_gai
Flux limitation per segment.
type(fluxline), dimension(:), allocatable fluxlinedata_flusec2
type(bief_obj), target thetac
Current direction (deg trigo)
type(bief_obj), target flbortra
Flux at the boundaries for tracer.
integer, pointer nptfrx
Maximum number number of boundary points.
type(bief_obj), target fw
Quadratic friction coefficient (waves)
type(bief_obj), target clt
Boundary conditions for tracer, u and v (modified litbor, liubor,livbor)
type(bief_obj), target it2
double precision, dimension(:), allocatable time_rfc
type(bief_obj), target qsyc
type(bief_obj), target parthe
Partheniades erosion coefficient (for output variables)
double precision, dimension(:), allocatable soldis
Prescribed solid discharges.
type(bief_obj), target ratios
Sediment composition : sand (for output variables)
type(bief_obj), target am2_s
type(bief_obj), target qu
X component of the flow rate.
type(bief_obj), target salfa_cl
Sinus of the angle between mean flow q and transport qs.
double precision, dimension(:,:), allocatable, target conc_mud
Mud concentration for each bed layer, at each point. It is variable in space for the active layer (la...
type(bief_obj), target elay
Active layer thickness.
type(bief_obj), target maskb
Mask.
double precision, dimension(:,:), allocatable mass_mud_tot
Surface total mass of mud (kg/m2), for ilayer,ipoin.
double precision, dimension(:), allocatable qe_moy
type(bief_obj), target te2
type(bief_obj), target w1
Working arrays.
type(bief_obj), target limdif
Type of boundary conditions for diffusion.
double precision, dimension(:,:), allocatable bedload_b_flux
Bedload boundary flux for every class (kg/s): variable for mass balance.
integer, dimension(:), allocatable nseg
integer, dimension(:), allocatable num_isand_icla
type(bief_obj), target clv
double precision, dimension(:), allocatable volposc
type(bief_obj), target hw
Significant wave height.
type(bief_obj), target estrat
Active stratum thickness.
type(bief_obj), target coefpn
Correction of transport for sloping bed effect.
type(bief_obj), target hprop
Propagation height.
type(bief_obj), target advar
Block of differentiated vectors.
integer, dimension(:), allocatable num_imud_icla
Tables to switch from mud number to class number and from sand number to class number.
type(bief_obj), target flbor
Zf values on boundaries.
type(bief_obj), target radsec
Curve radius for secondary currents.
integer, pointer nelem
Number of elements in the mesh.
type(bief_obj), target tb
Blocks of working arrays.
character(len=3), dimension(:), allocatable type_sed
Type of sediment (co or nco)
type(bief_obj), target te1
Work arrays for elements.
type(bief_obj), target bfbor
type(bief_obj), target hn
Water depth.
type(bief_obj), target unorm
Norm of the mean flow velocity.
double precision, dimension(:), allocatable sed0
Initial values of suspended sediments.
type(bief_obj), target unsv2d
Inverse of integral of bases.
double precision, dimension(:), allocatable mass_sand_masked
Contains "hidden" sand mass if more than 30% of mud (because no bedload for sand in that case) ...
logical, dimension(:), allocatable okqgl
Used in function qgl_gaia!
type(bief_obj), target masktr
Mask.
double precision, dimension(:,:), allocatable fluer_pur_sand
type(bief_obj), target layconc
Layer concentration (for output variables)
type(bief_obj), target q2bor
Imposed solid transport at the boundary In m2/s, total, read in the boundary conditions file...
double precision, dimension(:,:), allocatable, target partheniades
Partheniades erosion coefficient, for each bed layer, for each point. It is variable in space for the...
type(bief_obj), target evol_mm
Total mass evolution due to sliding (maxslope) [kg/m2].
type(bief_obj), target zref
Reference elevation.
integer vsmtype
For the Continous Vertical Sorting MODEL.
double precision, dimension(:,:,:), allocatable, target ratio_mud
Ratio of mud to all muds, for imud,ilayer,ipoin.
double precision, dimension(:,:), allocatable, target m2t
conversion mass to thickness
type(bief_obj), target z
Free surface elevation.
type(bief_obj), target csratio
Ratio between bottom concentration and average concentration.
double precision, dimension(:), allocatable presed
Prescribed values of suspended sediments.
type(bief_obj), target qsclyc
type(bief_obj), target qv
Y component of the flow rate.
double precision, dimension(:), allocatable volposs
type(bief_obj), target msktmp
Mask.
type(bief_obj), target mudb
Mass of mud in bedload added in suspension.
double precision, dimension(:), allocatable mass_mud_activ_layer
Mass of mud in active layer (array defined for temporary work in some subroutines) ...
type(bief_obj), target mu
Skin friction correction factor for bed roughness: Ratio between shear stress due skin friction and t...
type(bief_obj), target mask
Block of masks.
type(bief_obj), target qs_c
Solid discharge (bedload)
integer, dimension(:), allocatable schadvsed
Scheme for advection of suspended sediments.
double precision, dimension(:,:,:), allocatable, target ratio_sand
Ratio of sand to all sands, for isand,ilayer,ipoin.
type(bief_obj), target prive
Block of private vectors.
double precision, dimension(:,:,:), allocatable flux_pos_sand_activ_layer
Positive sand transfer to apply in case we have an active layer.
double precision, dimension(:), allocatable dnusedv
Coefficient for vertical diffusion of suspended sediments in 3D.
type(bief_obj), target qsy
integer, dimension(:), allocatable num_icla_imud
Tables to switch from class number to mud number and from class number to sand number.
double precision, dimension(:), allocatable qer_mud
type(bief_obj), target mpm_aray
Meyer Peter Mueller factor.
type(bief_obj), target mass_m
Mud mass (for output variables)
type(bief_obj), target qsclxc
Bedload transport rate in the x and y direction for a sediment class [kg*(m-1*s-1)].
double precision, dimension(:,:), allocatable toce_sand
Critical erosion shear stress of the sand, for each sand, for each point.
double precision, dimension(:,:), allocatable volflux_flusec2
type(bief_obj), target taup
Shear stress modified by skin friction.
type(bief_obj), target cumbe
Cumulated bed evolution over time [m].
type(bief_obj), target ifamas
Like IFABOR but ignoring masked elements.
type(bief_obj), target qscl_c
Bedload transport rate for a sediment class [kg*(m-1*s-1)].
type(bief_obj), target tw
Mean wave period.
double precision, dimension(:,:,:), allocatable flux_pos_mud_activ_layer
Positive mud transfer to apply in case we have an active layer.
integer, dimension(:), allocatable num_transf
Layer of mud of active layer into which mud will be added.
type(bief_obj), target vconv_gai
logical, dimension(:), allocatable okcgl
Used in function cgl_gaia.
type(bief_obj), target ebor
Imposed bed evolution at the boundary [m].
type(bief_obj), target cstaeq
Sediment equilibrium concentration.
integer, dimension(:), allocatable num_icla_isand
type(bief_obj), target varsor
Block of variables for output.
double precision, dimension(:,:), allocatable flux_mass_mud_tot
Same as FLUX_MASS_MUD but summed on the mud layer.
type(bief_obj), target cf
Quadratic friction coefficient.
type(bief_obj), target zf
Bottom elevation.
type(bief_obj), target tobw
Wave induced shear stress.
type(bief_mesh), target mesh
Mesh structure.
double precision, dimension(:), allocatable, target mpa2t
conversion mass per area to thickness
type(bief_obj), target mass_s
Sand mass (for output variables)
double precision, dimension(:), allocatable mcumucla
cumulated mass through bedload boundary per class per time step (kg/s)
type(bief_obj), target fludp
Deposition flux.
type(bief_obj), target liebor
Type of boundary conditions on bed evolution.
type(bief_obj), target hiding
Hiding factor correction.
double precision, dimension(:,:,:), allocatable flux_mass_sand
Flux of sand transfert inbetween layers.
type(bief_obj), target fluer
Erosion flux.
type(bief_obj), target tobcw_mean
Mean of total current + wave shear stress.
double precision, dimension(:), allocatable volnegs
integer, pointer mxptvs
Maximum number of neighbouring points.
integer, pointer nelmax
Maximum number of elements in the mesh.
double precision, dimension(:), allocatable time
double precision, dimension(:), allocatable volpos
type(slvcfg), dimension(:), allocatable slvsed
Solver for the diffusion of sediment.
double precision, dimension(:,:,:), allocatable, target avail
Sediment fraction for each layer, class, point.
type(bief_obj), target ksr
Ripple bed roughness.
double precision, dimension(:,:), allocatable sedsce
Values of suspended sediments at the sources.
double precision, dimension(:), allocatable sumbedload_b_flux
Sum over classes of bedload boundary flux or cumulated bedload:
type(bief_obj), target ksp
Bed skin roughness.
subroutine deall_gaia
Definition: deall_gaia.f:4
double precision, dimension(:,:,:), allocatable toce_mix
type(bief_obj), target tocemud
Critical erosion shear stress of the mud (for output variables)
type(bief_obj), target f_mudb
Flux of mud in bedload added in suspension.
integer, dimension(:,:,:), allocatable liste
type(bief_obj), target it1
Integer working arrays.
type(bief_obj), target unladm
Mean diameter of under-layer.
double precision, dimension(:,:), allocatable mass_mix_tot
Surface total mass of sediments (kg/m2), for ilayer,ipoin.
integer nadvar
Number of differentiating arrays, and those with a given name.
type(bief_obj), target lihbor
Type of boundary conditions for h.
integer, pointer npoin
Number of 2d points in the mesh.
double precision, dimension(:,:), allocatable flux_neg_mud_activ_layer
Negative mud transfer to apply in case we have an active layer.
type(bief_obj), target limpro
Type of boundary conditions for propagation.
double precision, dimension(:,:), allocatable, target ivide
Void index of bed layers.
type(bief_obj), target v2dpar
Integral of bases in parallel.
type(bief_obj), target boundary_colour
Last column of the boundary condition file.
type(bief_obj), target emax
Maximum evolution to be defined by the user for automatic.
integer, dimension(:), allocatable ctrlsc
Array containing the global number of the points in the control sections.
integer, dimension(:), allocatable optadv_sed
Advection scheme options for suspended sediments.
double precision, dimension(:), allocatable conc_mud_activ_layer
Concentation of mud in active layer (array defined for temporary work in some subroutines) ...
double precision, dimension(:), allocatable sanfra
sand fraction for Wilcock & Crowe transport formula
type(bief_obj), target flbcla
Fluxes at boundary for every class.
Definition: bief.f:3
type(bief_obj), target mtransfer
Mass transfer for consolidation between layers (for output variables)
double precision, dimension(:,:), allocatable flux_flusec2
type(bief_obj), target maskpt
Mask on points.