The TELEMAC-MASCARET system  trunk
declarations_artemis.f
Go to the documentation of this file.
1 ! ***************************
3 ! ***************************
4 !
5 !
6 !***********************************************************************
7 ! ARTEMIS V8P2
8 !***********************************************************************
9 !
10 !brief DECLARATION OF PRINICIPAL ARTEMIS VARIABLES
11 !
12 !history J-M HERVOUET (LNH)
13 !+
14 !+
15 !+ LINKED TO BIEF 5.0
16 !
17 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
18 !+ 13/07/2010
19 !+ V6P0
20 !+ Translation of French comments within the FORTRAN sources into
21 !+ English comments
22 !
23 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
24 !+ 21/08/2010
25 !+ V6P0
26 !+ Creation of DOXYGEN tags for automated documentation and
27 !+ cross-referencing of the FORTRAN sources
28 !
29 !history C.PEYRARD(EDF)
30 !+ 2011
31 !+ V6P1
32 !+ VARIABLES FOR RAPIDLY VARYING TOPOGRAPHY ADDED
33 !
34 !history C.PEYRARD(EDF)
35 !+ 2012
36 !+ V6P2
37 !+ VARIABLES FOR INCIDENT POTENTIAL ADDED
38 !
39 !history C.PEYRARD(EDF)
40 !+ 2013
41 !+ V6P3
42 !+ VARIABLES FOR CURRENT ADDED
43 !
44 !history C.PEYRARD(EDF)
45 !+ 2014
46 !+ V7P0
47 !+ VARIABLES FOR AUTOMATIC TETAP CALCULATION ADDED
48 !+ VARIABLES FOR AUTOMATIC PHASE CALCULATION ADDED
49 !
50 !history N.DURAND (HRW)
51 !+ November 2016
52 !+ V7P2
53 !+ New variables introduced for use in LECLIM call : ALFAPS and HBS
54 !+ These and TETAPS store parameters read from the cli file (RP also
55 !+ read in but does not require a temporary storage variable for it)
56 !
57 !history N.DURAND (HRW)
58 !+ November 2016
59 !+ V7P2
60 !+ Addition of new keywords in the steering file re: animation of the
61 !+ free surface. ANIMFS and ARTAMP declared
62 !
63 !history N.DURAND (HRW)
64 !+ August 2017
65 !+ V7P3
66 !+ 1. Consolidation with CHAINTWC, meaning that it is now a choice of integers
67 !+ and no longer a logical
68 !+ 2. NDTWC (now NDIR) and NFTWC (now NF) are no longer read in the
69 !+ steering file but directly extracted from the spectrum file
70 !+ 3. Declaration of new structure type: SPECTRUM
71 !+ 4. PI (and related variables) now defined here and in ARTEMIS_CONSTANTS
72 !+ 5. Addition of new keywords in the steering file for nesting option 2:
73 !+ X_SFREF, Y_SFREF, ART_FILES(WACRES)%NAME, and related variables
74 !+ 6. Introduction of BDALE, block holding spatially varying DALEs for
75 !+ incident boundary nodes and for nesting option 2
76 !+ 7. Addition of IPTFR_REF, the node number for the reference spectrum
77 !
78 !history N.DURAND (HRW)
79 !+ December 2018
80 !+ V8P0
81 !+ Addition of new keywords in the steering file for nesting option 2:
82 !+ ART_FILES(WACLQD)%NAME
83 !
84 !history N.DURAND (HRW)
85 !+ January 2019
86 !+ V8P0
87 !+ TYPE SPECTRUM MIGRATED TO BIEF_DEF
88 !
89 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
90 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
91 !
92  USE bief_def
93 !
94 ! NOTE: THIS MODULE IS ORGANISED IN 10 PARTS
95 !
96 ! 1) VECTORS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
97 ! 2) MATRICES (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
98 ! 3) BLOCKS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
99 ! 4) INTEGERS
100 ! 5) LOGICAL VALUES
101 ! 6) REALS
102 ! 7) STRINGS
103 ! 8) SLVCFG STRUCTURES
104 ! 9) MESH STRUCTURE
105 ! 10) ALIASES
106 !
107 !-----------------------------------------------------------------------
108 !
109 ! 1) VECTORS
110 !
111 !-----------------------------------------------------------------------
112 !
113 !brief REAL PART OF WAVE POTENTIAL
114 ! partie reelle du potentiel
115  TYPE(bief_obj), TARGET :: phir
116 !brief IMAGINARY PART OF WAVE POTENTIAL
117 ! partie imaginaire du potentiel
118  TYPE(bief_obj), TARGET :: phii
119 !brief WATER DEPTH AT REST
120 ! hauteur d'eau au repos
121  TYPE(bief_obj), TARGET :: h
122 !brief (MEAN) WAVE NUMBER
123 ! nombre d'onde
124  TYPE(bief_obj), TARGET :: k
125 !brief (MEAN) PHASE CELERITY
126 ! vitesse de phase
127  TYPE(bief_obj), TARGET :: c
128 !brief (MEAN) GROUP CELERITY
129 ! vitesse de groupe
130  TYPE(bief_obj), TARGET :: cg
131 !brief SIGNIFICANT WAVE HEIGHT (REGULAR MODE)
132 ! hauteur de la houle
133  TYPE(bief_obj), TARGET :: hho
134 !brief WAVE PHASE (REGULAR MODE)
135 ! phase de la houle
136  TYPE(bief_obj), TARGET :: phas
137 !brief SURFACE WAVE VELOCITY COMPONENT
138 ! vitesse en surface (a t=0)
139  TYPE(bief_obj), TARGET :: u0
140 !brief SURFACE WAVE VELOCITY COMPONENT
141 ! vitesse en surface (a t=0)
142  TYPE(bief_obj), TARGET :: v0
143 !brief MEAN COSINE OF WAVE DIRECTION
144 ! moyennes des cosinus de la direction de houle
145  TYPE(bief_obj), TARGET :: mcos
146 !brief MEAN SINE OF WAVE DIRECTION
147 ! moyennes des sinus de la direction de houle
148  TYPE(bief_obj), TARGET :: msin
149 !brief WAVE INCIDENCE (OR DIRECTION)
150 ! incidence de la houle
151  TYPE(bief_obj), TARGET :: inci
152 !brief FREE SURFACE ELEVATION
153 ! cote de la surface libre
154  TYPE(bief_obj), TARGET :: s
155 !brief BOTTOM ELEVATION
156 ! cote du fond
157  TYPE(bief_obj), TARGET :: zf
158 !brief FRICTION FACTOR
159 ! coefficient de frottement (variable en espace)
160  TYPE(bief_obj), TARGET :: fw
161 !brief WAVE HEIGHT (RANDOM WAVE)
162 ! hauteur de la houle aleatoire
163  TYPE(bief_obj), TARGET :: hale
164 !brief WAVE PERIODS ARRAY (RANDOM MODE)
165 ! tableau des periodes de discretisation du spectre pour un calcul en houle aleatoire multidirectionnelle
166  TYPE(bief_obj), TARGET :: pale
167 !brief REFLEXION COEFFICIENT
168 ! coefficient de reflexion des parois
169  TYPE(bief_obj), TARGET :: rp
170 !brief ANGLE OF WAVE ATTACK (FROM NORMAL AXIS)
171 ! angle d'attaque de la houle sur les limites - pas seulement les parois (compte par rapport a a la normale , inclut dans [0;90])
172  TYPE(bief_obj), TARGET :: tetap
173 !brief DEPHASING CAUSED BY THE WALLS
174 ! dephasage induit par la paroi entre l'onde reflechie et l'onde incidente (si alfap est positif, l'onde reflechie est en retard)
175  TYPE(bief_obj), TARGET :: alfap
176 !brief INCIDENT WAVE HEIGHT AT THE BOUNDARY
177 ! hauteur de la houle aux frontieres ouvertes
178  TYPE(bief_obj), TARGET :: hb
179 !brief INCIDENT WAVE DIRECTION AT THE BOUNDARY (FROM X AXIS)
180 ! angle d'attaque de la houle aux frontieres ouvertes (compte par rapport a l'axe des x dans le sens direct)
181  TYPE(bief_obj), TARGET :: tetab
182 !brief REAL PART OF INCIDENT WAVE AT THE BOUNDARY
183 ! partie reelle du potentiel impose au bord (dirichlet)
184  TYPE(bief_obj), TARGET :: phirb
185 !brief IMAGINARY PART OF INCIDENT WAVE AT THE BOUNDARY
186 ! partie imaginaire du potentiel impose au bord (dirichlet)
187  TYPE(bief_obj), TARGET :: phiib
188 
189 !brief REAL PART OF INCIDENT POTENTIAL AT THE BOUNDARY
190 ! partie reelle du potentiel incident au bord
191  TYPE(bief_obj), TARGET :: prb
192 !brief IMMAGINARY PART OF INCIDENT POTENTIAL AT THE BOUNDARY
193 ! partie imaginaire du potentiel incident au bord
194  TYPE(bief_obj), TARGET :: pib
195 !brief REAL PART OF GRADIENT X COMPONENT OF INCIDENT POTENTIAL AT THE BOUNDARY
196 ! partie reele du gradient en X du potentiel incident au bord
197  TYPE(bief_obj), TARGET :: ddxprb
198 !brief REAL PART OF GRADIENT Y COMPONENT OF INCIDENT POTENTIAL AT THE BOUNDARY
199 ! partie reele du gradient en Y du potentiel incident au bord
200  TYPE(bief_obj), TARGET :: ddyprb
201 !brief IMMAGINARY PART OF GRADIENT X COMPONENT OF INCIDENT POTENTIAL AT THE BOUNDARY
202 ! partie imaginaire du gradient en X du potentiel incident au bord
203  TYPE(bief_obj), TARGET :: ddxpib
204 !brief IMMAGINARY PART OF GRADIENT Y COMPONENT OF INCIDENT POTENTIAL AT THE BOUNDARY
205 ! partie imaginaire du gradient en Y du potentiel incident au bord
206  TYPE(bief_obj), TARGET :: ddypib
207 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
208 ! coefficient pour les conditions aux limites
209  TYPE(bief_obj), TARGET :: aphi1b
210 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
211 ! coefficient pour les conditions aux limites
212  TYPE(bief_obj), TARGET :: bphi1b
213 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
214 ! coefficient pour les conditions aux limites
215  TYPE(bief_obj), TARGET :: cphi1b
216 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
217 ! coefficient pour les conditions aux limites
218  TYPE(bief_obj), TARGET :: dphi1b
219 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
220 ! coefficient pour les conditions aux limites
221  TYPE(bief_obj), TARGET :: cgrx1b
222  TYPE(bief_obj), TARGET :: cgry1b
223  TYPE(bief_obj), TARGET :: dgrx1b
224  TYPE(bief_obj), TARGET :: dgry1b
225 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
226 ! coefficient pour les conditions aux limites
227  TYPE(bief_obj), TARGET :: aphi2b
228 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
229 ! coefficient pour les conditions aux limites
230  TYPE(bief_obj), TARGET :: bphi2b
231 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
232 ! coefficient pour les conditions aux limites
233  TYPE(bief_obj), TARGET :: cphi2b
234 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
235 ! coefficient pour les conditions aux limites
236  TYPE(bief_obj), TARGET :: dphi2b
237 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
238 ! coefficient pour les conditions aux limites
239  TYPE(bief_obj), TARGET :: aphi3b
240 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
241 ! coefficient pour les conditions aux limites
242  TYPE(bief_obj), TARGET :: bphi3b
243 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
244 ! coefficient pour les conditions aux limites
245  TYPE(bief_obj), TARGET :: cphi3b
246 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
247 ! coefficient pour les conditions aux limites
248  TYPE(bief_obj), TARGET :: dphi3b
249 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
250 ! coefficient pour les conditions aux limites
251  TYPE(bief_obj), TARGET :: aphi4b
252 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
253 ! coefficient pour les conditions aux limites
254  TYPE(bief_obj), TARGET :: bphi4b
255 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
256 ! coefficient pour les conditions aux limites
257  TYPE(bief_obj), TARGET :: cphi4b
258 !brief COEFFICIENT FOR BOUNDARY CONDITIONS
259 ! coefficient pour les conditions aux limites
260  TYPE(bief_obj), TARGET :: dphi4b
261 !brief STORAGE OF ANGLE OF WAVE ATTACK (FROM NORMAL AXIS) FOR LECLIM
262 ! angle d'attaque de la houle sur les limites - pas seulement les parois (compte par rapport a a la normale , inclut dans [0;90])
263  TYPE(bief_obj), TARGET :: tetaps
264 !brief STORAGE OF DEPHASING CAUSED BY THE WALLS FOR LECLIM
265 ! dephasage induit par la paroi entre l'onde reflechie et l'onde incidente (si alfap est positif, l'onde reflechie est en retard)
266  TYPE(bief_obj), TARGET :: alfaps
267 !brief STORAGE OF INCIDENT WAVE HEIGHT AT THE BOUNDARY FOR LECLIM
268 ! hauteur de la houle aux frontieres ouvertes
269  TYPE(bief_obj), TARGET :: hbs
270 !brief ANGLE OF WAVE ATTACK - IN LOOP STORAGE (FROM NORMAL AXIS)
271 ! angle d'attaque de la houle sur les limites - pas seulement les parois (compte par rapport a a la normale , inclut dans [0;90])
272  TYPE(bief_obj), TARGET :: tetapm
273 !
274 !brief WORKING ARRAY
275 ! tableau de travail
276  TYPE(bief_obj), TARGET :: w1
277 !brief INTEGER WORKING ARRAY
278 !
279  TYPE(bief_obj), TARGET :: it1
280 !brief INTEGER WORKING ARRAY
281 !
282  TYPE(bief_obj), TARGET :: it2
283 !brief INTEGER WORKING ARRAY
284 !
285  TYPE(bief_obj), TARGET :: it3
286 !brief VOID STRUCTURE
287 !
288  TYPE(bief_obj), TARGET :: sbid
289 !brief RIGHT MEMBER OF SYSTEM TO BE SOLVED
290 !
291  TYPE(bief_obj), TARGET :: cv1
292 !brief RIGHT MEMBER OF SYSTEM TO BE SOLVED
293 !
294  TYPE(bief_obj), TARGET :: cv2
295 !brief WAVE DISSIPATION QUANTITY
296 ! tableau pour la dissipation
297  TYPE(bief_obj), TARGET :: mu
298 !brief WAVE DISSIPATION QUANTITY
299 ! tableau pour la dissipation
300  TYPE(bief_obj), TARGET :: mu2
301 !brief WAVE DISSIPATION QUANTITY
302 ! tableau pour la dissipation
303  TYPE(bief_obj), TARGET :: qb
304 !brief WAVE DISSIPATION QUANTITY
305 ! tableau pour la dissipation
306  TYPE(bief_obj), TARGET :: hmu
307 !brief WAVE DISSIPATION QUANTITY
308 ! tableau pour la dissipation
309  TYPE(bief_obj), TARGET :: hmuanc
310 !brief RADIATION STRESSES QUANTITY
311 ! tableau pour les contraintes de radiation
312  TYPE(bief_obj), TARGET :: sxx
313 !brief RADIATION STRESSES QUANTITY
314 ! tableau pour les contraintes de radiation
315  TYPE(bief_obj), TARGET :: sxy
316 !brief RADIATION STRESSES QUANTITY
317 ! tableau pour les contraintes de radiation
318  TYPE(bief_obj), TARGET :: syy
319 !brief RADIATION STRESSES QUANTITY
320 ! tableau pour les contraintes de radiation
321  TYPE(bief_obj), TARGET :: fx
322 !brief RADIATION STRESSES QUANTITY
323 ! tableau pour les contraintes de radiation
324  TYPE(bief_obj), TARGET :: fy
325 !brief MEAN WAVE PERIOD
326 ! periode moyenne issue du moment d'ordre 1
327  TYPE(bief_obj), TARGET :: t01
328 !brief MEAN WAVE PERIOD
329 ! periode moyenne issue du moment d'ordre 2
330  TYPE(bief_obj), TARGET :: t02
331 !brief MEAN WAVE PERIOD
332 ! periode moyenne issue du moment d'ordre 1
333  TYPE(bief_obj), TARGET :: tm
334 !brief WAVE DIRECTIONS AT THE BOUNDARY (RANDOM MODE)
335 !
336  TYPE(bief_obj), TARGET :: dale
337 !brief WAVE DIRECTIONS AT THE BOUNDARY (RANDOM MODE)
338 ! these vary spatially along the boundary
339  TYPE(bief_obj), TARGET :: bdale
340 !brief PERIODS ASSOCIATED TO WAVE DIRECTIONS AT THE BOUNDARY (RANDOM MODE)
341 !
342  TYPE(bief_obj), TARGET :: pdale
343 !brief BOUNDARY CONDITION TYPE
344 ! type de conditions aux limites sur u
345  TYPE(bief_obj), TARGET :: liubor
346 !brief BOUNDARY CONDITION TYPE
347 ! type de conditions aux limites sur v
348  TYPE(bief_obj), TARGET :: livbor
349 !brief BOUNDARY CONDITION TYPE
350 ! type de conditions aux limites sur h
351  TYPE(bief_obj), TARGET :: lihbor
352 !brief
353 !
354 
355 ! 'COLOUR' OF BOUNDARY NODES (TAKEN IN BOUNDARY CONDITIONS FILE)
356  TYPE(bief_obj), TARGET :: boundary_colour
357 !brief
358 ! 'BIDON INTEGER'
359  TYPE(bief_obj), TARGET :: itb1
360 !brief
361 ! 'BIDON REEL
362  TYPE(bief_obj), TARGET :: tb1
363 !brief
364 
366  TYPE(bief_obj), TARGET :: numliq
367 !brief
368 !
369  TYPE(bief_obj), TARGET :: lidir
370 !brief
371 !
372  TYPE(bief_obj), TARGET :: maskel
373 !brief MASKS FOR BOUNDARY NODES, CORRESPONDS TO INCIDENT WAVES (KINC)
374 ! masque pour les points de bord
375  TYPE(bief_obj), TARGET :: mask1
376 !brief MASKS FOR BOUNDARY NODES, CORRESPONDS TO FREE EXIT (KSORT)
377 ! masque pour les points de bord
378  TYPE(bief_obj), TARGET :: mask2
379 !brief MASKS FOR BOUNDARY NODES, CORRESPONDS TO SOLID BOUNDARY (KLOG)
380 ! masque pour les points de bord
381  TYPE(bief_obj), TARGET :: mask3
382 !brief MASKS FOR BOUNDARY NODES, CORRESPONDS TO IMPOSED WAVES (KENT)
383 ! masque pour les points de bord
384  TYPE(bief_obj), TARGET :: mask4
385 !brief MASKS FOR BOUNDARY NODES, CORRESPONDS TO INCIDENT POTENTIAL (KPOT)
386 ! masque pour les points de bord
387  TYPE(bief_obj), TARGET :: mask5
388 
389 
390 !brief FLOW
391 ! CURRENT VELOCITY IN X IDRECTION
392  TYPE(bief_obj), TARGET :: uc
393 !brief FLOW
394 ! CURRENT VELOCITY IN Y IDRECTION
395  TYPE(bief_obj), TARGET :: vc
396 !brief RELATIVE ANGULAR FREQUENCY
397 ! RELATIVE PULSATION
398  TYPE(bief_obj), TARGET :: wr
399 !brief table for wave-current interaction
400 ! WAVE VECTOR COMPONENT X
401  TYPE(bief_obj), TARGET :: kancx
402 !brief table for wave-current interaction
403 ! WAVE VECTOR COMPONENT Y
404  TYPE(bief_obj), TARGET :: kancy
405 !brief table for mean omega on the domain (from T01)
406 ! MEAN OMEGA = 2 PI / T01
407  TYPE(bief_obj), TARGET :: omegam
408 !brief table for bottom velocity in random waves
409 ! BOTTOM VELOCITY (RANDOM SEAS)
410  TYPE(bief_obj), TARGET :: ueb
411 !
412 !-----------------------------------------------------------------------
413 !
414 ! 2) MATRICES
415 !
416 !-----------------------------------------------------------------------
417 !
418 !brief MATRICE FOR SYSTEM SOLVING
419 !
420  TYPE(bief_obj), TARGET :: am1
421 !brief MATRICE FOR SYSTEM SOLVING
422 !
423  TYPE(bief_obj), TARGET :: am2
424 !brief MATRICE FOR SYSTEM SOLVING
425 !
426  TYPE(bief_obj), TARGET :: am3
427 !brief MATRICE FOR SYSTEM SOLVING
428 !
429  TYPE(bief_obj), TARGET :: bm1
430 !brief MATRICE FOR SYSTEM SOLVING
431 !
432  TYPE(bief_obj), TARGET :: bm2
433 !brief MATRICE FOR SYSTEM SOLVING
434 !
435  TYPE(bief_obj), TARGET :: mbor
436 !
437 !-----------------------------------------------------------------------
438 !
439 ! 3) BLOCKS
440 !
441 !-----------------------------------------------------------------------
442 !
443 !
444 !brief BLOCK OF POTENTIAL VECTORS
445 !
446  TYPE(bief_obj), TARGET :: phib
447 !brief BLOCK OF WORKING ARRAYS
448 !
449  TYPE(bief_obj), TARGET :: tb
450 !brief BLOCK OF WORKING ARRAYS
451 !
452  TYPE(bief_obj), TARGET :: tbbd
453 !brief BLOCK OF PRIVATE VECTORS
454 ! tableaux reserves a l'utilisateur
455  TYPE(bief_obj), TARGET :: prive
456 !brief BLOCK OF MATRICES
457 !
458  TYPE(bief_obj), TARGET :: mat
459 !brief BLOCK OF UNKNOWN VECTORS
460 !
461  TYPE(bief_obj), TARGET :: unk
462 !brief BLOCK OF RIGHT HAND SIDE VECTORS IN SOLVING SYSTEM
463 !
464  TYPE(bief_obj), TARGET :: rhs
465 !brief BLOCK OF VARIABLES FOR OUTPUT
466 !
467  TYPE(bief_obj), TARGET :: varsor
468 !brief Number of variables in varsor
469  INTEGER :: nvar_art
470 !brief BLOCK OF VARIABLES FOR ANIMATION OUTPUT
471 !
472  TYPE(bief_obj), TARGET :: varnim
473 !
474 !-----------------------------------------------------------------------
475 !
476 ! 4) INTEGERS
477 !
478 !-----------------------------------------------------------------------
479 !
480 !brief
481 ! maximum de variables de sortie
482  INTEGER, PARAMETER :: maxvar = 100
483 !brief
484 ! maximum de frontieres liquides
485  INTEGER, PARAMETER :: maxfro = 3000
486 !brief OPEN BOUNDARY WITH INCIDENT POTENTIAL
487 !
488  INTEGER, PARAMETER :: kpot = 7
489 !brief DEBUGGER
490 !
491  INTEGER, TARGET :: debug
492 !brief ORIGIN COORDINATE
493 ! coordonnee de l'origine
494  INTEGER i_orig
495 !brief ORIGIN COORDINATE
496 ! coordonnee de l'origine
497  INTEGER j_orig
498 !brief GRAPHIC PRINTOUT PERIOD
499 ! periode de sortie graphique
500  INTEGER leoprd
501 !brief LISTING PRINTOUT PERIOD
502 ! periode de sortie listing
503  INTEGER lisprd
504 !brief MAXIMUM NUMBER OF ITERATIONS FOR SOLVER
505 ! maximum d'iterations pour le solveur
506  INTEGER nitmax
507 !brief GEOMETRY FILE STANDARD
508 ! standard du fichier de geometrie
509 ! INTEGER STDGEO
510 !brief RESULTS FILE STANDARD
511 ! standard du fichier des resultats
512 ! INTEGER STDRES
513 !brief SOLVER OPTION
514 ! option du solveur
515  INTEGER isolve(2)
516 !brief BOTTOM TOPOGRAPHY SMOOTHINGS
517 ! nombre de lissages du fond
518  INTEGER lisfon
519 !brief DISCRETISATION IN SPACE
520 ! discretisation en espace
521 ! INTEGER DISESP
522 !brief NUMBER OF DISCRETISED PERIODS
523 ! nombre de periodes de discretisation du spectre de houle
524  INTEGER npale
525 !brief NUMBER OF DISCRETISED DIRECTIONS
526 ! nombre de directions de discretisation du spectre de houle
527  INTEGER ndale
528 !brief MATRIX STORAGE
529 ! stockage des matrices
530  INTEGER optass
531 !brief BREAKING LAW
532 ! formulation du deferlement
533  INTEGER ibreak
534 !brief MAXIMUM OF SUB-ITERATIONS
535 ! maximum de sous-iterations
536  INTEGER nitdis
537 !brief VECTOR LENGTH
538 ! longueur du vecteur
539  INTEGER lvmac
540 !brief LAW OF BOTTOM FRICTION
541 ! loi de frottement sur le fond
542 ! INTEGER KFROT
543 !brief BOTTOM FRICTION LAW
544 ! formulation du frottement de fond
545  INTEGER formfr
546 !brief HYDRAULIC REGIME TYPE
547 ! type du regime hydraulique
548  INTEGER regido
549 !brief MATRIX-VECTOR PRODUCT
550 ! produit matrice-vecteur
551  INTEGER produc
552 !
553 !brief NUMBER OF PRIVATE ARRAYS, NUMBER OF PRIVATE ARRAYS WITH GIVEN NAME
554 ! nombre de tableaux prives et tableaux prives avec noms
555  INTEGER npriv,n_names_priv
556 !brief
557 !
558  INTEGER ptinig
559 !brief
560 !
561  INTEGER ptinil
562 !brief
563 ! type d'element
564  INTEGER ielm
565 !brief
566 !
567  INTEGER ielm0
568 !brief
569 ! type d'element de bord
570  INTEGER ielmb
571 !brief
572 !
573  INTEGER ielmb0
574 !brief ORIGINAL DATE OF TIME
575 ! date de l'origine des temps
576  INTEGER mardat(3)
577 !brief ORIGINAL HOUR OF TIME
578 ! heure de l'origine des temps
579  INTEGER martim(3)
580 !brief
581 !
582  INTEGER nfrliq
583 !brief
584 ! prise en compte des effets de pente/courbure
585  INTEGER ipentco
586 !brief
587  INTEGER lper
588 !brief
589  INTEGER ldir
590 !brief MAX NUMBER OF ITERATION ON TETAP CALCULATION
591 ! nombre max d'iteration pour calcul auto du TETAP
592  INTEGER nittp
593 !brief NESTING WITHIN TOMAWAC OUTER MODEL
594 ! si different de 0, ARTEMIS est integre dans un modele TOMAWAC global
595  INTEGER chaintwc
596 !brief
597 ! nombre de spectres
598  INTEGER nspec
599 !brief
600 ! nombre de directions definissant le(s) spectre(s)
601  INTEGER ndir
602 !brief
603 ! nombre de frequences definissant le(s) spectre(s)
604  INTEGER nf
605 !brief
606 ! numero de reference TOMAWAC pour le spectre en freq
607  INTEGER n_sfref
608 !brief
609 ! noeud frontiere de reference pour le spectre en freq
610  INTEGER iptfr_ref
611 !
612 !-----------------------------------------------------------------------
613 !
614 ! 5) LOGICAL VALUES
615 !
616 !-----------------------------------------------------------------------
617 !
618 !brief LISTING PRINTOUT
619 ! si oui, sortie listing
620  LOGICAL listin
621 !brief
622 !
623  LOGICAL infogr
624 !brief PERIOD SCANNING
625 ! si oui, balayage en periodes
626  LOGICAL balaye
627 !brief MONODIRECTIONAL RANDOM WAVE
628 ! si oui, houle aleatoire monodirectionnelle
629  LOGICAL alemon
630 !brief MULTIDIRECTIONAL RANDOM WAVE
631 ! si oui, houle aleatoire multidirectionnelle
632  LOGICAL alemul
633 !brief
634 !
635  LOGICAL msk
636 !brief
637 !
638  LOGICAL spheri
639 !brief BREAKING
640 ! si oui, deferlement
641  LOGICAL deferl
642 !brief FRICTION
643 ! si oui, frottement
644  LOGICAL frotte
645 !brief FRICTION FACTOR IMPOSED
646 ! si oui, facteur de frottement impose
647  LOGICAL entfw
648 !brief HYDRAULIC REGIME IMPOSED
649 ! si oui, regime hydraulique impose
650  LOGICAL entreg
651 !brief SKIN ROUGHNESS ONLY
652 ! si oui, rugosite de peau seule
653  LOGICAL entrug
654 !brief WAVE HEIGHTS SMOOTHING
655 ! si oui, lissage des hauteurs de houle
656  LOGICAL lishou
657 !brief
658 !
659  LOGICAL sorleo(maxvar)
660 !brief
661 !
662  LOGICAL sorimp(maxvar)
663 !brief
664 ! a logical array holding variable output for phase and amplitude file
665  LOGICAL sornim(maxvar)
666 !brief VALIDATION
667 ! si oui, validation
668  LOGICAL valid
669 !brief COURANT
670 ! Yes = current taken into account (defined in condih.f)
671  LOGICAL courant
672 !!brief AUTOMATIC ANGLES
673 ! si oui,calcul automatique des angles de sortie TETAP
674  LOGICAL langauto
675 !!brief AUTOMATIC PHASES ON INCIDENT WAVE BOUNDARY
676 ! si oui, calcul automatique des phases sur la frontière de type onde incidente (KINC)
677  LOGICAL lphaseauto
678 !brief ANIMATION
679 ! si oui, generates the phase and amplitude file
680  LOGICAL animfs
681 !
682 !-----------------------------------------------------------------------
683 !
684 ! 6) REALS
685 !
686 !-----------------------------------------------------------------------
687 !
688 !brief GRAVITY ACCELERATION
689 ! acceleration de la pesanteur
690  DOUBLE PRECISION grav
691 !brief MINIMUM VALUE FOR H
692 ! valeur minimum de h
693 ! DOUBLE PRECISION HMIN
694 !brief WAVE PERIOD
695 ! periode de la houle en cours de calcul
696  DOUBLE PRECISION per
697 !brief ANGULAR FREQUENCY
698 ! pulsation de la houle
699  DOUBLE PRECISION omega
700 !brief DIRECTION OF WAVE PROPAGATION
701 ! direction principale de propagation de la houle
702  DOUBLE PRECISION tetah
703 !brief INITIAL WATER LEVEL
704 ! cote initiale
705  DOUBLE PRECISION cotini
706 !brief INITIAL DEPTH
707 ! hauteur initiale
708  DOUBLE PRECISION hautin
709 !brief BEGINNING PERIOD FOR PERIOD SCANNING
710 ! periode de debut pour le balayage en periode
711  DOUBLE PRECISION perdeb
712 !brief ENDING PERIOD FOR PERIOD SCANNING
713 ! periode de fin pour le balayage en periode
714  DOUBLE PRECISION perfin
715 !brief STEP FOR PERIOD SCANNING
716 ! pas pour le balayage en periode
717  DOUBLE PRECISION perpas
718 !brief PEAK PERIOD
719 ! periode de pic
720  DOUBLE PRECISION perpic
721 !brief GAMMA
722 ! gamma
723  DOUBLE PRECISION gamma
724 !brief MINIMUM ANGLE OF PROPAGATION
725 ! valeur minimum de l'angle de propagation
726  DOUBLE PRECISION tetmin
727 !brief MAXIMUM ANGLE OF PROPAGATION
728 ! valeur maximum de l'angle de propagation
729  DOUBLE PRECISION tetmax
730 !brief S EXPONENT
731 ! exposant s dans la formule du spectre
732  DOUBLE PRECISION expos
733 !brief
734 !
735 ! DOUBLE PRECISION RELAX
736 !brief FRICTION COEFFICIENT
737 ! coefficient de frottement
738  DOUBLE PRECISION ffon
739 !brief SUB-ITERATIONS ACCURACY
740 ! precision sur les sous-iterations
741  DOUBLE PRECISION epsdis
742 !brief DISSIPATION RELAXATION
743 ! relaxation sur la dissipation
744  DOUBLE PRECISION reldis
745 !brief ALPHA
746 ! alpha
747  DOUBLE PRECISION alfabj
748 !brief GAMMAS
749 ! gammas
750  DOUBLE PRECISION gammas
751 !brief
752 !
753  DOUBLE PRECISION kdally
754 !brief
755 !
756  DOUBLE PRECISION gdally
757 !brief FLUID KINEMATIC VISCOSITY
758 ! viscosite cinematique du fluide
759  DOUBLE PRECISION visco
760 !brief DIAMETER90
761 ! diametre90
762  DOUBLE PRECISION diam90
763 !brief DIAMETER50
764 ! diametre50
765  DOUBLE PRECISION diam50
766 !brief SEDIMENT SPECIFIC WEIGHT
767 ! masse volumique du sediment
768  DOUBLE PRECISION mvsed
769 !brief FLUID SPECIFIC MASS
770 ! masse volumique du fluide
771  DOUBLE PRECISION mveau
772 !brief FRICTION FACTOR
773 ! coefficient de frottement constant impose
774  DOUBLE PRECISION fwcoef
775 !brief RIPPLES COEFFICIENT
776 ! coefficient de rides
777  DOUBLE PRECISION ricoef
778 !brief MINIMUM SPECTRAL PERIOD
779 ! periode minimum du spectre
780  DOUBLE PRECISION pmin
781 !brief MAXIMUM SPECTRAL PERIOD
782 ! periode maximum du spectre
783  DOUBLE PRECISION pmax
784 !brief REFERENCE WATER DEPTH FOR AUTOMATIC PHASE
785 ! Profondeur d'eau de référence pour le calcul automatique des phases
786  DOUBLE PRECISION depref
787 !brief REFERENCE WAVE NUMBER FOR AUTOMATIC PHASE
788 ! Nombre d'onde de référence pour le calcul automatique des phases
789  DOUBLE PRECISION kphref
790 !brief X COORDINATE OF REFERENCE POINT FOR PHASE DEFINITION
791 ! coordonnee en x du point de reference pour les phases
792  DOUBLE PRECISION x_phref
793 !brief Y COORDINATE OF REFERENCE POINT FOR PHASE DEFINITION
794 ! coordonnee en y du point de reference pour les phases
795  DOUBLE PRECISION y_phref
796 !brief SUB-ITERATIONS ACCURACY FOR CURRENT(WAVE VECTOR DIRECTION)
797 ! precision sur les sous-iterations pour le courant (direction vecteur d'onde)
798  DOUBLE PRECISION epsdir
799 !brief SUB-ITERATIONS ACCURACY FOR TETAP
800 ! precision sur les sous-iterations sur le TETAP
801  DOUBLE PRECISION epstp
802 !brief RELAXATION COEFFICIENT FOR TETAP
803 ! coefficient de relaxation pour calcul automatic de TETAP
804  DOUBLE PRECISION reltp
805 !brief TIME AT WHICH TOMAWAC SPECTRUM IS TAKEN
806 !
807  DOUBLE PRECISION tpstwc
808 !brief Significant wave height corresponding to global
809 ! energy contained in TOMAWAC spectrum
810  DOUBLE PRECISION hscal
811 !brief ABSCISSAE FOR THE REFERENCE F SPECTRUM
812 ! x du point de reference pour le spectre en freq
813  DOUBLE PRECISION x_sfref
814 !brief ORDINATES FOR THE REFERENCE F SPECTRUM
815 ! y du point de reference pour le spectre en freq
816  DOUBLE PRECISION y_sfref
817 !
818  DOUBLE PRECISION fp,gam,delta
819  DOUBLE PRECISION expo
820 !
821 !brief ARTEMIS CONSTANTS (INITIALISED INTO ARTEMIS_CONSTANTS)
822 !
823  DOUBLE PRECISION :: pi,deupi,pisur2,usdpi,raddeg,degrad
824 !
825 !-----------------------------------------------------------------------
826 !
827 ! 7) STRINGS
828 !
829 !-----------------------------------------------------------------------
830 !
831 !brief TITLE
832 ! titre de l'etude
833  CHARACTER(LEN=72) titcas
834 !brief
835 !
836  CHARACTER(LEN=72) vardes
837 !brief VARIABLES TO BE PRINTED
838 ! variables a imprimer
839  CHARACTER(LEN=72) varimp
840 !brief INITIAL CONDITIONS
841 ! conditions initiales
842  CHARACTER(LEN=72) cdtini
843 !brief GEOMETRY FILE BINARY
844 ! binaire du fichier de geometrie
845 ! CHARACTER(LEN=3) BINGEO
846 !brief RESULTS FILE BINARY
847 ! binaire du fichier des resultats
848 ! CHARACTER(LEN=3) BINRES
849 !brief
850 !
851  CHARACTER(LEN=20) equa
852 !brief
853 !
854  CHARACTER(LEN=32) varcla(10)
855 !brief
856 !
857  CHARACTER(LEN=32) texte(maxvar)
858 !brief
859 !
860  CHARACTER(LEN=32) textpr(maxvar)
861 !brief
862 ! name of variables in phase and amplitude file
863  CHARACTER(LEN=32) textanim(maxvar)
864 !
865 ! NAMES OF PRIVATE ARRAYS (GIVEN BY USER)
866 !
867  CHARACTER(LEN=32) names_prive(4)
868 !
869 !-----------------------------------------------------------------------
870 !
871 ! 8) SLVCFG STRUCTURES
872 !
873 !-----------------------------------------------------------------------
874 !
875 !brief SLVCFG STRUCTURE
876 !
877  TYPE(slvcfg) :: slvart
878 !
879 !-----------------------------------------------------------------------
880 !
881 ! 9) MESH STRUCTURE
882 !
883 !-----------------------------------------------------------------------
884 !
885 !brief MESH STRUCTURE
886 !
887  TYPE(bief_mesh),TARGET :: mesh
888 !
889 !-----------------------------------------------------------------------
890 !
891 ! 10) ALIASES
892 !
893 !-----------------------------------------------------------------------
894 !
895 ! DECLARATION OF POINTERS FOR ALIASES.
896 ! TARGETS ARE DEFINED IN POINT_ARTEMIS
897 !
898 ! ALIASES FOR WORKING VECTORS IN TB AND TBBD
899 !
900 !brief WORKING VECTOR IN TB
901 !
902  TYPE(bief_obj),POINTER :: t1
903 !brief WORKING VECTOR IN TB
904 !
905  TYPE(bief_obj),POINTER :: t2
906 !brief WORKING VECTOR IN TB
907 !
908  TYPE(bief_obj),POINTER :: t3
909 !brief WORKING VECTOR IN TB
910 !
911  TYPE(bief_obj),POINTER :: t4
912 !brief WORKING VECTOR IN TB
913 !
914  TYPE(bief_obj),POINTER :: t5
915 !brief WORKING VECTOR IN TB
916 !
917  TYPE(bief_obj),POINTER :: t6
918 !brief WORKING VECTOR IN TB
919 !
920  TYPE(bief_obj),POINTER :: t7
921 !brief WORKING VECTOR IN TB
922 !
923  TYPE(bief_obj),POINTER :: t8
924 !brief WORKING VECTOR IN TB
925 !
926  TYPE(bief_obj),POINTER :: t9
927 !brief WORKING VECTOR IN TB
928 !
929  TYPE(bief_obj),POINTER :: t10
930 !brief WORKING VECTOR IN TB
931 !
932  TYPE(bief_obj),POINTER :: t11
933 !brief WORKING VECTOR IN TB
934 !
935  TYPE(bief_obj),POINTER :: t12
936 !brief WORKING VECTOR IN TBBD
937 !
938  TYPE(bief_obj),POINTER :: t13
939 !brief WORKING VECTOR IN TB
940 !
941  TYPE(bief_obj),POINTER :: t14
942 !brief WORKING VECTOR IN TB
943 !
944  TYPE(bief_obj),POINTER :: t15
945 !brief WORKING VECTOR IN TB
946 !
947  TYPE(bief_obj),POINTER :: t16
948 !brief WORKING VECTOR IN TB
949 !
950  TYPE(bief_obj),POINTER :: tbd1
951 !brief WORKING VECTOR IN TBBD
952 !
953  TYPE(bief_obj),POINTER :: tbd2
954 !brief WORKING VECTOR IN TBBD
955 !
956  TYPE(bief_obj),POINTER :: tbd3
957 !brief WORKING VECTOR IN TBBD
958 !
959  TYPE(bief_obj),POINTER :: tbd4
960 !
961 ! USEFUL COMPONENTS IN STRUCTURE MESH
962 !
963 !brief
964 ! table de connectivite
965  TYPE(bief_obj), POINTER :: ikle
966 !brief
967 ! coordonnees des points du maillage
968  DOUBLE PRECISION, DIMENSION(:), POINTER :: x
969 !brief
970 ! coordonnees des points du maillage
971  DOUBLE PRECISION, DIMENSION(:), POINTER :: y
972 !brief
973 ! nombre d'elements du maillage
974  INTEGER, POINTER :: nelem
975 !brief
976 !
977  INTEGER, POINTER :: nelmax
978 !brief
979 ! nombre de points frontiere
980  INTEGER, POINTER :: nptfr
981 !brief
982 !
983  INTEGER, POINTER :: nptfrx
984 !brief
985 !
986  INTEGER, POINTER :: dim1
987 !brief
988 !
989  INTEGER, POINTER :: typelm
990 !brief
991 ! nombre de points du maillage
992  INTEGER, POINTER :: npoin
993 !brief
994 !
995  INTEGER, POINTER :: npmax
996 !brief
997 !
998  INTEGER, POINTER :: mxptvs
999 !brief
1001  INTEGER, POINTER :: mxelvs
1002 !brief
1004  INTEGER, POINTER :: lv
1005 !
1006 !-----------------------------------------------------------------------
1007 !
1008 ! 11) ART_FILES AND ASSOCIATED
1009 !
1010 !-----------------------------------------------------------------------
1011 !
1012 !brief
1014  INTEGER, PARAMETER :: maxlu_art = 46
1015 !brief NAME OF THE GEOMETRY FILE
1016 ! nom du fichier de geometrie
1017  INTEGER, TARGET :: artgeo
1018 !brief NAME OF THE STEERING FILE
1019 ! nom du fichier des parametres
1020  INTEGER :: artcas
1021 !brief NAME OF THE BOUNDARY CONDITIONS FILE
1022 ! nom du fichier des conditions aux limites
1023  INTEGER, TARGET :: artcli
1024 !brief NAME OF THE BOTTOM TOPOGRAPHY FILE
1025 ! nom du fichier des fonds
1026  INTEGER :: artfon
1027 !brief NAME OF THE RESULTS FILE
1028 ! nom du fichier des resultats
1029  INTEGER, TARGET :: artres
1030 !brief NAME OF THE BINARY RESULTS FILE
1031 ! nom du fichier des resultats binaire
1032  INTEGER :: artrbi
1033 !brief NAME OF THE FORMATTED RESULTS FILE
1034 ! nom du fichier des resultats formate
1035  INTEGER :: artrfo
1036 !brief NAME OF THE AMPLITUDE AND PHASE FILE
1037 ! nom du fichier contenant phases et amplitudes
1038  INTEGER :: artamp
1039 !brief NAME OF THE REFERENCE FILE
1040 ! nom du fichier de reference
1041  INTEGER :: artref
1042 !brief NAME OF THE BINARY DATA FILE 1
1043 ! nom du fichier de donnees binaire 1
1044  INTEGER :: artbi1
1045 !brief NAME OF THE BINARY DATA FILE 2
1046 ! nom du fichier de donnees binaire 2
1047  INTEGER :: artbi2
1048 !brief NAME OF THE FORMATTED DATA FILE 1
1049 ! nom du fichier de donnees formate 1
1050  INTEGER :: artfo1
1051 !brief NAME OF THE FORMATTED DATA FILE 2
1052 ! nom du fichier de donnees formate 2
1053  INTEGER :: artfo2
1054 !brief NAME OF THE TOMAWAC OUTER SPECTRAL FILE
1055 ! nom du fichier de spectre global
1056  INTEGER :: wacspe
1057 !brief NAME OF THE TOMAWAC OUTER RESULT FILE
1058 ! nom du fichier de resultats global
1059  INTEGER :: wacres
1060 !brief NAME OF THE TOMAWAC LIQUID BOUNDARY FILE
1061 ! nom du fichier de frontieres liquides global
1062  INTEGER :: waclqd
1063 !brief
1065  TYPE(bief_file), TARGET :: art_files(maxlu_art)
1066 !
1067  SAVE
1068 !
1069  END MODULE declarations_artemis
type(bief_obj), target am2
type(bief_obj), target unk
type(bief_obj), target hale
type(bief_obj), target lihbor
type(bief_obj), pointer tbd1
type(bief_obj), target prb
type(bief_obj), target am1
logical, dimension(maxvar) sorimp
type(bief_obj), target mat
character(len=32), dimension(maxvar) texte
type(bief_obj), target hmuanc
integer, parameter maxfro
type(bief_obj), target h
type(bief_obj), target phiib
type(bief_obj), target it2
type(bief_obj), pointer t14
type(bief_obj), pointer t6
type(bief_obj), target w1
type(bief_obj), target tbbd
type(bief_obj), target mask3
type(bief_obj), target mask4
logical, dimension(maxvar) sornim
type(bief_obj), target aphi2b
type(bief_obj), target dgrx1b
type(bief_obj), target mask1
type(bief_obj), target bphi1b
type(bief_obj), target omegam
type(bief_obj), target tetap
type(bief_obj), target dphi2b
type(bief_obj), target zf
type(bief_obj), target mbor
type(bief_obj), pointer t11
integer, dimension(3) martim
type(bief_obj), target msin
type(bief_obj), target hbs
type(bief_obj), target livbor
integer, dimension(2) isolve
integer, parameter maxlu_art
type(bief_obj), target maskel
type(bief_obj), pointer tbd2
type(bief_obj), target kancx
double precision, dimension(:), pointer y
type(bief_obj), target aphi3b
type(bief_obj), target tb
type(bief_obj), target cgry1b
type(bief_obj), target mask5
type(bief_obj), target aphi1b
type(bief_obj), target hho
type(bief_obj), target phii
type(bief_obj), target it3
type(bief_obj), target ueb
type(bief_obj), pointer t5
type(bief_obj), pointer t7
type(bief_obj), target cphi1b
type(bief_obj), target itb1
type(bief_obj), target t01
type(bief_obj), target phib
type(bief_obj), target v0
type(bief_obj), pointer t15
type(bief_obj), pointer t9
type(bief_obj), target dale
type(bief_obj), target pib
type(bief_obj), target bdale
character(len=32), dimension(10) varcla
type(bief_obj), target tetapm
type(bief_obj), target cv2
type(bief_obj), target varnim
integer, parameter maxvar
integer, parameter kpot
integer, dimension(3) mardat
type(bief_obj), pointer t2
logical, dimension(maxvar) sorleo
type(bief_obj), target k
type(bief_obj), target bm1
type(bief_obj), target it1
type(bief_obj), target phirb
type(bief_mesh), target mesh
character(len=32), dimension(maxvar) textanim
type(bief_obj), target t02
type(bief_obj), target mask2
type(bief_obj), target hb
type(bief_obj), target prive
type(bief_obj), target wr
type(bief_obj), target tetab
type(bief_obj), target aphi4b
type(bief_obj), target dphi1b
type(bief_obj), target ddxpib
type(bief_obj), target fy
type(bief_obj), target pale
type(bief_obj), target bphi2b
type(bief_obj), target fw
type(bief_obj), target boundary_colour
character(len=32), dimension(maxvar) textpr
type(bief_obj), target inci
type(bief_obj), pointer t8
type(bief_obj), target mcos
type(bief_obj), target alfaps
type(bief_obj), target dphi4b
type(bief_obj), pointer t10
type(bief_obj), target alfap
type(bief_obj), target hmu
type(bief_obj), pointer t4
type(bief_file), dimension(maxlu_art), target art_files
type(bief_obj), target numliq
type(bief_obj), target phir
type(bief_obj), target uc
type(bief_obj), target vc
type(bief_obj), target phas
double precision, dimension(:), pointer x
type(bief_obj), target rp
type(bief_obj), target ddxprb
type(bief_obj), target rhs
type(bief_obj), pointer tbd4
type(bief_obj), target sbid
character(len=72) titcas
type(bief_obj), target liubor
type(bief_obj), pointer ikle
type(bief_obj), target cg
type(bief_obj), pointer t12
character(len=32), dimension(4) names_prive
type(bief_obj), target qb
type(bief_obj), target cgrx1b
type(bief_obj), target cv1
type(bief_obj), pointer t3
type(bief_obj), target bm2
type(bief_obj), target kancy
type(bief_obj), target mu
type(bief_obj), target pdale
type(bief_obj), target am3
type(bief_obj), target bphi4b
type(bief_obj), target cphi2b
type(bief_obj), target cphi4b
type(bief_obj), pointer t1
type(bief_obj), target mu2
type(bief_obj), target ddyprb
type(bief_obj), target lidir
character(len=72) vardes
type(bief_obj), target c
type(bief_obj), target varsor
type(bief_obj), target sxx
type(bief_obj), target sxy
type(bief_obj), target tb1
type(bief_obj), target s
type(bief_obj), pointer t13
type(bief_obj), target dgry1b
character(len=72) cdtini
type(bief_obj), target u0
type(bief_obj), pointer t16
type(bief_obj), target ddypib
character(len=72) varimp
type(bief_obj), target syy
type(bief_obj), target fx
type(bief_obj), pointer tbd3
type(bief_obj), target cphi3b
type(bief_obj), target bphi3b
type(bief_obj), target tetaps
type(bief_obj), target dphi3b
type(bief_obj), target tm