microstructure/src/mstruct_parametres.cpp

#include "mstruct_parametres.h"
#include "mstruct_definition.h"
#include "mstruct_generateur_rsa.h"
#include "mg_export.h"
#include <string.h>
#include <fstream>
#include <cstdlib>


void PARAMETRES::ini_param_generateur_rsa(OT_PARAMETRES* params, int type_inclusion)
{
  params->ajouter("Type_generateur","0",OT_PARAMETRES::DOUBLE,"Type de generateur : RSA [0], DCR [1]");
  params->ajouter("Nom_mgcg_modele","VES",OT_PARAMETRES::STRING,"Nom du MG_CG_MODELE");
  params->ajouter("Nom_groupe_inclusion","Particule",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion (Ex: Particule)");  
  params->ajouter("Boite3D_distribution_Xmin","-0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Ymin","-0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Zmin","-0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Xmax","1.1",OT_PARAMETRES::DOUBLE,"Xmax boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Ymax","1.1",OT_PARAMETRES::DOUBLE,"Ymax boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Zmax","1.1",OT_PARAMETRES::DOUBLE,"Zmax boite3D de distribution d'inclusion");  
  params->ajouter("Nb_pas_X","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon X pout la creation d'un octree");
  params->ajouter("Nb_pas_Y","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Y pout la creation d'un octree");
  params->ajouter("Nb_pas_Z","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Z pout la creation d'un octree");  
  params->ajouter("Fraction_volumique_cible","0.1",OT_PARAMETRES::DOUBLE,"Fraction volumique cible d'inclusion [0.0, 1.0]");
  params->ajouter("Eps_fraction_volumique","0.001",OT_PARAMETRES::DOUBLE,"Epsilon fraction volumique");
  params->ajouter("Avec_intersections","1",OT_PARAMETRES::DOUBLE,"Avec ou sans intersections avec le bord du VES : Sans [0] Avec [1]");  
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::SPHERE)
  {    
    params->ajouter("Type_inclusion","0",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
    params->ajouter("Mu_rayon","0.09",OT_PARAMETRES::DOUBLE,"Moyenne du rayon de la sphere");
    params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon de la sphere");
    params->ajouter("Type_distribution_rayon","0",OT_PARAMETRES::DOUBLE,"Type de distribution du rayon : Fixe[0], Uniforme[1], Normale[2]");
  }
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    params->ajouter("Type_inclusion","1",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
    params->ajouter("Mu_rayon","0.05",OT_PARAMETRES::DOUBLE,"Moyenne du rayon du cylindre");
    params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon du cylindre");
    params->ajouter("Type_distribution_rayon","0",OT_PARAMETRES::DOUBLE,"Type de distribution du rayon : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_longueur","0.4",OT_PARAMETRES::DOUBLE,"Moyenne de la longueur du cylindre");
    params->ajouter("Sigma_longueur","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la longueur du cylindre");
    params->ajouter("Type_distribution_longueur","0",OT_PARAMETRES::DOUBLE,"Type de distribution de la longueur : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_theta","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle theta");
    params->ajouter("Sigma_theta","1.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle theta");
    params->ajouter("Type_distribution_theta","1",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle theta : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_phi","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle phi");
    params->ajouter("Sigma_phi","1.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle phi");
    params->ajouter("Type_distribution_phi","1",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle phi : Fixe[0], Uniforme[1], Normale[2]");
  }
  params->ajouter("Porosite","0",OT_PARAMETRES::DOUBLE,"Considerer l'inclusion comme etant une porosite (vide) [1]");
  params->ajouter("Multicouche","0,(0.01,0.01)",OT_PARAMETRES::DOUBLE,"Generer une inclusion avec une ou plusieurs couches [1,2,3,etc] et d'epaisseurs (0.01,0.02,etc)");
  params->ajouter("Distance_min_inter_volume","0.04",OT_PARAMETRES::DOUBLE,"Distance minimale qui separe deux inclusions");
  params->ajouter("Volume_min","0.0008",OT_PARAMETRES::DOUBLE,"Volume minimal admissible de l'inclusion");
  params->ajouter("Aire_min","0.001963",OT_PARAMETRES::DOUBLE,"Aire minimal admissible d'une face de l'inclusion");
  params->ajouter("Longueur_min","0.02",OT_PARAMETRES::DOUBLE,"Longueur minimal admissible d'une arete ");
  params->ajouter("Angle_min","0.7854",OT_PARAMETRES::DOUBLE,"Angle minimal admissible entre deux faces");
  params->ajouter("Nb_iteration_max","2000",OT_PARAMETRES::DOUBLE,"Nombre maximal d'iteration");
}

void PARAMETRES::ini_param_generateur_dcr(OT_PARAMETRES* params, int type_inclusion)
{
  params->ajouter("Type_generateur","1",OT_PARAMETRES::DOUBLE,"Type de generateur : RSA [0], DCR [1]");
  params->ajouter("Nom_mgcg_modele","VES",OT_PARAMETRES::STRING,"Nom du MG_CG_MODELE");
  params->ajouter("Nom_groupe_inclusion","Particule",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion (Ex: Particule)");  
  params->ajouter("Boite3D_distribution_Xmin","-0.25",OT_PARAMETRES::DOUBLE,"Xmin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Ymin","-0.25",OT_PARAMETRES::DOUBLE,"Ymin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Zmin","-0.25",OT_PARAMETRES::DOUBLE,"Zmin boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Xmax","1.25",OT_PARAMETRES::DOUBLE,"Xmax boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Ymax","1.25",OT_PARAMETRES::DOUBLE,"Ymax boite3D de distribution d'inclusion");
  params->ajouter("Boite3D_distribution_Zmax","1.25",OT_PARAMETRES::DOUBLE,"Zmax boite3D de distribution d'inclusion");  
  params->ajouter("Nb_pas_X","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon X pout la creation d'un octree");
  params->ajouter("Nb_pas_Y","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Y pout la creation d'un octree");
  params->ajouter("Nb_pas_Z","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Z pout la creation d'un octree");  
  params->ajouter("Fraction_volumique_cible","0.1",OT_PARAMETRES::DOUBLE,"Fraction volumique cible d'inclusion [0.0, 1.0]");
  params->ajouter("Eps_fraction_volumique","0.001",OT_PARAMETRES::DOUBLE,"Epsilon fraction volumique");
  params->ajouter("Avec_intersections","1",OT_PARAMETRES::DOUBLE,"Avec ou sans intersections avec le bord du VES : Sans [0] Avec [1]");  
  params->ajouter("Avec_interface_graphique","0",OT_PARAMETRES::DOUBLE,"Avec l'interface graphique : Sans [0] Avec [1]");
  params->ajouter("Facteur_reduction_vitesse","0.75",OT_PARAMETRES::DOUBLE,"Facteur de reduction de la vitesse des particules");
  params->ajouter("Eps_vitesse","0.0001",OT_PARAMETRES::DOUBLE,"Eps vitesse");
  params->ajouter("Friction","0.0",OT_PARAMETRES::DOUBLE,"Friction entre les particules");
  params->ajouter("Pas_temps_sim","0.01",OT_PARAMETRES::DOUBLE,"Pas de temps pour la simulation");
  params->ajouter("Temps_max_sim","60.0",OT_PARAMETRES::DOUBLE,"Temps maximal de la simulation");
  params->ajouter("Avec_post_RSA","0",OT_PARAMETRES::DOUBLE,"Effectuer un RSA pour tenter de remplir les vides au bord du VES : Sans [0] Avec [1]");
  params->ajouter("Nb_iteration_max","2000",OT_PARAMETRES::DOUBLE,"Nombre maximal d'iteration");
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::SPHERE)
  {    
    params->ajouter("Type_inclusion","0",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
    params->ajouter("Mu_rayon","0.09",OT_PARAMETRES::DOUBLE,"Moyenne du rayon de la sphere");
    params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon de la sphere");
    params->ajouter("Type_distribution_rayon","0",OT_PARAMETRES::DOUBLE,"Type de distribution du rayon : Fixe[0], Uniforme[1], Normale[2]");
  }
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    params->ajouter("Type_inclusion","1",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
    params->ajouter("Mu_rayon","0.05",OT_PARAMETRES::DOUBLE,"Moyenne du rayon du cylindre");
    params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon du cylindre");
    params->ajouter("Type_distribution_rayon","0",OT_PARAMETRES::DOUBLE,"Type de distribution du rayon : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_longueur","0.4",OT_PARAMETRES::DOUBLE,"Moyenne de la longueur du cylindre");
    params->ajouter("Sigma_longueur","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la longueur du cylindre");
    params->ajouter("Type_distribution_longueur","0",OT_PARAMETRES::DOUBLE,"Type de distribution de la longueur : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_theta","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle theta");
    params->ajouter("Sigma_theta","1.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle theta");
    params->ajouter("Type_distribution_theta","1",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle theta : Fixe[0], Uniforme[1], Normale[2]");
    params->ajouter("Mu_phi","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle phi");
    params->ajouter("Sigma_phi","1.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle phi");
    params->ajouter("Type_distribution_phi","1",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle phi : Fixe[0], Uniforme[1], Normale[2]");
  }
  params->ajouter("Porosite","0",OT_PARAMETRES::DOUBLE,"Considerer l'inclusion comme etant une porosite (vide) [1]");
  params->ajouter("Multicouche","0,(0.01,0.01)",OT_PARAMETRES::DOUBLE,"Generer une inclusion avec une ou plusieurs couches [1,2,3,etc] et d'epaisseurs (0.01,0.02,etc)");
  params->ajouter("Distance_min_inter_volume","0.04",OT_PARAMETRES::DOUBLE,"Distance minimale qui separe deux inclusions");
  params->ajouter("Volume_min","0.0008",OT_PARAMETRES::DOUBLE,"Volume minimal admissible de l'inclusion");
  params->ajouter("Aire_min","0.001963",OT_PARAMETRES::DOUBLE,"Aire minimal admissible d'une face de l'inclusion");
  params->ajouter("Longueur_min","0.02",OT_PARAMETRES::DOUBLE,"Longueur minimal admissible d'une arete ");
  params->ajouter("Angle_min","0.7854",OT_PARAMETRES::DOUBLE,"Angle minimal admissible entre deux faces");
}


void PARAMETRES::ini_param_materiau(OT_PARAMETRES* params,char* nom_groupe_forme,bool etude_mecanique,bool etude_thermique)
{
  params->ajouter("Nom_groupe_forme",nom_groupe_forme,OT_PARAMETRES::STRING,"Nom du groupe de forme (Ex: Particule, Matrice)");
  if(etude_mecanique) params->ajouter("Etude_mecanique","1",OT_PARAMETRES::DOUBLE,"Etude mecanique");
  else params->ajouter("Etude_mecanique","0",OT_PARAMETRES::DOUBLE,"Etude mecanique");   
  params->ajouter("Mu_E","69.0e9",OT_PARAMETRES::DOUBLE,"Moyenne du module d'elasticite");
  params->ajouter("Sigma_E","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du module d'elasticite");
  params->ajouter("Type_distribution_E","0",OT_PARAMETRES::DOUBLE,"Type de distribution du module d'elasticite : Fixe[0], Uniforme[1], Normale[2]");
  params->ajouter("Mu_nu","0.3",OT_PARAMETRES::DOUBLE,"Moyenne du coefficient de Poisson");
  params->ajouter("Sigma_nu","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du coefficient de Poisson");
  params->ajouter("Type_distribution_nu","0",OT_PARAMETRES::DOUBLE,"Type de distribution du coefficient de Poisson : Fixe[0], Uniforme[1], Normale[2]");  
  if(etude_thermique) params->ajouter("Etude_thermique","1",OT_PARAMETRES::DOUBLE,"Etude thermique");
  else params->ajouter("Etude_thermique","0",OT_PARAMETRES::DOUBLE,"Etude thermique");  
  params->ajouter("Mu_Lambda","0.14",OT_PARAMETRES::DOUBLE,"Moyenne de la conductivite thermique");
  params->ajouter("Sigma_Lambda","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la conductivite thermique");
  params->ajouter("Type_distribution_Lambda","0",OT_PARAMETRES::DOUBLE,"Type de distribution de la conductivite thermique : Fixe[0], Uniforme[1], Normale[2]");
  params->ajouter("Mu_Ro","1.0",OT_PARAMETRES::DOUBLE,"Moyenne de la masse volumique");
  params->ajouter("Sigma_Ro","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la masse volumique");
  params->ajouter("Type_distribution_Ro","0",OT_PARAMETRES::DOUBLE,"Type de distribution de la masse volumique : Fixe[0], Uniforme[1], Normale[2]");
  params->ajouter("Mu_Cp","1.0",OT_PARAMETRES::DOUBLE,"Moyenne de capacite thermique massique");
  params->ajouter("Sigma_Cp","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la capacite thermique massique");
  params->ajouter("Type_distribution_Cp","0",OT_PARAMETRES::DOUBLE,"Type de distribution de la capacite thermique massique : Fixe[0], Uniforme[1], Normale[2]");     
}

void PARAMETRES::ini_param_carte_taille(OT_PARAMETRES* params)
{
  params->ajouter("Nom_carte","carte",OT_PARAMETRES::STRING,"Nom de la carte de taille");
  params->ajouter("Ecart_nodal","0.08",OT_PARAMETRES::DOUBLE,"Ecart nodal");
  params->ajouter("Ecart_nodal_face_particule","0.04",OT_PARAMETRES::DOUBLE,"Ecart nodal des faces des particules");
  params->ajouter("Fechantillonnage","20",OT_PARAMETRES::DOUBLE,"Nombre d'échantillon par entité");
  params->ajouter("Nb_cellule_direction","50",OT_PARAMETRES::DOUBLE,"Nombre de cellules de la FEM_SOLUTION dans les directions x,y,z");
  params->ajouter("Nom_fem_solution","solcarte",OT_PARAMETRES::STRING,"Nom de la FEM_SOLUTION");
  params->ajouter("Nb_couche_min","2",OT_PARAMETRES::DOUBLE,"Nombre de couches entre chaque solide");
  params->ajouter("Nb_pas","32",OT_PARAMETRES::DOUBLE,"Nombre de pas de balayage");
  params->ajouter("Facteur_augmentation","25",OT_PARAMETRES::DOUBLE,"Facteur d'augmentation de l'ecart nodal d'une cellule a l'autre");
}

void PARAMETRES::ini_param_mailleur(OT_PARAMETRES* params)
{
  params->ajouter("Niveau","3",OT_PARAMETRES::DOUBLE,"Niveau du maillage : 1D [1], 2D [2], 3D [3]");
  params->ajouter("Niveau_opt_2d","5",OT_PARAMETRES::DOUBLE,"Limite supérieure des mailles à optimiser");
  params->ajouter("Niveau_opt_3d","2",OT_PARAMETRES::DOUBLE,"Limite supérieure des mailles à optimiser");
  params->ajouter("Priorite_metrique","0.65",OT_PARAMETRES::DOUBLE,"Valeur de la prise en compte de la métrique versus la qualité");
}

void PARAMETRES::ini_param_mailleur_fem(OT_PARAMETRES* params)
{
  params->ajouter("Degre","2",OT_PARAMETRES::DOUBLE,"Degre du maillage a creer");
  params->ajouter("Analyse","0",OT_PARAMETRES::DOUBLE,"Analyse de la qualité du maillage après le maillage : Analyse [1]");
  params->ajouter("Optimisation_num_noeud","0",OT_PARAMETRES::DOUBLE,"Optimisation de la numerotation des noeuds : Optimisation [1]");
}

void PARAMETRES::ini_param_etude(OT_PARAMETRES* params, int type_etude)
{
  params->ajouter("Fichier_param_aster","param_aster.txt",OT_PARAMETRES::STRING,"Nom du fichier de parametres aster");
  params->ajouter("Type_etude",type_etude,OT_PARAMETRES::DOUBLE,"Type d'etude : Mecanique [0], Thermique [1]");
  params->ajouter("Code_resu","11111111",OT_PARAMETRES::STRING,"Code de resolution");
  if(type_etude==MSTRUCT::TYPE_ETUDE::MECANIQUE)
  {
    params->ajouter("Type_Chargement","0",OT_PARAMETRES::DOUBLE,"Type de chargement : Spherique [0], Deviatorique [1]");
    params->ajouter("Type_CL","0",OT_PARAMETRES::DOUBLE,"Type de conditions aux limites : CLDH [0], CLCH[1]");
    params->ajouter("Valeur_CL","0.0001",OT_PARAMETRES::DOUBLE,"Valeur de la condition aux limites");
  }  
  else if(type_etude==MSTRUCT::TYPE_ETUDE::THERMIQUE)
  {
    params->ajouter("Type_CL","2",OT_PARAMETRES::DOUBLE,"Type de conditions aux limites : Gradient de temperature[2], Flux de chaleur[3]");
    params->ajouter("Direction","0",OT_PARAMETRES::DOUBLE,"Direction du gradient ou flux : X[0], Y[1], Z[2]");
    params->ajouter("Valeur_CL","1.0",OT_PARAMETRES::DOUBLE,"Valeur de la condition aux limites");
  }
}

void PARAMETRES::ini_param_post_traitement(OT_PARAMETRES* params, int type_post_traitement, int nb_couche_erosion, double epaisseur_couche_erosion)
{
  bool erosion;
  if(nb_couche_erosion>0) erosion=true;
  else erosion=false;
  if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP)
  {
    params->ajouter("Identifiant","Champ_SIGMA",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Num_solution",2,OT_PARAMETRES::DOUBLE,"Numero de la FEM_SOLUTION solution (0,1,2,...)");
    params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
    params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
    params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION)
  {
    params->ajouter("Identifiant","Orientation_Particule",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Nom_groupe_forme","Particule",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion.");
    params->ajouter("Avec_fem_maillage","1",OT_PARAMETRES::DOUBLE,"Analyser le volume des particules avec le FEM_MAILLAGE");
    params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
    params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CAO)
  {
    params->ajouter("Identifiant","Cao",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CAO,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [NULL] prend le ves au complet");
    params->ajouter("Largeur_colonne_distribution_nb_volume","1.0",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_nb_forme","1.0",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_volume","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_fraction_volumique","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_volume_forme","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_PROPRIETE_MASSIQUE)
  {
    params->ajouter("Identifiant","Prop_mass",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_PROPRIETE_MASSIQUE,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Nom_groupe_forme","Particule",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [NULL] prend le ves au complet");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_MG)
  {
    params->ajouter("Identifiant","Maillage_MG",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_MG,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
    params->ajouter("Largeur_colonne_distribution_nb_element_2D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_nb_element_3D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_qualite_2D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_qualite_3D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_taille_2D","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_taille_3D","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_volume","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_fraction_volumique","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
    params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_FEM)
  {
    params->ajouter("Identifiant","Maillage_FEM",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_FEM,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
    params->ajouter("Largeur_colonne_distribution_nb_element_2D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_nb_element_3D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_jacobien_2D_min","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_jacobien_2D_max","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_jacobien_3D_min","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_jacobien_3D_max","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_distortion_2D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_distortion_3D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_volume","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_fraction_volumique","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
    params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
    params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_TEMPS)
  {
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_TEMPS,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Largeur_colonne_temps_geometrie","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_materiau","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_carte","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_mg_maillage","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_fem_maillage","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_etude","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_temps_calcul","1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_MODULES_ELASTICITE)
  {
    params->ajouter("Identifiant","Modules_elasticite",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_MODULES_ELASTICITE,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_epsilon","Champ_EPSILON",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ epsilon");
    params->ajouter("Identifiant_sigma","Champ_SIGMA",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ sigma");
    params->ajouter("Largeur_colonne_distribution_module_Kapp","1000000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_module_Gapp","1000000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_module_Eapp","1000000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
    params->ajouter("Largeur_colonne_distribution_module_Nuapp","0.00001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE)
  {
    params->ajouter("Identifiant","Champ_EPSILON_normalise",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_champ","Champ_EPSILON",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ epsilon");
    params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART)
  {
    params->ajouter("Identifiant","Champ_EPSILON_ecart",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_champ","Champ_EPSILON",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ");
    params->ajouter("Identifiant_champ_compare","Champ_EPSILON_compare",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ a comparer");
    params->ajouter("Num_couche","0",OT_PARAMETRES::DOUBLE,"Numero de la couche a comparer");
    params->ajouter("Type_ecart","0",OT_PARAMETRES::STRING,"Ecart relatif [0], ecart absolue [1], ecart quadratique [2]");
    params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT)
  {
    params->ajouter("Identifiant","Champ_EPSILON_ecart",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_champ","Champ_EPSILON",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ");
    params->ajouter("Type_chargement","0",OT_PARAMETRES::STRING,"Spherique [0], deviatorique [1]");
    params->ajouter("Type_ecart","0",OT_PARAMETRES::STRING,"Ecart relatif [0], ecart absolue [1], ecart quadratique [2]");
    params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION_ECART)
  {
    params->ajouter("Identifiant","Orientation_Particule_ecart",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION_ECART,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_orientation","Orientation_Particule",OT_PARAMETRES::STRING,"Identifiant de l'analyse de l'orientation");
    params->ajouter("Type_ecart","0",OT_PARAMETRES::STRING,"Ecart relatif [0], ecart absolue [1], ecart quadratique [2]");
    params->ajouter("a_11","0.33333333333",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
    params->ajouter("a_22","0.33333333333",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
    params->ajouter("a_33","0.33333333333",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
    params->ajouter("a_12","0.0",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
    params->ajouter("a_23","0.0",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
    params->ajouter("a_13","0.0",OT_PARAMETRES::DOUBLE,"Valeur du tenseur d'orientation");
  }
  else if(type_post_traitement==MSTRUCT::TYPE_POST_TRAITEMENT::POST_CONDUCTIVITE_THERMIQUE)
  {
    params->ajouter("Identifiant","Conductivite_thermique",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
    params->ajouter("Type_post_traitement",MSTRUCT::TYPE_POST_TRAITEMENT::POST_CONDUCTIVITE_THERMIQUE,OT_PARAMETRES::DOUBLE,"Type de post_traitement");
    params->ajouter("Identifiant_champ_temp","Champ_TEMP",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ");
    params->ajouter("Identifiant_champ_flux","Champ_FLUX",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ");
    params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
  }
  if(erosion)
  {
    params->ajouter("Analyse_erosion","1",OT_PARAMETRES::DOUBLE,"Analyse d'erosion [1], sans erosion [0]");
    params->ajouter("Nb_couche",nb_couche_erosion+1,OT_PARAMETRES::DOUBLE,"Nombre de couche d'erosion");
    params->ajouter("Epaisseur_couche",epaisseur_couche_erosion,OT_PARAMETRES::DOUBLE,"Epaisseur d'une couche d'erosion");
  }
  if(!erosion)
  {
    params->ajouter("Analyse_erosion","0",OT_PARAMETRES::DOUBLE,"Analyse d'erosion [1], sans erosion [0]");
  }
}


int PARAMETRES::generer_fichier_parametres(bool etude_mecanique,
                                           bool etude_thermique,
                                           int type_generateur,
                                           int type_inclusion,
                                           bool fem_maill_lin,
                                           bool fem_maill_quad,
                                           int nb_couche_erosion,
                                           double epaisseur_couche_erosion,
                                           fonction_affiche *fonc_affiche)
{
  bool erosion;
  if(nb_couche_erosion>0) erosion=true;
  else erosion=false;
  int mkdir_err=0;
  std::string sys_command = "mkdir parametre";
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/geometrie");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres/geometrie ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/materiau");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/materiau ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/carte");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres/carte ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/maillage");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/maillage ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/fem_maillage");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/fem_maillage ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/etude");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/etude ! ***" << std::endl;
    return FAIL;
  }
  sys_command.clear(); sys_command.append("mkdir -p "); sys_command.append("parametre/post_traitement");
  mkdir_err = system(sys_command.c_str());
  if(mkdir_err !=0)
  {
    std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/post_traitement ! ***" << std::endl;
    return FAIL;
  }
  std::ofstream f;
  
  if(type_generateur==MSTRUCT::TYPE_GENERATEUR::RSA)
  {
    OT_PARAMETRES params_geo_1;
    ini_param_generateur_rsa(&params_geo_1,type_inclusion);
    sys_command.clear(); sys_command.append("parametre/geometrie/param_geo_particule.txt");
    params_geo_1.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  else if(type_generateur==MSTRUCT::TYPE_GENERATEUR::DCR)
  {
    OT_PARAMETRES params_geo_1;
    ini_param_generateur_dcr(&params_geo_1,type_inclusion);
    sys_command.clear(); sys_command.append("parametre/geometrie/param_geo_particule.txt");
    params_geo_1.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }  
  sys_command.clear(); sys_command.append("parametre/geometrie/lstparam_geo.txt");
  f.open(sys_command.c_str(),std::ios::out);
  f << "// Fichier de parametres GEOMETRIE" << std::endl;  
  f << "param_geo_particule.txt" << std::endl;
  f.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_mat_particule;
  ini_param_materiau(&params_mat_particule,(char*)"Particule",etude_mecanique,etude_thermique);
  if(etude_mecanique)
  {
    params_mat_particule.change_valeur("Mu_E",1.0e9);
    params_mat_particule.change_valeur("Mu_nu",0.1);
  }
  if(etude_thermique)
  {
    params_mat_particule.change_valeur("Mu_Lambda",0.07);
    params_mat_particule.change_valeur("Mu_Ro",1000);
    params_mat_particule.change_valeur("Mu_Cp",1000);
  }
  sys_command.clear(); sys_command.append("parametre/materiau/param_mat_particule.txt");
  params_mat_particule.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_mat_matrice;
  ini_param_materiau(&params_mat_matrice,(char*)"Matrice",etude_mecanique,etude_thermique);
  if(etude_mecanique)
  {
    params_mat_matrice.change_valeur("Mu_E",6.0e9);
    params_mat_matrice.change_valeur("Mu_nu",0.15);
  }
  if(etude_thermique)
  {
    params_mat_particule.change_valeur("Mu_Lambda",0.38);
    params_mat_particule.change_valeur("Mu_Ro",1000);
    params_mat_particule.change_valeur("Mu_Cp",1000);
  }  
  sys_command.clear(); sys_command.append("parametre/materiau/param_mat_matrice.txt");
  params_mat_matrice.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  sys_command.clear(); sys_command.append("parametre/materiau/lstparam_mat.txt");
  f.open(sys_command.c_str(),std::ios::out);
  f << "// Fichier de parametres MATERIAU" << std::endl;
  f << "param_mat_particule.txt" << std::endl;
  f << "param_mat_matrice.txt" << std::endl;
  f.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  
  OT_PARAMETRES params_carte;
  ini_param_carte_taille(&params_carte);
  sys_command.clear(); sys_command.append("parametre/carte/param_carte.txt");
  params_carte.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_maill;
  ini_param_mailleur(&params_maill);
  sys_command.clear(); sys_command.append("parametre/maillage/param_maill.txt");
  params_maill.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
   
  if(fem_maill_lin)
  {
    OT_PARAMETRES params_fem_maillage_lin;
    ini_param_mailleur_fem(&params_fem_maillage_lin);
    params_fem_maillage_lin.change_valeur("Degre",1);
    sys_command.clear(); sys_command.append("parametre/fem_maillage/param_fem_maill_lin.txt");
    params_fem_maillage_lin.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  
  if(fem_maill_quad)
  {
    OT_PARAMETRES params_fem_maillage_quad;
    ini_param_mailleur_fem(&params_fem_maillage_quad);
    params_fem_maillage_quad.change_valeur("Degre",2);
    sys_command.clear(); sys_command.append("parametre/fem_maillage/param_fem_maill_quad.txt");
    params_fem_maillage_quad.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  
  if(etude_mecanique)
  {
    OT_PARAMETRES params_etude_CLDH_sph;
    ini_param_etude(&params_etude_CLDH_sph,MSTRUCT::TYPE_ETUDE::MECANIQUE);
    params_etude_CLDH_sph.change_valeur("Type_Chargement",MSTRUCT::TYPE_CHARGEMENT::SPHERIQUE);
    params_etude_CLDH_sph.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::DEFORMATION_HOMOGENE);
    params_etude_CLDH_sph.change_valeur("Valeur_CL",0.0001);
    sys_command.clear(); sys_command.append("parametre/etude/param_etude_CLDH_sph.txt");
    params_etude_CLDH_sph.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLCH_sph;
    ini_param_etude(&params_etude_CLCH_sph,MSTRUCT::TYPE_ETUDE::MECANIQUE);
    params_etude_CLCH_sph.change_valeur("Type_Chargement",MSTRUCT::TYPE_CHARGEMENT::SPHERIQUE);
    params_etude_CLCH_sph.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::CONTRAINTE_HOMOGENE);
    params_etude_CLCH_sph.change_valeur("Valeur_CL",1000);
    sys_command.clear(); sys_command.append("parametre/etude/param_etude_CLCH_sph.txt");
    params_etude_CLCH_sph.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLDH_dev;
    ini_param_etude(&params_etude_CLDH_dev,MSTRUCT::TYPE_ETUDE::MECANIQUE);
    params_etude_CLDH_dev.change_valeur("Type_Chargement",MSTRUCT::TYPE_CHARGEMENT::DEVIATORIQUE);
    params_etude_CLDH_dev.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::DEFORMATION_HOMOGENE);
    params_etude_CLDH_dev.change_valeur("Valeur_CL",0.0001);
    sys_command.clear(); sys_command.append("parametre/etude/param_etude_CLDH_dev.txt");
    params_etude_CLDH_dev.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLCH_dev;
    ini_param_etude(&params_etude_CLCH_dev,MSTRUCT::TYPE_ETUDE::MECANIQUE);
    params_etude_CLCH_dev.change_valeur("Type_Chargement",MSTRUCT::TYPE_CHARGEMENT::DEVIATORIQUE);
    params_etude_CLCH_dev.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::CONTRAINTE_HOMOGENE);
    params_etude_CLCH_dev.change_valeur("Valeur_CL",1000);
    sys_command.clear(); sys_command.append("parametre/etude/param_etude_CLCH_dev.txt");
    params_etude_CLCH_dev.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  if(etude_thermique)
  {
    OT_PARAMETRES params_etude_CLG_x;
    ini_param_etude(&params_etude_CLG_x,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLG_x.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::GRADIENT_TEMPERATURE);
    params_etude_CLG_x.change_valeur("Valeur_CL",1.0);
    params_etude_CLG_x.change_valeur("Direction",0);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLG_x.txt");
    params_etude_CLG_x.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLG_y;
    ini_param_etude(&params_etude_CLG_y,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLG_y.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::GRADIENT_TEMPERATURE);
    params_etude_CLG_y.change_valeur("Valeur_CL",1.0);
    params_etude_CLG_y.change_valeur("Direction",1);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLG_y.txt");
    params_etude_CLG_y.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLG_z;
    ini_param_etude(&params_etude_CLG_z,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLG_z.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::GRADIENT_TEMPERATURE);
    params_etude_CLG_z.change_valeur("Valeur_CL",1.0);
    params_etude_CLG_z.change_valeur("Direction",2);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLG_z.txt");
    params_etude_CLG_z.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLF_x;
    ini_param_etude(&params_etude_CLF_x,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLF_x.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::FLUX_CHALEUR);
    params_etude_CLF_x.change_valeur("Valeur_CL",1.0);
    params_etude_CLF_x.change_valeur("Direction",0);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLF_x.txt");
    params_etude_CLF_x.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLF_y;
    ini_param_etude(&params_etude_CLF_y,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLF_y.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::FLUX_CHALEUR);
    params_etude_CLF_y.change_valeur("Valeur_CL",1.0);
    params_etude_CLF_y.change_valeur("Direction",1);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLF_y.txt");
    params_etude_CLF_y.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_etude_CLF_z;
    ini_param_etude(&params_etude_CLF_z,MSTRUCT::TYPE_ETUDE::THERMIQUE);
    params_etude_CLF_z.change_valeur("Type_CL",MSTRUCT::TYPE_CONDITIONS_LIMITES::FLUX_CHALEUR);
    params_etude_CLF_z.change_valeur("Valeur_CL",1.0);
    params_etude_CLF_z.change_valeur("Direction",2);
    sys_command.clear(); sys_command.append("parametre/etude/params_etude_CLF_z.txt");
    params_etude_CLF_z.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  
  MG_EXPORT exp;
  OT_PARAMETRES *params_aster=exp.get_param_aster();
  params_aster->change_valeur("Memoire",1024);
  params_aster->change_valeur("Temps_max",250000);
  params_aster->change_valeur("Nb_CPU",8);
  params_aster->change_valeur("Noeud_ele",2);
  params_aster->change_valeur("Base_num_mail",36);
  params_aster->change_valeur("Pas",0);
  sys_command.clear(); sys_command.append("parametre/etude/param_aster.txt");
  params_aster->enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_post_cao_all;
  ini_param_post_traitement(&params_post_cao_all,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CAO);
  params_post_cao_all.change_valeur("Identifiant","Cao_ALL");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_cao_all.txt");
  params_post_cao_all.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  OT_PARAMETRES params_post_cao_particule;
  ini_param_post_traitement(&params_post_cao_particule,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CAO);
  params_post_cao_particule.change_valeur("Identifiant","Cao_Particule");
  params_post_cao_particule.change_valeur("Nom_groupe_forme","Particule");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_cao_particule.txt");
  params_post_cao_particule.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  OT_PARAMETRES params_post_cao_matrice;
  ini_param_post_traitement(&params_post_cao_matrice,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CAO);
  params_post_cao_matrice.change_valeur("Identifiant","Cao_Matrice");
  params_post_cao_matrice.change_valeur("Nom_groupe_forme","Matrice");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_cao_matrice.txt");
  params_post_cao_matrice.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  OT_PARAMETRES params_post_prop_mass_particule;
  ini_param_post_traitement(&params_post_prop_mass_particule,MSTRUCT::TYPE_POST_TRAITEMENT::POST_PROPRIETE_MASSIQUE);
  params_post_prop_mass_particule.change_valeur("Identifiant","Prop_mass_Particule");
  params_post_prop_mass_particule.change_valeur("Nom_groupe_forme","Particule");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_prop_mass_particule.txt");
  params_post_prop_mass_particule.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    OT_PARAMETRES params_post_orientation_cao;
    ini_param_post_traitement(&params_post_orientation_cao,MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION);
    params_post_orientation_cao.change_valeur("Identifiant","Orientation_Particule");
    params_post_orientation_cao.change_valeur("Avec_fem_maillage","0");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_orientation_cao.txt");
    params_post_orientation_cao.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_orientation_ecart;
    ini_param_post_traitement(&params_post_orientation_ecart,MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION_ECART);
    params_post_orientation_ecart.change_valeur("Identifiant","Orientation_Particule_ecart");
    params_post_orientation_ecart.change_valeur("Identifiant_orientation","Orientation_Particule");
    params_post_orientation_ecart.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_ABSOLUE);
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_orientation_ecart.txt");
    params_post_orientation_ecart.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }  
  
  sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_geometrie.txt");
  f.open(sys_command.c_str(),std::ios::out);
  f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
  f << "params_post_cao_all.txt" << std::endl;
  f << "params_post_cao_particule.txt" << std::endl;
  f << "params_post_cao_matrice.txt" << std::endl;
  f << "params_post_prop_mass_particule.txt" << std::endl;
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    f << "params_post_orientation_cao.txt" << std::endl;
    f << "params_post_orientation_ecart.txt" << std::endl;
  }
  
  f.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_post_mg_maillage_all;
  ini_param_post_traitement(&params_post_mg_maillage_all,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_MG);
  params_post_mg_maillage_all.change_valeur("Identifiant","MG_maillage_ALL");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_mg_maillage_all.txt");
  params_post_mg_maillage_all.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_post_mg_maillage_particule;
  ini_param_post_traitement(&params_post_mg_maillage_particule,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_MG);
  params_post_mg_maillage_particule.change_valeur("Identifiant","MG_maillage_Particule");
  params_post_mg_maillage_particule.change_valeur("Nom_groupe_forme","Particule");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_mg_maillage_particule.txt");
  params_post_mg_maillage_particule.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
  OT_PARAMETRES params_post_mg_maillage_matrice;
  ini_param_post_traitement(&params_post_mg_maillage_matrice,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_MG);
  params_post_mg_maillage_matrice.change_valeur("Identifiant","MG_maillage_Matrice");
  params_post_mg_maillage_matrice.change_valeur("Nom_groupe_forme","Matrice");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_mg_maillage_matrice.txt");
  params_post_mg_maillage_matrice.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  
  sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_mg_maillage.txt");
  f.open(sys_command.c_str(),std::ios::out); 
  f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
  f << "params_post_mg_maillage_all.txt" << std::endl;
  f << "params_post_mg_maillage_particule.txt" << std::endl;
  f << "params_post_mg_maillage_matrice.txt" << std::endl;  
  f.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_post_fem_maillage_all;
  ini_param_post_traitement(&params_post_fem_maillage_all,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_FEM,nb_couche_erosion,epaisseur_couche_erosion);
  params_post_fem_maillage_all.change_valeur("Identifiant","FEM_maillage_ALL");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_fem_maillage_all.txt");
  params_post_fem_maillage_all.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  OT_PARAMETRES params_post_fem_maillage_particule;
  ini_param_post_traitement(&params_post_fem_maillage_particule,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_FEM,nb_couche_erosion,epaisseur_couche_erosion);
  params_post_fem_maillage_particule.change_valeur("Identifiant","FEM_maillage_Particule");
  params_post_fem_maillage_particule.change_valeur("Nom_groupe_forme","Particule");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_fem_maillage_particule.txt");
  params_post_fem_maillage_particule.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
  OT_PARAMETRES params_post_fem_maillage_matrice;
  ini_param_post_traitement(&params_post_fem_maillage_matrice,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MAILLAGE_FEM,nb_couche_erosion,epaisseur_couche_erosion);
  params_post_fem_maillage_matrice.change_valeur("Identifiant","FEM_maillage_Matrice");
  params_post_fem_maillage_matrice.change_valeur("Nom_groupe_forme","Matrice");
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_fem_maillage_matrice.txt");
  params_post_fem_maillage_matrice.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    OT_PARAMETRES params_post_orientation_ponderee;
    ini_param_post_traitement(&params_post_orientation_ponderee,MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_orientation_ponderee.change_valeur("Identifiant","Orientation_Particule_ponderee");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_orientation_ponderee.txt");
    params_post_orientation_ponderee.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_orientation_ponderee_ecart;
    ini_param_post_traitement(&params_post_orientation_ponderee_ecart,MSTRUCT::TYPE_POST_TRAITEMENT::POST_ORIENTATION_ECART,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_orientation_ponderee_ecart.change_valeur("Identifiant","Orientation_Particule_ponderee_ecart");
    params_post_orientation_ponderee_ecart.change_valeur("Identifiant_orientation","Orientation_Particule_ponderee");
    params_post_orientation_ponderee_ecart.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_ABSOLUE);
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_orientation_ponderee_ecart.txt");
    params_post_orientation_ponderee_ecart.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_fem_maillage.txt");
  f.open(sys_command.c_str(),std::ios::out);  
  f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
  f << "params_post_fem_maillage_all.txt" << std::endl;
  f << "params_post_fem_maillage_particule.txt" << std::endl;  
  f << "params_post_fem_maillage_matrice.txt" << std::endl;
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    f << "params_post_orientation_ponderee.txt" << std::endl; 
    f << "params_post_orientation_ponderee_ecart.txt" << std::endl; 
  }
  f.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str()); 

  if(etude_mecanique)
  {
    OT_PARAMETRES params_post_champ_sigma;
    ini_param_post_traitement(&params_post_champ_sigma,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_sigma.change_valeur("Identifiant","Champ_SIGMA");
    params_post_champ_sigma.change_valeur("Num_solution",2);
    params_post_champ_sigma.change_valeur("Largeur_colonne_distribution","0.1");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_sigma.txt");
    params_post_champ_sigma.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_sigma_normalise;
    ini_param_post_traitement(&params_post_champ_sigma_normalise,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_sigma_normalise.change_valeur("Identifiant","Champ_SIGMA_normalise");
    params_post_champ_sigma_normalise.change_valeur("Identifiant_champ","Champ_SIGMA");
    params_post_champ_sigma_normalise.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_sigma_normalise.txt");
    params_post_champ_sigma_normalise.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());   
    
    OT_PARAMETRES params_post_champ_epsilon;
    ini_param_post_traitement(&params_post_champ_epsilon,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_epsilon.change_valeur("Identifiant","Champ_EPSILON");
    params_post_champ_epsilon.change_valeur("Num_solution",1);
    params_post_champ_epsilon.change_valeur("Largeur_colonne_distribution","1.0e-11");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_epsilon.txt");
    params_post_champ_epsilon.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_epsilon_normalise;
    ini_param_post_traitement(&params_post_champ_epsilon_normalise,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_epsilon_normalise.change_valeur("Identifiant","Champ_EPSILON_normalise");
    params_post_champ_epsilon_normalise.change_valeur("Identifiant_champ","Champ_EPSILON");
    params_post_champ_epsilon_normalise.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_epsilon_normalise.txt");
    params_post_champ_epsilon_normalise.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_energie;
    ini_param_post_traitement(&params_post_champ_energie,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_energie.change_valeur("Identifiant","Champ_ENERGIE");
    params_post_champ_energie.change_valeur("Num_solution",5);
    params_post_champ_energie.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_energie.txt");
    params_post_champ_energie.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_energie_normalise;
    ini_param_post_traitement(&params_post_champ_energie_normalise,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_energie_normalise.change_valeur("Identifiant","Champ_ENERGIE_normalise");
    params_post_champ_energie_normalise.change_valeur("Identifiant_champ","Champ_ENERGIE");
    params_post_champ_energie_normalise.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_energie_normalise.txt");
    params_post_champ_energie_normalise.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());    
    
    OT_PARAMETRES params_post_champ_energie_ecart;
    ini_param_post_traitement(&params_post_champ_energie_ecart,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_energie_ecart.change_valeur("Identifiant","Champ_ENERGIE_ecart");
    params_post_champ_energie_ecart.change_valeur("Identifiant_champ","Champ_ENERGIE");
    params_post_champ_energie_ecart.change_valeur("Identifiant_champ_compare","Champ_ENERGIE");
    params_post_champ_energie_ecart.change_valeur("Num_couche","0");
    params_post_champ_energie_ecart.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_RELATIF_ABSOLUE);
    params_post_champ_energie_ecart.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_energie_ecart.txt");
    params_post_champ_energie_ecart.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());    
    
    sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_champ_meca.txt");
    f.open(sys_command.c_str(),std::ios::out); 
    f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
    f << "params_post_champ_sigma.txt" << std::endl;
    f << "params_post_champ_sigma_normalise.txt" << std::endl;
    f << "params_post_champ_epsilon.txt" << std::endl;
    f << "params_post_champ_epsilon_normalise.txt" << std::endl;
    f << "params_post_champ_energie.txt" << std::endl;
    f << "params_post_champ_energie_normalise.txt" << std::endl;
    f << "params_post_champ_energie_ecart.txt" << std::endl;
    f.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_sigma_ecart_chargement_sph;
    ini_param_post_traitement(&params_post_champ_sigma_ecart_chargement_sph,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_sigma_ecart_chargement_sph.change_valeur("Identifiant","Champ_SIGMA_ecart_chargement_sph");
    params_post_champ_sigma_ecart_chargement_sph.change_valeur("Identifiant_champ","Champ_SIGMA");
    params_post_champ_sigma_ecart_chargement_sph.change_valeur("Type_chargement","0");
    params_post_champ_sigma_ecart_chargement_sph.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_RELATIF_ABSOLUE);
    params_post_champ_sigma_ecart_chargement_sph.change_valeur("Largeur_colonne_distribution","0.1");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_sigma_ecart_chargement_sph.txt");
    params_post_champ_sigma_ecart_chargement_sph.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_epsilon_ecart_chargement_sph;
    ini_param_post_traitement(&params_post_champ_epsilon_ecart_chargement_sph,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_epsilon_ecart_chargement_sph.change_valeur("Identifiant","Champ_EPSILON_ecart_chargement_sph");
    params_post_champ_epsilon_ecart_chargement_sph.change_valeur("Identifiant_champ","Champ_EPSILON");
    params_post_champ_epsilon_ecart_chargement_sph.change_valeur("Type_chargement","0");
    params_post_champ_epsilon_ecart_chargement_sph.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_RELATIF_ABSOLUE);
    params_post_champ_epsilon_ecart_chargement_sph.change_valeur("Largeur_colonne_distribution","0.1");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_epsilon_ecart_chargement_sph.txt");
    params_post_champ_epsilon_ecart_chargement_sph.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt");
    f.open(sys_command.c_str(),std::ios::out); 
    f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
    f << "params_post_champ_sigma_ecart_chargement_sph.txt" << std::endl;
    f << "params_post_champ_epsilon_ecart_chargement_sph.txt" << std::endl;
    f.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_sigma_ecart_chargement_dev;
    ini_param_post_traitement(&params_post_champ_sigma_ecart_chargement_dev,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_sigma_ecart_chargement_dev.change_valeur("Identifiant","Champ_SIGMA_ecart_chargement_dev");
    params_post_champ_sigma_ecart_chargement_dev.change_valeur("Identifiant_champ","Champ_SIGMA");
    params_post_champ_sigma_ecart_chargement_dev.change_valeur("Type_chargement","1");
    params_post_champ_sigma_ecart_chargement_dev.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_RELATIF_ABSOLUE);
    params_post_champ_sigma_ecart_chargement_dev.change_valeur("Largeur_colonne_distribution","0.1");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_sigma_ecart_chargement_dev.txt");
    params_post_champ_sigma_ecart_chargement_dev.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());  
    
    OT_PARAMETRES params_post_champ_epsilon_ecart_chargement_dev;
    ini_param_post_traitement(&params_post_champ_epsilon_ecart_chargement_dev,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_ECART_CHARGEMENT,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_epsilon_ecart_chargement_dev.change_valeur("Identifiant","Champ_EPSILON_ecart_chargement_dev");
    params_post_champ_epsilon_ecart_chargement_dev.change_valeur("Identifiant_champ","Champ_EPSILON");
    params_post_champ_epsilon_ecart_chargement_dev.change_valeur("Type_chargement","1");
    params_post_champ_epsilon_ecart_chargement_dev.change_valeur("Type_ecart",MSTRUCT::TYPE_COMPARAISON::ECART_RELATIF_ABSOLUE);
    params_post_champ_epsilon_ecart_chargement_dev.change_valeur("Largeur_colonne_distribution","0.1");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_epsilon_ecart_chargement_dev.txt");
    params_post_champ_epsilon_ecart_chargement_dev.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt");
    f.open(sys_command.c_str(),std::ios::out); 
    f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
    f << "params_post_champ_sigma_ecart_chargement_dev.txt" << std::endl;
    f << "params_post_champ_epsilon_ecart_chargement_dev.txt" << std::endl;
    f.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_modules_elasticite;
    ini_param_post_traitement(&params_post_modules_elasticite,MSTRUCT::TYPE_POST_TRAITEMENT::POST_MODULES_ELASTICITE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_modules_elasticite.change_valeur("Identifiant_epsilon","Champ_EPSILON");
    params_post_modules_elasticite.change_valeur("Identifiant_sigma","Champ_SIGMA");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_modules_elasticite.txt");
    params_post_modules_elasticite.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());     
  }
  if(etude_thermique)
  {
    OT_PARAMETRES params_post_champ_temp;
    ini_param_post_traitement(&params_post_champ_temp,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_temp.change_valeur("Identifiant","Champ_TEMP");
    params_post_champ_temp.change_valeur("Num_solution",0);
    params_post_champ_temp.change_valeur("Largeur_colonne_distribution","0.01");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_temp.txt");
    params_post_champ_temp.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_temp_normalise;
    ini_param_post_traitement(&params_post_champ_temp_normalise,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_temp_normalise.change_valeur("Identifiant","Champ_TEMP_normalise");
    params_post_champ_temp_normalise.change_valeur("Identifiant_champ","Champ_TEMP");
    params_post_champ_temp_normalise.change_valeur("Num_solution",0);
    params_post_champ_temp_normalise.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_temp_normalise.txt");
    params_post_champ_temp_normalise.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_flux;
    ini_param_post_traitement(&params_post_champ_flux,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_flux.change_valeur("Identifiant","Champ_FLUX");
    params_post_champ_flux.change_valeur("Num_solution",1);
    params_post_champ_flux.change_valeur("Largeur_colonne_distribution","0.01");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_flux.txt");
    params_post_champ_flux.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_champ_flux_normalise;
    ini_param_post_traitement(&params_post_champ_flux_normalise,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CHAMP_NORMALISE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_flux_normalise.change_valeur("Identifiant","Champ_FLUX_normalise");
    params_post_champ_flux_normalise.change_valeur("Identifiant_champ","Champ_FLUX");
    params_post_champ_flux_normalise.change_valeur("Num_solution",0);
    params_post_champ_flux_normalise.change_valeur("Largeur_colonne_distribution","1e-06");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_champ_flux_normalise.txt");
    params_post_champ_flux_normalise.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
    
    OT_PARAMETRES params_post_conductivite_thermique;
    ini_param_post_traitement(&params_post_conductivite_thermique,MSTRUCT::TYPE_POST_TRAITEMENT::POST_CONDUCTIVITE_THERMIQUE,nb_couche_erosion,epaisseur_couche_erosion);
    params_post_champ_flux_normalise.change_valeur("Identifiant","Conductivite_thermique");
    params_post_conductivite_thermique.change_valeur("Identifiant_champ_temp","Champ_TEMP");
    params_post_conductivite_thermique.change_valeur("Identifiant_champ_flux","Champ_FLUX");
    sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_conductivite_thermique.txt");
    params_post_conductivite_thermique.enregistrer((char*)sys_command.c_str());
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str()); 
    
    sys_command.clear(); sys_command.append("parametre/post_traitement/lstparam_post_champ_ther.txt");
    f.open(sys_command.c_str(),std::ios::out); 
    f << "// Fichier de parametres POST_TRAITEMENT" << std::endl;
    f << "params_post_champ_temp.txt" << std::endl;
    f << "params_post_champ_temp_normalise.txt" << std::endl;
    f << "params_post_champ_flux.txt" << std::endl;
    f << "params_post_champ_flux_normalise.txt" << std::endl;
    f << "params_post_conductivite_thermique.txt" << std::endl;
    f.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  }
  OT_PARAMETRES params_post_temps;
  ini_param_post_traitement(&params_post_temps,MSTRUCT::TYPE_POST_TRAITEMENT::POST_TEMPS,false);
  sys_command.clear(); sys_command.append("parametre/post_traitement/params_post_temps.txt");
  params_post_temps.enregistrer((char*)sys_command.c_str());
  if(fonc_affiche!=NULL) fonc_affiche((char*)sys_command.c_str());
  
  return OK;
}

int PARAMETRES::generer_script_ves(bool etude_mecanique,
                                   bool etude_thermique,
                                   int type_inclusion,
                                   bool fem_maill_lin,
                                   bool fem_maill_quad,
                                   fonction_affiche *fonc_affiche)
{
 
    std::ofstream ofstrm((char*)"generer_nb_ves",std::ios::out|std::ios::trunc);
    ofstrm.precision(16);
    ofstrm.setf(std::ios::showpoint);
    ofstrm << "#!/bin/bash" << std::endl;
    ofstrm << "INDEX_DEPART=$1" << std::endl;
    ofstrm << "INDEX_FIN=$2" << std::endl;
    ofstrm << "CHEMINCOMPILE=$3" << std::endl;
    ofstrm << "if [ -z ${INDEX_DEPART// } ]" << std::endl;
    ofstrm << "then" << std::endl;
    ofstrm << "  echo \"Erreur : Numero du ves de depart manquant\"" << std::endl;
    ofstrm << "  exit" << std::endl;
    ofstrm << "fi" << std::endl;
    ofstrm << "if [ -z ${INDEX_FIN// } ]" << std::endl;
    ofstrm << "then" << std::endl;
    ofstrm << "  echo \"Erreur : Numero du ves de fin manquant\"" << std::endl;
    ofstrm << "  exit" << std::endl;
    ofstrm << "fi" << std::endl;
    ofstrm << "if [ -z ${CHEMINCOMPILE// } ]" << std::endl;
    ofstrm << "then" << std::endl;
    ofstrm << "  if [ -z $PATHMAGICEXEAPP ]" << std::endl;
    ofstrm << "  then" << std::endl;
    ofstrm << "    echo \"Erreur : Chemin vers exe/app manquant\"" << std::endl;
    ofstrm << "    exit" << std::endl;
    ofstrm << "  else" << std::endl;
    ofstrm << "    CHEMINCOMPILE=$PATHMAGICEXEAPP" << std::endl;
    ofstrm << "  fi" << std::endl;    
    ofstrm << "fi" << std::endl;
    ofstrm << "function executer" << std::endl;
    ofstrm << "{" << std::endl;
    ofstrm << "OPERATION_REUSSIE=1;" << std::endl;
    ofstrm << "$*" << std::endl;
    ofstrm << "VAR=$?" << std::endl;
    ofstrm << "if (( $VAR )) ; then " << std::endl;
    ofstrm << "    OPERATION_REUSSIE=0; " << std::endl;
    ofstrm << "    STAT=$VAR;" << std::endl;
    ofstrm << "fi" << std::endl;
    ofstrm << "}" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "echo \"*********************************************************\"" << std::endl;
    ofstrm << "echo \"Generation des microstructures $INDEX_DEPART - $INDEX_FIN\"" << std::endl;
    ofstrm << "echo \"*********************************************************\"" << std::endl;
    ofstrm << "STAT=$?" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "if [ $INDEX_DEPART -eq 1 ]" << std::endl;
    ofstrm << "then" << std::endl;
    ofstrm << "mkdir etude" << std::endl;
    ofstrm << "fi" << std::endl;
    ofstrm << "cd etude" << std::endl;
    ofstrm << "for (( k=$INDEX_DEPART; k<=$INDEX_FIN; k++ ))" << std::endl;
    ofstrm << "do" << std::endl;
    ofstrm << "  echo \"Microstructure #$k\"" << std::endl;
    ofstrm << "  mkdir e$k" << std::endl;
    ofstrm << "  cd e$k" << std::endl;
    ofstrm << "  OPERATION_REUSSIE=0" << std::endl;
    ofstrm << "  until [ $OPERATION_REUSSIE -eq 1 ]; do" << std::endl;
    ofstrm << "  executer ../../ves $CHEMINCOMPILE" << std::endl;
    ofstrm << "  done" << std::endl;
    ofstrm << "  DOSSIER_ACTUEL=$(pwd)" << std::endl;       
    ofstrm << "  cd .." << std::endl;
    ofstrm << "done" << std::endl;
    ofstrm << "cd .." << std::endl;
    ofstrm << "echo \"**********************\"" << std::endl;
    ofstrm << "echo \"        FIN           \"" << std::endl;
    ofstrm << "echo \"**********************\"" << std::endl;
    ofstrm << "exit $STAT" << std::endl;
    ofstrm.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)"generer_nb_ves");
    system((char*)"chmod 744 generer_nb_ves");
    
    ofstrm.open((char*)"ves",std::ios::out|std::ios::trunc);
    ofstrm.precision(16);
    ofstrm.setf(std::ios::showpoint);    
    ofstrm << "#!/bin/bash" << std::endl;
    ofstrm << "CHEMINCOMPILE=$1" << std::endl;
    ofstrm << "function executer" << std::endl;
    ofstrm << "{" << std::endl;
    ofstrm << "$*" << std::endl;
    ofstrm << "VAR=$?" << std::endl;
    ofstrm << "if (( $VAR )) ; then " << std::endl;
    ofstrm << "    STATUSTXT=\"ERREUR dans les tests\"; " << std::endl;
    ofstrm << "    STAT=$VAR;" << std::endl;
    ofstrm << "    exit $STAT" << std::endl;
    ofstrm << "fi" << std::endl;
    ofstrm << "}" << std::endl;
    ofstrm << "STATUSTXT=\"Test sans ERREUR\"" << std::endl;
    ofstrm << "echo \"*******************************\"" << std::endl;
    ofstrm << "echo \"Generation d'une microstructure\"" << std::endl;
    ofstrm << "echo \"*******************************\"" << std::endl;
    ofstrm << "STAT=$?" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "#Geometrie" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_geo -ves_out ves_geo.ves -lstparam ../../parametre/geometrie/lstparam_geo.txt" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/mgoperation/mgoperation.exe -brep_out -in ves_geo.magic -out ves_brep.brep" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_geo.ves -ves_out ves_geo_post.ves -lstparam ../../parametre/post_traitement/lstparam_post_geometrie.txt" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "#Materiau" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_mat -ves_in ves_geo.ves -ves_out ves_mat.ves -lstparam ../../parametre/materiau/lstparam_mat.txt" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "#Carte" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_carte -ves_in ves_mat.ves -ves_out ves_carte.ves -param ../../parametre/carte/param_carte.txt" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "#MG Maillage" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_mg_maill -ves_in ves_carte.ves -ves_out ves_mg_maill.ves -param ../../parametre/maillage/param_maill.txt -carte carte.ctt" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_mg_maill.ves -ves_out ves_mg_maill_post.ves -lstparam ../../parametre/post_traitement/lstparam_post_mg_maillage.txt" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/mgoperation/mgoperation.exe -gmsh -in ves_mg_maill.magic" << std::endl;
    ofstrm << "" << std::endl;
    ofstrm << "#Fem Maillage" << std::endl;
    if(fem_maill_lin)
    {
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_fem_maill -ves_in ves_mg_maill.ves -ves_out ves_fem_maill_lin.ves -param ../../parametre/fem_maillage/param_fem_maill_lin.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_fem_maill_lin.ves -ves_out ves_fem_maill_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_fem_maillage.txt" << std::endl;
    }
    if(fem_maill_quad)
    {
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_fem_maill -ves_in ves_mg_maill.ves -ves_out ves_fem_maill_quad.ves -param ../../parametre/fem_maillage/param_fem_maill_quad.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_fem_maill_quad.ves -ves_out ves_fem_maill_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_fem_maillage.txt" << std::endl;
    }    
    ofstrm << "" << std::endl;
    ofstrm << "#Etude" << std::endl;
    if(fem_maill_lin)
    {
      if(etude_mecanique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLDH_sph_lin.ves -param ../../parametre/etude/param_etude_CLDH_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLCH_sph_lin.ves -param ../../parametre/etude/param_etude_CLCH_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLDH_dev_lin.ves -param ../../parametre/etude/param_etude_CLDH_dev.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLCH_dev_lin.ves -param ../../parametre/etude/param_etude_CLCH_dev.txt" << std::endl;
      }
      if(etude_thermique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLG_x_lin.ves -param ../../parametre/etude/params_etude_CLG_x.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLG_y_lin.ves -param ../../parametre/etude/params_etude_CLG_y.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLG_z_lin.ves -param ../../parametre/etude/params_etude_CLG_z.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLF_x_lin.ves -param ../../parametre/etude/params_etude_CLF_x.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLF_y_lin.ves -param ../../parametre/etude/params_etude_CLF_y.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_lin.ves -ves_out ves_etude_CLF_z_lin.ves -param ../../parametre/etude/params_etude_CLF_z.txt" << std::endl;
      }
    }
    if(fem_maill_quad)
    {
      if(etude_mecanique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLDH_sph_quad.ves -param ../../parametre/etude/param_etude_CLDH_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLCH_sph_quad.ves -param ../../parametre/etude/param_etude_CLCH_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLDH_dev_quad.ves -param ../../parametre/etude/param_etude_CLDH_dev.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLCH_dev_quad.ves -param ../../parametre/etude/param_etude_CLCH_dev.txt" << std::endl;
      }
      if(etude_thermique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLG_x_quad.ves -param ../../parametre/etude/params_etude_CLG_x.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLG_y_quad.ves -param ../../parametre/etude/params_etude_CLG_y.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLG_z_quad.ves -param ../../parametre/etude/params_etude_CLG_z.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLF_x_quad.ves -param ../../parametre/etude/params_etude_CLF_x.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLF_y_quad.ves -param ../../parametre/etude/params_etude_CLF_y.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill_quad.ves -ves_out ves_etude_CLF_z_quad.ves -param ../../parametre/etude/params_etude_CLF_z.txt" << std::endl;
        
      }
    }    
    ofstrm << "" << std::endl;
    ofstrm << "#Calcul" << std::endl;
    if(fem_maill_lin)
    {
      if(etude_mecanique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLDH_sph_lin.ves -ves_out ves_calcul_CLDH_sph_lin.ves -param ../../parametre/etude/param_etude_CLDH_sph.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_sph_lin.ves -ves_out ves_calcul_CLDH_sph_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_sph_post_lin.ves -ves_out ves_calcul_CLDH_sph_ecart_chargement_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLCH_sph_lin.ves -ves_out ves_calcul_CLCH_sph_lin.ves -param ../../parametre/etude/param_etude_CLCH_sph.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_sph_lin.ves -ves_out ves_calcul_CLCH_sph_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_sph_post_lin.ves -ves_out ves_calcul_CLCH_sph_ecart_chargement_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLDH_dev_lin.ves -ves_out ves_calcul_CLDH_dev_lin.ves -param ../../parametre/etude/param_etude_CLDH_dev.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_dev_lin.ves -ves_out ves_calcul_CLDH_dev_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_dev_post_lin.ves -ves_out ves_calcul_CLDH_dev_ecart_chargement_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLCH_dev_lin.ves -ves_out ves_calcul_CLCH_dev_lin.ves -param ../../parametre/etude/param_etude_CLCH_dev.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_dev_lin.ves -ves_out ves_calcul_CLCH_dev_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_dev_post_lin.ves -ves_out ves_calcul_CLCH_dev_ecart_chargement_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt" << std::endl;     
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_sph ves_calcul_CLDH_sph_post_lin.ves -ves_dev ves_calcul_CLDH_dev_post_lin.ves -ves_out ves_CLDH_modules_elasticite_post_lin.ves -param ../../parametre/post_traitement/params_post_modules_elasticite.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_sph ves_calcul_CLCH_sph_post_lin.ves -ves_dev ves_calcul_CLCH_dev_post_lin.ves -ves_out ves_CLCH_modules_elasticite_post_lin.ves -param ../../parametre/post_traitement/params_post_modules_elasticite.txt" << std::endl;
      }
      if(etude_thermique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_x_lin.ves -ves_out ves_calcul_CLG_x_lin.ves -param ../../parametre/etude/params_etude_CLG_x.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_x_lin.ves -ves_out ves_calcul_CLG_x_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_y_lin.ves -ves_out ves_calcul_CLG_y_lin.ves -param ../../parametre/etude/params_etude_CLG_y.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_y_lin.ves -ves_out ves_calcul_CLG_y_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_z_lin.ves -ves_out ves_calcul_CLG_z_lin.ves -param ../../parametre/etude/params_etude_CLG_z.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_z_lin.ves -ves_out ves_calcul_CLG_z_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_x_lin.ves -ves_out ves_calcul_CLF_x_lin.ves -param ../../parametre/etude/params_etude_CLF_x.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_x_lin.ves -ves_out ves_calcul_CLF_x_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_y_lin.ves -ves_out ves_calcul_CLF_y_lin.ves -param ../../parametre/etude/params_etude_CLF_y.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_y_lin.ves -ves_out ves_calcul_CLF_y_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_z_lin.ves -ves_out ves_calcul_CLF_z_lin.ves -param ../../parametre/etude/params_etude_CLF_z.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_z_lin.ves -ves_out ves_calcul_CLF_z_post_lin.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_lambda_x ves_calcul_CLG_x_post_lin.ves "
                                                                                                 "-ves_lambda_y ves_calcul_CLG_y_post_lin.ves "
                                                                                                 "-ves_lambda_z ves_calcul_CLG_z_post_lin.ves "
                                                                                                 "-ves_out ves_CLG_conductivite_thermique_post_lin.ves "
                                                                                                 "-param ../../parametre/post_traitement/params_post_conductivite_thermique.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_lambda_x ves_calcul_CLF_x_post_lin.ves "
                                                                                                 "-ves_lambda_y ves_calcul_CLF_y_post_lin.ves "
                                                                                                 "-ves_lambda_z ves_calcul_CLF_z_post_lin.ves "
                                                                                                 "-ves_out ves_CLF_conductivite_thermique_post_lin.ves "
                                                                                                 "-param ../../parametre/post_traitement/params_post_conductivite_thermique.txt" << std::endl;                                                                                       
      }
    }
    
    if(fem_maill_quad)
    {
      if(etude_mecanique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLDH_sph_quad.ves -ves_out ves_calcul_CLDH_sph_quad.ves -param ../../parametre/etude/param_etude_CLDH_sph.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_sph_quad.ves -ves_out ves_calcul_CLDH_sph_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_sph_post_quad.ves -ves_out ves_calcul_CLDH_sph_ecart_chargement_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLCH_sph_quad.ves -ves_out ves_calcul_CLCH_sph_quad.ves -param ../../parametre/etude/param_etude_CLCH_sph.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_sph_quad.ves -ves_out ves_calcul_CLCH_sph_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_sph_post_quad.ves -ves_out ves_calcul_CLCH_sph_ecart_chargement_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLDH_dev_quad.ves -ves_out ves_calcul_CLDH_dev_quad.ves -param ../../parametre/etude/param_etude_CLDH_dev.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_dev_quad.ves -ves_out ves_calcul_CLDH_dev_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLDH_dev_post_quad.ves -ves_out ves_calcul_CLDH_dev_ecart_chargement_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLCH_dev_quad.ves -ves_out ves_calcul_CLCH_dev_quad.ves -param ../../parametre/etude/param_etude_CLCH_dev.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_dev_quad.ves -ves_out ves_calcul_CLCH_dev_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca.txt" << std::endl;      
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLCH_dev_post_quad.ves -ves_out ves_calcul_CLCH_dev_ecart_chargement_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt" << std::endl;     
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_sph ves_calcul_CLDH_sph_post_quad.ves -ves_dev ves_calcul_CLDH_dev_post_quad.ves -ves_out ves_CLDH_modules_elasticite_post_quad.ves -param ../../parametre/post_traitement/params_post_modules_elasticite.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_sph ves_calcul_CLCH_sph_post_quad.ves -ves_dev ves_calcul_CLCH_dev_post_quad.ves -ves_out ves_CLCH_modules_elasticite_post_quad.ves -param ../../parametre/post_traitement/params_post_modules_elasticite.txt" << std::endl;
      }
      if(etude_thermique)
      {
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_x_quad.ves -ves_out ves_calcul_CLG_x_quad.ves -param ../../parametre/etude/params_etude_CLG_x.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_x_quad.ves -ves_out ves_calcul_CLG_x_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_y_quad.ves -ves_out ves_calcul_CLG_y_quad.ves -param ../../parametre/etude/params_etude_CLG_y.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_y_quad.ves -ves_out ves_calcul_CLG_y_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLG_z_quad.ves -ves_out ves_calcul_CLG_z_quad.ves -param ../../parametre/etude/params_etude_CLG_z.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLG_z_quad.ves -ves_out ves_calcul_CLG_z_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_x_quad.ves -ves_out ves_calcul_CLF_x_quad.ves -param ../../parametre/etude/params_etude_CLF_x.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_x_quad.ves -ves_out ves_calcul_CLF_x_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_y_quad.ves -ves_out ves_calcul_CLF_y_quad.ves -param ../../parametre/etude/params_etude_CLF_y.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_y_quad.ves -ves_out ves_calcul_CLF_y_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_calcul -ves_in ves_etude_CLF_z_quad.ves -ves_out ves_calcul_CLF_z_quad.ves -param ../../parametre/etude/params_etude_CLF_z.txt -paramaster ../../parametre/etude/param_aster.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_in ves_calcul_CLF_z_quad.ves -ves_out ves_calcul_CLF_z_post_quad.ves -lstparam ../../parametre/post_traitement/lstparam_post_champ_ther.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_lambda_x ves_calcul_CLG_x_post_quad.ves "
                                                                                                 "-ves_lambda_y ves_calcul_CLG_y_post_quad.ves "
                                                                                                 "-ves_lambda_z ves_calcul_CLG_z_post_quad.ves "
                                                                                                 "-ves_out ves_CLG_conductivite_thermique_post_quad.ves "
                                                                                                 "-param ../../parametre/post_traitement/params_post_conductivite_thermique.txt" << std::endl;
        ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_post_traitement -ves_lambda_x ves_calcul_CLF_x_post_quad.ves "
                                                                                                 "-ves_lambda_y ves_calcul_CLF_y_post_quad.ves "
                                                                                                 "-ves_lambda_z ves_calcul_CLF_z_post_quad.ves "
                                                                                                 "-ves_out ves_CLF_conductivite_thermique_post_quad.ves "
                                                                                                 "-param ../../parametre/post_traitement/params_post_conductivite_thermique.txt" << std::endl;
      }
    }    
    ofstrm << "" << std::endl;
    ofstrm << "rm *.msh" << std::endl;
    ofstrm << "rm *.sol" << std::endl;
    ofstrm << "rm *.resu" << std::endl;
    ofstrm << "rm *.comm" << std::endl;
    ofstrm << "rm *.mail" << std::endl;
    ofstrm << "rm *.retour" << std::endl;
    ofstrm << "rm *.export" << std::endl;
    ofstrm << "rm *.log" << std::endl;
    ofstrm << "rm *.mess" << std::endl;
    ofstrm << "#find . ! -name '*.ves' -type f -exec rm -f {} +" << std::endl;
    ofstrm << "cd .." << std::endl;
    ofstrm << "echo \"**********************\"" << std::endl;
    ofstrm << "echo $STATUSTXT" << std::endl;
    ofstrm << "echo \"**********************\"" << std::endl;
    ofstrm << "exit $STAT" << std::endl;
    ofstrm.close();
    if(fonc_affiche!=NULL) fonc_affiche((char*)"ves");
    system((char*)"chmod 744 ves");
  
  return OK;
}

int PARAMETRES::generer_script_post_traitement(bool etude_mecanique,
                                               bool etude_thermique,
                                               int type_inclusion,
                                               bool fem_maill_lin,
                                               bool fem_maill_quad,
                                               double epaisseur_couche_erosion,
                                               fonction_affiche *fonc_affiche)
{
  std::ofstream ofstrm((char*)"post_traitement",std::ios::out|std::ios::trunc);
  ofstrm.precision(16);
  ofstrm.setf(std::ios::showpoint);
  ofstrm << "#!/bin/bash" << std::endl;
  ofstrm << "CHEMINCOMPILE=$1" << std::endl;
  ofstrm << "if [ -z ${CHEMINCOMPILE// } ]" << std::endl;
  ofstrm << "then" << std::endl;
  ofstrm << "  if [ -z $PATHMAGICEXEAPP ]" << std::endl;
  ofstrm << "  then" << std::endl;
  ofstrm << "    echo \"Erreur : Chemin vers exe/app manquant\"" << std::endl;
  ofstrm << "    exit" << std::endl;
  ofstrm << "  else" << std::endl;
  ofstrm << "    CHEMINCOMPILE=$PATHMAGICEXEAPP" << std::endl;
  ofstrm << "  fi" << std::endl;    
  ofstrm << "fi" << std::endl;
  ofstrm << "function executer" << std::endl;
  ofstrm << "{" << std::endl;
  ofstrm << "$*" << std::endl;
  ofstrm << "VAR=$?" << std::endl;
  ofstrm << "if (( $VAR )) ; then " << std::endl;
  ofstrm << "    STATUSTXT=\"ERREUR dans les tests\"; " << std::endl;
  ofstrm << "    STAT=$VAR;" << std::endl;
  ofstrm << "    exit $STAT" << std::endl;
  ofstrm << "fi" << std::endl;
  ofstrm << "}" << std::endl;
  ofstrm << "STATUSTXT=\"Test sans ERREUR\"" << std::endl;
  ofstrm << "echo \"*******************************\"" << std::endl;
  ofstrm << "echo \"Generation d'une microstructure\"" << std::endl;
  ofstrm << "echo \"*******************************\"" << std::endl;
  ofstrm << "STAT=$?" << std::endl;
  ofstrm << "" << std::endl;
  ofstrm << "mkdir resultats" << std::endl;
  ofstrm << "" << std::endl;
  ofstrm << "#Geometrie" << std::endl;
  ofstrm << "mkdir -p resultats/cao" << std::endl;
  ofstrm << "find etude/ -name \"ves_geo_post.ves\" | sort > tmp_lst_ves_geo_post.txt" << std::endl;
  ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_geo_post.txt -lstparam parametre/post_traitement/lstparam_post_geometrie.txt -dossier resultats/cao" << std::endl;
  ofstrm << "rm tmp_lst_ves_geo_post.txt" << std::endl;
  ofstrm << "" << std::endl;
  ofstrm << "#MG Maillage" << std::endl;
  ofstrm << "mkdir -p resultats/maillage_mg" << std::endl;
  ofstrm << "find etude/ -name \"ves_mg_maill_post.ves\" | sort > tmp_lst_ves_mg_maill_post.txt" << std::endl;
  ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_mg_maill_post.txt -lstparam parametre/post_traitement/lstparam_post_mg_maillage.txt -dossier resultats/maillage_mg" << std::endl;
  ofstrm << "rm tmp_lst_ves_mg_maill_post.txt" << std::endl;
  ofstrm << "" << std::endl;
  ofstrm << "#Fem Maillage" << std::endl;
  if(fem_maill_lin)
  {
    ofstrm << "mkdir -p resultats/maillage_fem/lin" << std::endl;
    ofstrm << "find etude/ -name \"ves_fem_maill_post_lin.ves\" | sort  > tmp_lst_ves_fem_maill_post_lin.txt" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_fem_maill_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_fem_maillage.txt -dossier resultats/maillage_fem/lin" << std::endl;
    ofstrm << "rm tmp_lst_ves_fem_maill_post_lin.txt" << std::endl;
  }
  if(fem_maill_quad)
  {
    ofstrm << "mkdir -p resultats/maillage_fem/quad" << std::endl;
    ofstrm << "find etude/ -name \"ves_fem_maill_post_quad.ves\" | sort  > tmp_lst_ves_fem_maill_post_quad.txt" << std::endl;
    ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_fem_maill_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_fem_maillage.txt -dossier resultats/maillage_fem/quad" << std::endl;
    ofstrm << "rm tmp_lst_ves_fem_maill_post_quad.txt" << std::endl;
  }
  ofstrm << "" << std::endl;
  ofstrm << "#Calculs" << std::endl;
  if(fem_maill_lin)
  {
    if(etude_mecanique)
    {
      ofstrm << "mkdir -p resultats/champ/DH/sph/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/DH/dev/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/CH/sph/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/CH/dev/lin" << std::endl;
      ofstrm << "mkdir -p resultats/modules_elasticite/DH/lin" << std::endl;
      ofstrm << "mkdir -p resultats/modules_elasticite/CH/lin" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_sph_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLDH_sph_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_sph_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLCH_sph_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_dev_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLDH_dev_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_dev_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLCH_dev_post_lin.txt" << std::endl;    
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_sph_ecart_chargement_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_sph_ecart_chargement_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_dev_ecart_chargement_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_dev_ecart_chargement_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_lin.txt" << std::endl;    
      ofstrm << "find etude/ -name \"ves_CLDH_modules_elasticite_post_lin.ves\" | sort > tmp_lst_ves_CLDH_modules_elasticite_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_CLCH_modules_elasticite_post_lin.ves\" | sort > tmp_lst_ves_CLCH_modules_elasticite_post_lin.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_sph_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/DH/sph/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_sph_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/CH/sph/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_dev_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/DH/dev/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_dev_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/CH/dev/lin" << std::endl;    
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt -dossier resultats/champ/DH/sph/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt -dossier resultats/champ/CH/sph/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt -dossier resultats/champ/DH/dev/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt -dossier resultats/champ/CH/dev/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLDH_modules_elasticite_post_lin.txt -param parametre/post_traitement/params_post_modules_elasticite.txt -dossier resultats/modules_elasticite/DH/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLCH_modules_elasticite_post_lin.txt -param parametre/post_traitement/params_post_modules_elasticite.txt -dossier resultats/modules_elasticite/CH/lin" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_sph_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_sph_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_dev_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_dev_post_lin.txt" << std::endl;    
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLDH_modules_elasticite_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLCH_modules_elasticite_post_lin.txt" << std::endl;
    }
    if(etude_thermique)
    {
      ofstrm << "mkdir -p resultats/champ/FH/F_x/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/FH/F_y/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/FH/F_z/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_x/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_y/lin" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_z/lin" << std::endl;
      ofstrm << "mkdir -p resultats/conductivite_thermique/FH/lin" << std::endl;
      ofstrm << "mkdir -p resultats/conductivite_thermique/GH/lin" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_x_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLF_x_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_y_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLF_y_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_z_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLF_z_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_x_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLG_x_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_y_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLG_y_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_z_post_lin.ves\" | sort > tmp_lst_ves_calcul_CLG_z_post_lin.txt" << std::endl;      
      ofstrm << "find etude/ -name \"ves_CLF_conductivite_thermique_post_lin.ves\" | sort > tmp_lst_ves_CLF_conductivite_thermique_post_lin.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_CLG_conductivite_thermique_post_lin.ves\" | sort > tmp_lst_ves_CLG_conductivite_thermique_post_lin.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_x_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_x/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_y_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_y/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_z_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_z/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_x_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_x/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_y_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_y/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_z_post_lin.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_z/lin" << std::endl;      
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLF_conductivite_thermique_post_lin.txt -param parametre/post_traitement/params_post_conductivite_thermique.txt -dossier resultats/conductivite_thermique/FH/lin" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLG_conductivite_thermique_post_lin.txt -param parametre/post_traitement/params_post_conductivite_thermique.txt -dossier resultats/conductivite_thermique/GH/lin" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_x_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_y_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_z_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_x_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_y_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_z_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLF_conductivite_thermique_post_lin.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLG_conductivite_thermique_post_lin.txt" << std::endl; 
    }
  }
  if(fem_maill_quad)
  {
    if(etude_mecanique)
    {
      ofstrm << "mkdir -p resultats/champ/DH/sph/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/DH/dev/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/CH/sph/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/CH/dev/quad" << std::endl;
      ofstrm << "mkdir -p resultats/modules_elasticite/DH/quad" << std::endl;
      ofstrm << "mkdir -p resultats/modules_elasticite/CH/quad" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_sph_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLDH_sph_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_sph_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLCH_sph_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_dev_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLDH_dev_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_dev_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLCH_dev_post_quad.txt" << std::endl;    
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_sph_ecart_chargement_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_sph_ecart_chargement_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLDH_dev_ecart_chargement_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLCH_dev_ecart_chargement_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_quad.txt" << std::endl;    
      ofstrm << "find etude/ -name \"ves_CLDH_modules_elasticite_post_quad.ves\" | sort > tmp_lst_ves_CLDH_modules_elasticite_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_CLCH_modules_elasticite_post_quad.ves\" | sort > tmp_lst_ves_CLCH_modules_elasticite_post_quad.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_sph_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/DH/sph/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_sph_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/CH/sph/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_dev_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/DH/dev/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_dev_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca.txt -dossier resultats/champ/CH/dev/quad" << std::endl;    
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt -dossier resultats/champ/DH/sph/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_sph.txt -dossier resultats/champ/CH/sph/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt -dossier resultats/champ/DH/dev/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_meca_ecart_chargement_dev.txt -dossier resultats/champ/CH/dev/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLDH_modules_elasticite_post_quad.txt -param parametre/post_traitement/params_post_modules_elasticite.txt -dossier resultats/modules_elasticite/DH/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLCH_modules_elasticite_post_quad.txt -param parametre/post_traitement/params_post_modules_elasticite.txt -dossier resultats/modules_elasticite/CH/quad" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_sph_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_sph_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_dev_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_dev_post_quad.txt" << std::endl;    
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_sph_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_sph_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLDH_dev_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLCH_dev_ecart_chargement_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLDH_modules_elasticite_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLCH_modules_elasticite_post_quad.txt" << std::endl;
    }
    if(etude_thermique)
    {
      ofstrm << "mkdir -p resultats/champ/FH/F_x/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/FH/F_y/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/FH/F_z/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_x/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_y/quad" << std::endl;
      ofstrm << "mkdir -p resultats/champ/GH/G_z/quad" << std::endl;
      ofstrm << "mkdir -p resultats/conductivite_thermique/FH/quad" << std::endl;
      ofstrm << "mkdir -p resultats/conductivite_thermique/GH/quad" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_x_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLF_x_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_y_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLF_y_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLF_z_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLF_z_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_x_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLG_x_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_y_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLG_y_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_calcul_CLG_z_post_quad.ves\" | sort > tmp_lst_ves_calcul_CLG_z_post_quad.txt" << std::endl;      
      ofstrm << "find etude/ -name \"ves_CLF_conductivite_thermique_post_quad.ves\" | sort > tmp_lst_ves_CLF_conductivite_thermique_post_quad.txt" << std::endl;
      ofstrm << "find etude/ -name \"ves_CLG_conductivite_thermique_post_quad.ves\" | sort > tmp_lst_ves_CLG_conductivite_thermique_post_quad.txt" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_x_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_x/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_y_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_y/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLF_z_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/FH/F_z/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_x_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_x/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_y_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_y/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_calcul_CLG_z_post_quad.txt -lstparam parametre/post_traitement/lstparam_post_champ_ther.txt -dossier resultats/champ/GH/G_z/quad" << std::endl;      
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLF_conductivite_thermique_post_quad.txt -param parametre/post_traitement/params_post_conductivite_thermique.txt -dossier resultats/conductivite_thermique/FH/quad" << std::endl;
      ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -cumuler_post_traitement -lstves tmp_lst_ves_CLG_conductivite_thermique_post_quad.txt -param parametre/post_traitement/params_post_conductivite_thermique.txt -dossier resultats/conductivite_thermique/GH/quad" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_x_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_y_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLF_z_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_x_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_y_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_calcul_CLG_z_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLF_conductivite_thermique_post_quad.txt" << std::endl;
      ofstrm << "rm tmp_lst_ves_CLG_conductivite_thermique_post_quad.txt" << std::endl; 
    }
  }
  ofstrm << "" << std::endl;
  ofstrm << "./produire_figures 0" << std::endl;
  ofstrm << "echo \"**********************\"" << std::endl;
  ofstrm << "echo \"        FIN           \"" << std::endl;
  ofstrm << "echo \"**********************\"" << std::endl;
  ofstrm << "exit $STAT" << std::endl;
  ofstrm.close();
  if(fonc_affiche!=NULL) fonc_affiche((char*)"post_traitement");
  system((char*)"chmod 744 post_traitement");
  
  
  ofstrm.open((char*)"produire_figures",std::ios::out|std::ios::trunc);
  ofstrm.precision(16);
  ofstrm.setf(std::ios::showpoint);
  ofstrm << "#!/bin/bash" << std::endl;
  ofstrm << "COUCHE_REF=\"$1\"" << std::endl;
  ofstrm << "#Data pour graph" << std::endl;
  ofstrm << "PWD=pwd" << std::endl;
  ofstrm << "mkdir -p resultats/data_graph resultats/graph" << std::endl;
  if(type_inclusion==MSTRUCT::TYPE_INCLUSION::CYLINDRE)
  {
    ofstrm << "#Orientation" << std::endl;
    ofstrm << "awk 'NR>=2 {print $1,($2+$3+$4+$5+$6+$7)}'  resultats/cao/cumul_Orientation_Particule_ecart.txt > resultats/data_graph/cumul_Err_abs_ori.txt" << std::endl;
    ofstrm << "sed -i '1s/^/Nb_ves somme_ecart_absolu\\n/' resultats/data_graph/cumul_Err_abs_ori.txt" << std::endl;
    if(fem_maill_lin)
    {
      ofstrm << "awk 'NR>=2 {print $1,$2,($2+$3+$4+$5+$6+$7)}'  resultats/maillage_fem/lin/cumul_erosion_Orientation_Particule_ponderee_ecart.txt > resultats/data_graph/cumul_erosion_Err_abs_ori_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves somme_ecart_absolu\\n/' resultats/data_graph/cumul_erosion_Err_abs_ori_lin.txt" << std::endl;
    }
    if(fem_maill_quad)
    {
      ofstrm << "awk 'NR>=2 {print $1,$2,($2+$3+$4+$5+$6+$7)}'  resultats/maillage_fem/quad/cumul_erosion_Orientation_Particule_ponderee_ecart.txt > resultats/data_graph/cumul_erosion_Err_abs_ori_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves somme_ecart_absolu\\n/' resultats/data_graph/cumul_erosion_Err_abs_ori_quad.txt" << std::endl;
    }
  }  
  if(fem_maill_lin)
  {
    ofstrm << "awk 'NR>=2 {print $1,$2,$7,$8}'  resultats/maillage_fem/lin/cumul_erosion_FEM_maillage_Particule.txt > resultats/data_graph/tmp_erosion_fem_maill_lin" << std::endl;
    ofstrm << "awk 'NR>=2 {print $1,$3}'  resultats/maillage_fem/lin/c_$COUCHE_REF/liste_FEM_maillage_Particule.txt > resultats/data_graph/tmp_liste_fem_maill_lin" << std::endl;
    if(etude_mecanique)
    {
      ofstrm << "#Module elasticite" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4,$7,$8,$11,$12,$15,$16}'  resultats/modules_elasticite/CH/lin/cumul_erosion_Modules_elasticite.txt  > resultats/data_graph/tmp_erosion_Modules_app_CH_lin" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4,$7,$8,$11,$12,$15,$16}'  resultats/modules_elasticite/DH/lin/cumul_erosion_Modules_elasticite.txt  > resultats/data_graph/tmp_erosion_Modules_app_DH_lin" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_lin resultats/data_graph/tmp_erosion_Modules_app_CH_lin > resultats/data_graph/erosion_Modules_app_CH_lin.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_lin resultats/data_graph/tmp_erosion_Modules_app_DH_lin > resultats/data_graph/erosion_Modules_app_DH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Eapp(moy ±) Nuapp(moy ±) Gapp(moy ±) Kapp(moy ±)\\n/' resultats/data_graph/erosion_Modules_app_CH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Eapp(moy ±) Nuapp(moy ±) Gapp(moy ±) Kapp(moy ±)\\n/' resultats/data_graph/erosion_Modules_app_DH_lin.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_erosion_Modules_app_CH_lin resultats/data_graph/tmp_erosion_Modules_app_DH_lin" << std::endl;    
      ofstrm << "awk 'NR>=2 {print $2,$3,$4,$5}'  resultats/modules_elasticite/CH/lin/c_$COUCHE_REF/liste_Modules_elasticite.txt > resultats/data_graph/tmp_Modules_app_CH_lin" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2,$3,$4,$5}'  resultats/modules_elasticite/DH/lin/c_$COUCHE_REF/liste_Modules_elasticite.txt > resultats/data_graph/tmp_Modules_app_DH_lin" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_lin resultats/data_graph/tmp_Modules_app_CH_lin > resultats/data_graph/liste_Modules_app_CH_lin.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_lin resultats/data_graph/tmp_Modules_app_DH_lin > resultats/data_graph/liste_Modules_app_DH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Eapp Nuapp Gapp Kapp\\n/' resultats/data_graph/liste_Modules_app_CH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Eapp Nuapp Gapp Kapp\\n/' resultats/data_graph/liste_Modules_app_DH_lin.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_Modules_app_CH_lin resultats/data_graph/tmp_Modules_app_DH_lin" << std::endl;
    }
    if(etude_thermique)
    {
      ofstrm << "#Conductivite thermique" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4}'  resultats/conductivite_thermique/FH/lin/cumul_erosion_Conductivite_thermique.txt  > resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_lin" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4}'  resultats/conductivite_thermique/GH/lin/cumul_erosion_Conductivite_thermique.txt  > resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_lin" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_lin resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_lin > resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_lin resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_lin > resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Lambda(moy ±)\\n/' resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Lambda(moy ±)\\n/' resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_lin resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_lin" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2}'  resultats/conductivite_thermique/FH/lin/c_$COUCHE_REF/liste_Conductivite_thermique.txt > resultats/data_graph/tmp_Conductivite_thermique_app_FH_lin" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2}'  resultats/conductivite_thermique/GH/lin/c_$COUCHE_REF/liste_Conductivite_thermique.txt > resultats/data_graph/tmp_Conductivite_thermique_app_GH_lin" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_lin resultats/data_graph/tmp_Conductivite_thermique_app_FH_lin > resultats/data_graph/liste_Conductivite_thermique_app_FH_lin.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_lin resultats/data_graph/tmp_Conductivite_thermique_app_GH_lin > resultats/data_graph/liste_Conductivite_thermique_app_GH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Lambda_app\\n/' resultats/data_graph/liste_Conductivite_thermique_app_FH_lin.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Lambda_app\\n/' resultats/data_graph/liste_Conductivite_thermique_app_GH_lin.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_Conductivite_thermique_app_FH_lin resultats/data_graph/tmp_Conductivite_thermique_app_GH_lin" << std::endl;
    }
    ofstrm << "rm resultats/data_graph/tmp_erosion_fem_maill_lin resultats/data_graph/tmp_liste_fem_maill_lin" << std::endl;
  }
  if(fem_maill_quad)
  {
    ofstrm << "awk 'NR>=2 {print $1,$2,$7,$8}'  resultats/maillage_fem/quad/cumul_erosion_FEM_maillage_Particule.txt > resultats/data_graph/tmp_erosion_fem_maill_quad" << std::endl;
    ofstrm << "awk 'NR>=2 {print $1,$3}'  resultats/maillage_fem/quad/c_$COUCHE_REF/liste_FEM_maillage_Particule.txt > resultats/data_graph/tmp_liste_fem_maill_quad" << std::endl;
    if(etude_mecanique)
    {
      ofstrm << "#Module elasticite" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4,$7,$8,$11,$12,$15,$16}'  resultats/modules_elasticite/CH/quad/cumul_erosion_Modules_elasticite.txt  > resultats/data_graph/tmp_erosion_Modules_app_CH_quad" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4,$7,$8,$11,$12,$15,$16}'  resultats/modules_elasticite/DH/quad/cumul_erosion_Modules_elasticite.txt  > resultats/data_graph/tmp_erosion_Modules_app_DH_quad" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_quad resultats/data_graph/tmp_erosion_Modules_app_CH_quad > resultats/data_graph/erosion_Modules_app_CH_quad.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_quad resultats/data_graph/tmp_erosion_Modules_app_DH_quad > resultats/data_graph/erosion_Modules_app_DH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Eapp(moy ±) Nuapp(moy ±) Gapp(moy ±) Kapp(moy ±)\\n/' resultats/data_graph/erosion_Modules_app_CH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Eapp(moy ±) Nuapp(moy ±) Gapp(moy ±) Kapp(moy ±)\\n/' resultats/data_graph/erosion_Modules_app_DH_quad.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_erosion_Modules_app_CH_quad resultats/data_graph/tmp_erosion_Modules_app_DH_quad" << std::endl;    
      ofstrm << "awk 'NR>=2 {print $2,$3,$4,$5}'  resultats/modules_elasticite/CH/quad/c_$COUCHE_REF/liste_Modules_elasticite.txt > resultats/data_graph/tmp_Modules_app_CH_quad" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2,$3,$4,$5}'  resultats/modules_elasticite/DH/quad/c_$COUCHE_REF/liste_Modules_elasticite.txt > resultats/data_graph/tmp_Modules_app_DH_quad" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_quad resultats/data_graph/tmp_Modules_app_CH_quad > resultats/data_graph/liste_Modules_app_CH_quad.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_quad resultats/data_graph/tmp_Modules_app_DH_quad > resultats/data_graph/liste_Modules_app_DH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Eapp Nuapp Gapp Kapp\\n/' resultats/data_graph/liste_Modules_app_CH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Eapp Nuapp Gapp Kapp\\n/' resultats/data_graph/liste_Modules_app_DH_quad.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_Modules_app_CH_quad resultats/data_graph/tmp_Modules_app_DH_quad" << std::endl;
    }
    if(etude_thermique)
    {
      ofstrm << "#Conductivite thermique" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4}'  resultats/conductivite_thermique/FH/quad/cumul_erosion_Conductivite_thermique.txt  > resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_quad" << std::endl;
      ofstrm << "awk 'NR>=2 {print $3,$4}'  resultats/conductivite_thermique/GH/quad/cumul_erosion_Conductivite_thermique.txt  > resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_quad" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_quad resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_quad > resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_erosion_fem_maill_quad resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_quad > resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Lambda(moy ±)\\n/' resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/Couche Nb_ves Fvol(moy ±) Lambda(moy ±)\\n/' resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_erosion_Conductivite_thermique_app_FH_quad resultats/data_graph/tmp_erosion_Conductivite_thermique_app_GH_quad" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2}'  resultats/conductivite_thermique/FH/quad/c_$COUCHE_REF/liste_Conductivite_thermique.txt > resultats/data_graph/tmp_Conductivite_thermique_app_FH_quad" << std::endl;
      ofstrm << "awk 'NR>=2 {print $2}'  resultats/conductivite_thermique/GH/quad/c_$COUCHE_REF/liste_Conductivite_thermique.txt > resultats/data_graph/tmp_Conductivite_thermique_app_GH_quad" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_quad resultats/data_graph/tmp_Conductivite_thermique_app_FH_quad > resultats/data_graph/liste_Conductivite_thermique_app_FH_quad.txt" << std::endl;
      ofstrm << "paste -d ' ' resultats/data_graph/tmp_liste_fem_maill_quad resultats/data_graph/tmp_Conductivite_thermique_app_GH_quad > resultats/data_graph/liste_Conductivite_thermique_app_GH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Lambda_app\\n/' resultats/data_graph/liste_Conductivite_thermique_app_FH_quad.txt" << std::endl;
      ofstrm << "sed -i '1s/^/# Fvol Lambda_app\\n/' resultats/data_graph/liste_Conductivite_thermique_app_GH_quad.txt" << std::endl;
      ofstrm << "rm resultats/data_graph/tmp_Conductivite_thermique_app_FH_quad resultats/data_graph/tmp_Conductivite_thermique_app_GH_quad" << std::endl;
    }
    ofstrm << "rm resultats/data_graph/tmp_erosion_fem_maill_quad resultats/data_graph/tmp_liste_fem_maill_quad" << std::endl;
  }
  ofstrm << "gnuplot fig.gnu" << std::endl;
  if(fonc_affiche!=NULL) fonc_affiche((char*)"produire_figures");
  system((char*)"chmod 744 produire_figures");
  ofstrm.close();

  
  ofstrm.open((char*)"fig.gnu",std::ios::out|std::ios::trunc);
  ofstrm.precision(16);
  ofstrm.setf(std::ios::showpoint);
  ofstrm << "" << std::endl;
  ofstrm << "set terminal postscript eps size 4,2.25 enhanced color font 'Helvetica,14'" << std::endl;
  
  if(etude_mecanique)
  {
    ofstrm << "" << std::endl;
    ofstrm << "set output \"resultats/graph/nuage_Eapp.eps\"" << std::endl;
    ofstrm << "set xrange [0:50]" << std::endl;
    ofstrm << "set title \"Module de Young (GPa)\"" << std::endl;
    ofstrm << "set xlabel \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ylabel \"Module de Young (GPa)\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Modules_app_DH_lin.txt\" u ($2*100):($3/(1000000000)) pt 4 ps 2 lc rgb 'blue' lw 1 title \"DH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Modules_app_CH_lin.txt\" u ($2*100):($3/(1000000000)) pt 4 ps 2 lc rgb 'red' lw 1 title \"CH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Modules_app_DH_quad.txt\" u ($2*100):($3/(1000000000)) pt 5 ps 2 lc rgb 'blue' lw 1 title \"DH quad\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Modules_app_CH_quad.txt\" u ($2*100):($3/(1000000000)) pt 5 ps 2 lc rgb 'red' lw 1 title \"DH quad\",\\" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Modules_app_DH_lin.txt\" u ($2*100):($3/(1000000000)) pt 4 ps 2 lc rgb 'blue' lw 1 title \"DH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Modules_app_CH_lin.txt\" u ($2*100):($3/(1000000000)) pt 4 ps 2 lc rgb 'red' lw 1 title \"CH lin\",\\" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Modules_app_DH_quad.txt\" u ($2*100):($3/(1000000000)) pt 5 ps 2 lc rgb 'blue' lw 1 title \"DH quad\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Modules_app_CH_quad.txt\" u ($2*100):($3/(1000000000)) pt 5 ps 2 lc rgb 'red' lw 1 title \"DH quad\",\\" << std::endl;
    }
    ofstrm << "\"resE.dat\" u($1*100):($2) with line lw 2 lc rgb 'black' lt 2 title \"E_{Reuss}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($3) with line lw 2 lc rgb 'grey' lt 2 title \"E_{HS-}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($7) with line lw 2 lc rgb 'grey' lt 2 dashtype 2 title \"E_{HS+}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($8) with line lw 2 lc rgb 'black' lt 2 dashtype 2 title \"E_{Voigt}\"" << std::endl;
    ofstrm << "#pause -1" << std::endl;    
    
    ofstrm << "" << std::endl;
    ofstrm << "set output \"resultats/graph/cumul_Eapp.eps\"" << std::endl;
    ofstrm << "set title \"Module de Young (GPa)\"" << std::endl;
    ofstrm << "set xlabel \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ylabel \"Module de Young (GPa)\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_DH_lin.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"DH lin\",\\" << std::endl;    
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_CH_lin.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"CH lin\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_DH_quad.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"DH quad\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_CH_quad.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"DH quad\",\\" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_DH_lin.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"DH lin\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_CH_lin.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"CH lin\",\\" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_DH_quad.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"DH quad\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Modules_app_CH_quad.txt)\" u ($3*100):($5/(1000000000)):($4*100):($6/(1000000000)) with xyerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"CH quad\",\\" << std::endl;
    }
    ofstrm << "\"resE.dat\" u($1*100):($2) with line lw 2 lc rgb 'black' lt 2 title \"E_{Reuss}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($3) with line lw 2 lc rgb 'grey' lt 2 title \"E_{HS-}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($7) with line lw 2 lc rgb 'grey' lt 2 dashtype 2 title \"E_{HS+}\" ,\\" << std::endl;
    ofstrm << "\"resE.dat\" u($1*100):($8) with line lw 2 lc rgb 'black' lt 2 dashtype 2 title \"E_{Voigt}\"" << std::endl;
    ofstrm << "#pause -1" << std::endl;    
    
    ofstrm << "set output \"resultats/graph/cumul_erosion_Eapp.eps\"" << std::endl;
    ofstrm << "set xrange [0:0.5]" << std::endl;
    ofstrm << "set y2range[0:100]" << std::endl;
    ofstrm << "set title \"Module de Young (GPa)\"" << std::endl;
    ofstrm << "set xlabel \"Distance d'erosion\"" << std::endl;
    ofstrm << "set ylabel \"Module de Young (GPa)\"" << std::endl;
    ofstrm << "set y2label \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ytics nomirror" << std::endl;
    ofstrm << "set y2tics" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Modules_app_DH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'blue' lw 1 title \"DH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_DH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'red' lw 1 title \"CH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_DH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'blue' lw 1 title \"DH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_DH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'red' lw 1 title \"CH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH quad\" axes x1y1 ,\\" << std::endl;  
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'black' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'black'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Modules_app_DH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'blue' lw 1 title \"DH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_DH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'red' lw 1 title \"CH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2" << std::endl;
    }
    if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Modules_app_DH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'blue' lw 1 title \"DH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_DH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)) with lines lc rgb 'red' lw 1 title \"CH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5/(1000000000)):($6/(1000000000)) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH quad\" axes x1y1 ,\\" << std::endl;  
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'black' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Modules_app_CH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'black'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2" << std::endl;
    }
    ofstrm << "unset y2label" << std::endl;
    ofstrm << "unset y2tics" << std::endl;
    ofstrm << "#pause -1" << std::endl;  
    
    ofstrm << "set output \"resultats/graph/cumul_erosion_ecart_energie.eps\"" << std::endl;
    ofstrm << "set xrange [0:0.5]" << std::endl;
    ofstrm << "set title \"Ecart d'energie relatif absolue (%)\"" << std::endl;
    ofstrm << "set xlabel \"Distance d'erosion\"" << std::endl;
    ofstrm << "set ylabel \"Ecart d'energie relatif absolue (%)\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH dev quad\",\\" << std::endl;    
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH sph quad\"" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH dev lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH sph lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH sph lin\"" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH dev quad\",\\" << std::endl;    
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"erreur CH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH dev quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH sph quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_ENERGIE_ecart.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"erreur DH sph quad\"" << std::endl;
    }
    ofstrm << "#pause -1" << std::endl;
    
    
    ofstrm << "set output \"resultats/graph/cumul_erosion_ecart_chargement_dev.eps\"" << std::endl;
    ofstrm << "set xrange [0:0.5]" << std::endl;
    ofstrm << "set title \"Ecart de chargement relatif absolue (%)\"" << std::endl;
    ofstrm << "set xlabel \"Distance d'erosion\"" << std::endl;
    ofstrm << "set ylabel \"Ecart de chargement relatif absolue (%)\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA quad\",\\" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA lin\"" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/dev/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_dev.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA quad\"" << std::endl;
    }    
    ofstrm << "#pause -1" << std::endl;
    
    
    ofstrm << "set output \"resultats/graph/cumul_erosion_ecart_chargement_sph.eps\"" << std::endl;
    ofstrm << "set xrange [0:0.5]" << std::endl;
    ofstrm << "set title \"Ecart de chargement relatif absolue (%)\"" << std::endl;
    ofstrm << "set xlabel \"Distance d'erosion\"" << std::endl;
    ofstrm << "set ylabel \"Ecart de chargement relatif absolue (%)\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA quad\",\\" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA lin\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/lin/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA lin\"" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 1 lw 1 title \"CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'red' dt 2 lw 1 title \"CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/CH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur CH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 1 lw 1 title \"DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_EPSILON_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH EPSILON quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0) with lines lc rgb 'blue' dt 2 lw 1 title \"DH SIGMA quad\",\\" << std::endl;
      ofstrm << "     \"resultats/champ/DH/sph/quad/cumul_erosion_Champ_SIGMA_ecart_chargement_sph.txt\" u ($1*0.02):($3*100.0):($4*100.0) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur DH SIGMA quad\"" << std::endl;
    }    
    ofstrm << "#pause -1" << std::endl;    
  }
  if(etude_thermique)
  {
    ofstrm << "" << std::endl;
    ofstrm << "set output \"resultats/graph/nuage_Lambda_app.eps\"" << std::endl;
    ofstrm << "set xrange [0:50]" << std::endl;
    ofstrm << "set title \"Conductivite thermique (W/(m*K))\"" << std::endl;
    ofstrm << "set xlabel \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ylabel \"Conductivite thermique (W/(m*K))\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Conductivite_thermique_app_FH_lin.txt\" u ($2*100):($3) pt 4 ps 2 lc rgb 'blue' lw 1 title \"FH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Conductivite_thermique_app_GH_lin.txt\" u ($2*100):($3) pt 4 ps 2 lc rgb 'red' lw 1 title \"GH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Conductivite_thermique_app_FH_quad.txt\" u ($2*100):($3) pt 5 ps 2 lc rgb 'blue' lw 1 title \"FH quad\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Conductivite_thermique_app_GH_quad.txt\" u ($2*100):($3) pt 5 ps 2 lc rgb 'red' lw 1 title \"GH quad\"" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Conductivite_thermique_app_FH_lin.txt\" u ($2*100):($3) pt 4 ps 2 lc rgb 'blue' lw 1 title \"FH lin\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Conductivite_thermique_app_GH_lin.txt\" u ($2*100):($3) pt 4 ps 2 lc rgb 'red' lw 1 title \"GH lin\"" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/liste_Conductivite_thermique_app_FH_quad.txt\" u ($2*100):($3) pt 5 ps 2 lc rgb 'blue' lw 1 title \"FH quad\",\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/liste_Conductivite_thermique_app_GH_quad.txt\" u ($2*100):($3) pt 5 ps 2 lc rgb 'red' lw 1 title \"GH quad\"" << std::endl;
    }
    ofstrm << "#pause -1" << std::endl;
    
    ofstrm << "" << std::endl;
    ofstrm << "set output \"resultats/graph/cumul_Lambda_app.eps\"" << std::endl;
    ofstrm << "set title \"Conductivite thermique (W/(m*K))\"" << std::endl;
    ofstrm << "set xlabel \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ylabel \"Conductivite thermique (W/(m*K))\"" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"FH lin\",\\" << std::endl;    
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"GH lin\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"FH quad\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"GH quad\"" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"FH lin\",\\" << std::endl;    
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"GH lin\"" << std::endl;
    }
    else if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 5 ps 2 lc rgb 'blue' lw 1 title \"FH quad\",\\" << std::endl;
      ofstrm << "     \"<(sed -n 2p resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt)\" u ($3*100):($5):($4*100):($6) with xyerrorbars pt 5 ps 2 lc rgb 'red' lw 1 title \"GH quad\"" << std::endl;
    }
    ofstrm << "#pause -1" << std::endl;
    
    
    ofstrm << "set output \"resultats/graph/cumul_erosion_Lambda_app.eps\"" << std::endl;
    ofstrm << "set xrange [0:0.5]" << std::endl;
    ofstrm << "set y2range[0:100]" << std::endl;
    ofstrm << "set title \"Conductivite thermique (W/(m*K))\"" << std::endl;
    ofstrm << "set xlabel \"Distance d'erosion\"" << std::endl;
    ofstrm << "set ylabel \"Conductivite thermique (W/(m*K))\"" << std::endl;
    ofstrm << "set y2label \"Fraction volumique (%)\"" << std::endl;
    ofstrm << "set ytics nomirror" << std::endl;
    ofstrm << "set y2tics" << std::endl;
    if(fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'blue' lw 1 title \"FH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur FH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'red' lw 1 title \"GH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur GH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'blue' lw 1 title \"FH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur FH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'red' lw 1 title \"GH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur GH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2" << std::endl;
    }
    else if(fem_maill_lin && !fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'blue' lw 1 title \"FH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur FH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'red' lw 1 title \"GH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur GH lin\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_lin.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2 " << std::endl;
    }
    if(!fem_maill_lin && fem_maill_quad)
    {
      ofstrm << "plot \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'blue' lw 1 title \"FH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_FH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'blue' lw 1 title \"erreur FH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5) with lines lc rgb 'red' lw 1 title \"GH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($5):($6) with yerrorbars pt 4 ps 2 lc rgb 'red' lw 1 title \"erreur GH quad\" axes x1y1 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100) with linespoint lc rgb 'grey' lt 0 lw 1 pt 15 ps 1 title \"Frac. Vol.\" axes x1y2 ,\\" << std::endl;
      ofstrm << "     \"resultats/data_graph/erosion_Conductivite_thermique_app_GH_quad.txt\" u ($1*" << epaisseur_couche_erosion << "):($3*100):($4*100) with yerrorbars lc rgb 'grey'  pt 15 ps 1 title \"erreur Frac. Vol.\" axes x1y2" << std::endl;
    }
    ofstrm << "unset y2label" << std::endl;
    ofstrm << "unset y2tics" << std::endl;
    ofstrm << "#pause -1" << std::endl;  
    
  }
  ofstrm.close();  
  if(fonc_affiche!=NULL) fonc_affiche((char*)"fig_Eapp.gnu");
}


int PARAMETRES::importer_fichier_parametres(char* fichier, OT_PARAMETRES* params)
{
  return params->lire(fichier);
}

int PARAMETRES::importer_fichier_liste_parametres(char* fichier, std::vector< OT_PARAMETRES* >& vector_params)
{
  std::string str_fichier = fichier;
  std::size_t found = str_fichier.rfind((char*)"/");
  std::string chemin;
  if(found!=-1)
  {
    chemin = str_fichier.substr(0,found+1);
  }
  else chemin = "";
  FILE* in=fopen(fichier,"rt");
  char ligne[4000];
  fgets(ligne,4000,in);
  while(!feof(in))
  {
    char *aide = strstr(ligne,"//");
    if(aide==NULL)
    {
      char sous_fichier[500];
      sscanf(ligne,"%s",sous_fichier);
      if(strlen(sous_fichier)==0) continue;
      OT_PARAMETRES *param = new OT_PARAMETRES;
      std::string str_sous_fichier = chemin;
      str_sous_fichier.append(sous_fichier);
      if(param->lire((char*)str_sous_fichier.c_str())==FAIL) return FAIL;
      vector_params.push_back(param);
    }
    fgets(ligne,4000,in);
  }
  return OK;
}
Revision:	971
Committed:	Thu Sep 20 19:04:41 2018 UTC (6 years, 7 months ago) by couturad
File size:	189342 byte(s)
Log Message:	MG_EXPORT: Changement du facteur de correction dans l'application de la conductivite thermique (corr=1) MICROSTRUCTURE: -> Ajout de l'etude thermique (Calcul de la conductivite thermique apparente) -> Ajout d'une analyse des proprietes massiques (Centre de masse, moment d'inertie)