Télécharger la plaquette Cast3M


fr en

Télécharger eau3d.eso

Retour à la liste

Numérotation des lignes : 

eau3d

  1. C EAU3D     SOURCE    FD218221  24/02/07    21:15:10     11834          
  2.        subroutine eau3d(xmat,nmat,vnmat,nstype,var0,varf,nvari,vnvari,
  3.      # ndtens,nbelas,teta1,teta2,tetaref,mfr,keau)
  4.  
  5. c      calcul des coeff hydrique pour mc3d
  6.  
  7. c      declarations externes
  8.        implicit integer(i-n)
  9.        implicit real*8(a-h,o-z)
  10.  
  11.        integer nmat,nstype,ndtens,nvari,nbelas
  12.        integer vnmat(nstype)
  13.        integer vnvari(nstype,2)
  14.        real*8 xmat(nmat)
  15.        real*8 var0(nvari),varf(nvari)
  16.        real*8 teta1,teta2,tetaref,mfr
  17.        real*8 pwcap,bwcap,swcap,keau
  18.  
  19. c      nrj activation pour l eau
  20.        real*8 nrjw 
  21.        parameter (nrjw=17000.d0)
  22.        real*8 biotw,poro2,vw2,xnsat,mvgn
  23.        real*8 pw0,Mwat,dphiw,bw0,bw1
  24.        integer numpw,numbw,numsw
  25.  
  26.  
  27.        call xmat3d(poro2,xmat,nmat,vnmat,nstype,nbelas,1,2)
  28.        call xmat3d(vw2,xmat,nmat,vnmat,nstype,nbelas,1,3)
  29.        call xmat3d(biotw,xmat,nmat,vnmat,nstype,nbelas,1,4)
  30.        call xmat3d(xnsat,xmat,nmat,vnmat,nstype,nbelas,1,5)
  31.        call xmat3d(mvgn,xmat,nmat,vnmat,nstype,nbelas,1,6)
  32.        call xmat3d(tkvg,xmat,nmat,vnmat,nstype,nbelas,1,7)
  33.  
  34.        call xvar3d(pw0,numpw,var0,nvari,vnvari,nstype,ndtens,1,3)
  35.        call xvar3d(bw0,numbw,var0,nvari,vnvari,nstype,ndtens,1,4)
  36.        call xvar3d(sw0,numsw,var0,nvari,vnvari,nstype,ndtens,1,5)
  37.  
  38.        call watr3d(mfr,biotw,poro2,vw2,xnsat,mvgn,pwcap,bwcap,swcap,
  39.      # teta2,tetaref,tkvg,pw0,Mwat,dphiw,bw0,bw1)
  40.  
  41.        varf(numpw)=pwcap
  42.        varf(numbw)=bwcap
  43.        varf(numsw)=swcap
  44.  
  45.        keau=swcap
  46. c       print*,'ds eau3d keau',keau
  47.  
  48.        return
  49.        end
  50.  
  51.  
  52.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales