Bonjour,
Je veux discrétiser un problème en deux dimensions avec des conditions aux bords ( neumann + dirichlet ) :
la température est imposée si y < ou égale à la hauteur sur deux ( H/2)
et un pour y > H/ 2, le flux est nul.
Ma question est la suivante : lors du passage aux indices, comment on va discrétiser cette condition ?
J'ai procédé comme suit : y>H/2 ==> j>(n-1)/2, mais je ne sais pas comment l'implémenter sur fortran.
Merci
-----