init_zero.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\init_zero.f
00002 !
00098                      SUBROUTINE INIT_ZERO
00099 !                    ********************
00100 !
00101 !
00102 !***********************************************************************
00103 ! SISYPHE   V6P1                                   21/07/2011
00104 !***********************************************************************
00105 !
00106 !
00107 !
00108 !
00109 !
00110 !
00111 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00112 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00113 !
00114       USE DECLARATIONS_SISYPHE
00115       USE BIEF
00116 !
00117       IMPLICIT NONE
00118       INTEGER LNG,LU
00119       COMMON/INFO/LNG,LU
00120 !
00121 !-----------------------------------------------------------------------
00122 !
00123       INTEGER I
00124 !
00125 !-----------------------------------------------------------------------
00126 !
00127 !========================================================================
00128 !                         INITIALISES
00129 ! =======================================================================
00130 !
00131 !---- THE WORKING ARRAYS
00132 !
00133       IF(NPRIV > 0) CALL OS ('X=0     ', X=PRIVE)
00134       CALL OS('X=0     ', X=T12   )
00135       CALL OS('X=0     ', X=COEFPN)
00136 !
00137 !---- THE SEDIMENT VARIABLES :
00138 !
00139       CALL OS('X=0     ', X=QS)
00140       CALL OS('X=0     ', X=QSX)
00141       CALL OS('X=0     ', X=QSY)
00142       CALL OS('X=0     ', X=QSCLXC )
00143       CALL OS('X=0     ', X=QSCLYC )
00144       CALL OS('X=0     ', X=QSCLXS )
00145       CALL OS('X=0     ', X=QSCLYS )
00146 !
00147 ! 7 FOLLOWING LINES ADDED BY JMH 22/04/2005
00148 ! PROVISIONAL INITIALISATION FOR FIRST OUTPUT IN RESULTS FILE
00149 !
00150       CALL OS('X=0     ', X=QSCL )  ! BLOCK OF SIZE NSICLA
00151       CALL OS('X=0     ', X=QS_S )
00152       CALL OS('X=0     ', X=QSXS )
00153       CALL OS('X=0     ', X=QSYS )
00154       CALL OS('X=0     ', X=QS_C )
00155       CALL OS('X=0     ', X=QSXC )
00156       CALL OS('X=0     ', X=QSYC )
00157 !
00158 ! FOLLOWING LINE ADDED BY JMH 04/05/2005
00159 ! PROBABLY USEFUL ONLY IF(CHARR) AND WITH FINITE ELEMENTS
00160 !
00161       CALL OS('X=0     ', X=ZFCL_C )
00162 !
00163 !---- THE DEPOSITION MASSES FOR EVERY CLASS IN SUSPENSION
00164 !
00165       DO I=1,NSICLA
00166         MASDEP(I)=0.D0
00167       ENDDO
00168 !
00169       CALL OS('X=0     ', X=E    )
00170       CALL OS('X=0     ', X=ESOMT)
00171       CALL OS('X=0     ', X=CS   )
00172 !
00173 !---- THE HYDRODYNAMIC VARIABLES :
00174 !
00175       CALL OS('X=0     ', X=QU )
00176       CALL OS('X=0     ', X=QV )
00177       CALL OS('X=0     ', X=U2D )
00178       CALL OS('X=0     ', X=V2D )
00179       CALL OS('X=0     ', X=HN )
00180       CALL OS('X=0     ', X=Q  )
00181       CALL OS('X=0     ', X=TOB)
00182 !
00183 !---- THE WAVE PARAMETERS IF NEED BE
00184 !
00185 !     ALL INITIALISATIONS OF THE WAVES ARE TO BE REMOVED
00186 !     WHEN ALL CHECKS WILL BE DONE
00187 !     SEE BEDLOAD_BAILARD, DIBWAT, BIJKER AND SOULSBY
00188 !
00189 !
00190 !     FW=0.3 CORRESPONDS TO NO WAVES, 0 WOULD DO A LOG(0)
00191       CALL OS('X=C     ', X=FW ,C=0.3D0   )   !
00192       CALL OS('X=0     ', X=HW    )   !
00193       CALL OS('X=0     ', X=TW    )   !
00194       CALL OS('X=C     ', X=THETAW, C=90.D0)  !
00195       CALL OS('X=0     ', X=UW    )   !
00196       CALL OS('X=0     ', X=TOBW)     !
00197 !
00198 !-----------------------------------------------------------------------
00199 !
00200       RETURN
00201       END

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