calctm.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\artemis\calctm.f
00002 !
00086                      SUBROUTINE CALCTM
00087 !                    *****************
00088 !
00089 !
00090 !***********************************************************************
00091 ! ARTEMIS   V6P1                                   21/08/2010
00092 !***********************************************************************
00093 !
00094 !
00095 !
00096 !
00097 !
00098 !
00099 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00100 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00101 !
00102       USE BIEF
00103       USE DECLARATIONS_TELEMAC
00104       USE DECLARATIONS_ARTEMIS
00105 !
00106       IMPLICIT NONE
00107       INTEGER LNG,LU
00108       COMMON/INFO/LNG,LU
00109 !
00110       DOUBLE PRECISION PONDER,RADDEG
00111 !
00112       DOUBLE PRECISION BID
00113 !
00114       INTRINSIC SQRT, ATAN2, DMOD, ABS, COS, SIN
00115 !
00116 !-----------------------------------------------------------------------
00117 !
00118 ! STRUCTURES
00119 !
00120 !
00121 !-----------------------------------------------------------------------
00122 !
00123       PONDER = 1.D0
00124 !      PONDER = 1.D0/DBLE(NDALE*NPALE)
00125       RADDEG = 180.D0/3.141592654D0
00126 !
00127 !=======================================================================
00128 ! COMPUTES M0 MOMENTUM AND STORES IT IN T2
00129 !=======================================================================
00130 !
00131       CALL OS( 'X=YZ    ', X=T1 , Y=HHO  , Z=HHO )
00132       CALL OS( 'X=CX    ', T1 , SBID  , SBID  , PONDER )
00133       CALL OS( 'X=Y+Z   ', T2 , HALE , T1   , BID )
00134 !
00135 !=======================================================================
00136 ! T01 = M0 / M1 (adverage period from adverae frequency)
00137 !=======================================================================
00138 !
00139       CALL OS( 'X=Y     ', T3 , T01 , SBID  , BID )
00140       CALL OS( 'X=Y/Z   ', T01, T2  , T3   , BID )
00141 !
00142 !=======================================================================
00143 ! T02 = SQRT( M0 / M2 ) (~ 0 - crossing period Tz)
00144 !=======================================================================
00145 !
00146       CALL OS( 'X=Y     ', T3 , T02 , SBID , BID )
00147       CALL OS( 'X=Y/Z   ', T1 , T2  , T3  , BID )
00148       CALL OS( 'X=SQR(Y)', T02, T1  , SBID , BID )
00149 !
00150 !=======================================================================
00151 ! TM =  MT1 / M0 (adverage period)
00152 !=======================================================================
00153 !
00154       CALL OS( 'X=Y     ', T3 , TM , SBID , BID )
00155       CALL OS( 'X=Y/Z   ', TM , T3 , T2  , BID )
00156 !
00157 !
00158 !=======================================================================
00159 ! MEAN DIRECTION: INCI
00160 !=======================================================================
00161 !
00162       CALL OS( 'X=A(Y,Z)',INCI , MSIN, MCOS , BID )
00163 !
00164       RETURN
00165       END

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