MINICERN/mathlib/gen/d/epdloc.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:18  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE EPDLOC(IR,XM,XBPR,NBPR,UBPR,HL,HR,UL,UR,LAMDA,DX)
  11 C
  12 C     LOCATION OF MESH-POINT ROW-WISE OR COLUMN-WISE.
  13 C
  14       DIMENSION XBPR(200,20),UBPR(200,20),NBPR(200)
  15       HL=1.0
  16       HR=1.0
  17       UL=0.0
  18       UR=0.0
  19       LAMDA=0
  20       N=NBPR(IR)
  21       IF(N)30,30,3
  22     3 J=1
  23     4 IF(XM-XBPR(IR,J)) 10,10,5
  24     5 IF(J-N)6,8,8
  25     6 J=J+1
  26       GO TO 4
  27     8 HR=1.0
  28       J=J+1
  29       GO TO 16
  30    10 HR=MIN((XBPR(IR,J)-XM)/DX,1.0)
  31       IF(HR-1.0)12,13,12
  32    12 UR=UBPR(IR,J)
  33    13 IF(J-1)16,15,16
  34    15 HL=1.0
  35       GO TO 20
  36    16 HL=MIN((XM-XBPR(IR,J-1))/DX,1.0)
  37       IF(HL-1.0)18,20,18
  38    18 UL=UBPR(IR,J-1)
  39    20 IF((J/2)*2-J)24,22,24
  40    22 LAMDA=2
  41       RETURN
  42    24 IF(HL-1.0)28,26,28
  43    26 IF(HR-1.0)28,30,28
  44    28 LAMDA=1
  45    30 RETURN
  46       END