Hi Bernd, thanks for the time you spend with my questions it helps me to move on. I have done what you suggested to me and it worked well, however I still have few errors : tutorialproblem_coupled.hh:166:16: erreur: expected initializer before 'ScalarField' tutorialproblem_coupled.hh:199:23: erreur: 'Kxx' was not declared in this scope tutorialproblem_coupled.hh:206:52: erreur: 'Kxx' was not declared in this scope I checked with heterogeneousproblem.hh but the declaration of the variables are identical to my current file. tutorialproblem_coupled.hh:198:102: error: cannot convert 'const Dune::FieldMatrix<double, 3, 3>' to 'Dumux::TutorialProblemCoupled<Dumux::Properties::TTag::TutorialProblemCoupled>::Scalar' in initialization If I understand well this one, there is a dimension problem for a variable (my case is declared in 3D in the heterogeneousprolem.hh) ? Thanks for the support. Have a nice day, Cheers, Romain.
________________________________ De : dumux-boun...@listserv.uni-stuttgart.de [mailto:dumux-boun...@listserv.uni-stuttgart.de] De la part de Bernd Flemisch Envoyé : vendredi 18 janvier 2013 16:44 À : DuMuX User Forum Objet : Re: [DuMuX] Print out permeability (with Paraview) via addOutputVtkFields method Hi Romain, please have a look at the compiler messages: tutorialproblem_coupled.hh:173:14: erreur: déclaration ne déclarant rien du tout tutorialproblem_coupled.hh:173:33: erreur: 'fvElemGeom' was not declared in this scope tutorialproblem_coupled.hh:174:14: erreur: 'VolumeVariables' was not declared in this scope tutorialproblem_coupled.hh:174:31: erreur: 'volVars' was not declared in this scope The first semicolon in each of the lines FVElementGeometry; fvElemGeom; VolumeVariables; volVars; is too much. tutorialproblem_coupled.hh:176:14: erreur: 'ElementIterator' was not declared in this scope tutorialproblem_coupled.hh:176:30: erreur: expected ';' before 'elemIt' tutorialproblem_coupled.hh:177:30: erreur: expected ';' before 'elemEndIt' tutorialproblem_coupled.hh:178:21: erreur: 'elemIt' was not declared in this scope tutorialproblem_coupled.hh:178:31: erreur: 'elemEndIt' was not declared in this scope You have to declare the type "ElementIterator." You can copy the typedef from heterogeneousproblem.hh. tutorialproblem_coupled.hh:186:38: erreur: 'numVerts' was not declared in this scope You have to uncomment the line defining numVerts (184?). tutorialproblem_coupled.hh:198:22: erreur: 'Kzz' was not declared in this scope You have to comment this line since you only want to provide Kxx. Kind regards Bernd .. 1- But I don't how to declare variables and where to do it. Because in the " "test/boxmodels/co2/heterogeneousproblem.hh.", variables are only declared in the routine and it compiles but not when I copy it in "tutorialproblem_coupled.hh". 2- If I understood well it is from intrinsicPermeability that my permeability will come from "tutorialspatialparams.hh" ? Thanks to the community, Cheers, Romain. ps : the addOutputVtkFields : void addOutputVtkFields() { typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField; unsigned numVertices = this->gridView().size(dim); // //create required scalar fields ScalarField *Kxx = this->resultWriter().allocateManagedBuffer(numVertices); // ScalarField *boxVolume = this->resultWriter().allocateManagedBuffer(numVertices); // // // // (*boxVolume) = 0; // // //Fill the scalar fields with values unsigned numElements = this->gridView().size(0); ScalarField *rank = this->resultWriter().allocateManagedBuffer(numElements); // FVElementGeometry; fvElemGeom; VolumeVariables; volVars; ElementIterator elemIt = this->gridView().template begin<0>(); ElementIterator elemEndIt = this->gridView().template end<0>(); for (; elemIt != elemEndIt; ++elemIt) { int idx = this->elementMapper().map(*elemIt); (*rank)[idx] = this->gridView().comm().rank(); fvElemGeom.update(this->gridView(), *elemIt); // // int numVerts = elemIt->template count<dim> (); // for (int i = 0; i < numVerts; ++i) { int globalIdx = this->vertexMapper().map(*elemIt, i, dim); volVars.update(this->model().curSol()[globalIdx], *this, *elemIt, fvElemGeom, i, false); // (*boxVolume)[globalIdx] += fvElemGeom.subContVol[i].volume; Scalar perm = this->spatialParams().intrinsicPermeability(*elemIt, fvElemGeom, i); (*Kxx)[globalIdx] = perm; (*Kzz)[globalIdx] = perm[dim-1][dim-1]; } } //pass the scalar fields to the vtkwriter // this->resultWriter().attachVertexData(*boxVolume, "boxVolume"); this->resultWriter().attachVertexData(*Kxx, "Kxx"); }; __________________________ Avant d'imprimer, pensez à l'environnement ! Please consider the environment before printing ! Ce message et toutes ses pièces jointes sont confidentiels et établis à l'intention exclusive de ses destinataires. Toute utilisation non conforme à sa destination, toute diffusion ou toute publication, totale ou partielle, est interdite, sauf autorisation expresse. IFP Energies nouvelles décline toute responsabilité au titre de ce message. This message and any attachments are confidential and intended solely for the addressees. Any unauthorised use or dissemination is prohibited. IFP Energies nouvelles should not be liable for this message. __________________________ _______________________________________________ Dumux mailing list Dumux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux -- _____________________________________________________________________ Bernd Flemisch phone: +49 711 685 69162 IWS, Universität Stuttgart fax: +49 711 685 60430 Pfaffenwaldring 61 email: be...@iws.uni-stuttgart.de D-70569 Stuttgart url: www.hydrosys.uni-stuttgart.de _____________________________________________________________________ __________________________ Avant d'imprimer, pensez à l'environnement ! Please consider the environment before printing ! Ce message et toutes ses pièces jointes sont confidentiels et établis à l'intention exclusive de ses destinataires. Toute utilisation non conforme à sa destination, toute diffusion ou toute publication, totale ou partielle, est interdite, sauf autorisation expresse. IFP Energies nouvelles décline toute responsabilité au titre de ce message. This message and any attachments are confidential and intended solely for the addressees. Any unauthorised use or dissemination is prohibited. IFP Energies nouvelles should not be liable for this message. __________________________
_______________________________________________ Dumux mailing list Dumux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
