Esope source files

   1 : C DELSKM    SOURCE    CHAT      05/01/12    22:41:23     5004           
   2 : 
   3 :       SUBROUTINE DETERMINA_LS_KMA(SIG,NESCRI,KMA,IZONE)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer kma,nescri,izone
   6 :       real*8  sig(*),xi1,xj2,valor
   7 :       real*8               nn1,nn2,eta0,aa1,aa2,bb1,bb2,
   8 :      .                     pia,pib,sigy0,pcc,rrr
   9 :       common /poder2_data/ nn1,nn2,eta0,aa1,aa2,bb1,bb2,
  10 :      .                     pia,pib,sigy0,pcc,rrr
  10 :       call Invari_J2(sig,3,xj2)
  10 :       call Invari_I1(sig,3,xi1)
  11 :       if (xi1.lt.pia) then 
  12 :          kma=2
  13 :          valor=rrr**2*xj2+(xi1-pia)**2-(pia-pib)**2
  14 :          valor=xj2+((xi1-pia)/rrr)**2-aa2*sigy0**2/3.D0
  15 :          if (nescri.eq.1) 
  16 :      .   write(199,*)' CAP ',xi1,xj2,
  17 :      .   pia,pib,rrr,((xi1-pia)/rrr)**2,aa2,valor
  18 :          izone = 0
  19 :       else
  20 :          kma=3
  21 :          valor=dsqrt(2.D0*xj2)+bb1*xi1/3.D0-bb2
  22 :          if (nescri.eq.1) 
  23 :      .   write(199,*)' CONE',xi1,pia,valor,kma
  24 :          izone = 1
  25 :       endif
  26 :       return
  27 :       end
  28 :  
  29 :  
  30 :