intfop1.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "task_c.inc"
#include "param_c.inc"
#include "parit_c.inc"
#include "impl1_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	intfop1 (output, ipari, x, a, icodt, fsav, wa, v, ms, dt2t, neltst, ityptst, itab, stifn, npc, tf, fskyi, isky, vr, fcont, in, igrsurf, bufsf, fncont, ftcont, icontact, rcontact, num_imp, ns_imp, ne_imp, nt_imp, sensor_tab, intbuf_tab, h3d_data, nsensor)

Function/Subroutine Documentation

intfop1()

subroutine intfop1	(	type(output_), intent(inout)	output,
		integer, dimension(npari,*)	ipari,
			x,
			a,
		integer, dimension(*)	icodt,
			fsav,
			wa,
			v,
			ms,
			dt2t,
		integer	neltst,
		integer	ityptst,
		integer, dimension(*)	itab,
			stifn,
		integer, dimension(*)	npc,
			tf,
			fskyi,
		integer, dimension(*)	isky,
			vr,
			fcont,
			in,
		type (surf_), dimension(nsurf)	igrsurf,
			bufsf,
			fncont,
			ftcont,
		integer, dimension(*)	icontact,
			rcontact,
		integer, dimension(*)	num_imp,
		integer, dimension(*)	ns_imp,
		integer, dimension(*)	ne_imp,
		integer	nt_imp,
		type (sensor_str_), dimension(nsensor)	sensor_tab,
		type(intbuf_struct_), dimension(*)	intbuf_tab,
		type(h3d_database)	h3d_data,
		integer, intent(in)	nsensor )

Definition at line 40 of file intfop1.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE output_mod, ONLY : output_
      USE intbufdef_mod
      USE h3d_mod
      USE groupdef_mod
      USE sensor_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "task_c.inc"
#include      "param_c.inc"
#include      "parit_c.inc"
#include      "impl1_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
      INTEGER ,INTENT(IN) :: NSENSOR
      INTEGER NELTST,ITYPTST 
      INTEGER NUM_IMP(*),NS_IMP(*),NE_IMP(*),NT_IMP
      INTEGER IPARI(NPARI,*), ICODT(*),ITAB(*), 
     .        NPC(*), ISKY(*), 
     .        ICONTACT(*)
C     REAL
      my_real dt2t, 
     .   x(*), a(*), fsav(nthvki,*) , wa(*), v(*),ms(*),
     .   stifn(*), tf(*),fskyi(lsky,nfskyi),vr(3,*),fcont(3,*),in(*),
     .   bufsf(*), fncont(3,*),ftcont(3,*),rcontact(*)
 
      TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
      TYPE(H3D_DATABASE) :: H3D_DATA
      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF
      TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) :: SENSOR_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N, IAD,IDUM, ISENS, NTY
      my_real ts
C-----------------------------------------------
      IF (impl_s==1) iad = 1
C
      DO 200 n=1,ninter
       nty   =ipari(7,n)
       IF(nty==14.OR.nty==15) i7kglo = 1
       IF(ispmd/=0) GOTO 200
       IF(nty==3.OR.nty==4.OR.nty==5
     .    .OR.nty==14.OR.nty==15)THEN
#include "lockon.inc"
        IF(ipari(16,n)==ncycle)THEN
#include "lockoff.inc"
        ELSE
         ipari(16,n) = ncycle
#include "lockoff.inc"
C 
         isens = 0
         IF(nty == 5 ) isens = ipari(64,n) ! INTERFACE SENSOR NUMBER
 
         IF (isens > 0) THEN    ! IF INTERFACE IS ACTIVATED BY SENSOR 
            ts = sensor_tab(isens)%TSTART
         ELSE
            ts = tt
         ENDIF
C
         IF(n2d==0)THEN
          IF (nty==14) THEN
           CALL i14cmp(output,ipari(1,n),intbuf_tab(n),igrsurf ,bufsf ,
     2              a ,x ,v ,wa ,fsav(1,n) ,
     3              in ,stifn ,fcont ,ms ,fskyi , 
     4              isky ,npc ,tf ,h3d_data)
          ELSEIF (nty==15) THEN
           CALL i15cmp(output,ipari(1,n),intbuf_tab(n),
     2            igrsurf,bufsf,a, x, v,
     3            wa, fsav(1,n),in, stifn, fcont,
     4            ms, fskyi,isky, npc, tf,
     5            dt2t, neltst,ityptst, itab,h3d_data)
          ELSE
           IF(impl_s/=1)THEN
              IF(tt>=ts) THEN
                CALL intvo3(
     1        ipari(1,n) ,x,a    ,
     2        icodt      ,fsav(1,n)    ,wa           ,v ,ms       ,
     3        fskyi      ,isky         ,fcont   ,fncont ,ftcont   ,
     4        icontact   ,rcontact     ,idum    ,idum   ,idum     ,
     5          stifn      ,itab,intbuf_tab(n),h3d_data)
             ENDIF
           ELSE
             num_imp(n) =0
              IF(tt>=ts) THEN
                CALL intvo3(
     1          ipari(1,n) ,x ,a        ,
     2          icodt      ,fsav(1,n)    ,wa           ,v ,ms       ,
     3          fskyi      ,isky         ,fcont   ,fncont ,ftcont   ,
     4          icontact   ,rcontact,num_imp(n),ns_imp(iad),ne_imp(iad),
     5          stifn      ,itab,intbuf_tab(n),h3d_data)
              ENDIF
             iad=iad+num_imp(n)
           END IF
          END IF
         ELSE
          IF(tt>=ts) THEN
            CALL intvo2(
     1        ipari(1,n) ,intbuf_tab(n),x   ,a ,
     2        icodt      ,fsav(1,n)    ,fskyi        ,isky   ,
     3        fcont      ,fncont       , ftcont     ,h3d_data )
          ENDIF
         ENDIF
        ENDIF
       ELSEIF(nty==6)THEN
#include "lockon.inc"
         IF(ipari(16,n)==ncycle)THEN
#include "lockoff.inc"
         ELSE
           ipari(16,n) = ncycle
#include "lockoff.inc"
           CALL i6main(
     1        ipari(1,n),x ,a        ,
     2        icodt     ,npc          ,tf           ,v ,fsav(1,n),
     3        dt2t      ,neltst       ,ityptst      ,fskyi ,isky ,
     4        fcont     ,ms           ,intbuf_tab(n),h3d_data )
         ENDIF
       ELSEIF(nty==7.OR.nty==22)THEN
       ELSEIF(nty==10)THEN
       ELSEIF(nty==11)THEN
       ENDIF
  200 CONTINUE
C
#include "lockon.inc"
      
      IF (impl_s==1) nt_imp=iad-1
#include "lockoff.inc"
 
C      
      RETURN