module12.f

油田化学驱模拟的经典 fortran源代码

C
        MODULE MODULE12
        DOUBLE PRECISION   AKWC7,AKWS7,AKM7,AK7,PT7 
     &                    ,AKWC8,AKWS8,AKM8,AK8,PT8
        DOUBLE PRECISION   EPSALC
        INTEGER     NALMAX
        INTEGER     NALK1,NALK2,NALK3,NALK4,NALK5,NALK6
        INTEGER     IMES
        DOUBLE PRECISION   HBNS70,HBNC70,HBNS71,HBNC71,HBNS72
     &                    ,HBNC72,HBNS80
        DOUBLE PRECISION   HBNC80,HBNS81,HBNC81,HBNS82,HBNC82
        INTEGER     IBKIN,IBNONB,NBC,NBS,NBCNOB,NBIOEQ,NRLIM,NMET
        INTEGER     NBCAQ,NBCNAQ,IBINAQ,IBFNAQ,IBIAQ,NAPTOT
        DOUBLE PRECISION   EPSBIO
        DOUBLE PRECISION   CMIN
        DOUBLE PRECISION   DENBLK,DIAMP
        INTEGER     NARTOT
        DOUBLE PRECISION   SUBMAX,EAMAX,RMTMAX,CEXIT
        INTEGER     IEND
        INTEGER     LIML,LIMU
        INTEGER     IDXYZ
        DOUBLE PRECISION   ZERO,ONE,ONEM,ONEM4,ONEM5,ONEM6,ONEM7 
     &                    ,ONEM8,ONEM9
        DOUBLE PRECISION   ONEM10,ONEM12,ONEP12,ONEM50  
     &                    ,ONEP50,ONEM5M,PONEM,ONE199,PRCSN
        DOUBLE PRECISION   PIE,F1P8
        INTEGER     IREACT,ICAP,ICWM,IMODE,IBIO
        DOUBLE PRECISION   DCLIM,DTMIN,DTMAX,CNMIN,CNMAX,TMAX,AVGIT
        INTEGER     ITSTEP,ISTOP
        DOUBLE PRECISION   D111,THETAX,THETAY
        INTEGER     IDEPTH
        DOUBLE PRECISION   OPSK7O,OPSK7S,OPSK8O,OPSK8S,BETA6
     &                    ,BETA7,BETA8
        DOUBLE PRECISION   WATLOG
        DOUBLE PRECISION   AG1,AG2,CRG,AGK,BGK
        DOUBLE PRECISION   A15D,B15D
        DOUBLE PRECISION   CRNAK,HNAK,A14D,B14D,C160
        INTEGER     ICREX
        DOUBLE PRECISION   AK1,AK2,SCR,X4,X14,X16,WM4
        DOUBLE PRECISION   AK1T,AK2T,X13
        INTEGER     KGOPT
        INTEGER     NH,NNA,NCA,NMG,NCARB,NACD,NALU,NSILI,NOXY      
        DOUBLE PRECISION   EPXLO,SOSTAR,CSTAR,RFMAX,SHRTN,VELGR
        DOUBLE PRECISION   S4RC,P4RC,E4C,T44
        DOUBLE PRECISION   SGI
        INTEGER     IGAS
        DOUBLE PRECISION   VIS4,VSLOPG,XIFTG,XIFTGW,S2RC4
        INTEGER     NGP1,NGP2,NGP3,NGP4,NG,NGP5
        DOUBLE PRECISION   EPHI1
        DOUBLE PRECISION   TQLOS
        DOUBLE PRECISION   TCONO,TCONU
        INTEGER     IHLOS,IANAL
        DOUBLE PRECISION   CRTC,CVSPR,DENS,CUMHI,CUMHP,TEMPI
        INTEGER     IENG
        DOUBLE PRECISION   DENO,DENU,CVSPO,CVSPU
        DOUBLE PRECISION   TEMREF
        DOUBLE PRECISION   HBNT0,HBNT1,HBNT2,CSET
        DOUBLE PRECISION   C2PLC,C2PRC,B1,B2,B3,F1,F2,F3 
     &                    ,G11,G12,G13,G21,G22,G23
        DOUBLE PRECISION   XIFTW,CHUH,AHUH
        INTEGER     IFT,IHAND
        DOUBLE PRECISION   ALPHA,BETA
        INTEGER     IJK1,NP
        DOUBLE PRECISION   A,B,G,F,XR,XL
        INTEGER     IFX
        DOUBLE PRECISION   BRK,CRK,VIS1,VIS2,AP1,AP2,AP3,GAMMAC
     &                    ,EPHI4,EPHI3
        DOUBLE PRECISION   GAMHF,SSLOPE,POWN,CSE1,BETAP,TSTAND
        INTEGER     IPOLYM
        INTEGER     IERROR,IONCNV,NEGSAT
        DOUBLE PRECISION   T11,T22,T33
        DOUBLE PRECISION   P1RC,P2RC,P3RC
        DOUBLE PRECISION   S1RC,S2RC,S3RC,E13C,E23C,E31C
        INTEGER     ITRAP
        INTEGER     IPERM
        INTEGER     IPRESS
        DOUBLE PRECISION   DT,CURANT
        INTEGER     NXM1,NXNY,NBL,NBLW
        DOUBLE PRECISION   DEN1,DEN2,DEN23,DEN3,DEN7,DEN8
        INTEGER     IDEN
        INTEGER     IOW
        DOUBLE PRECISION   AD31,AD32,B3D,FADS,REFK
        INTEGER     IADS1,IADK
        DOUBLE PRECISION   AD41,AD42,B4D
        DOUBLE PRECISION   PSTAND,COMPR
        INTEGER     IHYST,IPARK 
        DOUBLE PRECISION   QV,XKC,XKS,EQW
        INTEGER     IXYZ
        INTEGER     IP1,IP2
        DOUBLE PRECISION   VHS,VHT,TDLN,CTD
        INTEGER     IPRFL
        INTEGER     IVARI,IVARP
        DOUBLE PRECISION   TEMPOB,TEMPUB
        INTEGER     MTWM1,MTAM1,MGLM1
        LOGICAL LTW,LTA,LGC,LACD,LGL
        INTEGER     MGCM1,NGC
        LOGICAL LCAL2
        INTEGER     MMOM1
        LOGICAL  LMO
        DOUBLE PRECISION   SSC,EMIN,SDS,BDS,SS0,BS0,SS1
     &                    ,BS1,SS2,BS2
        LOGICAL  LCOMP
        INTEGER     IOVIS
        INTEGER     IFGHBN
        DOUBLE PRECISION   CUMPR1,CUMHI1,WRHPV,WRPRF,RSTC
        INTEGER     ICOPSM,ICUMTM
        INTEGER     IADS,ICKL,IVEL,IVIS,IPER,ICNM,IRKF,IPHSE
        INTEGER     ICSE
        INTEGER     ITS,ISTATE,IOUTGMS,IGMSS,IGMSE,ITUNIT
        INTEGER     IPPRES,IPSAT,IPCTOT,IPBIO,IPCAP,IPGEL,IPALK
        INTEGER     IPTEMP,IPOBS,NOBS
        DOUBLE PRECISION   FOREC,FOREC1,RESPV
        DOUBLE PRECISION   PRESS1,PINIT,HINIT
        DOUBLE PRECISION   RESPVI,RESPVM
        INTEGER     NELET,NFLD,NSLD,NSORB,NACAT,NIAQ,NEX
        INTEGER     NFLASH,NSLEL,NIND,NIND1
        INTEGER     ICHRGE,NSURF1,NSURF2,IRSPS
        INTEGER     INFLD,INSORB,ININD1,ININD,INFLSH
        DOUBLE PRECISION   CAC2I,ACIDIS
        DOUBLE PRECISION   C50,C60
        DOUBLE PRECISION   C1I,C2I
        INTEGER     NONSLD,NRSLD,NITER,IDAMP
        DOUBLE PRECISION   CONFCT,CANION,CSURF,CSURFA,WTRINV,FVBPV
        DOUBLE PRECISION   EQWPS,CSELP,CSEUP,CMEQ2
        INTEGER     NACTEM,NINDT
        DOUBLE PRECISION   PHC,PHT,PHT1,HPHAD
        INTEGER     IPHAD  
        INTEGER     NALK7
        INTEGER     NALK8,NHFD,NCRFD
        INTEGER     ITCUM    
        DOUBLE PRECISION   T,TINJ,WHPV
        DOUBLE PRECISION   PRF
        INTEGER     ICNT,IINJ,INEC,IRST
        DOUBLE PRECISION   DCMAX
        INTEGER     IDISPC,ICOORD,ITC,IUNIT
        INTEGER     NBLHV
        INTEGER     ITREAC,NRT
        DOUBLE PRECISION   C5INI
        DOUBLE PRECISION   CUMPV,CSEL7,CSEL8,CSEU7,CSEU8           
        INTEGER     ISHIFT,NNFLA
        DOUBLE PRECISION   CUMPRO,CUMHIS,TIRST,RSTCON 
        DOUBLE PRECISION   EPSME
        INTEGER     IALC
        DOUBLE PRECISION   SWI
        DOUBLE PRECISION   DCMIN
        INTEGER     ICOUNT
        INTEGER     ISTB
        INTEGER     NWREL
        INTEGER     IRO
	DOUBLE PRECISION   WSOL
        INTEGER     IMASS,ISOL,ICOR
	DOUBLE PRECISION   EQB0,EQB1,EQB2,EQB3
        INTEGER     IBOUND,IBL,IBR,IBT,IZONE
        DOUBLE PRECISION   ENTP,ALPDS,ALPWS,VNDS,VNWS,PNEG,CLAMDA
        INTEGER     IADSO,IREV
        INTEGER     IBPP
        INTEGER     IBTIM,IMAUTO,IMTVAR,NBTS 
        DOUBLE PRECISION   BTMAX,BTMIN,BIOTIM,BIORME,EFMIN,DAMX,BTSAVG
        INTEGER      NXZ,NXZP,NWBC1,NPHASE1,NNOB,NBEQ
        INTEGER IOUHS,IOUTOT
        INTEGER IHYSTP,IFOAMP,INONEQ
        INTEGER ICWI
