read_pch_file.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	read_pch_file (flag, matrix, itag, matrix_add, cpt_stiff, cpt_mass, itabm1, pch_file, id, titr)

Function/Subroutine Documentation

read_pch_file()

subroutine read_pch_file	(	integer	flag,
			matrix,
		integer, dimension(*)	itag,
		integer, dimension(4,*)	matrix_add,
		integer	cpt_stiff,
		integer	cpt_mass,
		integer, dimension(*)	itabm1,
		character	pch_file,
		integer	id,
		character(len=nchartitle)	titr )

Definition at line 35 of file read_pch_file.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE inoutfile_mod
      USE message_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ID,ITABM1(*),ITAG(*),MATRIX_ADD(4,*), CPT_STIFF,CPT_MASS,FLAG
      my_real matrix(*)
      CHARACTER PCH_FILE*100
      CHARACTER(LEN=NCHARTITLE) :: TITR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,REF,IFO,TIN,TOUT,POLAR,NCOL,DUM,ID_NOD1,ID_NOD2,J,ISTOP
      my_real rell
      CHARACTER DEBUT*8,NAME*8,MATRIX_TYPE*16,NWLINE*100,MESS*40
 
      INTEGER :: LEN_TMP_NAME,I_ERR
      CHARACTER(len=2148) :: TMP_NAME
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER USR2SYS
      DATA mess/'FLEXIBLE BODY DEF - READ OF PCH FILE '/
C=======================================================================
C
C     FLAG = 0 -> PREREADING FOR NODES DETECTION AND NUMBERING (SETRFXBYON.F)
C     FLAG = 1 -> READING OF NODES AND MATRIX (LECFXB1.F)
C
      ref = 71
      i_err = 0  
C
      len_tmp_name = infile_name_len+len_trim(pch_file)
      tmp_name=infile_name(1:infile_name_len)//pch_file(1:len_trim(pch_file))
      OPEN(unit=ref,file=tmp_name(1:len_tmp_name),
     .      access='SEQUENTIAL',form='FORMATTED',status='UNKNOWN')  
C
      CALL fxrline(ref,nwline,id,titr)
      READ(nwline,'((A,A,5I8,I16))') debut,name,dum,ifo,tin,tout,polar,ncol
C
      IF ((ifo == 6).AND.(tin < 3)) THEN
C
C-------  Symmetric matrix - real (no imaginary part) -----------
C
C       Reading of a new block --
        DO WHILE (i_err == 0)
C
          istop=0
          DO WHILE (istop==0)
            READ(ref,'(A)',END=100) nwline
            IF (nwline(1:5)=='DMIG ') THEN
              READ(nwline,'((A,A,5I8,I16))') debut,name,dum,ifo,tin,tout,polar,ncol
            ELSEIF ((nwline(1:1)/='#').AND.((nwline(1:1)/='$')).AND.((len_trim(nwline)/=0))) THEN
              istop=1
            ENDIF
          ENDDO
C
          READ(nwline,'((A,A,2I16))') debut,matrix_type,id_nod1,i
          id_nod1=usr2sys(id_nod1,itabm1,mess,id)
          IF (id_nod1 /= 0) THEN
            itag(id_nod1) = 1
          ELSE
            RETURN
          ENDIF
C
          CALL fxrline(ref,nwline,id,titr)
C
          IF (adjustl(matrix_type) == 'KAAX') THEN
            DO WHILE (nwline(1:1) == '*')
              READ(nwline,'((A,I16,I16,F16.3))') debut,id_nod2,j,rell
              cpt_stiff = cpt_stiff + 1
              id_nod2=usr2sys(id_nod2,itabm1,mess,id)
              IF (id_nod2 /= 0) THEN
                itag(id_nod2) = 1
                IF (flag == 1) THEN
                  matrix(cpt_stiff-1) = rell
                  matrix_add(1,cpt_stiff-1) = id_nod1   
                  matrix_add(2,cpt_stiff-1) = id_nod2
                  matrix_add(3,cpt_stiff-1) = i 
                  matrix_add(4,cpt_stiff-1) = j 
                ENDIF
              ELSE
                RETURN
              ENDIF   
              READ(ref,'(A)',iostat=i_err,END=100) nwline
            END DO
C
          ELSEIF (adjustl(matrix_type) == 'MAAX') THEN
            DO WHILE (nwline(1:1) == '*')
              READ(nwline,'((A,I16,I16,F16.3))') debut,id_nod2,j,rell
              cpt_mass = cpt_mass + 1
              id_nod2=usr2sys(id_nod2,itabm1,mess,id)
              IF (id_nod2 /= 0) THEN
                IF (flag == 1) THEN
                  matrix(cpt_mass-1) = rell
                  matrix_add(1,cpt_mass-1) = id_nod1   
                  matrix_add(2,cpt_mass-1) = id_nod2
                  matrix_add(3,cpt_mass-1) = i 
                  matrix_add(4,cpt_mass-1) = j 
                ENDIF
              ELSE
                RETURN
              ENDIF   
              READ(ref,'(A)',iostat=i_err,END=100) nwline
            END DO
C
          ELSEIF ((adjustl(matrix_type)=='baax.OR.')(ADJUSTL(MATRIX_TYPE)=='pax.OR.')ADJUSTL(MATRIX_TYPE)=='k4aax') THEN
C--         Input not compatible
            I_ERR = 1
            IF (FLAG == 0) THEN
              CALL ANCMSG(MSGID=1738,
     .                  MSGTYPE=MSGERROR,
     .                  ANMODE=ANINFO_BLIND_1,
     .                  I1=ID,
     .                  C1=TITR,
     .                  C2=ADJUSTL(MATRIX_TYPE))
            ENDIF            
C
          ENDIF      
C
          BACKSPACE(REF)
C
        END DO
C
100     CONTINUE
C
      ENDIF
C
      CLOSE(UNIT=REF,STATUS='keep')                            
C 
      RETURN     
C