00001 #ifndef RICHARDSCONTAM_H
00002 #define RICHARDSCONTAM_H
00003
00004 #include "genfile.h"
00005 #include "richards.h"
00006
00007 class richardscontam
00008 {
00009 public:
00010 richardscontam();
00011 ~richardscontam();
00012
00013 void read(XFILE *in);
00014 void print(FILE *out);
00015
00016 double cc_value (double h);
00017
00018 void matcond (matrix &d,long ri,long ci,long ipp);
00019 void matcap (double &c,long ri,long ci,long ipp);
00020
00021 void matcond1d (matrix &d,long ri,long ci,long ipp);
00022 void matcond2d (matrix &d,long ri,long ci,long ipp);
00023 void matcond3d (matrix &d,long ri,long ci,long ipp);
00024
00025 void matcond2 (matrix &d,long ri,long ci,long ipp);
00026 void matcond1d_2 (matrix &d,long ri,long ci,long ipp);
00027 void matcond2d_2 (matrix &d,long ri,long ci,long ipp);
00028 void matcond3d_2 (matrix &d,long ri,long ci,long ipp);
00029
00030
00031
00032 double alpha;
00033 double n;
00034 double m;
00035
00036
00037 long dim;
00038
00039
00040 double kksxx,kksxy,kksxz,kksyy,kksyz,kkszz;
00041
00042 double thetas;
00043
00044 double thetar;
00045
00046 double storage;
00047
00048
00049 richards rich;
00050 };
00051
00052 #endif