rho_cp = 8869000

DEBUT(identifier=u'1:1',
      LANG='EN')

mesh = LIRE_MAILLAGE(identifier=u'2:1',
                     UNITE=2)

model0 = AFFE_MODELE(identifier=u'3:1',
                     AFFE=_F(MODELISATION=('3D', ),
                             PHENOMENE='THERMIQUE',
                             TOUT='OUI'),
                     MAILLAGE=mesh)

mater = DEFI_MATERIAU(identifier=u'4:1',
                      THER=_F(LAMBDA=237.0,
                              RHO_CP=rho_cp))

fieldmat = AFFE_MATERIAU(identifier=u'5:1',
                         AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh)

temp_f = DEFI_FONCTION(identifier=u'6:1',
                       NOM_PARA='INST',
                       VALE=(0.0, 20.0, 100.0, 214.9, 600.0, 290.0))

listr = DEFI_LIST_REEL(identifier=u'7:1',
                       DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=600.0,
                                     NOMBRE=20))

Hin = DEFI_FONCTION(identifier=u'8:1',
                    NOM_PARA='INST',
                    VALE=(0.0, 1.0, 100.0, 1.0, 600.0, 1.0))

load = AFFE_CHAR_THER_F(identifier=u'9:1',
                        ECHANGE=_F(COEF_H=Hin,
                                   GROUP_MA=('sokumen', ),
                                   TEMP_EXT=Hin),
                        MODELE=model0,
                        RAYONNEMENT=_F(EPSILON=temp_f,
                                       SIGMA=temp_f,
                                       TEMP_EXT=temp_f,
                                       TOUT='OUI'),
                        TEMP_IMPO=_F(GROUP_MA=('netugen', ),
                                     TEMP=temp_f))

load2 = AFFE_CHAR_THER(identifier=u'10:1',
                       MODELE=model0,
                       TEMP_IMPO=_F(GROUP_MA=('taikitemp', ),
                                    TEMP=20.0))

resther = THER_LINEAIRE(identifier=u'11:1',
                        CHAM_MATER=fieldmat,
                        ETAT_INIT=_F(VALE=20.0),
                        EXCIT=(_F(CHARGE=load),
                               _F(CHARGE=load2)),
                        INCREMENT=_F(LIST_INST=listr),
                        MODELE=model0)

IMPR_RESU(identifier=u'12:1',
          RESU=_F(RESULTAT=resther),
          UNITE=3)

FIN(identifier=u'13:1',
    )
