marast.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\marast.f
00002 !
00064                      SUBROUTINE MARAST
00065 !                    *****************
00066 !
00067      &(MARDAT,MARTIM,PHI0,NPOIN,AT,FU1,FV1,X,SINLAT,COSLAT,GRAV)
00068 !
00069 !***********************************************************************
00070 ! TELEMAC2D   V6P2                                  21/08/2010
00071 !***********************************************************************
00072 !
00073 !
00074 !
00075 !
00076 !
00077 !!
00078 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00079 !| AT             |-->| TIME IN SECONDS
00080 !| COSLAT         |-->| COSINUS OF LATITUDE (IN SPHERICAL COORDINATES)
00081 !| FU1            |<--| X-COMPONENT OF TIDE GENERATING FORCE
00082 !| FV1            |<--| Y-COMPONENT OF TIDE GENERATING FORCE
00083 !| GRAV           |-->| GRAVITY
00084 !| MARDAT         |-->| DATE (YEAR, MONTH,DAY)
00085 !| MARTIM         |-->| TIME (HOUR, MINUTEUTE,SECOND)
00086 !| NPOIN          |-->| NUMBER OF POINTS
00087 !| PHI0           |-->| LONGITUDE OF ORIGIN POINT
00088 !| SINLAT         |-->| SINUS OF DE LATITUDE (IN SPHERICAL COORDINATES)
00089 !| X              |-->| ABSCISSAE OF POINTS IN THE MESH
00090 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00091 !
00092       USE INTERFACE_TELEMAC2D, EX_MARAST => MARAST
00093 !
00094       IMPLICIT NONE
00095       INTEGER LNG,LU
00096       COMMON/INFO/LNG,LU
00097 !
00098 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00099 !
00100       INTEGER, INTENT(IN) :: MARDAT(3),MARTIM(3),NPOIN
00101       DOUBLE PRECISION, INTENT(INOUT) :: FU1(NPOIN),FV1(NPOIN)
00102       DOUBLE PRECISION, INTENT(IN) :: COSLAT(NPOIN),SINLAT(NPOIN)
00103       DOUBLE PRECISION, INTENT(IN) :: X(NPOIN),GRAV,AT,PHI0
00104 !
00105 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00106 !
00107       INTEGER YEAR,MONTH,DAY,HOUR,MINUTE,SEC,I
00108 !
00109       DOUBLE PRECISION ARL,ARS,DL,DS,AL,AS
00110       DOUBLE PRECISION RT,LONG,AHL,AHS,MLT,MST,LONG0
00111       DOUBLE PRECISION F0L,F0S,FXL,FYL,FXS,FYS
00112       DOUBLE PRECISION TLOC,TSLOC,TLOC1
00113       DOUBLE PRECISION K2,H2
00114 !
00115       EXTERNAL TSLOC
00116 !
00117       INTRINSIC ACOS,COS,SIN
00118 !
00119 !-----------------------------------------------------------------------
00120 !
00121 !     SPLITS ARRAYS MARDAT AND MARTIM
00122 !
00123       YEAR  = MARDAT(1)
00124       MONTH = MARDAT(2)
00125       DAY   = MARDAT(3)
00126       HOUR  = MARTIM(1)
00127       MINUTE   = MARTIM(2)
00128       SEC   = MARTIM(3)
00129 !
00130 ! REMINUTEDER : HOUR IN UNIVERSAL TIME |
00131 ! GENERAL REMARK : T=TERRE, L=LUNE , S=SOLEIL
00132 !
00133 ! LONG0: REFERENCE LONGITUDE IN RADIAN (0,2PI)
00134 !
00135       LONG0=PHI0*ACOS(-1.D0)/180.D0
00136 !
00137 ! CALLS THE MAIN FUNCTION COMPUTING THE LUNAR AND SOLAR ANGLES
00138 !
00139       CALL ASTRO(YEAR,MONTH,DAY,HOUR,MINUTE,SEC,AT,ARL,ARS,DL,DS,AL,AS)
00140 !
00141 ! RT: EARTH RADIUS IN M
00142 !
00143       RT   = 6378000.D0
00144 !
00145 ! MASS RATIO MOON/EARTH
00146 !
00147       MLT  = 1.D0 / 81.53D0
00148 !
00149 ! MASS RATIO SUN/EARTH
00150 !
00151       MST  = 331954.D0
00152 !
00153 ! AMPLITUDE OF THE FORCE INDUCED BY :
00154 !
00155 !     - THE MOON
00156 !
00157       F0L  = GRAV * MLT * ARL**2
00158 !
00159 !     - THE SUN
00160 !
00161       F0S  = GRAV * MST * ARS**2
00162 !
00163 ! SIDEREAL TIME
00164 !
00165       TLOC1 = TSLOC(YEAR,MONTH,DAY,HOUR,MINUTE,SEC,AT)
00166 !
00167       DO I=1,NPOIN
00168 !
00169 ! LONGITUDE OF THE CONSIDERED NODE
00170 !
00171         LONG = X(I)/RT+LONG0
00172 !
00173 ! LOCAL SIDEREAL TIME
00174 !
00175         TLOC = TLOC1 + LONG
00176 !
00177 ! TIME ANGLE OF THE MOON
00178 !
00179         AHL  = TLOC - AL
00180 !
00181 ! TIME ANGLE OF THE SUN
00182 !
00183         AHS  = TLOC - AS
00184 !
00185 ! FORCE INDUCED BY THE ASTRONOMICAL POTENTIAL ONLY
00186 !
00187 !    FXL : FORCE ALONG X FOR THE MOON
00188 !     Y  : ALONG Y
00189 !     S  : SAME THING FOR THE SUN
00190 !
00191         FXL  = F0L * COS(DL) * SIN(AHL) *
00192      &         ( ( 1.D0-2*ARL*(SINLAT(I)*SIN(DL)+
00193      &           COSLAT(I)*COS(DL)*COS(AHL))+ARL*ARL )**(-1.5D0) -1.D0 )
00194 !
00195         FXS  = F0S * COS(DS) * SIN(AHS) *
00196      &         ( ( 1.D0-2*ARS*(SINLAT(I)*SIN(DS)+
00197      &           COSLAT(I)*COS(DS)*COS(AHS))+ARS*ARS )**(-1.5D0) -1.D0 )
00198 !
00199         FYL  = F0L*(COSLAT(I)*SIN(DL)-SINLAT(I)*COS(DL)*COS(AHL))*
00200      &         ( ( 1.D0-2*ARL*(SINLAT(I)*SIN(DL)+
00201      &           COSLAT(I)*COS(DL)*COS(AHL))+ARL*ARL )**(-1.5D0) -1.D0 )
00202 !
00203         FYS  = F0S*( COSLAT(I)*SIN(DS)-SINLAT(I)*COS(DS)*COS(AHS))*
00204      &         ( ( 1.D0-2*ARS*(SINLAT(I)*SIN(DS)+
00205      &           COSLAT(I)*COS(DS)*COS(AHS))+ARS*ARS )**(-1.5D0) -1.D0 )
00206 !
00207 ! TAKES INTO ACCOUNT :
00208 !
00209 !    - THE TERRESTRIAL TIDE (LOVE NUMBER H2)
00210 !
00211         H2=0.61D0
00212 !
00213 !    - THE STATIC PERTUBATIONS (LOVE NUMBER K2)
00214 !
00215         K2=0.30D0
00216 !
00217 ! MISSES :
00218 !
00219 !    - DYNAMIC PERTUBATION OF AUTO-ATTRACTION
00220 !    - DYNAMIC PERTUBATION OF THE EFFECTS OF LOADS
00221 !
00222 ! FINAL FORCE
00223 !
00224         FU1(I)=FU1(I)+(1.D0+K2-H2)*(FXL+FXS)
00225         FV1(I)=FV1(I)+(1.D0+K2-H2)*(FYL+FYS)
00226 !
00227       ENDDO ! I
00228 !
00229 !-----------------------------------------------------------------------
00230 !
00231       RETURN
00232       END

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