Numérotation des lignes :

clmbda
C CLMBDA    SOURCE    AM        12/12/20    21:15:01     7631        SUBROUTINE CLMBDA(D1,D2,D3,LAMBDA,YOUNG,NU,ALPH,COMP)      IMPLICIT INTEGER(I-N)      implicit real*8 (a-h,o-z)        REAL*8  D1,D2,D3,ALPH        REAL*8  YOUNG,NU        REAL*8  S0302,S0320,S2100        REAL*8  LAMBDA(6)        REAL*8  PI        REAL*8  X1,sinpixs2,sinpiys2,sinpizs2        LOGICAL COMP************************************************************************************************************************************* CALCUL DE LA MATRICE LAMDA DANS LE CAS OU IL N'Y A PAS DE DELTA**     (PAS DE DECOMPOSITION DE SIGMA-N EN SIGMA-D + SIGMA-V)     *************************************************************************************************************************************CC INITIALISATION DES VARIABLESC        PI=4.D0*ATAN(1.D0)        IF (COMP) THEN          AF=1.D0-(0.5*(D2+D3))**ALPH          BF=1.D0-(0.5*(D3+D1))**ALPH          CF=1.D0-(0.5*(D1+D2))**ALPH        ELSE          AF=1.D0-D1          BF=1.D0-D2          CF=1.D0-D3        END IFC        print*,'dans callambda: alpha=',alphC        print*,'af,bf,cf',af,bf,cfC        print*,'D3',D3,'d1',d1,'D2',D2CC ON A DIVISE EN HAUT ET EN BAS DE X1/DET PAR AF*BF*CFC -> X1 VAUT 1 (CF FIJ)C        X1=AF*BF*CFCC CALCUL DES INTEGRALES SUR 1/4 DE SPHEREC       sinpixs2=sin(pi*AF/2.D0)       sinpiys2=sin(pi*BF/2.D0)       sinpizs2=sin(pi*CF/2.D0)      S0302=     & +  .49587724947700D+02*sinpixs2    *sinpiys2    *sinpizs2     & -.11826019614403D+03  *sinpixs2    *sinpiys2     *sinpizs2**2     & +  .15888740080185D+03*sinpixs2    *sinpiys2     *sinpizs2**3     & -.12090181599425D+03  *sinpixs2    *sinpiys2     *sinpizs2**4     & +  .48134780378626D+02*sinpixs2    *sinpiys2     *sinpizs2**5     & -.76097170495877D+01  *sinpixs2    *sinpiys2     *sinpizs2**6     & -.42982186153999D+02  *sinpixs2    *sinpiys2**2 *sinpizs2     & +  .92334205518749D+02*sinpixs2    *sinpiys2**2  *sinpizs2**2     & -.94595653059631D+02  *sinpixs2    *sinpiys2**2  *sinpizs2**3     & +  .46508202980390D+02*sinpixs2    *sinpiys2**2  *sinpizs2**4     & -.85557270110063D+01  *sinpixs2    *sinpiys2**2  *sinpizs2**5     & -.59795082278790D+01  *sinpixs2    *sinpiys2**3 *sinpizs2     & -.26820923709862D+02  *sinpixs2    *sinpiys2**3  *sinpizs2**2     & +  .23967186388666D+02*sinpixs2    *sinpiys2**3  *sinpizs2**3     & -.68917781168937D+01  *sinpixs2    *sinpiys2**3  *sinpizs2**4     & +  .49597006251687D+02*sinpixs2    *sinpiys2**4 *sinpizs2     & -.38087148378371D+01  *sinpixs2    *sinpiys2**4  *sinpizs2**2      S0302=S0302     & -.62685158801033D+00  *sinpixs2    *sinpiys2**4  *sinpizs2**3     & -.45758829115125D+02  *sinpixs2    *sinpiys2**5 *sinpizs2     & +  .21517303738659D+01*sinpixs2    *sinpiys2**5  *sinpizs2**2     & +  .15441307857503D+02*sinpixs2    *sinpiys2**6 *sinpizs2     & -.11825524477172D+03  *sinpixs2**2 *sinpiys2    *sinpizs2     & +  .16747889813112D+03*sinpixs2**2 *sinpiys2     *sinpizs2**2     & -.13676189436854D+03  *sinpixs2**2 *sinpiys2     *sinpizs2**3     & +  .59638830553943D+02*sinpixs2**2 *sinpiys2     *sinpizs2**4     & -.10868512220306D+02  *sinpixs2**2 *sinpiys2     *sinpizs2**5     & +  .92341016512039D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2     & -.82914555064696D+02  *sinpixs2**2 *sinpiys2**2  *sinpizs2**2     & +  .41590292781164D+02*sinpixs2**2 *sinpiys2**2  *sinpizs2**3     & -.83827514727064D+01  *sinpixs2**2 *sinpiys2**2  *sinpizs2**4     & -.26813764702232D+02  *sinpixs2**2 *sinpiys2**3 *sinpizs2     & +  .14967505013368D+02*sinpixs2**2 *sinpiys2**3  *sinpizs2**2     & -.43826806253910D+01  *sinpixs2**2 *sinpiys2**3  *sinpizs2**3     & -.38031226905467D+01  *sinpixs2**2 *sinpiys2**4 *sinpizs2     & +  .10826472143751D+01*sinpixs2**2 *sinpiys2**4  *sinpizs2**2      S0302=S0302     & +  .21539434629142D+01*sinpixs2**2 *sinpiys2**5 *sinpizs2     & +  .15885885518183D+03*sinpixs2**3 *sinpiys2    *sinpizs2     & -.13675927885166D+03  *sinpixs2**3 *sinpiys2     *sinpizs2**2     & +  .62883150455112D+02*sinpixs2**3 *sinpiys2     *sinpizs2**3     & -.11888439526808D+02  *sinpixs2**3 *sinpiys2     *sinpizs2**4     & -.94627040869973D+02  *sinpixs2**3 *sinpiys2**2 *sinpizs2     & +  .41591649499595D+02*sinpixs2**3 *sinpiys2**2  *sinpizs2**2     & -.81213735252988D+01  *sinpixs2**3 *sinpiys2**2  *sinpizs2**3     & +  .23945386833297D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2     & -.43826106229617D+01  *sinpixs2**3 *sinpiys2**3  *sinpizs2**2     & -.63409891017392D+00  *sinpixs2**3 *sinpiys2**4 *sinpizs2     & -.12083164855111D+03  *sinpixs2**4 *sinpiys2    *sinpizs2     & +  .59635568525663D+02*sinpixs2**4 *sinpiys2     *sinpizs2**2     & -.11888504501967D+02  *sinpixs2**4 *sinpiys2     *sinpizs2**3     & +  .46560572564192D+02*sinpixs2**4 *sinpiys2**2 *sinpizs2     & -.83833965641336D+01  *sinpixs2**4 *sinpiys2**2  *sinpizs2**2     & -.68748680983310D+01  *sinpixs2**4 *sinpiys2**3 *sinpizs2     & +  .48058063063930D+02*sinpixs2**5 *sinpiys2    *sinpizs2      S0302=S0302     & -.10867120963554D+02  *sinpixs2**5 *sinpiys2     *sinpizs2**2     & -.85838947053063D+01  *sinpixs2**5 *sinpiys2**2 *sinpizs2     & -.75792226760522D+01  *sinpixs2**6 *sinpiys2    *sinpizs2CCC      S0320=     & +  .49587727033586D+02*sinpixs2    *sinpiys2    *sinpizs2     & -.11826018579995D+03  *sinpixs2    *sinpiys2     *sinpizs2**2     & +  .15888737682738D+03*sinpixs2    *sinpiys2     *sinpizs2**3     & -.12090180560884D+03  *sinpixs2    *sinpiys2     *sinpizs2**4     & +  .48134782941498D+02*sinpixs2    *sinpiys2     *sinpizs2**5     & -.76097182547511D+01  *sinpixs2    *sinpiys2     *sinpizs2**6     & -.11825524141371D+03  *sinpixs2    *sinpiys2**2 *sinpizs2     & +  .16747884821826D+03*sinpixs2    *sinpiys2**2  *sinpizs2**2     & -.13676182945011D+03  *sinpixs2    *sinpiys2**2  *sinpizs2**3     & +  .59638796128496D+02*sinpixs2    *sinpiys2**2  *sinpizs2**4     & -.10868505584411D+02  *sinpixs2    *sinpiys2**2  *sinpizs2**5     & +  .15885885515815D+03*sinpixs2    *sinpiys2**3 *sinpizs2     & -.13675922155796D+03  *sinpixs2    *sinpiys2**3  *sinpizs2**2     & +  .62883118221758D+02*sinpixs2    *sinpiys2**3  *sinpizs2**3     & -.11888433966211D+02  *sinpixs2    *sinpiys2**3  *sinpizs2**4     & -.12083167330411D+03  *sinpixs2    *sinpiys2**4 *sinpizs2     & +  .59635549771322D+02*sinpixs2    *sinpiys2**4  *sinpizs2**2     & -.11888500424114D+02  *sinpixs2    *sinpiys2**4  *sinpizs2**3      S0320=S0320     & +  .48058086050282D+02*sinpixs2    *sinpiys2**5 *sinpizs2     & -.10867121049526D+02  *sinpixs2    *sinpiys2**5  *sinpizs2**2     & -.75792275647182D+01  *sinpixs2    *sinpiys2**6 *sinpizs2     & -.42982213136941D+02  *sinpixs2**2 *sinpiys2    *sinpizs2     & +  .92334241631458D+02*sinpixs2**2 *sinpiys2     *sinpizs2**2     & -.94595651357197D+02  *sinpixs2**2 *sinpiys2     *sinpizs2**3     & +  .46508203600717D+02*sinpixs2**2 *sinpiys2     *sinpizs2**4     & -.85557322153729D+01  *sinpixs2**2 *sinpiys2     *sinpizs2**5     & +  .92341059860355D+02*sinpixs2**2 *sinpiys2**2 *sinpizs2     & -.82914586469419D+02  *sinpixs2**2 *sinpiys2**2  *sinpizs2**2     & +  .41590274527138D+02*sinpixs2**2 *sinpiys2**2  *sinpizs2**3     & -.83827399930577D+01  *sinpixs2**2 *sinpiys2**2  *sinpizs2**4     & -.94627049453367D+02  *sinpixs2**2 *sinpiys2**3 *sinpizs2     & +  .41591614676660D+02*sinpixs2**2 *sinpiys2**3  *sinpizs2**2     & -.81213580886453D+01  *sinpixs2**2 *sinpiys2**3  *sinpizs2**3     & +  .46560584004622D+02*sinpixs2**2 *sinpiys2**4 *sinpizs2     & -.83833812204707D+01  *sinpixs2**2 *sinpiys2**4  *sinpizs2**2     & -.85839039617849D+01  *sinpixs2**2 *sinpiys2**5 *sinpizs2      S0320=S0320     & -.59794809896161D+01  *sinpixs2**3 *sinpiys2    *sinpizs2     & -.26820979413895D+02  *sinpixs2**3 *sinpiys2     *sinpizs2**2     & +  .23967192097949D+02*sinpixs2**3 *sinpiys2     *sinpizs2**3     & -.68917762822240D+01  *sinpixs2**3 *sinpiys2     *sinpizs2**4     & -.26813831452726D+02  *sinpixs2**3 *sinpiys2**2 *sinpizs2     & +  .14967600520974D+02*sinpixs2**3 *sinpiys2**2  *sinpizs2**2     & -.43826992963629D+01  *sinpixs2**3 *sinpiys2**2  *sinpizs2**3     & +  .23945404982568D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2     & -.43826231635356D+01  *sinpixs2**3 *sinpiys2**3  *sinpizs2**2     & -.68748689102157D+01  *sinpixs2**3 *sinpiys2**4 *sinpizs2     & +  .49597019950936D+02*sinpixs2**4 *sinpiys2    *sinpizs2     & -.38087036456114D+01  *sinpixs2**4 *sinpiys2     *sinpizs2**2     & -.62684827232429D+00  *sinpixs2**4 *sinpiys2     *sinpizs2**3     & -.38031088909489D+01  *sinpixs2**4 *sinpiys2**2 *sinpizs2     & +  .10826222088093D+01*sinpixs2**4 *sinpiys2**2  *sinpizs2**2     & -.63410176286672D+00  *sinpixs2**4 *sinpiys2**3 *sinpizs2     & -.45758849889839D+02  *sinpixs2**5 *sinpiys2    *sinpizs2     & +  .21517325427401D+01*sinpixs2**5 *sinpiys2     *sinpizs2**2      S0320=S0320     & +  .21539481683039D+01*sinpixs2**5 *sinpiys2**2 *sinpizs2     & +  .15441312720879D+02*sinpixs2**6 *sinpiys2    *sinpizs2CCC      S2100=     & +  .49587467973979D+02*sinpixs2    *sinpiys2    *sinpizs2     & -.42980769227454D+02  *sinpixs2    *sinpiys2     *sinpizs2**2     & -.59805404737685D+01  *sinpixs2    *sinpiys2     *sinpizs2**3     & +  .49599317855294D+02*sinpixs2    *sinpiys2     *sinpizs2**4     & -.45759239038724D+02  *sinpixs2    *sinpiys2     *sinpizs2**5     & +  .15441895510618D+02*sinpixs2    *sinpiys2     *sinpizs2**6     & -.11825825078992D+03  *sinpixs2    *sinpiys2**2 *sinpizs2     & +  .92332492764995D+02*sinpixs2    *sinpiys2**2  *sinpizs2**2     & -.26819401505354D+02  *sinpixs2    *sinpiys2**2  *sinpizs2**3     & -.38091597509388D+01  *sinpixs2    *sinpiys2**2  *sinpizs2**4     & +  .21520684276939D+01*sinpixs2    *sinpiys2**2  *sinpizs2**5     & +  .15888385621676D+03*sinpixs2    *sinpiys2**3 *sinpizs2     & -.94594703123626D+02  *sinpixs2    *sinpiys2**3  *sinpizs2**2     & +  .23966025562355D+02*sinpixs2    *sinpiys2**3  *sinpizs2**3     & -.62703651861082D+00  *sinpixs2    *sinpiys2**3  *sinpizs2**4     & -.12089720898763D+03  *sinpixs2    *sinpiys2**4 *sinpizs2     & +  .46508836417526D+02*sinpixs2    *sinpiys2**4  *sinpizs2**2     & -.68912765229417D+01  *sinpixs2    *sinpiys2**4  *sinpizs2**3      S2100=S2100     & +  .48129967460623D+02*sinpixs2    *sinpiys2**5 *sinpizs2     & -.85566605377921D+01  *sinpixs2    *sinpiys2**5  *sinpizs2**2     & -.76074803611026D+01  *sinpixs2    *sinpiys2**6 *sinpizs2     & -.11825824332720D+03  *sinpixs2**2 *sinpiys2    *sinpizs2     & +  .92332489273672D+02*sinpixs2**2 *sinpiys2     *sinpizs2**2     & -.26819394292098D+02  *sinpixs2**2 *sinpiys2     *sinpizs2**3     & -.38091608415685D+01  *sinpixs2**2 *sinpiys2     *sinpizs2**4     & +  .21520663179297D+01*sinpixs2**2 *sinpiys2     *sinpizs2**5     & +  .16747419796953D+03*sinpixs2**2 *sinpiys2**2 *sinpizs2     & -.82911722127338D+02  *sinpixs2**2 *sinpiys2**2  *sinpizs2**2     & +  .14966176247075D+02*sinpixs2**2 *sinpiys2**2  *sinpizs2**3     & +  .10825591719184D+01*sinpixs2**2 *sinpiys2**2  *sinpizs2**4     & -.13675642532919D+03  *sinpixs2**2 *sinpiys2**3 *sinpizs2     & +  .41589140819666D+02*sinpixs2**2 *sinpiys2**3  *sinpizs2**2     & -.43821050567095D+01  *sinpixs2**2 *sinpiys2**3  *sinpizs2**3     & +  .59639580135436D+02*sinpixs2**2 *sinpiys2**4 *sinpizs2     & -.83812695652628D+01  *sinpixs2**2 *sinpiys2**4  *sinpizs2**2     & -.10870415781963D+02  *sinpixs2**2 *sinpiys2**5 *sinpizs2      S2100=S2100     & +  .15888383027537D+03*sinpixs2**3 *sinpiys2    *sinpizs2     & -.94594699500019D+02  *sinpixs2**3 *sinpiys2     *sinpizs2**2     & +  .23966016553017D+02*sinpixs2**3 *sinpiys2     *sinpizs2**3     & -.62703218218355D+00  *sinpixs2**3 *sinpiys2     *sinpizs2**4     & -.13675640659286D+03  *sinpixs2**3 *sinpiys2**2 *sinpizs2     & +  .41589154124408D+02*sinpixs2**3 *sinpiys2**2  *sinpizs2**2     & -.43821101691482D+01  *sinpixs2**3 *sinpiys2**2  *sinpizs2**3     & +  .62875012679184D+02*sinpixs2**3 *sinpiys2**3 *sinpizs2     & -.81231387874529D+01  *sinpixs2**3 *sinpiys2**3  *sinpizs2**2     & -.11886187812260D+02  *sinpixs2**3 *sinpiys2**4 *sinpizs2     & -.12089717335754D+03  *sinpixs2**4 *sinpiys2    *sinpizs2     & +  .46508832417159D+02*sinpixs2**4 *sinpiys2     *sinpizs2**2     & -.68912742079210D+01  *sinpixs2**4 *sinpiys2     *sinpizs2**3     & +  .59639553850800D+02*sinpixs2**4 *sinpiys2**2 *sinpizs2     & -.83812725667904D+01  *sinpixs2**4 *sinpiys2**2  *sinpizs2**2     & -.11886185734126D+02  *sinpixs2**4 *sinpiys2**3 *sinpizs2     & +  .48129948086467D+02*sinpixs2**5 *sinpiys2    *sinpizs2     & -.85566592717808D+01  *sinpixs2**5 *sinpiys2     *sinpizs2**2      S2100=S2100     & -.10870405546200D+02  *sinpixs2**5 *sinpiys2**2 *sinpizs2     & -.76074772239863D+01  *sinpixs2**6 *sinpiys2    *sinpizs2CCCCC DEFINITION DE LA MATRICE DE RIGIDITE LAMBDAC ON CALCULE LES INTEGRALES ENTRE 0 ET PI/2 POUR TETA ET 0 ET PI POUR PHIC -> IL FAUT MULTIPLIER PAR 4 LES INTEGRALESC          LAMBDA(1)=S0320          LAMBDA(2)=S0302          LAMBDA(3)=S2100          LAMBDA(4)=(S0320+S0302)/2.D0          LAMBDA(5)=(S0320+S2100)/2.D0          LAMBDA(6)=(S0302+S2100)/2.D0C       WRITE (6,200)  LAMBDAC       WRITE (6,100)  S0320,S0302,S2100,S0300C       WRITE (6,100)  S2340,S0322,S2300,S0302,S2100,S0300C       WRITE (6,300) D1,D2,D3C200    FORMAT ('LAMBDA = ',/,6 ( ' |',6 (E13.6,' , '),' |',/),//)C100    FORMAT ('SOM = ',/,6 ( ' |',6 (E11.4,' , '),' |',/),//)C300    FORMAT ('D1=',E11.4,'D2=',E11.4,'D3=',E11.4)        RETURN        END     

© Cast3M 2003 - Tous droits réservés.
Mentions légales