The TELEMAC-MASCARET system  trunk
couuti.f
Go to the documentation of this file.
1 ! *****************
2  SUBROUTINE couuti
3 ! *****************
4 !
5  &(ncou,fmtcou)
6 !
7 !***********************************************************************
8 ! TOMAWAC V6P1 14/06/2011
9 !***********************************************************************
10 !
11 !brief READS IN THE CURRENTS USING A USER-DEFINED FORMAT.
12 !
13 !note DURING THE FIRST PASS THE USER MUST IDENTIFY THE TIMES TC1 AND TC2
14 !+ WHICH SURROUND THE FIRST TIME STEP. NEXT, USING THE ARRAYS
15 !+ XRELC,YRELC,UR,VR OR DIRECTLY FROM THE CURRENT FILE, THE USER
16 !+ MAY HAVE TO INTERPOLATE THE CURRENTS READ FROM THE FILE INTO THE
17 !+ ARRAYS UC1,VC1 UC2,VC2.
18 !+
19 !+ INTERPOLATION SUBROUTINE FASP :
20 !+
21 !+ CALL FASP(X,Y,UC1,NPOIN,XRELC,YRELC,UR,NP,NBOR,MESH%KP1BOR%I,NPTFR,0.D0)
22 !+
23 !+ CALL FASP(X,Y,VC1,NPOIN,XRELC,YRELC,VR,NP,NBOR,MESH%KP1BOR%I,NPTFR,0.D0)
24 !+
25 !+ THE CODE WILL INTERPOLATE THE CURRENT AUTOMATICALLY BETWEEN THESE
26 !+ 2 TIME STEPS.
27 !+
28 !+ THE OTHER PASSES OCCUR WHEN A NEW RECORD IS REQUIRED (AT>TC2).
29 !+ IN THIS CASE TC2,UC2,VC2 ONLY ARE TO BE COMPUTED.
30 !
31 !history F.MARCOS (LNH)
32 !+ 30/08/95
33 !+ V1P0
34 !+
35 !
36 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
37 !+ 13/07/2010
38 !+ V6P0
39 !+ Translation of French comments within the FORTRAN sources into
40 !+ English comments
41 !
42 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
43 !+ 21/08/2010
44 !+ V6P0
45 !+ Creation of DOXYGEN tags for automated documentation and
46 !+ cross-referencing of the FORTRAN sources
47 !
48 !history G.MATTAROLO (EDF - LNHE)
49 !+ 14/06/2011
50 !+ V6P1
51 !+ Translation of French names of the variables in argument
52 !
53 !history J-M HERVOUET (EDF-LNHE)
54 !+ 23/11/20012
55 !+ V6P3
56 !+ XRELC, YRELC, UR, VR removed, must be declared locally
57 !+ if necessary.
58 !
59 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
60 !| AT |-->| COMPUTATION TIME
61 !| FMTCOU |-->| CURRENTS FILE BINARY FORMAT
62 !| DDC |-->| DATE OF COMPUTATION BEGINNING
63 !| NBOR |-->| GLOBAL NUMBER OF BOUNDARY POINTS
64 !| NCOU |-->| LOGICAL UNIT NUMBER OF THE CURRENTS FILE
65 !| NPMAX |-->| MAXIMUM NUMBER OF POINTS THAT CAN BE READ
66 !| NPOIN |-->| NUMBER OF POINTS IN 2D MESH
67 !| NPTFR |-->| NUMBER OF BOUNDARY POINTS
68 !| TC1 |-->| TIME T1 IN THE CURRENT FILE
69 !| TC2 |-->| TIME T2 IN THE CURRENT FILE
70 !| UC1,VC1 |<->| CURRENTS VALUES AT TIME T1 IN THE MESH
71 !| UC2,VC2 |<->| CURRENTS VALUES AT TIME T2 IN THE MESH
72 !| X |-->| ABSCISSAE OF POINTS IN THE MESH
73 !| Y |-->| ORDINATES OF POINTS IN THE MESH
74 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
75 !
76  USE interface_tomawac, ex_couuti => couuti
77 ! USE DECLARATIONS_TOMAWAC, ONLY : NPTFR, NBOR
78 !
80  IMPLICIT NONE
81 !
82 !
83 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
84 !
85  INTEGER, INTENT(IN) :: NCOU
86  CHARACTER(LEN=8), INTENT(IN) :: FMTCOU
87 !
88 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
89 !
90  CALL user_current(ncou,fmtcou)
91 !
92 !-----------------------------------------------------------------------
93 !
94  RETURN
95  END
subroutine user_current(NCOU, FMTCOU)
Definition: user_current.f:7
subroutine couuti(NCOU, FMTCOU)
Definition: couuti.f:7