The TELEMAC-MASCARET system  trunk
init_zero_gaia.f
Go to the documentation of this file.
1 ! ********************
2  SUBROUTINE init_zero_gaia
3 ! ********************
4 !
5 !
6 !***********************************************************************
7 ! GAIA
8 !***********************************************************************
9 !
11 !
12 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
14 !
16  USE bief
17 !
19  IMPLICIT NONE
20 !
21 !-----------------------------------------------------------------------
22 !
23  INTEGER I,K
24 !
25 !-----------------------------------------------------------------------
26 !
27 !========================================================================
28 ! INITIALISES
29 ! =======================================================================
30 !
31 !---- THE WORKING ARRAYS
32 !
33  IF(npriv > 0) CALL os ('X=0 ', x=prive)
34  CALL os('X=0 ', x=t12 )
35  CALL os('X=0 ', x=coefpn)
36 !
37 !---- THE SEDIMENT VARIABLES :
38 !
39  CALL os('X=0 ', x=qs)
40  CALL os('X=0 ', x=qsx)
41  CALL os('X=0 ', x=qsy)
42  CALL os('X=0 ', x=qsclxc )
43  CALL os('X=0 ', x=qsclyc )
44  CALL os('X=0 ', x=qscl_c)
45  CALL os('X=0 ', x=fluer)
46  CALL os('X=0 ', x=fludp)
47  CALL os('X=0 ', x=fludpt)
48 !
49 ! 7 FOLLOWING LINES ADDED BY JMH 22/04/2005
50 ! PROVISIONAL INITIALISATION FOR FIRST OUTPUT IN RESULTS FILE
51 !
52  CALL os('X=0 ', x=qscl ) ! BLOCK OF SIZE NSICLA
53  CALL os('X=0 ', x=qs_c )
54  CALL os('X=0 ', x=qsxc )
55  CALL os('X=0 ', x=qsyc )
56 !
57  CALL os('X=0 ', x=mudb )
58  CALL os('X=0 ', x=f_mudb )
59 !
60  CALL os('X=0 ', x=evcl_mb )
61  CALL os('X=0 ', x=flbcla )
62 !
63 !---- BED EVOLUTION, CUMULATED BED EVOLUTION
64 ! IN KG/M2 AND IN M
65 !
66  CALL os('X=0 ', x=e )
67  CALL os('X=0 ', x=esomt)
68  CALL os('X=0 ', x=cumbe)
69 !
70 !---- THE HYDRODYNAMIC VARIABLES :
71 !
72  CALL os('X=0 ', x=qu )
73  CALL os('X=0 ', x=qv )
74  CALL os('X=0 ', x=u2d )
75  CALL os('X=0 ', x=v2d )
76  CALL os('X=0 ', x=hn )
77  CALL os('X=0 ', x=q )
78  CALL os('X=0 ', x=tob)
79 !
80 !---- THE WAVE PARAMETERS IF NEED BE
81 !
82 ! ALL INITIALISATIONS OF THE WAVES ARE TO BE REMOVED
83 ! WHEN ALL CHECKS WILL BE DONE
84 ! SEE BEDLOAD_BAILARD_GAIA, DIBWAT, BIJKER AND SOULSBY
85 !
86 !
87 ! FW=0.3 CORRESPONDS TO NO WAVES, 0 WOULD DO A LOG(0)
88  CALL os('X=C ', x=fw ,c=0.3d0 ) !
89  CALL os('X=0 ', x=hw ) !
90  CALL os('X=0 ', x=tw ) !
91  CALL os('X=C ', x=thetaw, c=0.d0) !
92  CALL os('X=C ', x=thetac, c=0.d0) !
93  CALL os('X=0 ', x=uw ) !
94  CALL os('X=0 ', x=tobw) !
95 !
96  IF(nmud.EQ.0) THEN
97  DO i=1,npoin
98  DO k=1,nomblay
99  conc_mud(k,i)=0.d0
100  ENDDO
101  ENDDO
102  ENDIF
103 !-----------------------------------------------------------------------
104 !
105  RETURN
106  END
type(bief_obj), target esomt
Cumulated bed evolution over time [kg/m2].
type(bief_obj), target uw
Orbital wave velocity.
type(bief_obj), target qs
Total solid discharge (bedload+suspension)
type(bief_obj), target e
Evolution of the bed mass at each point for each time step [kg/m2].
type(bief_obj), target fludpt
Deposition flux for implicitation.
type(bief_obj), target u2d
Components of depth-averaged velocity.
type(bief_obj), target tob
Bed shear stress [n/m2].
type(bief_obj), target q
Flow rate.
type(bief_obj), target v2d
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.
double precision, dimension(:), pointer x
2d coordinates of the mesh
type(bief_obj), target qscl
Total transport rate for a sediment class : bedload+suspended load.
type(bief_obj), target qsxc
Solid discharge (bedload), along x and y.
type(bief_obj), target qsx
Solid discharge, along x and y.
type(bief_obj), target thetaw
Wave direction (deg wrt ox axis) !!!!!some say oy axis!!!!!
type(bief_obj), target thetac
Current direction (deg trigo)
type(bief_obj), target fw
Quadratic friction coefficient (waves)
type(bief_obj), target qsyc
type(bief_obj), target qu
X component of the flow rate.
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 hw
Significant wave height.
type(bief_obj), target coefpn
Correction of transport for sloping bed effect.
type(bief_obj), target hn
Water depth.
type(bief_obj), target qsclyc
type(bief_obj), target qv
Y component of the flow rate.
type(bief_obj), target mudb
Mass of mud in bedload added in suspension.
integer nmud
Total number of muds.
type(bief_obj), target qs_c
Solid discharge (bedload)
type(bief_obj), target prive
Block of private vectors.
type(bief_obj), target qsy
type(bief_obj), target qsclxc
Bedload transport rate in the x and y direction for a sediment class [kg*(m-1*s-1)].
subroutine init_zero_gaia
Definition: init_zero_gaia.f:4
type(bief_obj), target cumbe
Cumulated bed evolution over time [m].
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.
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
type(bief_obj), target tobw
Wave induced shear stress.
type(bief_obj), target fludp
Deposition flux.
type(bief_obj), target fluer
Erosion flux.
type(bief_obj), pointer t12
integer npriv
Number of private arrays, number of private arrays with given name.
type(bief_obj), target f_mudb
Flux of mud in bedload added in suspension.
integer, pointer npoin
Number of 2d points in the mesh.
type(bief_obj), target flbcla
Fluxes at boundary for every class.
Definition: bief.f:3