pmat3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	pmat3 (jft, jlt, pm, geo, ym, g, a1, b1, b2, b3, al, ali, sh1, sh2, mtn, mat, pid)

Function/Subroutine Documentation

pmat3()

subroutine pmat3	(	integer	jft,
		integer	jlt,
			pm,
			geo,
			ym,
			g,
			a1,
			b1,
			b2,
			b3,
			al,
			ali,
			sh1,
			sh2,
		integer	mtn,
		integer, dimension(*)	mat,
		integer, dimension(*)	pid )

Definition at line 28 of file pmat3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT  ,MTN  
      INTEGER MAT(*), PID(*)
C     REAL
      my_real
     .   geo(npropg,*), pm(npropm,*), 
     .    ym(*),g(*),al(*),ali(*),sh1(*),sh2(*),
     .    a1(*),b1(*),b2(*),b3(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX,IPID,J
C     REAL
      my_real
     .   shf(mvsiz),sh,yma2,ali2,sh10,sh0,sh20
C------------- ---------------------------------
C
      DO i=jft,jlt
       g(i)   =pm(22,mat(i))
       ym(i)  =pm(20,mat(i))
       a1(i)  =geo(1,pid(i))
       b1(i)  =geo(2,pid(i))
       b2(i)  =geo(18,pid(i))
       b3(i)  =geo(4,pid(i))
       shf(i) =geo(37,pid(i))
      ENDDO
C
      DO i=jft,jlt
       sh=five_over_6*g(i)*a1(i)
       ali2 =ali(i)*ali(i) 
       yma2=twelve*ym(i)*ali2
       sh10=yma2*b1(i)
       sh20=yma2*b2(i)
       sh0=(one-shf(i))*sh
       sh1(i)=sh0*sh10/(sh+sh10) + shf(i)*sh10
       sh2(i)=sh0*sh20/(sh+sh20) + shf(i)*sh20
      ENDDO
C
      RETURN