C
      DOUBLE PRECISION, ALLOCATABLE ::  X7OLD(:),X8OLD(:)
      DOUBLE PRECISION, ALLOCATABLE ::  D(:,:),ALPHAL(:),ALPHAT(:)
      DOUBLE PRECISION, ALLOCATABLE ::  AREAXZ(:),HX(:),HZ(:),HXL(:)
     & ,HZL(:),HXB(:),HZB(:)
      DOUBLE PRECISION   DBACN(14)
      DOUBLE PRECISION, ALLOCATABLE ::  CB(:,:,:),BIOMIN(:,:),BIOCUM(:)
      DOUBLE PRECISION, ALLOCATABLE ::  AKA(:),AKN(:,:),AKS(:) 
     * ,BRMAX(:),BRMAXB(:) 
     * ,BSIHB(:,:),CBIOMN(:),COLMAS(:) 
     * ,COLNUM(:),COLSA(:),ENDOG(:) 
     * ,ENDOGB(:),FEA(:),FN(:,:) 
     * ,FP(:,:),FPABIO(:,:),RABIO(:),RCOL(:) 
     * ,TCOL(:),VCOL(:),YXS(:)
      INTEGER, ALLOCATABLE ::  ICSUB(:,:) 
     * ,IDMET(:,:,:),IPABIO(:,:) 
     * ,IRABIO(:),NCOMPS(:),NIHB(:),NNUT(:) 
     * ,NPABIO(:),NPROD(:)
      INTEGER, ALLOCATABLE ::  IMSUB(:),IMEA(:),IMBS(:) 
     * ,IHB(:,:),IPR(:,:),INUT(:,:) 
     * ,IKCB(:),IBIOC(:),IBS(:) 
      DOUBLE PRECISION, ALLOCATABLE ::  BMTC(:),SC(:)
      DOUBLE PRECISION, ALLOCATABLE ::  CD(:,:,:),SD(:,:),SF(:,:) 
     * ,CM(:,:),CTF(:,:),FFL(:),FFH(:),SLP(:)
      INTEGER, ALLOCATABLE ::  ICAPFL(:)
      DOUBLE PRECISION, ALLOCATABLE ::  PR(:),POLD(:),PORC(:)
     * ,CORFP(:),CORFF(:)
      DOUBLE PRECISION, ALLOCATABLE ::  AKR(:),FRC(:),FRP(:)