[updated 3.Mar.2000]
*
SUBROUTINE FIPRIM(LAMBDA)
*
* Auteur : F.Mahe 23 avril 1992
* Derniere modification (F.Mahe) le 21 juin 1995
*
*
* NATURE : Calcule et donne les coefficients FIPRIMm pour l'assemblage
* de la matrice E.F.L. :
* FIPRIMm = sqrt(LAMBDAm-LAMBDA)
* ou LAMBDA est la valeur propre calculee a l'iteration precdente
* et LAMBDAm la valeur propre du m-ieme mode du guide plan borne.
* Ils sont ranges dans le tableau associe a chaque frontiere
* elements finis localises.
*
*
* DEMARCHE SUIVIE: Ici, le cas est simple et on donne directement le contenu
* des tableaux definissant les fonctions, sans avoir besoin
* de faire des calculs prealables.
*
* PARAMETRES EN ENTREE:
* LAMBDA Valeur propre initiale adim. ou de l'iteration precedente:
* LAMBDA = K^2*Np^2 + VALPRO
* LAMBDAm Valeurs propres adim. du guide plan borne:
* LAMBDAm = K^2*Np^2 + GAMMAm
* elles sont passees par le tableau "melina" 'TABMOD.
*
* PARAMETRES EN SORTIE:
* Le sous-tableau des -FIPRIMm dans les tableaux associes aux frontieres
* E.F.L. est rempli.
*
*
*
INCLUDE 'ALLOC'
INCLUDE 'CONTEX'
*
CHARACTER ERCODE*120
REAL LAMBDA
INTEGER MCDOMA,MCDOMB,MCDOMC,NBFDEV,NBTFON
INTEGER M,MCTMOD
*
COMMON /FORMAH/ERCODE
*
CALL TBAR1(ERCODE,'TABMOD',1,MCTMOD)
*
CALL KLFDEV('PSI',NBFDEV,NBTFON)
LGSTAB=NBFDEV*NBTFON
*
* Recuperation des adresses des tableaux associes aux frontieres E.F.L.
*
CALL TBAR3(ERCODE,'TBAS3a',1,MCASCA,'TBAS3b',1,MCASCB,
* 'TBAS3c',1,MCASCC)
* Calcul des adresses des sous-tableaux reels dans RST.
MCDOMA=2*(MCASCA+4+3*LGSTAB)-1
MCDOMB=2*(MCASCB+4+3*LGSTAB)-1
MCDOMC=2*(MCASCC+4+3*LGSTAB)-1
* Calcul et affectation des FIPRIMm
*
DO 100 M=1,NBFDEV
RRR=SQRT(RST(MCTMOD+M-1)-LAMBDA)
RST(MCDOMA+M-1)=RRR
RST(MCDOMB+M-1)=RRR
RST(MCDOMC+M-1)=RRR
100 CONTINUE
*
RETURN
*
END
fiprim est appelé dans (6 procédures)