Télécharger fluxbc.eso

Retour à la liste

Numérotation des lignes :

  1. C FLUXBC SOURCE CHAT 05/01/13 00:05:13 5004
  2.  
  3. FUNCTION FLUXBC (RX1,ZX1,RX2,ZX2,DSURF,NE,ELIP)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. REAL*8 RX1,RX2,ZX1,ZX2,DSURF,AS1,F,DH,AX
  7. REAL*8 ELI1,ELI2,ELI
  8. INTEGER NF,NE
  9. DIMENSION ELIP(101)
  10.  
  11. C WRITE(6,*) '##FluxBC##'
  12. IF (RX2.EQ.0) THEN
  13. c WRITE(6,*) 'ERREUR il faut un rayon d inducteur > 0 '
  14. call erreur(959)
  15. else
  16.  
  17. IF (RX1.EQ.0) THEN
  18. FLUXBC=0.
  19. ELSE
  20. F=8.54
  21. DH=DSURF/F
  22. AS1= (RX1*RX1)+(RX2*RX2)+((ZX1-ZX2)**2)
  23. as1 = as1 + dh
  24. * IF ((AS1-(2*RX1*RX2)).LE.DH) THEN
  25. * AS1=(2*RX1*RX2)+DH
  26. * END IF
  27. AX= RX1*RX2/AS1
  28. NF=INT (2*AX*NE)
  29. ELI1=ELIP(NF+1)*(NF+1-(2*AX*NE))
  30. ELI2=ELIP(NF+2)*((2*AX*NE)-NF)
  31. ELI=ELI1+ELI2
  32. ELI=ELI-(0.707107*(LOG(1-(4*AX*AX))))+(0.3127313*AX*AX)
  33. FLUXBC=ELI*(AS1**0.5)/10.
  34. RETURN
  35. END IF
  36. end if
  37. return
  38. END
  39.  
  40.  
  41.  

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