spload3.F File Reference

#include "implicit_f.inc"
#include "sphcom.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	spload3 (spbuf, wa, rhoa, rhon, deltax, offg, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, mumax, lft, llt, nft)

Function/Subroutine Documentation

spload3()

subroutine spload3	(		spbuf,
			wa,
			rhoa,
			rhon,
			deltax,
			offg,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			mumax,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt,
		integer, intent(inout)	nft )

Definition at line 28 of file spload3.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      "sphcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: NFT
      my_real
     .   spbuf(nspbuf,*) ,wa(kwasph,*)  ,
     .   rhoa(*) ,rhon(*) ,deltax(*) ,offg(*) ,
     .   dxx(*), dxy(*), dxz(*),
     .   dyx(*), dyy(*), dyz(*),
     .   dzx(*), dzy(*), dzz(*),
     .   mumax(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N,I
C     REAL
      my_real
     .        off_l
C-----------------------------------------------
      DO 100  i=lft,llt
         n=nft+i
         rhon(i)  =spbuf(2,n)
 100  CONTINUE
C
      off_l = zero
      DO i=lft,llt
        off_l  = min(off_l,offg(i))
      ENDDO
C-----
      IF(off_l>=zero)THEN
      DO 200  i=lft,llt
         n=nft+i
         dxx(i)   =wa(1,n)
         dyy(i)   =wa(2,n)
         dzz(i)   =wa(3,n)
         dxy(i)   =wa(4,n)
         dyz(i)   =wa(5,n)
         dxz(i)   =wa(6,n)
         dyx(i)   =wa(7,n)
         dzy(i)   =wa(8,n)
         dzx(i)   =wa(9,n)
         rhoa(i)  =wa(10,n)
C        for stability time step computation:
         deltax(i)=wa(11,n)
         mumax(i) =wa(12,n)
 200  CONTINUE
      ELSE
      DO 300  i=lft,llt
       n=nft+i
       rhoa(i)  =wa(10,n)
       IF(offg(i)>=zero)THEN
         dxx(i)   =wa(1,n)
         dyy(i)   =wa(2,n)
         dzz(i)   =wa(3,n)
         dxy(i)   =wa(4,n)
         dyz(i)   =wa(5,n)
         dxz(i)   =wa(6,n)
         dyx(i)   =wa(7,n)
         dzy(i)   =wa(8,n)
         dzx(i)   =wa(9,n)
C        for stability time step computation:
         deltax(i)=wa(11,n)
         mumax(i) =wa(12,n)
       ELSE
         dxx(i)   =zero
         dyy(i)   =zero
         dzz(i)   =zero
         dxy(i)   =zero
         dyz(i)   =zero
         dxz(i)   =zero
         dyx(i)   =zero
         dzy(i)   =zero
         dzx(i)   =zero   
C        for stability time step computation:
         deltax(i)=ep20
         mumax(i) =zero
       ENDIF
 300  CONTINUE
      ENDIF
C----------------------------------
      RETURN