calcmn.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\artemis\calcmn.f
00002 !
00078                      SUBROUTINE CALCMN
00079 !                    *****************
00080 !
00081 !
00082 !***********************************************************************
00083 ! ARTEMIS   V6P1                                   21/08/2010
00084 !***********************************************************************
00085 !
00086 !
00087 !
00088 !
00089 !
00090 !
00091 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00092 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00093 !
00094       USE BIEF
00095       USE DECLARATIONS_TELEMAC
00096       USE DECLARATIONS_ARTEMIS
00097 !
00098       IMPLICIT NONE
00099       INTEGER LNG,LU
00100       COMMON/INFO/LNG,LU
00101 !
00102       DOUBLE PRECISION FREQ, FREQ2, PONDER
00103 !
00104       DOUBLE PRECISION BID
00105 !
00106       INTRINSIC SQRT, ATAN2, DMOD, ABS, COS, SIN
00107 !
00108 !-----------------------------------------------------------------------
00109 !
00110 ! STRUCTURES
00111 !
00112 !-----------------------------------------------------------------------
00113 !
00114       FREQ = 1.D0/PER
00115       FREQ2 = FREQ * FREQ
00116 !      PONDER = 1.D0/DBLE(NDALE*NPALE)
00117       PONDER= 1D0
00118 !
00119 !=======================================================================
00120 ! M1 = INTEGRAL OF ( F * S(F) * DF )
00121 !=======================================================================
00122 !
00123       CALL OS( 'X=YZ    ', X=T1 , Y=HHO  , Z=HHO )
00124       CALL OS( 'X=CY    ', T2, T1  , SBID , FREQ )
00125       CALL OS( 'X=CX    ', T2, SBID , SBID , PONDER )
00126       CALL OS( 'X=X+Y   ', T01, T2 , SBID , 1.D0 )
00127 !
00128 !=======================================================================
00129 ! M2 = INTEGRAL OF ( F**2 * S(F) * DF )
00130 !=======================================================================
00131 !
00132       CALL OS( 'X=CY    ', T2, T1  , SBID , FREQ2 )
00133       CALL OS( 'X=CX    ', T2, SBID , SBID , PONDER )
00134       CALL OS( 'X=X+Y   ', T02, T2 , SBID , 1.D0 )
00135 !
00136 !=======================================================================
00137 ! MT1 = INTEGRAL OF ( T * S(F) * DF )
00138 !=======================================================================
00139 !
00140       CALL OS( 'X=CY    ', T2 , T1  , SBID , PER )
00141       CALL OS( 'X=CX    ', T2 , SBID , SBID , PONDER )
00142       CALL OS( 'X=X+Y   ', TM , T2  , SBID , BID )
00143 !
00144 !=======================================================================
00145 ! MCOS = INTEGRAL OF ( COS(INCI) * S(F) * DF )
00146 !=======================================================================
00147 !
00148       CALL OS( 'X=COS(Y)',  T2 , INCI , SBID , BID )
00149       CALL OS( 'X=CXY   ',  T2 , T1   , SBID , PONDER )
00150       CALL OS( 'X=X+Y   ', MCOS, T2   , SBID , BID )
00151 !
00152 !=======================================================================
00153 ! MSIN = INTEGRAL OF ( SIN(INCI) * S(F) * DF )
00154 !=======================================================================
00155 !
00156       CALL OS( 'X=SIN(Y)',  T2 , INCI , SBID , BID )
00157       CALL OS( 'X=CXY   ',  T2 , T1   , SBID , PONDER )
00158       CALL OS( 'X=X+Y   ', MSIN, T2   , SBID , BID )
00159 !
00160       RETURN
00161       END

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