/*************************************************************************************************************************************/
/* */
/* 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