cbapinchthk.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cbapinchthk (jft, jlt, nplat, iplat, dt1c, thk, thk0, ezzpg)

Function/Subroutine Documentation

cbapinchthk()

subroutine cbapinchthk	(	integer	jft,
		integer	jlt,
		integer	nplat,
		integer, dimension(*)	iplat,
			dt1c,
			thk,
			thk0,
			ezzpg )

Definition at line 30 of file cbapinchthk.F.

C-------------------------------------------------------------------------------
C   M o d u l e s
C-------------------------------------------------------------------------------
      USE pinchtype_mod
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   D u m m y   A r g u m e n t s
C-------------------------------------------------------------------------------
      INTEGER JFT, JLT, NPLAT,IPLAT(*)
C     REAL
      my_real
     .   dt1c(*), thk(*), thk0(*), ezzpg(mvsiz,4)
C-------------------------------------------------------------------------------
C   L o c a l   V a r i a b l e s
C-------------------------------------------------------------------------------
      INTEGER I,J,EP
      my_real thkn(mvsiz,4), thkavg
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
 
#include "vectorize.inc"
      DO i=jft,jlt 
        ep=iplat(i)
        thkn(ep,1) = thk(ep) * (one + ezzpg(ep,1))
        thkn(ep,2) = thk(ep) * (one + ezzpg(ep,2))
        thkn(ep,3) = thk(ep) * (one + ezzpg(ep,3))
        thkn(ep,4) = thk(ep) * (one + ezzpg(ep,4))
        thkavg = fourth*(thkn(ep,1)+thkn(ep,2)+thkn(ep,3)+thkn(ep,4))
        thk(ep) = thkavg
        thk0(ep) = thkavg      
      ENDDO       
      RETURN