/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M O D E L E   T R I D I M E N S I O N N E L   D E   C O T E S   A V E C   E R O S I O N  :                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xci/CoastL_3D.11$K' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20011208091400).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        I N T E R F A C E   ' listG '  :                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        :Debut_listG:                                                                                                              */
/*        :Fin_listG:                                                                                                                */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D I R E C T I V E S   S P E C I F I Q U E S   D E   C O M P I L A T I O N  :                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
@define   PRAGMA_CL_____MODULE_NON_OPTIMISABLE

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        F I C H I E R S   D ' I N C L U D E S  :                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  INCLUDES_BASE
#include  image_image_DI_ALBUM_EXT

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xci/sequence.01.I"
                                        /* ATTENTION, on definit ainsi le symbole 'DERNIERE_IMAGE' qui ne sert a rien ici, puisque   */
                                        /* c'est en effet 'Zmax' qui joue ce role...                                                 */

#define   EDITER_LES_PARAMETRES_DU_NOYAU                                                                                                \
                    FAUX                                                                                                                \
                                        /* Faut-il editer les parametres du noyau et sortir immediatement ('VRAI') ou exploiter      */ \
                                        /* le modele d'erosion de cotes tridimensionnel ('FAUX') ?                                   */
#define   LISTER_LE_CONTENU_DU_NOYAU                                                                                                    \
                    FAUX                                                                                                                \
                                        /* Faut-il lister le contenu du noyau ('VRAI') ou pas ('FAUX') ?                             */

#define   IMPOSER_LES_EXTREMA_DES_CONDITIONS_INITIALES                                                                                  \
                    FAUX
#define   MINIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES                                                                                  \
                    NIVEAU_DE_LA_TERRE_DANS_Aerode__CoastLines_3D
#define   MAXIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES                                                                                  \
                    NIVEAU_DE_LA___MER_DANS_Aerode__CoastLines_3D
                                        /* Faut-il imposer ('VRAI') ou pas ('FAUX') les extrema des conditions initiales, et si      */
                                        /* 'IL_FAUT(imposer_les_extrema_des_conditions_initiales)', valeurs implicites ?             */
#define   EDITER_LES_EXTREMA_DES_CONDITIONS_INITIALES                                                                                   \
                    FAUX                                                                                                                \
                                        /* Faut-il editer les extrema des des conditions initiales, ('VRAI') ou pas ('FAUX') ?       */

#define   NOMBRE_D_ITERATIONS                                                                                                           \
                    UN                                                                                                                  \
                                        /* Nombre d'iterations.                                                                      */
#define   ALLER_PLUTOT_D_ALTERNANCES_EN_ALTERNANCES                                                                                     \
                    FAUX                                                                                                                \
                                        /* Effectuer un nombre d'iterations donne a l'avance ('FAUX') ou bien aller d'alternance     */ \
                                        /* ("erosion-affaiblissement") en alternance ('VRAI') ?                                      */

#define   DECREMENTER_LE_MARQUEUR_DES_POINTS_ERODES_APRES_L_EROSION                                                                     \
                    VRAI                                                                                                                \
                                        /* Pour decrementer 'Aerode__CoastLines_3D_____marqueur_des_points_erodes' apres chaque      */ \
                                        /* iteration ('VRAI') ou avant ('FAUX'). En decrementant avant, on garantit ainsi que le     */ \
                                        /* niveau de marquage de la mer est strictement superieur a celui de la cote. Ceci a ete     */ \
                                        /* introduit le 20040909133231...                                                            */
#define   DECREMENT_DU_MARQUEUR_DES_POINTS_ERODES                                                                                       \
                    FZERO                                                                                                               \
                                        /* Decrement de 'Aerode__CoastLines_3D_____marqueur_des_points_erodes' avant ou apres        */ \
                                        /* chaque iteration...                                                                       */

#define   TAUX_D_AFFAIBLISSEMENT_DE_LA_DURETE_DE_LA_TERRE                                                                               \
                    GRO5(FRA10(FRA10(FRA10(FU))))                                                                                       \
                                        /* Taux d'affaiblissement de la durete de la terre.                                          */
#define   FACTEUR_D_ANTI_CORROSITE_DE_LA_MER                                                                                            \
                    GRO5(FRA10(FU))                                                                                                     \
                                        /* Facteur d'anti-corrosite de la mer.                                                       */

#define   EDITER_LES_ALTERNANCES_D_EROSION_ET_D_AFFAIBLISSEMENT_DE_LA_TERRE                                                             \
                    FAUX                                                                                                                \
                                        /* Faut-il editer les alternances d'erosion et d'affaiblissement de la terre ('VRAI') ou     */ \
                                        /* pas ('FAUX') ?                                                                            */
#define   ALTERNER_L_EROSION_ET_L_AFFAIBLISSEMENT_DE_LA_TERRE                                                                           \
                    FAUX                                                                                                                \
                                        /* Faut-il alterner l'erosion et l'affaiblissement de la terre ('VRAI') ou bien faire les    */ \
                                        /* deux en meme temps ('FAUX') ?                                                             */
#define   FAIRE_EROSION_ET_AFFAIBLISSEMENT_DE_LA_TERRE_EN_CAS_D_ALTERNANCE                                                              \
                    FAUX                                                                                                                \
                                        /* Si 'IL_FAUT(alterner_l_erosion_et_l_affaiblissement_de_la_terre)' peut-on, au moment de   */ \
                                        /* l'alternance, faire simultanement de l'erosion et de l'affaiblissement de la terre        */ \
                                        /* ('VRAI') ou pas ('FAUX') ?                                                                */

#define   EDITER_LA_SURFACE_DE_LA_COTE_COURANTE                                                                                         \
                    FAUX                                                                                                                \
                                        /* Faut-il editer la surface de la cote courante ('VRAI') ou pas ('FAUX') ?                  */

