Salve a tutti, vorrei sapere da voi se il codice scritto va bene per quello che devo fare.il mio problema è il seguente: ho una geometria tridimensionale della quale mi interessa il piano y;z, ho una distribuzione a campana di dati, mi interessa che ad ogni valore di z sia assegnato il valore f(z).il codice è il seguente
DEFINE_PROFILE(flusso_radiativo, thread, position)
{
real x[ND_ND]; /* this will hold the position vector */
real z;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z = x[2];
F_PROFILE(f, thread, position) = 149529572.*1./(sqrt((2.*3.141593*(0.051*0.051))))*exp((-(z-0.2)*(z-0.2))/(2.*(0.051*0.051)));
}
end_f_loop(f, thread)
}
l'andamento di z e di y=f(z) è ne''allegato
Allegati:
10891_0b6eac46f34ba264711c90ed31386219.jpg