Esope source files

   1 : C RES2LV    SOURCE    CHAT      05/01/13    02:51:15     5004           
   2 : 
   3 :       SUBROUTINE RESIDUO_2LEVELS(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM,LAM)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n,i,j,ndims,ndimv,nmodel,ndimx
   6 :       real*8  x(n),r(n),void,sig(4),vecm(4),xtri(ndimx),lam
   7 :       integer             augla
   8 :       real*8                    c
   9 :       common /auglagrang1/ augla
  10 :       common /auglagrang2/ c
  11 :       void=0.D0
  12 :       ndimv=1
  13 :       ndims=3
  13 :       call der_enerelas_dpral(x,sig,nmodel)
  14 :       if (ndimx.eq.3) then
  14 :        call yieldd(sig,ndims,void,ndimv,r(n),nmodel)
  14 :        call vflsig(sig,ndims,void,ndimv,vecm,nmodel)
  15 :       else if (ndimx.eq.4) then 
  15 :        call yieldd(sig,ndims,x(ndimx),ndimv,r(n),nmodel)
  15 :        call vflsig(sig,ndims,x(ndimx),ndimv,vecm,nmodel)
  15 :        call vflvar(sig,ndims,x(ndimx),ndimv,vecm(ndimx),nmodel)
  16 :       endif
  17 :       do i=1,ndimx
  18 :          r(i)=x(i)-xtri(i)+dabs(lam+c*r(n))*vecm(i)
  19 :       enddo
  20 :       return
  21 :       end 
  22 :  
  23 :  
  24 :