#define   CONVERTIR_EN_IMAGES_STANDARDS                                                                                                 \
                    VRAI                                                                                                                \
                                        /* Faut-il convertir les images Resultats en images standards ('VRAI') ou les laisser        */ \
                                        /* sous forme d'images flottantes ('FAUX') ?                                                 */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M A C R O S   U T I L E S  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   NOMBRE_DE_CHIFFRES_NECESSAIRES(dimension)                                                                                     \
                    ADD2(NOMBRE_DE_CHIFFRES_DECIMAUX(dimension),UN)                                                                     \
                                        /* Choix du nombre de chiffres necessaires pour editer une dimension. On notera qu'il y      */ \
                                        /* a deux 'TRPU(...)', le deuxieme etant destine au signe qui compte pour un caractere.      */ \
                                        /*                                                                                           */ \
                                        /* La procedure 'NOMBRE_DE_CHIFFRES_DECIMAUX(...)' fut introduite le 20051210183809. Les     */ \
                                        /* deux 'TRPU(...)' furent donc remplaces par 'un 'ADD2(...,UN)'...                          */

#define   DECREMENTATION_DU_MARQUEUR_DES_POINTS_ERODES(decrementer)                                                                     \
                    Bblock                                                                                                              \
                    Test(IL_FAUT(decrementer))                                                                                          \
                         Bblock                                                                                                         \
                         Test(IL_NE_FAUT_PAS(Aerode__CoastLines_3D_____marquer_les_points_erodes_avec_la_force))                        \
                              Bblock                                                                                                    \
                              DECR(Aerode__CoastLines_3D_____marqueur_des_points_erodes                                                 \
                                  ,decrement_du_marqueur_des_points_erodes                                                              \
                                   );                                                                                                   \
                                        /* Gestion du marqueur des points erodes lorsqu'il est utile...                              */ \
                              Eblock                                                                                                    \
                         ATes                                                                                                           \
                              Bblock                                                                                                    \
                              Eblock                                                                                                    \
                         ETes                                                                                                           \
                         Eblock                                                                                                         \
                    ATes                                                                                                                \
                         Bblock                                                                                                         \
                         Eblock                                                                                                         \
                    ETes                                                                                                                \
                    Eblock                                                                                                              \
                                        /* Decrementation du marqueur des points erodes.                                             */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M O D E L E   T R I D I M E N S I O N N E L   D E   C O T E S   A V E C   E R O S I O N  :                                 */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
     DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE));
     DEFV(CHAR,INIC(POINTERc(nom_imageN),NOM_PIPE));
     DEFV(CHAR,INIC(POINTERc(nom_imageF),NOM_PIPE));

     DEFV(CHAR,INIC(POINTERc(nom_postfixe),NOM_UNDEF_VIDE));
                                        /* Nom d'un eventuel postfixe a placer derriere <nom_imageA><numero> (par exemple '$ROUGE'). */
     DEFV(Int,INIT(premiere_image,PREMIERE_IMAGE));
     DEFV(Int,INIT(premiere_imageN,PREMIERE_IMAGE));
                                        /* Numero de la premiere image,                                                              */
     DEFV(Int,INIT(pas_des_images,PAS_DES_IMAGES));
     DEFV(Int,INIT(pas_des_imagesN,PAS_DES_IMAGES));
                                        /* Pas de passage d'un numero d'image a une autre.                                           */
     DEFV(Int,INIT(nombre_de_chiffres,NOMBRE_DE_CHIFFRES));
     DEFV(Int,INIT(nombre_de_chiffresN,NOMBRE_DE_CHIFFRES));
                                        /* Nombre de chiffres codant le numero des images de la serie...                             */

     DEFV(Logical,INIT(editer_les_parametres_du_noyau,EDITER_LES_PARAMETRES_DU_NOYAU));
                                        /* Faut-il editer les parametres du noyau et sortir immediatement ('VRAI') ou exploiter      */
                                        /* le modele d'erosion de cotes tridimensionnel ('FAUX') ?                                   */
     DEFV(Logical,INIT(lister_le_contenu_du_noyau,LISTER_LE_CONTENU_DU_NOYAU));
                                        /* Faut-il lister le contenu du noyau ('VRAI') ou pas ('FAUX') ?                             */

     DEFV(Logical,INIT(imposer_les_extrema_des_conditions_initiales,IMPOSER_LES_EXTREMA_DES_CONDITIONS_INITIALES));
     DEFV(genere_Float,INIT(minimum_impose_pour_les_conditions_initiales,MINIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES));
     DEFV(genere_Float,INIT(maximum_impose_pour_les_conditions_initiales,MAXIMUM_IMPOSE_POUR_LES_CONDITIONS_INITIALES));
                                        /* Faut-il imposer ('VRAI') ou pas ('FAUX') les extrema des conditions initiales, et si      */
                                        /* 'IL_FAUT(imposer_les_extrema_des_conditions_initiales)', valeurs implicites ?             */
     DEFV(Logical,INIT(editer_les_extrema_des_conditions_initiales,EDITER_LES_EXTREMA_DES_CONDITIONS_INITIALES));
                                        /* Faut-il editer les extrema des des conditions initiales, ('VRAI') ou pas ('FAUX') ?       */

     DEFV(Positive,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
                                        /* Nombre d'iterations.                                                                      */
     DEFV(Logical,INIT(aller_plutot_d_alternances_en_alternances,ALLER_PLUTOT_D_ALTERNANCES_EN_ALTERNANCES));
                                        /* Effectuer un nombre d'iterations donne a l'avance ('FAUX') ou bien aller d'alternance     */
                                        /* ("erosion-affaiblissement") en alternance ('VRAI') ?                                      */

     DEFV(Logical,INIT(decrementer_le_marqueur_des_points_erodes_apres_l_erosion
                      ,DECREMENTER_LE_MARQUEUR_DES_POINTS_ERODES_APRES_L_EROSION
                       )
          );
                                        /* Pour decrementer 'Aerode__CoastLines_3D_____marqueur_des_points_erodes' apres chaque      */
                                        /* iteration ('VRAI') ou avant ('FAUX'). En decrementant avant, on garantit ainsi que le     */
                                        /* niveau de marquage de la mer est strictement superieur a celui de la cote. Ceci a ete     */
                                        /* introduit le 20040909133231...                                                            */
     DEFV(Float,INIT(decrement_du_marqueur_des_points_erodes,DECREMENT_DU_MARQUEUR_DES_POINTS_ERODES));
                                        /* Decrement de 'Aerode__CoastLines_3D_____marqueur_des_points_erodes' avant ou apres        */
                                        /* chaque iteration...                                                                       */

     DEFV(Float,INIT(taux_d_affaiblissement_de_la_durete_de_la_terre,TAUX_D_AFFAIBLISSEMENT_DE_LA_DURETE_DE_LA_TERRE));
                                        /* Taux d'affaiblissement de la durete de la terre.                                          */
     DEFV(Float,INIT(facteur_d_anti_corrosite_de_la_mer,FACTEUR_D_ANTI_CORROSITE_DE_LA_MER));
                                        /* Facteur d'anti-corrosite de la mer. L'experience montre qu'une valeur faible (par exemple */
                                        /* 0.1) pour l'anti-corrosite de la mer favorise l'existence d'iles et de cotes d'apparence  */
                                        /* plus tourmentees, alors qu'une valeur elevee (par exemple 0.9) favorise l'absence d'iles  */
                                        /* et produit des cotes plus "simples"...                                                    */

     DEFV(Logical,INIT(editer_les_alternances_d_erosion_et_d_affaiblissement_de_la_terre
                      ,EDITER_LES_ALTERNANCES_D_EROSION_ET_D_AFFAIBLISSEMENT_DE_LA_TERRE
                       )
          );
                                        /* Faut-il editer les alternances d'erosion et d'affaiblissement de la terre ('VRAI') ou     */
                                        /* pas ('FAUX') ?                                                                            */
     DEFV(Logical,INIT(alterner_l_erosion_et_l_affaiblissement_de_la_terre
                      ,ALTERNER_L_EROSION_ET_L_AFFAIBLISSEMENT_DE_LA_TERRE
                       )
          );
                                        /* Faut-il alterner l'erosion et l'affaiblissement de la terre ('VRAI') ou bien faire les    */
                                        /* deux en meme temps ('FAUX') ?                                                             */
     DEFV(Logical,INIT(faire_erosion_et_affaiblissement_de_la_terre_en_cas_d_alternance
                      ,FAIRE_EROSION_ET_AFFAIBLISSEMENT_DE_LA_TERRE_EN_CAS_D_ALTERNANCE
                       )
          );
                                        /* Si 'IL_FAUT(alterner_l_erosion_et_l_affaiblissement_de_la_terre)' peut-on, au moment de   */
                                        /* l'alternance, faire simultanement de l'erosion et de l'affaiblissement de la terre        */
                                        /* ('VRAI') ou pas ('FAUX') ?                                                                */

     DEFV(Logical,INIT(editer_la_surface_de_la_cote_courante,EDITER_LA_SURFACE_DE_LA_COTE_COURANTE));
                                        /* Faut-il editer la surface de la cote courante ('VRAI') ou pas ('FAUX') ?                  */

     DEFV(Logical,INIT(convertir_en_images_standards,CONVERTIR_EN_IMAGES_STANDARDS));
                                        /* Faut-il convertir les images Resultats en images standards ('VRAI') ou les laisser        */
                                        /* sous forme d'images flottantes ('FAUX') ?                                                 */
     /*..............................................................................................................................*/
     GET_ARGUMENTSi(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_C("imageA=""A=",nom_imageA);
                         GET_ARGUMENT_C("imageN=""N=",nom_imageN);
                         GET_ARGUMENT_C("imageR=""R=",nom_imageR);
                         GET_ARGUMENT_C("imageF=""F=",nom_imageF);
                         GET_ARGUMENT_C("postfixe=",nom_postfixe);
                         GET_ARGUMENT_I("premiere=",premiere_image);
                         GET_ARGUMENT_I("premiereN=",premiere_imageN);
                         GET_ARGUMENT_I("pas=",pas_des_images);
                         GET_ARGUMENT_I("pasN=",pas_des_imagesN);
                         GET_ARGUMENT_I("chiffres=",nombre_de_chiffres);
                         GET_ARGUMENT_I("chiffresN=",nombre_de_chiffresN);

                         GET_ARGUMENT_L("periodiser_X=""periodX=",Aerode__CoastLines_3D_____periodiser_X);
                         GET_ARGUMENT_L("periodiser_Y=""periodY=",Aerode__CoastLines_3D_____periodiser_Y);
                         GET_ARGUMENT_L("periodiser_Z=""periodZ=",Aerode__CoastLines_3D_____periodiser_Z);
                         GET_ARGUMENT_L("symetriser_X=""symetrX=",Aerode__CoastLines_3D_____symetriser_X);
                         GET_ARGUMENT_L("symetriser_Y=""symetrY=",Aerode__CoastLines_3D_____symetriser_Y);
                         GET_ARGUMENT_L("symetriser_Z=""symetrZ=",Aerode__CoastLines_3D_____symetriser_Z);
                         GET_ARGUMENT_L("prolonger_X=""prolX=",Aerode__CoastLines_3D_____prolonger_X);
                         GET_ARGUMENT_L("prolonger_Y=""prolY=",Aerode__CoastLines_3D_____prolonger_Y);
                         GET_ARGUMENT_L("prolonger_Z=""prolZ=",Aerode__CoastLines_3D_____prolonger_Z);
                         GET_ARGUMENT_F("niveau_hors_album=""niveau=",Aerode__CoastLines_3D_____niveau_hors_album);

                         GET_ARGUMENT_L("imposer_extrema=""extrema=",imposer_les_extrema_des_conditions_initiales);
                         GET_ARGUMENT_F("minimum=",minimum_impose_pour_les_conditions_initiales);
                         GET_ARGUMENT_F("maximum=",maximum_impose_pour_les_conditions_initiales);
                         GET_ARGUMENT_L("editer_extrema=",editer_les_extrema_des_conditions_initiales);

                         GET_ARGUMENT_I("ni=""iterations=",nombre_d_iterations);
                                        /* Le 20050623141247, "n=" et "N=" ont ete remplaces par "ni=" (double definition...).       */
                         GET_ARGUMENT_L("iterations_variables=",aller_plutot_d_alternances_en_alternances);

                         GET_ARGUMENT_L("noyau=",Aerode__CoastLines_3D_____utiliser_un_noyau);
                         GET_ARGUMENT_I("demi_dimension=""dd=",Aerode__CoastLines_3D_____demi_dimension_effective_du_noyau);
                         GET_ARGUMENT_L("editer_noyau=",editer_les_parametres_du_noyau);
                         GET_ARGUMENT_L("lister_noyau=",lister_le_contenu_du_noyau);

                         GET_ARGUMENT_F("sX__Y__Z=""sXYZ=",Aerode__CoastLines_3D_____ponderation_du_point_sX__Y__Z);
                         GET_ARGUMENT_F("_X_sY__Z=""XsYZ=",Aerode__CoastLines_3D_____ponderation_du_point__X_sY__Z);
                         GET_ARGUMENT_F("_X__Y_sZ=""XYsZ=",Aerode__CoastLines_3D_____ponderation_du_point__X__Y_sZ);
                         GET_ARGUMENT_F("pX__Y__Z=""pXYZ=",Aerode__CoastLines_3D_____ponderation_du_point_pX__Y__Z);
                         GET_ARGUMENT_F("_X_pY__Z=""XpYZ=",Aerode__CoastLines_3D_____ponderation_du_point__X_pY__Z);
                         GET_ARGUMENT_F("_X__Y_pZ=""XYpZ=",Aerode__CoastLines_3D_____ponderation_du_point__X__Y_pZ);

                         GET_ARGUMENT_I("dx=""dX=",Aerode__CoastLines_3D_____delta_X);
                         GET_ARGUMENT_I("dy=""dY=",Aerode__CoastLines_3D_____delta_Y);
                         GET_ARGUMENT_I("dz=""dZ=",Aerode__CoastLines_3D_____delta_Z);

                         GET_ARGUMENT_L("global=",Aerode__CoastLines_3D_____utiliser_la_surface_globale);
                         GET_ARGUMENT_I("ddx=""ddX="
                                       ,Aerode__CoastLines_3D_____demi_dimension_X_du_domaine_de_calcul_de_la_surface_locale
                                        );
                         GET_ARGUMENT_I("ddy=""ddY="
                                       ,Aerode__CoastLines_3D_____demi_dimension_Y_du_domaine_de_calcul_de_la_surface_locale
                                        );
                         GET_ARGUMENT_I("ddz=""ddZ="
                                       ,Aerode__CoastLines_3D_____demi_dimension_Z_du_domaine_de_calcul_de_la_surface_locale
                                        );
                         GET_ARGUMENT_F("fad=""ad="
                                       ,Aerode__CoastLines_3D_____agrandissement_du_domaine_de_calcul_de_la_surface_locale
                                        );

                         GET_ARGUMENT_F("amplificateur=",Aerode__CoastLines_3D_____amplificateur_de_la_force);
                         GET_ARGUMENT_L("force_locale=",Aerode__CoastLines_3D_____la_force_depend_du_niveau_local_de_la_mer);

                         GET_ARGUMENT_L("marquer_force=",Aerode__CoastLines_3D_____marquer_les_points_erodes_avec_la_force);
                         GET_ARGUMENT_F("points_erodes=""mpe=",Aerode__CoastLines_3D_____marqueur_des_points_erodes);
                         GET_ARGUMENT_L("decrementer_apres=""dap=",decrementer_le_marqueur_des_points_erodes_apres_l_erosion);
                         GET_ARGUMENT_N("decrementer_avant=""dav=",decrementer_le_marqueur_des_points_erodes_apres_l_erosion);
                         GET_ARGUMENT_F("dmpe=""decrement_points_erodes=",decrement_du_marqueur_des_points_erodes);

                         GET_ARGUMENT_F("terre=""durete=""affaiblissement=",taux_d_affaiblissement_de_la_durete_de_la_terre);
                         GET_ARGUMENT_F("mer=""anti_corrosite=",facteur_d_anti_corrosite_de_la_mer);

                         GET_ARGUMENT_L("editer_alternances=",editer_les_alternances_d_erosion_et_d_affaiblissement_de_la_terre);
                         GET_ARGUMENT_L("alterner=",alterner_l_erosion_et_l_affaiblissement_de_la_terre);
                         GET_ARGUMENT_L("simultane=",faire_erosion_et_affaiblissement_de_la_terre_en_cas_d_alternance);

                         GET_ARGUMENT_I("voisinsA="
                                       ,Aerode__CoastLines_3D_____nombre_de_voisins_mer_d_un_point_terre_amorcant_l_erosion
                                        );
                         GET_ARGUMENT_I("voisinsF="
                                       ,Aerode__CoastLines_3D_____nombre_minimal_de_voisins_mer_d_un_point_terre_forcant_l_erosion
                                        );
                         GET_ARGUMENT_F("Fvoisins=",Aerode__CoastLines_3D_____facteur_du_nombre_d_etats_de_type___mer);
                         GET_ARGUMENT_F("Tvoisins=",Aerode__CoastLines_3D_____translateur_du_nombre_d_etats_de_type___mer);
                         GET_ARGUMENT_F("epaisseur=",Aerode__CoastLines_3D_____facteur_d_epaisseur_de_la_cote);
                         GET_ARGUMENT_I("points=""convolution=",Aerode__CoastLines_3D_____nombre_de_points_du_noyau_de_convolution);

                         GET_ARGUMENT_L("editer_surface=",editer_la_surface_de_la_cote_courante);

                         GET_ARGUMENT_L("convertir=",convertir_en_images_standards);
                         )
                    );

     Test(IL_FAUT(alterner_l_erosion_et_l_affaiblissement_de_la_terre))
          Bblock
          Eblock
     ATes
          Bblock
          Test(IL_FAUT(aller_plutot_d_alternances_en_alternances))
               Bblock
               PRINT_ATTENTION("la progression d'alternance en alternance ne peut avoir lieu que s'il y a bien alternance");
               PRINT_ATTENTION("la progression d'alternance en alternance est donc inhibee");

               EGAL(aller_plutot_d_alternances_en_alternances,FAUX);
                                        /* Et oui, car sinon, cela bouclerait...                                                     */
               Eblock
          ATes
               Bblock
               Eblock
          ETes
          Eblock
     ETes

     Test(IL_FAUT(editer_les_parametres_du_noyau))
          Bblock
          CAL3(Prme6("Xmin=%d  Xmax=%d  Ymin=%d  Ymax=%d  Zmin=%d  Zmax=%d\n"
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                     )
               );
          CAL3(Prme4("XminR=%d XmaxR=%d YminR=%d YmaxR=%d\n"
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                     )
               );
          CALS(Fsauts_de_lignes(UN));

          CAL3(Prme4("XYmaxNe   %d %d %d %d\n"
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                    ,PINTE(Aerode__CoastLines_3D_____XYZmin),PINTE(Aerode__CoastLines_3D_____XYZmax)
                     )
               );
                                        /* ATTENTION : on ne devra pas etre surpris par les valeurs ainsi obtenues. Par exemple, on  */
                                        /* obtiendra :                                                                               */
                                        /*                                                                                           */
                                        /*                  XYmax     128 128                                                        */
                                        /*                                                                                           */
                                        /* avec la valeur 64 pour 'DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Aerode__CoastLines_3D', en  */
                                        /* se rappelant que l'on passe a 'DimNo_Aerode__CoastLines_3D' par un 'DOUP(...)'. Les       */
                                        /* extrema ainsi obtenus (128) different d'une unite de ceux fournis par 'Suq' (127) et      */
                                        /* les parites des deux dimensions respectives (129 et 128) sont elles-aussi differentes...  */
                                        /*                                                                                           */
                                        /* ATTENTION : passage de 'XYmax' a 'XYmaxNe' le 20020909112853.                             */
                                        /* ATTENTION : passage de 'XYmax' a 'XYmaxNe' le 20020909112152.                             */
                                        /*                                                                                           */
                                        /* Pour l'usage de 'PINTE(...)' voir 'v $xil/defi_K2$vv$DEF 20230311135805' qui a ete        */
                                        /* introduit pour ces impressions le 20230311144143...                                       */
          Eblock
     ATes
          Bblock
          BDEFV(albumF,Falbum_d_imagesA);
          BDEFV(albumF,Falbum_d_imagesR);
          BDEFV(album,album_d_imagesR);
                                        /* Definition des albums d'images Argument et Resultat...                                    */

          Test(IFNE_chaine(nom_imageN,NOM_PIPE))
               Bblock
               PUSH_DIMENSIONS_3D;
               SET_DIMENSIONS_3D_SANS_VALIDATION(Aerode__CoastLines_3D_____XYZmin,Aerode__CoastLines_3D_____XYZmax
                                                ,Aerode__CoastLines_3D_____XYZmin,Aerode__CoastLines_3D_____XYZmax
                                                ,Aerode__CoastLines_3D_____XYZmin,Aerode__CoastLines_3D_____XYZmax
                                                 );
                                        /* Le noyau va etre recupere comme une "petite" image dont les dimensions sont definies      */
                                        /* par {Aerode__CoastLines_3D_____XYZmin,Aerode__CoastLines_3D_____XYZmax} et ce             */
                                        /* independemment de 'Aerode__CoastLines_3D_____demi_dimension_effective_du_noyau'...        */

               CALi(dAloadF_album(PAerode__CoastLines_3D_____noyau
                                 ,nom_imageN,C_VIDE
                                 ,premiere_imageN,pas_des_imagesN,nombre_de_chiffresN
                                 ,FAUX
                                 ,UNDEF
                                 ,FAUX
                                  )
                    );
                                        /* Chargement de l'album du noyau...                                                         */
                                        /*                                                                                           */
                                        /* ATTENTION : on notera que l'album 'nom_imageN' doit avoir comme taille la taille          */
                                        /* maximale du noyau correspondant a 'v $xiii/di_album$DEF DimNo_Aerode__CoastLines_3D',     */
                                        /* et non pas celle qui correspond a la taille effective de ce meme noyau                    */
                                        /* 'v $xiii/di_album$FON Aerode__CoastLines_3D_____demi_dimension_effective_du_noyau'.       */
                                        /*                                                                                           */
                                        /* ATTENTION : on notera de plus que l'on passe par 'PAerode__CoastLines_3D_____noyau' car,  */
                                        /* en effet, 'Aerode__CoastLines_3D_____noyau' n'est pas un vrai 'albumF' puisqu'il n'est    */
                                        /* pas cree par une allocation dynamique et variable...                                      */
                                        /*                                                                                           */
                                        /* L'argument 'Inv' (='FAUX') a ete introduit le 20091125133113.                             */

               Test(PAS_D_ERREUR(CODE_ERREUR))
                    Bblock
                    EGAL(Aerode__CoastLines_3D_____initialiser_le_noyau,FAUX);
                                        /* Evidemment, l'initialisation est faite...                                                 */

                    Test(IL_FAUT(lister_le_contenu_du_noyau))
                         Bblock
                         begin_albumQ(DoIn,Aerode__CoastLines_3D_____XYZmin_effectif,Aerode__CoastLines_3D_____XYZmax_effectif,PasZ
                                     ,DoIn,Aerode__CoastLines_3D_____XYZmin_effectif,Aerode__CoastLines_3D_____XYZmax_effectif,PasY
                                     ,DoIn,Aerode__CoastLines_3D_____XYZmin_effectif,Aerode__CoastLines_3D_____XYZmax_effectif,PasX
                                      )
                              Bblock
                              CAL3(Prme7("noyau(%+0*d,%+0*d,%+0*d)=%+.^^^\n"
                                        ,NOMBRE_DE_CHIFFRES_NECESSAIRES(DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Aerode__CoastLines_3D),X
                                        ,NOMBRE_DE_CHIFFRES_NECESSAIRES(DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Aerode__CoastLines_3D),Y
                                        ,NOMBRE_DE_CHIFFRES_NECESSAIRES(DEMI_DIMENSION_MAXIMALE_DU_NOYAU_DANS_Aerode__CoastLines_3D),Z
                                        ,ACCES_NOYAU_DANS_Aerode__CoastLines_3D(X,Y,Z)
                                         )
                                   );
                                        /* Le 20060105153819, le format "16g" est passe a "^^g" pour plus de souplesse...            */
                                        /*                                                                                           */
                                        /* Le 20091123123111, le format "^^g" est passe a "^^^" pour plus de souplesse...            */
                              Eblock
                         end_albumQ(EDoI,EDoI,EDoI)
                         Eblock
                    ATes
                         Bblock
                         Eblock
                    ETes
                    Eblock
               ATes
                    Bblock
                    Test__CODE_ERREUR__ERREUR07;

                    INITIALISATION_EVENTUELLE_DU_NOYAU_DANS_Aerode__CoastLines_3D;
                                        /* Initialisation du noyau en cas d'erreur...                                                */
                    Eblock
               ETes

               PULL_DIMENSIONS_3D;
                                        /* Restauration des dimensions de travail...                                                 */
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          CALi(dAFinitialisation(Falbum_d_imagesR,FZERO));
                                        /* Initialisation de l'image Resultat flottante (dans tous les cas car elle est calculee     */
                                        /* meme si 'IL_FAUT(convertir_en_images_standards)').                                        */

          Test(IL_FAUT(convertir_en_images_standards))
               Bblock
               CALi(dAnoir(album_d_imagesR));
                                        /* Initialisation de l'album Resultat...                                                     */
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          CALi(dAloadF_album(Falbum_d_imagesA
                            ,nom_imageA
                            ,nom_postfixe
                            ,premiere_image
                            ,pas_des_images
                            ,nombre_de_chiffres
                            ,FAUX
                            ,UNDEF
                            ,FAUX
                             )
               );
                                        /* Chargement de l'album Argument...                                                         */
                                        /*                                                                                           */
                                        /* L'argument 'Inv' (='FAUX') a ete introduit le 20091125133113.                             */

          Test(PAS_D_ERREUR(CODE_ERREUR))
               Bblock
               DEFV(Logical,INIT(c_est_la_premiere_iteration,VRAI));
                                        /* Afin de pouvoir faire la difference entre la premiere iteration et les suivantes...       */

               DEFV(genere_Float,INIT(niveau_minimum_de_Falbum_d_imagesA,minimum_impose_pour_les_conditions_initiales));
               DEFV(genere_Float,INIT(niveau_maximum_de_Falbum_d_imagesA,maximum_impose_pour_les_conditions_initiales));
                                        /* Extrema du champ initial ('Falbum_d_imagesA').                                            */

               DEFV(Logical,INIT(iterer_l_erosion,VRAI));
                                        /* Le 20020213104540 l'iteration de l'erosion est devenue un 'Tant(...)' a la place d'un     */
                                        /* 'Repe(...)' pour permettre un nombre variable d'iterations...                             */
               DEFV(Int,INIT(compteur_des_iterations_d_erosion,PREMIERE_ITERATION_D_UN_Repe));
                                        /* Afin de simuler un 'Repe(...)' a nombre variable d'iterations...                          */

               Test(IL_NE_FAUT_PAS(imposer_les_extrema_des_conditions_initiales))
                    Bblock
                    CALS(AFnivo_extrema(Falbum_d_imagesA
                                       ,ADRESSE(niveau_minimum_de_Falbum_d_imagesA)
                                       ,ADRESSE(niveau_maximum_de_Falbum_d_imagesA)
                                        )
                         );
                                        /* Recherche des extrema de 'Falbum_d_imagesA'.                                              */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               Test(IL_FAUT(editer_les_extrema_des_conditions_initiales))
                    Bblock
                    CAL3(Prme1("minimum des conditions initiales=%+.^^^\n",niveau_minimum_de_Falbum_d_imagesA));
                    CAL3(Prme1("maximum des conditions initiales=%+.^^^\n",niveau_maximum_de_Falbum_d_imagesA));
                                        /* Edition des extrema de 'Falbum_d_imagesA'.                                                */
                                        /*                                                                                           */
                                        /* Le 20060105153819, le format "16g" est passe a "^^g" pour plus de souplesse...            */
                                        /*                                                                                           */
                                        /* Le 20091123123111, le format "^^g" est passe a "^^^" pour plus de souplesse...            */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               EGAL(Aerode__CoastLines_3D_____niveau_de_la___mer,niveau_maximum_de_Falbum_d_imagesA);
               EGAL(Aerode__CoastLines_3D_____niveau_de_la_terre,niveau_minimum_de_Falbum_d_imagesA);
                                        /* Mise a jour du modele...                                                                  */

               Test(IL_NE_FAUT_PAS(Aerode__CoastLines_3D_____marquer_les_points_erodes_avec_la_force))
                    Bblock
                    Test(IFOU(SONT_DE_MEME_SIGNE(NEUT(Aerode__CoastLines_3D_____marqueur_des_points_erodes)
                                                ,Aerode__CoastLines_3D_____niveau_de_la_terre
                                                 )
                             ,SONT_DE_MEME_SIGNE(SOUS(Aerode__CoastLines_3D_____marqueur_des_points_erodes
                                                     ,MUL2(PRED(nombre_d_iterations)
                                                          ,decrement_du_marqueur_des_points_erodes
                                                           )
                                                      )
                                                ,Aerode__CoastLines_3D_____niveau_de_la_terre
                                                 )
                              )
                         )
                         Bblock
                         PRINT_ATTENTION("il va y avoir confusion entre le marqueur des points erodes et le niveau de la terre");

                         CAL3(Prme1("niveau de la terre................=%+.^^^\n",Aerode__CoastLines_3D_____niveau_de_la_terre));
                                        /* Le 20060105153819, le format "16g" est passe a "^^g" pour plus de souplesse...            */
                                        /*                                                                                           */
                                        /* Le 20091123123111, le format "^^g" est passe a "^^^" pour plus de souplesse...            */
                         CAL3(Prme1("nombre total d'iterations.........=%d\n"
                                   ,nombre_d_iterations
                                    )
                              );
                         CAL3(Prme1("marqueur initial des points erodes=%+.^^^\n"
                                   ,NEUT(Aerode__CoastLines_3D_____marqueur_des_points_erodes)
                                    )
                              );
                                        /* Le 20060105153819, le format "16g" est passe a "^^g" pour plus de souplesse...            */
                                        /*                                                                                           */
                                        /* Le 20091123123111, le format "^^g" est passe a "^^^" pour plus de souplesse...            */
                         CAL3(Prme1("marqueur final des points erodes..=%+.^^^\n"
                                   ,SOUS(Aerode__CoastLines_3D_____marqueur_des_points_erodes
                                        ,MUL2(PRED(nombre_d_iterations)
                                             ,decrement_du_marqueur_des_points_erodes
                                              )
                                         )
                                    )
                              );
                                        /* Le 20060105153819, le format "16g" est passe a "^^g" pour plus de souplesse...            */
                                        /*                                                                                           */
                                        /* Le 20091123123111, le format "^^g" est passe a "^^^" pour plus de souplesse...            */
                         Eblock
                    ATes
                         Bblock
                         Eblock
                    ETes
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               Tant(IL_FAUT(iterer_l_erosion))
                    Bblock
                    Test(IFOU(EST_VRAI(c_est_la_premiere_iteration)
                             ,IFET(EST_FAUX(c_est_la_premiere_iteration)
                                  ,IZGT(Aerode__CoastLines_3D_____surface_de_la_cote_courante)
                                   )
                              )
                         )
                         Bblock
                         DEFV(Float,INIT(taux_d_affaiblissement_de_la_durete_de_la_terre_effectif
                                        ,taux_d_affaiblissement_de_la_durete_de_la_terre
                                         )
                              );
                         DEFV(Float,INIT(facteur_d_anti_corrosite_de_la_mer_effectif
                                        ,facteur_d_anti_corrosite_de_la_mer
                                         )
                              );
                                        /* Parametres respectifs de la "terre" et de la "mer" effectivement utilises...              */

                         Test(IL_FAUT(alterner_l_erosion_et_l_affaiblissement_de_la_terre))
                              Bblock
                              Test(EST_VRAI(c_est_la_premiere_iteration))
                                   Bblock
                                   EGAL(taux_d_affaiblissement_de_la_durete_de_la_terre_effectif,FZERO);
                                        /* Dans le cas de la premiere iteration, il est evidemment impossible de tester              */
                                        /* 'Aerode__CoastLines_3D_____nombre_de_points_erodes_a_l_iteration_courante'.               */
                                        /* La terre ne doit alors pas etre affaiblie...                                              */
                                   Eblock
                              ATes
                                   Bblock
                                   Test(IZEQ(Aerode__CoastLines_3D_____nombre_de_points_erodes_a_l_iteration_courante))
                                        Bblock
                                        /* Cas ou la mer n'a erode aucun point de la terre : la terre va etre affaiblie, avec ou     */
                                        /* pas erosion par la mer simultanement...                                                   */
                                        Test(IL_FAUT(editer_les_alternances_d_erosion_et_d_affaiblissement_de_la_terre))
                                             Bblock
                                             CAL3(Prme1("iteration %d : "
                                                       ,compteur_des_iterations_d_erosion
                                                        )
                                                  );
                                             CALS(FPrme0("alternance 'erosion-affaiblissement'\n"));
                                             Eblock
                                        ATes
                                             Bblock
                                             Eblock
                                        ETes

                                        Test(IL_FAUT(faire_erosion_et_affaiblissement_de_la_terre_en_cas_d_alternance))
                                             Bblock
                                        /* Erosion simultanee par la mer...                                                          */
                                             Eblock
                                        ATes
                                             Bblock
                                             EGAL(facteur_d_anti_corrosite_de_la_mer_effectif,FZERO);
                                        /* Pas d'erosion simultanee par la mer...                                                    */
                                             Eblock
                                        ETes

                                        Test(IL_FAUT(aller_plutot_d_alternances_en_alternances))
                                             Bblock
                                             EGAL(iterer_l_erosion,FAUX);
                                        /* Le nombre d'iterations est variable et limite par l'alternance "erosion-affaiblissement". */
                                             Eblock
                                        ATes
                                             Bblock
                                             Eblock
                                        ETes
                                        Eblock
                                   ATes
                                        Bblock
                                        /* Cas ou la mer a erode au moins un point de la terre :                                     */
                                        EGAL(taux_d_affaiblissement_de_la_durete_de_la_terre_effectif,FZERO);
                                        /* La terre ne subit alors pas d'affaiblissement, mais uniquement de l'erosion...            */
                                        Eblock
                                   ETes
                                   Eblock
                              ETes
                              Eblock
                         ATes
                              Bblock
                              Eblock
                         ETes

                         Test(IL_FAUT(editer_la_surface_de_la_cote_courante))
                              Bblock
                              CAL3(Prme1("iteration %d : "
                                        ,compteur_des_iterations_d_erosion
                                         )
                                   );
                              Eblock
                         ATes
                              Bblock
                              Eblock
                         ETes

