thrnur.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "task_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	thrnur (iad, nn, iadv, nvar, iparg, ithbuf, bufel, wa)

Function/Subroutine Documentation

thrnur()

subroutine thrnur	(	integer	iad,
		integer	nn,
		integer	iadv,
		integer	nvar,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(*)	ithbuf,
			bufel,
			wa )

Definition at line 28 of file thrnur.F.

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      "com01_c.inc"
#include      "task_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IAD,NN,IADV,NVAR,
     .   IPARG(NPARG,*),ITHBUF(*)
      my_real
     .   bufel(*),wa(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER II, I, N, IH, NG, ITY, MTE, NB1, NB2, NB3,
     .   NB4, NB5, NB6, NB10, NB11, NB12, NB13, K, M1, M2, M3, M4, M5,
     .   M6, N1, N2, N3, N4, N5, NB7, NB8, M11, M10, NB9,IP,L,
     .   NB2A, NB2B, NB4A, NB4B,  NB9A, NB9B, M8,NB14, NB15, NB16, NB17,
     .   NB10A, NB10B, NB12A, NB12B, NB18,NB8A, NB8B, OFFSET1,OFFSET2,
     .   LWA,NEL,NFT
      my_real
     .   wwa(100),xfn
      ii=0
      ih=iad
C   specifique spmd
      lwa = nn*nvar
      DO i = 1, lwa
        wa(i) = zero
      ENDDO
C decalage IH
      DO WHILE((ithbuf(ih+nn)/=ispmd).AND.(ih<iad+nn))
        ih = ih + 1
      ENDDO
      IF (ih>=iad+nn) RETURN 
C
      DO ng=1,ngroup
        ity=iparg(5,ng)
        IF(ity==50) THEN
          mte=iparg(1,ng)
          nel=iparg(2,ng)
          nft=iparg(3,ng)
          nb1=iparg(4,ng)
          nb2=nb1+nel*11
          nb3=nb2+nel*12
          nb4=nb3+nel
          nb5=nb4+nel*13
          nb6=nb5+nel*3
C
          DO i=1,nel
            n=i+nft
            k=ithbuf(ih)
            ip=ithbuf(ih+nn)
C
            IF (k==n)THEN
              ih=ih+1
C   traitement specifique spmd
 
C recherche du ii correct
              ii = ((ih-1) - iad)*nvar
              DO WHILE((ithbuf(ih+nn)/=ispmd).AND.(ih<iad+nn))
                ih = ih + 1
              ENDDO
C
              IF(ih>iad+nn) RETURN
C             IF (IP==ISPMD)THEN
C
              m1=nb2+3*i-3
              m2=nb5+3*i-3
              m3=nb6+3*i-3
              xfn=bufel(m1)*bufel(m2)+bufel(m1+1)*bufel(m2+1)
     .           +bufel(m1+2)*bufel(m2+2)
              wwa(1)=bufel(nb1+i)
              wwa(2)=bufel(nb3+i)
              wwa(3)=bufel(nb4+i)
              wwa(4)=xfn*bufel(m1)
              wwa(5)=xfn*bufel(m1+1)
              wwa(6)=xfn*bufel(m1+2)
              wwa(7)=bufel(m2)  -xfn*bufel(m1)
              wwa(8)=bufel(m2+1)-xfn*bufel(m1+1)
              wwa(9)=bufel(m2+2)-xfn*bufel(m1+2)
              wwa(10)=bufel(m3)
              wwa(11)=bufel(m3+1)
              wwa(12)=bufel(m3+2)
              DO k=1,12
                ii=ii+1
                wa(ii)=wwa(k)
              ENDDO
C             ENDIF
            ENDIF
          ENDDO 
        ENDIF
      ENDDO
C
      RETURN