c******************************************************* SUBROUTINE base_chim c routine subordonnée de diffus c Formation de la base pour la diffusion des éléments chimiques c la fonction d'interpolation est d'ordre m c le vecteur nodal a m point en x(1), x(n) et m-i aux singularités c s'il n'y a pas de singularité une base continue est crée c dans les ZC, base continue non dérivable c entrées c 0 <= i < m : ordre de continuité aux points de singularité c i = 0 discontinue c i = 1 dérivée première discontinue c i = 2 dérivée seconde discontinue c ................. c i = m-1 fonction continue partout c is(0:ns+1) : indices des abscisses des singularités c ns : nombre de singularités c Auteur: P. Morel, Département Cassiopée, O.C.A. c CESAM2k c------------------------------------------------------------------ INTEGER, ALLOCATABLE, DIMENSION(:) :: mult c------------------------------------------------------------------ 2000 FORMAT(8es10.3) c multiplicités ALLOCATE(mult(n_ch)) c continuité et discontinuité des dérivées 1-ières aux limites ZR/ZC mult=1 ; mult(idis(1:ndis))=MAX(1,m_ch-1) c construction du vecteur nodal mult(1)=m_ch ; mult(n_ch)=m_ch ; knotc=SUM(mult) DEALLOCATE(mct) ; ALLOCATE(mct(knotc)) CALL noeud(mc,mct,mult,n_ch,knotc) DEALLOCATE(mult) RETURN END SUBROUTINE base_chim