Hi everybody,
I'm trying to use Homard to refine the mesh in the case of a crack propagation but I'm having big trubles trying to putting it in the FOR cycle.
This is the code I'm working on
`DEBUT(LANG='FR')

#----------------------------------------------

Variables

#----------------------------------------------

NbFond = 3
Forza = 3500
Nodi = 32
FlexNoda = -Forza / Nodi
Nume_Fond = 1
nbcalc = 5

#----------------------------------------------

Mesh

#----------------------------------------------

mesh = LIRE_MAILLAGE(FORMAT='MED', UNITE=20)

mesh = MODI_MAILLAGE(reuse=mesh,
MAILLAGE=mesh,
ORIE_PEAU=_F(GROUP_MA_PEAU=('Fix_surface', )))

mesh2 = [None] * nbcalc
mesh2[0] = CREA_MAILLAGE(CREA_POI1=_F(GROUP_NO=('Load', ),
NOM_GROUP_MA=('assegnazione', )),
MAILLAGE=mesh)

#----------------------------------------------

Sane Model

#----------------------------------------------

Sane = AFFE_MODELE(AFFE=(_F(GROUP_MA=('FlexLine', 'Fix_surface', 'FaceRot', 'tutto'),
MODELISATION='3D',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('assegnazione', ),
MODELISATION='DIS_TR',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh2[0])

#----------------------------------------------

Material

#----------------------------------------------

mater = DEFI_MATERIAU(ECRO_LINE=F(D_SIGM_EPSI=60000.0,
SY=300.0),
ELAS=F(E=210000.0,
NU=0.3))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
TOUT='OUI'),
MAILLAGE=mesh2[0],
MODELE=Sane)

#----------------------------------------------

Functions

#----------------------------------------------
listr = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1.0,
PAS=0.1))

xLoad = DEFI_FONCTION(NOM_PARA='INST',
VALE_FONC=listr,
VALE_PARA=listr)

times = DEFI_LIST_INST(DEFI_LIST=F(LIST_INST=listr),
ECHEC=F(EVENEMENT='ERREUR',
SUBD_METHODE='MANUEL',
SUBD_NIVEAU=10,
SUBD_PAS=10),
METHODE='AUTO')

#----------------------------------------------

Initial Crack

#----------------------------------------------

FISS = [None] * nbcalc

FISS[0] = DEFI_FISS_XFEM(DEFI_FISS=_F(DTAN=(0.0, 0.0, -1.0),
FORM_FISS='DEMI_PLAN',
NORMALE=(1.0, 0.0, 0.0),
PFON=(54.0, 0.0, 4.95)),
INFO=1,
MAILLAGE=mesh2[0],
TYPE_DISCONTINUITE='FISSURE',
TYPE_ENRI_FOND='TOPOLOGIQUE')

#----------------------------------------------

Chargements

#----------------------------------------------

To apply displacements

#----------------------------------------------

Propagation sur nbcalc pas de temps

#----------------------------------------------

modFrac = [None] * nbcalc
Constr = [None] * nbcalc
rot = [None] * nbcalc
vert = [None] * nbcalc
elemprop = [None] * nbcalc
RESU = [None] * nbcalc
SIF = [None] * nbcalc
Propa = [None] * nbcalc
Ma_visu = [None] * nbcalc
Mo_visu = [None] * nbcalc
raff = [None] * nbcalc
Lobster = [None] * nbcalc
Res_visu = [None] * nbcalc
h = 0.0025
RI = 0.1
RS = 4 * h
RP = RS

for i in range(nbcalc):
modFrac = MODI_MODELE_XFEM(FISSURE=FISS,
INFO=1,
MODELE_IN=Sane)

elemprop[i] = AFFE_CARA_ELEM(DISCRET=_F(CARA='K_TR_D_N',
                                        GROUP_MA=('assegnazione', ),
                                        VALE=(0.001,) * 6),
                             MODELE=modFrac[i])

Constr[i] = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                       DY=0.0,
                                       DZ=0.0,
                                       GROUP_MA=('Fix_surface', )),
                           MODELE=modFrac[i])

rot[i] = AFFE_CHAR_MECA(DDL_IMPO=_F(DRX=0.056478,
                                    GROUP_NO=('Load', )),
                        LIAISON_SOLIDE=_F(GROUP_NO=('Tw2', )),
                        MODELE=modFrac[i])

vert[i] = AFFE_CHAR_MECA(DDL_IMPO=_F(DZ=-5,
                                     GROUP_MA=('FlexLine', )),
                         MODELE=modFrac[i])

RESU[i] = MECA_STATIQUE(CARA_ELEM=elemprop[i],
                        CHAM_MATER=fieldmat,
                        EXCIT=(_F(CHARGE=Constr[i],
                                  FONC_MULT=xLoad),
                               _F(CHARGE=vert[i],
                                  FONC_MULT=xLoad),
                               _F(CHARGE=rot[i],
                                  FONC_MULT=xLoad)),
                        INST=1,
                        MODELE=modFrac[i],
                        SOLVEUR=_F(METHODE='MUMPS'))

SIF[i] = CALC_G_XFEM(LISSAGE=_F(LISSAGE_G='LAGRANGE',
                                LISSAGE_THETA='LAGRANGE'),
                     OPTION='CALC_K_G',
                     RESULTAT=RESU[i],
                     THETA=_F(FISSURE=FISS[i],
                              NB_POINT_FOND=NbFond,
                              NUME_FOND=Nume_Fond,
                              R_INF=0.1,
                              R_SUP=1))

IMPR_TABLE(TABLE=SIF[i], UNITE=3)

if i != (nbcalc - 1):
    FISS[i + 1] = CO('FISS_%d' % (i + 1))
    Propa[i] = PROPA_FISS(COMP_LINE=_F(COEF_MULT_MAXI=1.0, COEF_MULT_MINI=0.0),
                          DA_MAX=1,
                          FISSURE=_F(FISS_ACTUELLE=FISS[i],
                                     FISS_PROPAGEE=FISS[i + 1],
                                     NB_POINT_FOND=(NbFond,),
                                     TABLE=SIF[i]),
                          INFO=2,
                          LOI_PROPA=_F(C=1.0, LOI='PARIS', M=1.0, MATER=mater),
                          METHODE_PROPA='GEOMETRIQUE',
                          MODELE=modFrac[i],
                          RAYON=2,
                          TEST_MAIL='NON',
                          ZONE_MAJ='TOUT')

    Ma_visu[i] = POST_MAIL_XFEM(INFO=2, MODELE=modFrac[i])

    Mo_visu[i] = AFFE_MODELE(AFFE=_F(GROUP_MA=('FlexLine', 'Fix_surface', 'FaceRot', 'tutto'),
                                     MODELISATION='3D',
                                     PHENOMENE='MECANIQUE'),
                             MAILLAGE=Ma_visu[i])

    raff[i] = RAFF_XFEM(FISSURE=(FISS[i],), TYPE='DISTANCE')

    mesh2[i + 1] = CO("mesh2_%d" % (i + 1))

    Lobster[i] = MACR_ADAP_MAIL(ADAPTATION='RAFFINEMENT',
                                CHAM_GD=raff[i],
                                CRIT_RAFF_PE=0.1,
                                NIVE_MAX=5,
                                QUALITE="OUI",
                                DIAMETRE="OUI",  # Adjusted for missing comma
                                INFO=4,
                                LANGUE='ENGLISH',
                                MAILLAGE_N=mesh2[i],
                                MAILLAGE_NP1=mesh2[i + 1])

    Res_visu[i] = POST_CHAM_XFEM(MODELE_VISU=Mo_visu[i], RESULTAT=RESU[i])

    IMPR_RESU(RESU=_F(MAILLAGE=Ma_visu[i], RESULTAT=Res_visu[i]), UNITE=2)

FIN()
`
Thank you very much for the help.