user_utimp_art.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE user_utimp_art
!                   *************************
!
!***********************************************************************
! ARTEMIS   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    ALMOST ALL THE COMPUTATION VARIABLES ARE AVAILABLE
!+             HERE TO WRITE OUT SPECIFIC OUTPUT, COMPUTE ANALYTICAL
!+             SOLUTIONS...
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE interface_artemis, ex_user_utimp_art=> user_utimp_art
      USE declarations_artemis
      USE bief
!
      USE declarations_special
      IMPLICIT NONE
!
!------------------------------------------------------------------
!
      INTEGER I
      TYPE(bief_obj) :: BID
!
!------------------------------------------------------------------
!
      IF(.false.) THEN
!     EXAMPLE : U1 AND V1
!              (HORIZONTAL VELOCITIES AT T/4)
!               ARE TRANSFERRED TO PRIVE(I,1) AND PRIVE(I,2)
        CALL vector(t2, '=' , 'GRADF          X' , ielm ,
     &              1.d0 , phii , bid , bid , bid , bid , bid ,
     &              mesh , msk , maskel )

        CALL vector(t3 , '=' , 'GRADF          Y' , ielm ,
     &              1.d0 , phii , bid , bid , bid , bid , bid ,
     &              mesh , msk , maskel )

        CALL vector(t1 , '=' , 'MASBAS          ' , ielm ,
     &              1.d0 , bid , bid , bid , bid , bid , bid ,
     &              mesh , msk , maskel )

        CALL os('X=Y/Z   ', x=t2, y=t2, z=t1)
        CALL os('X=Y/Z   ', x=t3, y=t3, z=t1)

        ! TODO: Correct the computation below they are wrong
        DO i = 1,npoin
          prive%ADR(1)%P%R(i) = omegam%R(i)*t2%R(i)
          prive%ADR(2)%P%R(i) = 2d0*3.1415d0/perpic*t3%R(i)
        ENDDO
      ENDIF
!
      RETURN
      END