00001 #ifndef CONSOLAWF2MAT_H
00002 #define CONSOLAWF2MAT_H
00003
00004 #include "genfile.h"
00005
00006 class con_awf2mat
00007 {
00008 public:
00009 con_awf2mat();
00010 ~con_awf2mat();
00011
00012 void read(XFILE *in);
00013 void print(FILE *out);
00014
00015 void matcond (matrix &d,long ri,long ci,long ipp);
00016 void matcap (double &c,long ri,long ci,long ipp);
00017
00018 void matcond1d (matrix &d,long ri,long ci,long ipp);
00019 void matcond2d (matrix &d,long ri,long ci,long ipp);
00020 void matcond3d (matrix &d,long ri,long ci,long ipp);
00021
00022 void rhs_volume (matrix &d,long ri,long ci,long ipp);
00023 void rhs1d1 (matrix &d,long ri,long ci,long ipp);
00024 void rhs2d1 (matrix &d,long ri,long ci,long ipp);
00025 void rhs3d1 (matrix &d,long ri,long ci,long ipp);
00026
00027 double get_sw(double pw, double pg);
00028 double get_cs(double pw, double pg);
00029 double get_se(double sw);
00030 double get_krl(double pw, double pg);
00031 double get_krg(double pw, double pg);
00032
00033 double get_kww(double pw, double pg);
00034 double get_kwg(double pw, double pg);
00035 double get_kgw(double pw, double pg);
00036 double get_kgg(double pw, double pg);
00037
00038 double get_capww(double pw, double pg);
00039 double get_capwg(double pw, double pg);
00040 double get_capgw(double pw, double pg);
00041 double get_capgg(double pw, double pg);
00042
00043 double get_fw1(double pw, double pg);
00044 double get_fg1(double pw, double pg);
00045
00046 double get_othervalue(long compother,double pw, double pg, long ipp);
00047 void print_othervalue_name(FILE *out,long compother);
00048 void values_correction (vector &nv, long ipp);
00049 void water_pressure_check(double &pw,double &pg,long ipp);
00050
00051
00052 private:
00053
00054 double alpha,ks,phi0,kw,rhow,muw0,kintr,mug;
00055
00056 };
00057
00058 #endif