The TELEMAC-MASCARET system  trunk
init_zero.f
Go to the documentation of this file.
1 ! ********************
2  SUBROUTINE init_zero
3 ! ********************
4 !
5 !
6 !***********************************************************************
7 ! SISYPHE V6P1 21/07/2011
8 !***********************************************************************
9 !
10 !brief INITIALISES VARIABLES.
11 !
12 !history JMH
13 !+ 04/06/2008
14 !+
15 !+ INITIALISATION OF MASDEP
16 !
17 !history C. VILLARET
18 !+
19 !+ V5P9
20 !+
21 !
22 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
23 !+ 13/07/2010
24 !+ V6P0
25 !+ Translation of French comments within the FORTRAN sources into
26 !+ English comments
27 !
28 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
29 !+ 21/08/2010
30 !+ V6P0
31 !+ Creation of DOXYGEN tags for automated documentation and
32 !+ cross-referencing of the FORTRAN sources
33 !
34 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
35 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
36 !
38  USE bief
39 !
41  IMPLICIT NONE
42 !
43 !-----------------------------------------------------------------------
44 !
45  INTEGER I
46 !
47 !-----------------------------------------------------------------------
48 !
49 !========================================================================
50 ! INITIALISES
51 ! =======================================================================
52 !
53 !---- THE WORKING ARRAYS
54 !
55  IF(npriv > 0) CALL os ('X=0 ', x=prive)
56  CALL os('X=0 ', x=t12 )
57  CALL os('X=0 ', x=coefpn)
58 !
59 !---- THE SEDIMENT VARIABLES :
60 !
61  CALL os('X=0 ', x=qs)
62  CALL os('X=0 ', x=qsx)
63  CALL os('X=0 ', x=qsy)
64  CALL os('X=0 ', x=qsclxc )
65  CALL os('X=0 ', x=qsclyc )
66  CALL os('X=0 ', x=qsclxs )
67  CALL os('X=0 ', x=qsclys )
68 !
69 ! PROVISIONAL INITIALISATION FOR FIRST OUTPUT IN RESULTS FILE
70 !
71  CALL os('X=0 ', x=qscl ) ! BLOCK OF SIZE NSICLA
72  CALL os('X=0 ', x=qs_s )
73  CALL os('X=0 ', x=qsxs )
74  CALL os('X=0 ', x=qsys )
75  CALL os('X=0 ', x=qs_c )
76  CALL os('X=0 ', x=qsxc )
77  CALL os('X=0 ', x=qsyc )
78 !
79 ! PROBABLY USEFUL ONLY IF(CHARR) AND WITH FINITE ELEMENTS
80 !
81  CALL os('X=0 ', x=zfcl_c )
82 !
83 !---- THE DEPOSITION MASSES FOR EVERY CLASS IN SUSPENSION
84 !
85  DO i=1,nsicla
86  masdep(i)=0.d0
87  ENDDO
88 !
89  CALL os('X=0 ', x=e )
90  CALL os('X=0 ', x=esomt)
91  CALL os('X=0 ', x=cs )
92 !
93 !---- THE HYDRODYNAMIC VARIABLES :
94 !
95  CALL os('X=0 ', x=qu )
96  CALL os('X=0 ', x=qv )
97  CALL os('X=0 ', x=u2d )
98  CALL os('X=0 ', x=v2d )
99  CALL os('X=0 ', x=hn )
100  CALL os('X=0 ', x=q )
101  CALL os('X=0 ', x=tob)
102 !
103 !---- THE WAVE PARAMETERS IF NEED BE
104 !
105 ! ALL INITIALISATIONS OF THE WAVES ARE TO BE REMOVED
106 ! WHEN ALL CHECKS WILL BE DONE
107 ! SEE BEDLOAD_BAILARD, DIBWAT, BIJKER AND SOULSBY
108 !
109 !
110 ! FW=0.3 CORRESPONDS TO NO WAVES, 0 WOULD DO A LOG(0)
111  CALL os('X=C ', x=fw ,c=0.3d0 ) !
112  CALL os('X=0 ', x=hw ) !
113  CALL os('X=0 ', x=tw ) !
114  CALL os('X=C ', x=thetaw, c=90.d0) !
115  CALL os('X=0 ', x=uw ) !
116  CALL os('X=0 ', x=tobw) !
117 !
118 !-----------------------------------------------------------------------
119 !
120  RETURN
121  END
type(bief_obj), target thetaw
type(bief_obj), target e
type(bief_obj), target zfcl_c
type(bief_obj), target qsy
type(bief_obj), target esomt
type(bief_obj), target qsxc
type(bief_obj), target fw
type(bief_obj), target coefpn
type(bief_obj), target qsys
type(bief_obj), target qsxs
type(bief_obj), target qsclxc
type(bief_obj), target tobw
type(bief_obj), target qsclxs
type(bief_obj), target cs
subroutine init_zero
Definition: init_zero.f:4
double precision, dimension(:), pointer x
type(bief_obj), target tw
type(bief_obj), target qs
type(bief_obj), target qs_c
type(bief_obj), target qsyc
type(bief_obj), target qsclyc
type(bief_obj), target qs_s
type(bief_obj), target v2d
type(bief_obj), target hw
type(bief_obj), target q
type(bief_obj), target qv
type(bief_obj), target qu
type(bief_obj), pointer t12
type(bief_obj), target qsclys
type(bief_obj), target prive
type(bief_obj), target hn
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
type(bief_obj), target qscl
type(bief_obj), target tob
double precision, dimension(nsiclm) masdep
type(bief_obj), target u2d
type(bief_obj), target uw
type(bief_obj), target qsx
Definition: bief.f:3