cm27in3.F File Reference

#include "implicit_f.inc"
#include "vect01_c.inc"
#include "param_c.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cm27in3 (elbuf_str, geo, igeo, pm, ipm, ix, nix, nlay, ir, is, imat)

Function/Subroutine Documentation

cm27in3()

subroutine cm27in3	(	type(elbuf_struct_), target	elbuf_str,
			geo,
		integer, dimension(npropgi,*)	igeo,
			pm,
		integer, dimension(npropmi,*)	ipm,
		integer, dimension(nix,*)	ix,
		integer	nix,
		integer	nlay,
		integer	ir,
		integer	is,
		integer	imat )

Definition at line 35 of file cm27in3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod            
      USE message_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C INITIALISE LES DIRECTIONS DE FISSURES
C INITIALISE LES EPAISSEURS ET LES MATERIAUX DES COUCHES
C
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
#include      "param_c.inc"
#include      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IGEO(NPROPGI,*),IX(NIX,*),NIX,IPM(NPROPMI,*),NLAY,IR,IS,IMAT
      my_real geo(npropg,*),pm(npropm,*)
      TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IPTHK,IPMAT,IPPOS,MID,MLN,IMID,IPID,PID
      INTEGER I,N,I1,I2,I3,J,IGTYP,II,L_DMG,ILAYER,IT,NPTT
      my_real thkl,pos0,dp
      CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
      my_real, DIMENSION(:), POINTER :: dir_dmg
C
      TYPE(BUF_LAY_) ,POINTER :: BUFLY
      TYPE(L_BUFEL_) ,POINTER :: LBUF
C======================================================================|
      IF(npt==0)THEN
        imid=ix(1,1)
        ipid = ix(nix-1,1)
        pid = igeo(1,ipid)
        mid = ipm(1,imid)
        mln = nint(pm(19,imid))
        CALL fretitl2(titr,igeo(npropgi-ltitr+1,pid),ltitr)
        CALL fretitl2(titr1,ipm(npropmi-ltitr+1,imid),ltitr)
        CALL ancmsg(msgid=23, anmode=aninfo, msgtype=msgerror, i1=pid, c1=titr, i2=mid, c2=titr1, i3=27)
      ENDIF
C
      igtyp=nint(geo(12,imat))
      IF (igtyp /= 51 .AND. igtyp /= 52) THEN
        DO n=1,npt
          ilayer = n
          IF (nlay > 1) THEN
            lbuf => elbuf_str%BUFLY(ilayer)%LBUF(ir,is,1)
          ELSE
            lbuf => elbuf_str%BUFLY(1)%LBUF(ir,is,ilayer)
          ENDIF
C
          l_dmg = elbuf_str%BUFLY(1)%L_DMG
          dir_dmg => lbuf%DMG(1:l_dmg*llt)
C
          DO i=lft,llt
            dir_dmg(i) = one
            dir_dmg(i+llt) = zero
          ENDDO
        ENDDO
      ELSEIF (igtyp == 51) THEN
        DO ilayer=1,nlay
          nptt = elbuf_str%BUFLY(ilayer)%NPTT
          DO it=1,nptt
            lbuf => elbuf_str%BUFLY(ilayer)%LBUF(ir,is,it)
C
            l_dmg = elbuf_str%BUFLY(ilayer)%L_DMG
            dir_dmg => lbuf%DMG(1:l_dmg*llt)
C
            DO i=lft,llt
              dir_dmg(i) = one
              dir_dmg(i+llt) = zero
            ENDDO
          ENDDO
        ENDDO
      ENDIF
C
      ipthk = 300
      ippos = 400
      ipmat = 100
      IF(igtyp/=11 .AND. igtyp/=17 . and. igtyp/=51) THEN
        thkl = one / npt
        pos0 =-half*(one + thkl)
        DO n=1,npt                                                
          i1=ippos+n                                              
          i2=ipthk+n                                              
          i3=ipmat+n                                              
          dp =  n*thkl
          DO i=lft,llt                                            
            geo(i1,imat) = pos0 + dp 
            geo(i2,imat) = thkl                              
          ENDDO                                                   
        ENDDO                                                     
      ENDIF
C-----------
      RETURN