mintp1_rt.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	mintp1_rt (xx, nx, yy, ny, zz, x, y, z, dzdx, dzdy)

Function/Subroutine Documentation

mintp1_rt()

subroutine mintp1_rt	(		xx,
		integer	nx,
			yy,
		integer	ny,
			zz,
			x,
			y,
			z,
			dzdx,
			dzdy )

Definition at line 30 of file mintp1_rt.F.

C=======================================================================
C
C                           |         |                |          |
C                         --+---------+--------*-------+----------+-- YY(IY+2)
C     INPUT :               |         |       Z3       |          |
C          XX(NX)           |         |                |          |
C          NX               |         |                |          |
C          YY(NY)           |         |                |          |
C          ZZ(NX,NY)        |     Z(IX,IY+1)     Z(IX+1,IY+1)     |
C          X              --+---------+--------*-------+----------+-- YY(IY+1)
C          Y                |         |       Z2       |          |
C                           |         |                |          |
C      OUTPUT:              |         |                |          |
C          Z                |         |      Z(X,Y)    |          |
C          DZDX             |         |        * - - - | - - - - -| - - - Y
C          DZDY             |         |        .       |          |
C                           |         |        .       |          |
C                           |     Z(IX,IY)     .   Z(IX+1,IY)     |
C                         --+---------+--------*-------+----------+-- YY(IY)
C                           |         |       Z1       |          |
C                           |         |        .       |          |
C                           |         |        .       |          |
C                           |         |       Z0       |          |
C                         --+---------+--------*-------+----------+-- YY(IY-1)
C                           |         |        .       |          |
C                       XX(IX-1)   XX(IX)      .   XX(IX+1)   XX(IX+2)
C                                              X
C
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER NX, NY
C     REAL
        my_real
     .     x, y, z, dzdx, dzdy
C     REAL
        my_real
     .     xx(nx), yy(ny), zz(nx,ny)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER IX, IY, IXM1, IXP1, IXP2, IYM1, IYP1, IYP2
C     REAL
        my_real
     .     xx0, xx1, xx2, xx3, zz00, zz10, zz20, zz30, zz01, zz11, zz21,
     .     zz31, zz02, zz12, zz22, zz32, zz03, zz13, zz23, zz33, z0,
     .     dzdx0, z1, dzdx1, z2, dzdx2, z3, dzdx3, d2zdxdy
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
        INTEGER MINDEX
C-----------------------------------------------
        ix = mindex(xx,1,nx,x)
        iy = mindex(yy,1,ny,y)
C
        ix = max(1,min(nx-1,ix))
        iy = max(1,min(ny-1,iy))
C
        ix = max(1,min(nx-1,ix))
        iy = max(1,min(ny-1,iy))
C
        ixm1 = max(1,ix-1)
        ixp1 = ix+1
        ixp2 = min(nx,ix+2)
C
        iym1 = max(1,iy-1)
        iyp1 = iy+1
        iyp2 = min(ny,iy+2)
C
        xx0 = xx(ixm1)
        xx1 = xx(ix)
        xx2 = xx(ixp1)
        xx3 = xx(ixp2)
C
        zz00 = zz(ixm1,iym1)
        zz10 = zz(ix  ,iym1)
        zz20 = zz(ixp1,iym1)
        zz30 = zz(ixp2,iym1)
C
        zz01 = zz(ixm1,iy)
        zz11 = zz(ix  ,iy)
        zz21 = zz(ixp1,iy)
        zz31 = zz(ixp2,iy)
C
        zz02 = zz(ixm1,iyp1)
        zz12 = zz(ix  ,iyp1)
        zz22 = zz(ixp1,iyp1)
        zz32 = zz(ixp2,iyp1)
C
        zz03 = zz(ixm1,iyp2)
        zz13 = zz(ix  ,iyp2)
        zz23 = zz(ixp1,iyp2)
        zz33 = zz(ixp2,iyp2)
C
        CALL minter1d_rat(
     +    xx0 ,xx1 ,xx2 ,xx3 ,zz00 ,zz10 ,zz20 ,zz30 ,x ,z0 ,dzdx0,ix,nx)
        CALL minter1d_rat(
     +    xx0 ,xx1 ,xx2 ,xx3 ,zz01 ,zz11 ,zz21 ,zz31 ,x ,z1 ,dzdx1,ix,nx)
        CALL minter1d_rat(
     +    xx0 ,xx1 ,xx2 ,xx3 ,zz02 ,zz12 ,zz22 ,zz32 ,x ,z2 ,dzdx2,ix,nx)
        CALL minter1d_rat(
     +    xx0 ,xx1 ,xx2 ,xx3 ,zz03 ,zz13 ,zz23 ,zz33 ,x ,z3 ,dzdx3,ix,nx)
C
        CALL minter1d_rat(
     +           yy(iym1) ,yy(iy) ,yy(iyp1) ,yy(iyp2),
     +           z0       ,z1     ,z2       ,z3      ,
     +           y        ,z      ,dzdy     ,iy      ,ny)
C
        CALL minter1d_rat(
     +           yy(iym1) ,yy(iy) ,yy(iyp1) ,yy(iyp2),
     +           dzdx0    ,dzdx1  ,dzdx2    ,dzdx3   ,
     +           y        ,dzdx   ,d2zdxdy  ,iy      ,ny)
C
        RETURN