//#include <alberta.h>
//#include "newalberta.hh"
//#include "header mit PROBLEM_DATA - Definition" muss vor Einbindung passiert sein
#ifndef _heat_problem_
#define _heat_problem_ 1
//extern const  PROB_DATA_DIFF* prob_data_heat;

extern const REAL Celsius;//=273.0

extern  REAL th0(const REAL_D x) ;           /* initial temperature */
//{
//  return(723.0/*+Celsius*/);
//}

extern  REAL gth_D(const REAL_D x, REAL t); /* external temp */
//{
//  return(20.0+Celsius);/*dies ist keine Dirichlet-Bedingung, sondern Aussentemperatur f die Robin-Bedi
////}


extern  REAL gth_N_3(const REAL_D x, REAL t, S_CHAR bound);  /* Robin coef   */
//{
//  if (bound == -1 || bound == -6)
//    return(0.1);
//  else
//    return(0.01);
//}



// const REAL (*LALt(const EL_INFO *el_info, const QUAD *quad,
//                         int iq, void *ud))[DIM+1];

//extern  OPERATOR_INFO operator_heat;


//extern  PROB_DATA_DIFF prob_data_stdheat ; 24.4

void init_problem_stdheat(const PROB_DATA_DIFF **prob_data_heat)
;
#endif
