init_compo_coh.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\init_compo_coh.f
00002 !
00033                      SUBROUTINE INIT_COMPO_COH
00034 !                    *************************
00035 !
00036      &(ES,CONC_VASE,CONC,NPOIN,NOMBLAY,NSICLA,AVAIL,AVA0)
00037 !
00038 !***********************************************************************
00039 ! SISYPHE   V6P2                                   21/07/2011
00040 !***********************************************************************
00041 !
00042 !
00043 !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00044 !| AVA0           |-->| VOLUME PERCENT
00045 !| AVAIL          |<->| VOLUME PERCENT OF EACH CLASS
00046 !| CONC           |<->| CONC OF EACH BED LAYER (KG/M3)
00047 !| CONC_VASE      |<->| MUD CONCENTRATION FOR EACH LAYER
00048 !| ES             |<->| LAYER THICKNESSES AS DOUBLE PRECISION
00049 !| NOMBLAY        |-->| NUMBER OF LAYERS FOR CONSOLIDATION
00050 !| NPOIN          |-->| NUMBER OF POINTS
00051 !| NSICLA         |-->| NUMBER OF SIZE CLASSES FOR BED MATERIALS
00052 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00053 !
00054       USE BIEF
00055       USE INTERFACE_SISYPHE, EX_INIT_COMPO_COH=> INIT_COMPO_COH
00056       USE DECLARATIONS_SISYPHE, ONLY : NLAYMAX
00057 !
00058       IMPLICIT NONE
00059       INTEGER LNG,LU
00060       COMMON/INFO/LNG,LU
00061 !
00062 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00063 !
00064       INTEGER, INTENT(IN)              :: NPOIN,NOMBLAY,NSICLA
00065       DOUBLE PRECISION, INTENT(INOUT)  :: ES(NPOIN,NOMBLAY)
00066       DOUBLE PRECISION, INTENT(IN)     :: CONC_VASE(NOMBLAY)
00067       DOUBLE PRECISION,  INTENT(INOUT) :: CONC(NPOIN,NOMBLAY)
00068       DOUBLE PRECISION, INTENT(INOUT)  :: AVAIL(NPOIN,NOMBLAY,NSICLA)
00069       DOUBLE PRECISION, INTENT(IN)     :: AVA0(NSICLA)
00070 !
00071 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00072 !
00073       DOUBLE PRECISION  EPAI_VASE(NLAYMAX),EPAI_SABLE(NLAYMAX)
00074 
00075       INTEGER I,J
00076 !
00077 !-----------------------------------------------------------------------
00078 !    EXAMPLE FOR NOMBLAY = 10
00079 !
00080 !     EPAI_VASE(1)=0.0525D0
00081 !     EPAI_VASE(2)=0.0385D0
00082 !     EPAI_VASE(3)=0.03995D0
00083 !     EPAI_VASE(4)=0.0437D0
00084 !     EPAI_VASE(5)=0.0517D0
00085 !     EPAI_VASE(6)=0.1259D0
00086 !     EPAI_VASE(7)=0.4889D0
00087 !     EPAI_VASE(8)=1.5071D0
00088 !     EPAI_VASE(9)=0.86410D0
00089 !     EPAI_VASE(9)=0.80D0
00090 !
00091 !
00092 !     HERE A CONSTANT
00093 !
00094       DO J= 1,NOMBLAY
00095         EPAI_VASE(J) = 0.1D0
00096         IF(NSICLA.GT.1) THEN
00097           EPAI_SABLE(J) = AVA0(1)/AVA0(2)*EPAI_VASE(J)
00098         ENDIF
00099       ENDDO
00100 !-----------------------------------------------------------------------
00101 !
00102 !     INITIALISING OF LAYER THICKNESS AND CONC
00103 !
00104 
00105 !     BY DEFAULT : UNIFORM BED COMPOSITION (KEY WORDS)
00106 !     V6P3: IT WILL BE POSSIBLE TO HAVE A SPATIAL DISTRIBUTION OF THE BED CONC
00107 !     V6P2: SO FAR THE MUD CONC IS CONSTANT PER LAYER
00108 !     si mixte: calculer aussi les AVAI!
00109 !
00110       DO I=1,NPOIN
00111         DO J= 1,NOMBLAY
00112 !
00113           CONC(I,J) = CONC_VASE(J)
00114           ES(I,J)   = EPAI_VASE(J)
00115 !
00116           IF(NSICLA.GT.1) THEN
00117             ES(I,J)= ES(I,J) + EPAI_SABLE(J)
00118             IF(ES(I,J).GE.1.D-6) THEN
00119 ! Class 1 is for sand, class 2 is mud
00120               AVAIL(I,J,1)= EPAI_SABLE(J)/ES(I,J)
00121               AVAIL(I,J,2)= EPAI_VASE(J)/ES(I,J)
00122             ELSE
00123               AVAIL(I,J,1)= AVA0(1)
00124               AVAIL(I,J,2)= AVA0(2)
00125             ENDIF
00126           ENDIF
00127 !
00128         ENDDO
00129       ENDDO
00130 !
00131 !-----------------------------------------------------------------------
00132 !
00133       RETURN
00134       END

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