calcue.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\artemis\calcue.f
00002 !
00076                      SUBROUTINE CALCUE
00077 !                    *****************
00078 !
00079 !
00080 !***********************************************************************
00081 ! ARTEMIS   V6P1                                   21/08/2010
00082 !***********************************************************************
00083 !
00084 !
00085 !
00086 !
00087 !
00088 !
00089 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00090 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00091 !
00092       USE BIEF
00093       USE DECLARATIONS_TELEMAC
00094       USE DECLARATIONS_ARTEMIS
00095 !
00096       IMPLICIT NONE
00097       INTEGER LNG,LU
00098       COMMON/INFO/LNG,LU
00099 !
00100       DOUBLE PRECISION BID
00101 !
00102       INTRINSIC ABS,EXP
00103 !
00104       CALL VECTOR(U0 , '=' , 'GRADF          X' , IELM ,
00105      &            1.D0 , PHIR , SBID , SBID , SBID , SBID , SBID ,
00106      &            MESH , MSK  , MASKEL )
00107 !
00108       CALL VECTOR(V0 , '=' , 'GRADF          Y' , IELM ,
00109      &            1.D0 , PHIR , C , C , C , C , C ,
00110      &            MESH , MSK  , MASKEL )
00111 !
00112 !
00113 !      THE OLD VARIABLE U1 IS STORED IN T1
00114 !
00115       CALL VECTOR(T1, '=' , 'GRADF          X' , IELM ,
00116      &            1.D0 , PHII , C , C , C , C , C ,
00117      &            MESH , MSK  , MASKEL )
00118 !
00119 !
00120 !      THE OLD VARIABLE V1 IS STORED IN T2
00121 !
00122       CALL VECTOR(T2, '=' , 'GRADF          Y' , IELM ,
00123      &            1.D0 , PHII , C , C , C , C , C ,
00124      &            MESH , MSK  , MASKEL )
00125 !
00126       CALL VECTOR(T4 , '=' , 'MASBAS          ' , IELM ,
00127      &            1.D0 , C , C , C , C , C , C ,
00128      &            MESH , MSK  , MASKEL )
00129 !
00130       CALL OS( 'X=Y/Z   ' , U0 , U0 , T4 , BID )
00131       CALL OS( 'X=Y/Z   ' , V0 , V0 , T4 , BID )
00132       CALL OS( 'X=Y/Z   ' , T1 , T1 , T4 , BID )
00133       CALL OS( 'X=Y/Z   ' , T2 , T2 , T4 , BID )
00134 !
00135 !--------------------------------------------------------------
00136 !             COMPUTES UE
00137 !--------------------------------------------------------------
00138 !
00139       CALL OS( 'X=C     ' , T4 , SBID , SBID , 0.D0 )
00140       CALL OS( 'X=YZ    ' , T4 , U0 , U0 , 0.D0 )
00141       CALL OS( 'X=X+YZ  ' , T4 , V0 , V0 , 0.D0 )
00142       CALL OS( 'X=X+YZ  ' , T4 , T1 , T1 , 0.D0 )
00143       CALL OS( 'X=X+YZ  ' , T4 , T2 , T2 , 0.D0 )
00144 !
00145       CALL OS( 'X=CX    ' , T4 , SBID , SBID  , 0.5D0 )
00146       CALL OS( 'X=SQR(Y)' , T1 , T4   , SBID  , BID   )
00147       CALL OS( 'X=CY    ' , T4 , T1   , SBID  , 1.2D0 )
00148 !
00149       RETURN
00150       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0