fredamp.F File Reference

#include "implicit_f.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fredamp (ikad, key0, kdamp, ndamp)

Function/Subroutine Documentation

fredamp()

subroutine fredamp	(	integer, dimension(0:*)	ikad,
		character, dimension(*)	key0,
		integer	kdamp,
		integer	ndamp )

Definition at line 34 of file fredamp.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
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 IKAD(0:*),KDAMP,NDAMP
      CHARACTER KEY0(*)*5
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER K, NN, ID, IKEY, IDG, IGR, NBC, KORTH
      my_real dampa, dampb
      
C=======================================================================
      k=0
      ikey = kdamp
      igr = 0
      IF (ndamp == 1) THEN
        READ(iusc1,rec=ikad(ikey)+k,fmt='(6X,I10,85X,I10)',err=9990)
     .                                                       id,nbc
 
        IF (id/=1.AND.nbc==1) THEN
          korth=0
          k=k+1
          CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
          READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
          WRITE(iin,'(I10,1P2G20.13,2I10)') id,dampa,dampb,igr,korth
        ELSEIF(id==1) THEN
          IF(nbc==2) THEN
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) igr
          ENDIF
          IF(nbc/=6)THEN
            korth=0
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(I10,1P2G20.13,2I10)') id,dampa,dampb,igr,korth
          ELSE
            korth=1
C           dampax, dampbx
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(I10,1P2G20.13,2I10)') id,dampa,dampb,igr,korth
C           dampay, dampby
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           dampaz, dampbz
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           damparx, dampbrx
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           dampary, dampbry
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           damparz, damprbz
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
          END IF
        ELSE
          GOTO 9990
        ENDIF
      ELSE
        DO nn=1,ndamp
          READ(iusc1,rec=ikad(ikey)+k,fmt='(6X,I10,85X,I10)',err=9990)
     .                                                       id,nbc
     
          IF(nbc/=6)THEN
            korth=0
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(I10,1P2G20.13,10X,I10)') id,dampa,dampb,korth
          ELSE
            korth=1
C           dampax, dampbx
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(I10,1P2G20.13,10X,I10)') id,dampa,dampb,korth
C           dampay, dampby
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           dampaz, dampbz
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           damparx, dampbrx
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           dampary, dampbry
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
C           damparz, damprbz
            k=k+1
            CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
            READ(iusc2,*,err=9990,END=9990) DAMPA,dampb
            WRITE(iin,'(1P2G20.13)') dampa,dampb
          END IF
          k=k+1
        ENDDO
      ENDIF
C---
      RETURN
 9990 CONTINUE
      CALL ancmsg(msgid=73,anmode=aninfo,
     .            c1=key0(ikey))
      CALL arret(0)