#define   decrementer_apres_l_erosion                                                                                                   \
                    decrementer_le_marqueur_des_points_erodes_apres_l_erosion                                                           \
                                        /* Pour raccourcir certaines lignes qui suivent...                                           */

                         DECREMENTATION_DU_MARQUEUR_DES_POINTS_ERODES(NOTL(decrementer_apres_l_erosion));

                         CALS(Aerode__CoastLines_3D(Falbum_d_imagesR
                                                   ,Falbum_d_imagesA
                                                   ,Aerode__CoastLines_3D_____albumA_convolue_initial
                                                   ,taux_d_affaiblissement_de_la_durete_de_la_terre_effectif
                                                   ,facteur_d_anti_corrosite_de_la_mer_effectif
                                                   ,editer_la_surface_de_la_cote_courante
                                                    )
                              );
                                        /* Et calcul de la generation suivante...                                                    */

                         CALS(AFmove(Falbum_d_imagesA,Falbum_d_imagesR));
                                        /* Preparation de l'iteration suivante sans optimiser par rapport a la derniere ou cette     */
                                        /* operation est inutile.                                                                    */

                         DECREMENTATION_DU_MARQUEUR_DES_POINTS_ERODES(NEUL(decrementer_apres_l_erosion));

#undef    decrementer_apres_l_erosion

                         Eblock
                    ATes
                         Bblock
                         Eblock
                    ETes

                    Test(EST_VRAI(c_est_la_premiere_iteration))
                         Bblock
                         EGAL(c_est_la_premiere_iteration,FAUX);
                                        /* Pour les iterations suivantes...                                                          */
                         Eblock
                    ATes
                         Bblock
                         Eblock
                    ETes

                    INCR(compteur_des_iterations_d_erosion,I);
                                        /* Simulation du 'Repe(...)' d'erosion...                                                    */

                    Test(IFET(IL_NE_FAUT_PAS(aller_plutot_d_alternances_en_alternances)
                             ,IFGT(compteur_des_iterations_d_erosion,nombre_d_iterations)
                              )
                         )
                         Bblock
                         EGAL(iterer_l_erosion,FAUX);
                                        /* On a fait le nombre d'iterations demandees...                                             */
                         Eblock
                    ATes
                         Bblock
                         Eblock
                    ETes
                    Eblock
               ETan

               Test(IL_FAUT(convertir_en_images_standards))
                    Bblock
                    CALS(Afloat_std(album_d_imagesR
                                   ,Falbum_d_imagesR
                                   ,niveau_minimum_de_Falbum_d_imagesA,niveau_maximum_de_Falbum_d_imagesA
                                    )
                         );
                                        /* Renormalisation avec forcage des extrema en utilisant ceux de l'image initiale            */
                                        /* 'Falbum_d_imagesA'.                                                                       */

                    CALi(dAupdate_album(nom_imageR
                                       ,nom_postfixe
                                       ,premiere_image
                                       ,pas_des_images
                                       ,album_d_imagesR
                                       ,nombre_de_chiffres
                                       ,FAUX
                                        )
                         );
                                        /* Generation de l'album Resultat...                                                         */
                                        /*                                                                                           */
                                        /* L'argument 'Inv' (='FAUX') a ete introduit le 20100624091304...                           */
                    Eblock
               ATes
                    Bblock
                    CALi(dAupdateF_album(nom_imageR
                                        ,nom_postfixe
                                        ,premiere_image
                                        ,pas_des_images
                                        ,Falbum_d_imagesR
                                        ,nombre_de_chiffres
                                        ,FAUX
                                         )
                         );
                                        /* Generation de l'album Resultat...                                                         */
                                        /*                                                                                           */
                                        /* L'argument 'Inv' (='FAUX') a ete introduit le 20100624091304...                           */
                    Eblock
               ETes

               Test(IFET(IFNE_chaine(nom_imageF,NOM_PIPE),IL_FAUT(convertir_en_images_standards)))
                    Bblock
                    CALi(dAupdateF_album(nom_imageF
                                        ,nom_postfixe
                                        ,premiere_image
                                        ,pas_des_images
                                        ,Falbum_d_imagesR
                                        ,nombre_de_chiffres
                                        ,FAUX
                                         )
                         );
                                        /* Dans le cas ou la sortie est du type "image standard", la possibilite est offerte de      */
                                        /* generer les conditions Finales en mode "non standard" afin qu'elles puissent servir       */
                                        /* eventuellement de conditions initiales dans un autre calcul...                            */
                                        /*                                                                                           */
                                        /* L'argument 'Inv' (='FAUX') a ete introduit le 20100624091304...                           */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          EDEFV(album,album_d_imagesR);
          EDEFV(albumF,Falbum_d_imagesR);
          EDEFV(albumF,Falbum_d_imagesA);
                                        /* Definition des albums d'images Argument et Resultat...                                    */
          Eblock
     ETes

     RETU_Commande;
     Eblock
ECommande



Copyright © Jean-François COLONNA, 2019-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2019-2024.