Esope source files

   1 : C RESDPR    SOURCE    CHAT      05/01/13    02:51:37     5004           
   2 : 
   3 :       SUBROUTINE RESIDUO_DPRAL(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM)
   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)
   7 :       integer             augla
   8 :       real*8                    c
   9 :       common /auglagrang1/ augla
  10 :       common /auglagrang2/ c
  10 :       call zzero(sig,4)
  10 :       call zzero(vecm,4)
  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(x(n)+c*r(n))*vecm(i)
  19 :       enddo
  20 : C
  21 : c VISCOSIDAD FICTICIA AÑADIDA
  22 : c      r(n)=r(n)-1.D-8
  23 : c
  24 :       return
  25 :       end 
  26 :  
  27 :  
  28 :