00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #ifndef __bedrockOnlySystem_hh
00020 #define __bedrockOnlySystem_hh
00021
00022 #include "columnSystem.hh"
00023 #include "NCVariable.hh"
00024
00026
00032 class bedrockOnlySystemCtx : public columnSystemCtx {
00033
00034 public:
00035 bedrockOnlySystemCtx(const NCConfigVariable &config, int my_Mbz);
00036 ~bedrockOnlySystemCtx();
00037
00038 PetscErrorCode initAllColumns(
00039 const PetscScalar my_dtTemp, const PetscScalar my_dzbEQ);
00040
00041 PetscErrorCode setBoundaryValuesThisColumn(
00042 const PetscScalar my_Tbedtop, const PetscScalar my_Ghf);
00043
00044 PetscErrorCode viewConstants(PetscViewer viewer, bool show_col_dependent);
00045
00046 PetscErrorCode solveThisColumn(PetscScalar **x);
00047 PetscScalar extractHeatFluxFromSoln(const PetscScalar *x);
00048
00049 public:
00050
00051 PetscScalar *Tb;
00052
00053 private:
00054 PetscInt Mbz;
00055 PetscScalar dtTemp, dzbEQ, bed_rho, bed_c, bed_k, bedK, bedR, Ghf, Tbedtop;
00056 };
00057
00058 #endif // ifndef __bedrockOnlySystem_hh
00059
00060