Thanks Nori for your customized file. I could make the compiler happy
after revising the following.
in fem_ele_std.h
public:
void Config();
in vtk.cpp
/* PCH Comment out just for compilation
case 3:
mat_value = mmp->FractureAperture(i_e);
break;
case 4:
mat_value = mmp->FractureContactAreaRatio(i_e);
break;
*/
And with Georg suggestion of permeability and PVD output option as you
suggested.
in .out file
$ELE_VALUES
PERMEABILITY
That does not make it work. It prints all zero instead of 1.e-11 I have
assigned. Just for your info and thanks. In fact, I prefer the PVD
option as you may notice it from our conversation here in the google group.
Cheers,
Chan-Hee
On Monday, April 1, 2013 7:52:08 PM UTC+9, Norihiro Watanabe wrote:
Chan-Hee,
Attached is my ad-hoc version of PVD output supporting some MMP and MFP
properties as element data (evaluated at the 1st gauss point). you
cannot compile it as it is because some MMP functions exist only in my
branch. but you can simply remove lines to fix errors. Viscosity is not
supported but you can easily add it.
best,
nori
On 04/01/2013 05:28 AM, Chan-Hee Park wrote:
> Thanks Marc and Georg,
>
> I have tried permeability output in vtk with the latest trunk. It
works
> with VTK output as Georg mentioned. I didn't work with VISOCITY
as an
> ELE_VALUES though. Meanwhile, "PVD" output did not work with
> PERMEABILITY, either. Currently, VTK works only for PERMEABILITY
by the
> following lines.
>
> // PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR? //
JTARON 2010
> else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)
> {
> vtk_file << "VECTORS permeability double " << "\n";
> for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)
> {
> MeshLib::CElem* ele = _mesh->ele_vector[j];
> CMediumProperties* MediaProp =
> mmp_vector[ele->GetPatchIndex()];
> for (size_t i = 0; i < 3; i++)
> vtk_file <<
> MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";
> vtk_file << "\n";
> }
>
> As for VISCOSITY, I guess this is a good starting point how I
revise the
> code. First, I have to choose whether I put it as node or element
value.
> Thanks.
>
> On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:
>
> Hoi,
> permeability output in vtk legacy files is supported:
> $ELE_VALUES
> PERMEABILITY
>
> Other properties need probably modifications of the source
code.
> Greetings, Georg
>
> Am Freitag, 29. M�rz 2013 01:29:35 UTC+1 schrieb Chan-Hee > Park:
>
> Hi guys,
>
> Just have a quick sit to print out material(i.e.,
permeability)
> or fluid properties (i.e., viscosity, density) of
elements(or
> nodes) in the output file. Anybody has an example
configured for
> this or should I code to do it?
>
> Cheers,
>
> Chan-Hee
>
> --
> You received this message because you are subscribed to the Google
> Groups "ogs-users" group.
> To unsubscribe from this group and stop receiving emails from it,
send
> an email to ogs-users+...@googlegroups.com <javascript:>.
> For more options, visit https://groups.google.com/groups/opt_out
<https://groups.google.com/groups/opt_out>\.
>
--
Dr.-Ing. Norihiro Watanabe
Department of Environmental Informatics (ENVINF)
Wissenschaftler
Helmholtz-Zentrum f�r Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstra�e 15 / 04318 Leipzig / Germany
norihiro...@ufz.de <javascript:> / http://www.ufz.de
Telefon +49 341 235 1806
Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Gesch�ftsf�hrer: Prof. Dr. Georg Teutsch
Administrativer Gesch�ftsf�hrer: Dr. Heike Gra�mann
--
You received this message because you are subscribed to the Google
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an email to ogs-users+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out\.