condin.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\condin.f
00002 !
00097                      SUBROUTINE CONDIN
00098 !                    *****************
00099 !
00100 !
00101 !***********************************************************************
00102 ! TELEMAC2D   V6P3                                   21/08/2010
00103 !***********************************************************************
00104 !
00105 !
00106 !
00107 !
00108 !
00109 !
00110 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00111 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00112 !
00113       USE BIEF
00114       USE DECLARATIONS_TELEMAC
00115       USE DECLARATIONS_TELEMAC2D
00116       USE TPXO
00117       USE OKADA
00118 !
00119       IMPLICIT NONE
00120       INTEGER LNG,LU
00121       COMMON/INFO/LNG,LU
00122 !
00123 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00124 !
00125 !
00126 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00127 !
00128       INTEGER ITRAC
00129 !
00130 !-----------------------------------------------------------------------
00131 !
00132 !   INITIALISES THE TIME
00133 !
00134       AT = 0.D0
00135 !
00136 !-----------------------------------------------------------------------
00137 !
00138 !   INITIALISES THE VELOCITIES: ZERO VELOCITIES
00139 !
00140       CALL OS('X=0     ',X=U)
00141       CALL OS('X=0     ',X=V)
00142 !
00143 !-----------------------------------------------------------------------
00144 !
00145 !   INITIALISES THE WATER DEPTH H
00146 !
00147       IF(CDTINI(1:10).EQ.'COTE NULLE'.OR.
00148      &   CDTINI(1:14).EQ.'ZERO ELEVATION') THEN
00149         CALL OS( 'X=0     ' , X=H )
00150         CALL OS( 'X=X-Y   ' , X=H , Y=ZF )
00151       ELSEIF(CDTINI(1:14).EQ.'COTE CONSTANTE'.OR.
00152      &       CDTINI(1:18).EQ.'CONSTANT ELEVATION') THEN
00153         CALL OS( 'X=C     ' , H , H  , H , COTINI )
00154         CALL OS( 'X=X-Y   ' , H , ZF , H , 0.D0   )
00155       ELSEIF(CDTINI(1:13).EQ.'HAUTEUR NULLE'.OR.
00156      &       CDTINI(1:10).EQ.'ZERO DEPTH') THEN
00157         CALL OS( 'X=C     ' , H , H  , H , 0.D0  )
00158       ELSEIF(CDTINI(1:17).EQ.'HAUTEUR CONSTANTE'.OR.
00159      &       CDTINI(1:14).EQ.'CONSTANT DEPTH') THEN
00160         CALL OS( 'X=C     ' , H , H  , H , HAUTIN )
00161       ELSEIF(CDTINI(1:25).EQ.'ALTIMETRIE SATELLITE TPXO'.OR.
00162      &       CDTINI(1:24).EQ.'TPXO SATELLITE ALTIMETRY') THEN
00163         CALL OS('X=-Y    ',X=H,Y=ZF)
00164         CALL CONDI_TPXO(NPOIN,MESH%NPTFR,MESH%NBOR%I,
00165      &                  X,Y,H%R,U%R,V%R,
00166      &                  LIHBOR%I,LIUBOR%I,KENT,KENTU,
00167      &                  GEOSYST,NUMZONE,LAMBD0,PHI0,
00168      &                  T2D_FILES,T2DBB1,T2DBB2,
00169      &                  MARDAT,MARTIM,INTMICON,MSL)
00170       ELSEIF(CDTINI(1:13).EQ.'PARTICULIERES'.OR.
00171      &       CDTINI(1:10).EQ.'PARTICULAR'.OR.
00172      &       CDTINI(1:07).EQ.'SPECIAL') THEN
00173 !
00174 !  TO BE MODIFIED BY USER IF SPECIAL INITIAL CONDITIONS
00175 !
00176         IF(LNG.EQ.1) WRITE(LU,10)
00177         IF(LNG.EQ.2) WRITE(LU,11)
00178 10      FORMAT(1X,'CONDIN : AVEC DES CONDITIONS INITIALES PARTICULIERES'
00179      &         ,/,'         VOUS DEVEZ MODIFIER CONDIN')
00180 11      FORMAT(1X,'CONDIN : WITH SPECIAL INITIAL CONDITIONS'
00181      &         ,/,'         YOU HAVE TO MODIFY CONDIN')
00182         CALL PLANTE(1)
00183         STOP
00184 !
00185 !  END OF CODE TO BE MODIFIED BY USER
00186 !
00187       ELSE
00188         IF(LNG.EQ.1) THEN
00189         WRITE(LU,*) 'CONDIN : CONDITION INITIALE NON PREVUE : ',CDTINI
00190         ENDIF
00191         IF(LNG.EQ.2) THEN
00192         WRITE(LU,*) 'CONDIN: INITIAL CONDITION UNKNOWN: ',CDTINI
00193         ENDIF
00194         CALL PLANTE(1)
00195         STOP
00196       ENDIF
00197 !
00198 !-----------------------------------------------------------------------
00199 !
00200 !   INITIALISES TSUNAMI DISPLACEMENT
00201 !
00202       IF(OPTTSUNAMI.EQ.1) THEN
00203         CALL CONDI_OKADA(NPOIN,X,Y,H%R,COETSUNAMI,LAMBD0,PHI0)
00204       ENDIF
00205 !
00206 !-----------------------------------------------------------------------
00207 !
00208 !   INITIALISES THE TRACERS
00209 !
00210       IF(NTRAC.GT.0) THEN
00211         DO ITRAC=1,NTRAC
00212           CALL OS('X=C     ',X=T%ADR(ITRAC)%P,C=TRAC0(ITRAC))
00213         ENDDO
00214       ENDIF
00215 !
00216 !-----------------------------------------------------------------------
00217 !
00218 ! INITIALISES THE VISCOSITY
00219 !
00220       CALL OS('X=C     ',X=VISC,C=PROPNU)
00221 !
00222 !-----------------------------------------------------------------------
00223 !
00224       RETURN
00225       END

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