The TELEMAC-MASCARET system  trunk
user_utimp_art.f
Go to the documentation of this file.
1 ! *************************
2  SUBROUTINE user_utimp_art
3 ! *************************
4 !
5 !***********************************************************************
6 ! ARTEMIS V6P1 21/08/2010
7 !***********************************************************************
8 !
9 !brief ALMOST ALL THE COMPUTATION VARIABLES ARE AVAILABLE
10 !+ HERE TO WRITE OUT SPECIFIC OUTPUT, COMPUTE ANALYTICAL
11 !+ SOLUTIONS...
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
15 !
16  USE interface_artemis, ex_user_utimp_art=> user_utimp_art
18  USE bief
19 !
21  IMPLICIT NONE
22 !
23 !------------------------------------------------------------------
24 !
25  INTEGER I
26  TYPE(bief_obj) :: BID
27 !
28 !------------------------------------------------------------------
29 !
30  IF(.false.) THEN
31 ! EXAMPLE : U1 AND V1
32 ! (HORIZONTAL VELOCITIES AT T/4)
33 ! ARE TRANSFERRED TO PRIVE(I,1) AND PRIVE(I,2)
34  CALL vector(t2, '=' , 'GRADF X' , ielm ,
35  & 1.d0 , phii , bid , bid , bid , bid , bid ,
36  & mesh , msk , maskel )
37 
38  CALL vector(t3 , '=' , 'GRADF Y' , ielm ,
39  & 1.d0 , phii , bid , bid , bid , bid , bid ,
40  & mesh , msk , maskel )
41 
42  CALL vector(t1 , '=' , 'MASBAS ' , ielm ,
43  & 1.d0 , bid , bid , bid , bid , bid , bid ,
44  & mesh , msk , maskel )
45 
46  CALL os('X=Y/Z ', x=t2, y=t2, z=t1)
47  CALL os('X=Y/Z ', x=t3, y=t3, z=t1)
48 
49  ! TODO: Correct the computation below they are wrong
50  DO i = 1,npoin
51  prive%ADR(1)%P%R(i) = omegam%R(i)*t2%R(i)
52  prive%ADR(2)%P%R(i) = 2d0*3.1415d0/perpic*t3%R(i)
53  ENDDO
54  ENDIF
55 !
56  RETURN
57  END
type(bief_obj), target omegam
type(bief_obj), target maskel
double precision, dimension(:), pointer y
type(bief_obj), target phii
type(bief_obj), pointer t2
type(bief_mesh), target mesh
type(bief_obj), target prive
subroutine user_utimp_art
Definition: user_utimp_art.f:4
double precision, dimension(:), pointer x
subroutine vector(VEC, OP, FORMUL, IELM1, XMUL, F, G, H, U, V, W, MESH, MSK, MASKEL, LEGO, ASSPAR)
Definition: vector.f:7
type(bief_obj), pointer t3
subroutine os(OP, X, Y, Z, C, IOPT, INFINI, ZERO)
Definition: os.f:7
type(bief_obj), pointer t1
Definition: bief.f:3