REAL FUNCTION cut4()
      include ?
      REAL C1,C2,C3,C4,C5,C6,C7,C8,C9,C10,C11,PROD
      REAL D1,D2,D3,D4,D5
      REAL E1,E2
      REAL F1,F2
      REAL cute
*
      C1=-0.0137*ph_tgth+0.0502-th_tgth
      C2=0.0502+th_tgth

      C3=-0.136*y_tgth+0.02518-0.991*ph_tgth
      C4=y_tgth*0.136+0.02518+0.991*ph_tgth

      C5=-0.975*dph+0.0396-0.220*th_tgth
      C6=dph+0.0409+0.068*th_tgth
      C9=-0.421*dph-0.907*th_tgth+0.0561
      
      C7=-dph*0.1776+0.984*ph_tgth+0.0237
      C8=-dph*0.1738+0.0236-0.985*ph_tgth

      C10=0.319* y_tgth-0.948*ph_tgth+0.0336
      C11=-0.319*y_tgth+0.948*ph_tgth+0.0336

      D1=PROD(C1,C2)
      D2=PROD(C3,C4)
      D3=PROD(C5,C6)
      D3=PROD(C9,D3)
      D4=PROD(C7,C8)
      D5=PROD(C10,C11)

      E1=PROD(D1,D2)
      E2=PROD(D3,D4)
      E2=PROD(E2,D5)
      
      F1=PROD(E1,E2)
      F2=cute(y_tgte,ph_tgte,th_tgte,dpe)
      cut4=PROD(F1,F2)
      END

      REAL FUNCTION cute(y_tgth,ph_tgth,th_tgth,dph)
      REAL C1,C2,C3,C4,C5,C6,C7,C8,C9,C10,C11,PROD
      REAL D1,D2,D3,D4,D5
      REAL E1,E2
      REAL y_tgth,ph_tgth,th_tgth,dph
*
      C1=-0.0137*ph_tgth+0.0502-th_tgth
      C2=0.0502+th_tgth

      C3=-0.136*y_tgth+0.02518-0.991*ph_tgth
      C4=y_tgth*0.136+0.02518+0.991*ph_tgth

      C5=-0.975*dph+0.0396-0.220*th_tgth
      C6=dph+0.0409+0.068*th_tgth
      C9=-0.421*dph-0.907*th_tgth+0.0561
      
      C7=-dph*0.1776+0.984*ph_tgth+0.0237
      C8=-dph*0.1738+0.0236-0.985*ph_tgth

      C10=0.319* y_tgth-0.948*ph_tgth+0.0336
      C11=-0.319*y_tgth+0.948*ph_tgth+0.0336

      D1=PROD(C1,C2)
      D2=PROD(C3,C4)
      D3=PROD(C5,C6)
      D3=PROD(C9,D3)
      D4=PROD(C7,C8)
      D5=PROD(C10,C11)
      
      E1=PROD(D1,D2)
      E2=PROD(D3,D4)
      E2=PROD(E2,D5)

      cute=PROD(E1,E2)
      END

      REAL FUNCTION PROD(X,Y)
      REAL X,Y
c      PROD=X+Y-SQRT(X*X+Y*Y+2*0.9*X*Y)
c      PROD=X+Y-SQRT(X*X+Y*Y)
      PROD=MIN(X,Y)
      RETURN
      END