cvsp_init_from_layers.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\cvsp_init_from_layers.f
00002 !
00056                      SUBROUTINE CVSP_INIT_FROM_LAYERS
00057 !                    ********************************
00058 !
00059 !
00060 !***********************************************************************
00061 ! SISYPHE V6P3                                   14/03/2013
00062 !***********************************************************************
00063 !
00064 !
00065 !
00066 !
00067 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00068 !|           |---|
00069 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00070 !
00071       USE DECLARATIONS_TELEMAC
00072       USE DECLARATIONS_SISYPHE
00073 !
00074       IMPLICIT NONE
00075 !
00076       INTEGER  I, J, K, M, L, MDISC, UBS
00077       DOUBLE PRECISION DEPTH
00078 !
00079 !-----------------------------------------------------------------------
00080 !
00081       ALLOCATE(PRO_D(NPOIN,PRO_MAX_MAX,NSICLA))
00082       ALLOCATE(PRO_F(NPOIN,PRO_MAX_MAX,NSICLA))
00083       ALLOCATE(PRO_MAX(NPOIN))
00084 !
00085 !-----------------------------------------------------------------------
00086 !
00087       DO J=1,NPOIN
00088         DEPTH = 0                    ! INIT DEPTH OF THE VSP
00089         PRO_MAX(J) =  2* NLAYER%I(J) ! 2 SECTION POINTS PER LAYER
00090         L = PRO_MAX(J)
00091 !
00092 !-----------------------------------------------------------------------
00093 ! WATER / BOTTOM
00094 !-----------------------------------------------------------------------
00095 !
00096         DO I=1,NSICLA
00097           PRO_D(J,L,I) = ZF%R(J)
00098           PRO_F(J,L,I) = AVAIL(J,1,I)
00099         ENDDO
00100 !
00101 !-----------------------------------------------------------------------
00102 ! SECTIONS
00103 !-----------------------------------------------------------------------
00104 !
00105         DO M=1,NLAYER%I(J)-1   !FOR THE UPPER 8 LAYERS
00106           DEPTH = DEPTH + ES(J,M)
00107           L = L - 1
00108           DO I=1,NSICLA
00109             PRO_D(J,L,I) = ZF%R(J) - DEPTH
00110             PRO_F(J,L,I) = AVAIL(J,M,I)
00111           ENDDO
00112           L = L - 1
00113           DO I=1,NSICLA
00114             PRO_D(J,L,I) = ZF%R(J) - DEPTH
00115             PRO_F(J,L,I) = AVAIL(J,M+1,I)
00116           ENDDO
00117         ENDDO
00118 !
00119 !-----------------------------------------------------------------------
00120 ! BOTTOM / RIGID BED
00121 !-----------------------------------------------------------------------
00122 !
00123         L = L - 1
00124         DO I=1,NSICLA
00125           PRO_D(J,L,I) = ZR%R(J)
00126           PRO_F(J,L,I) = AVAIL(J,NLAYER%I(J),I)
00127         ENDDO
00128 !
00129         CALL CVSP_COMPRESS_DP(J,1.D-5)
00130 !
00131       ENDDO
00132 !
00133 !-----------------------------------------------------------------------
00134 !
00135       RETURN
00136       END SUBROUTINE CVSP_INIT_FROM_LAYERS

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