Besoin d'aide en Fortran

[alvarogill]

Bonjour,

je voudrais que quelqu'un puisse m'aider à résoudre un exercice qui demande d'être programmé en langage FORTRAN.
Voici l'exercice:

j'ai une ligne divisée en 13 parties que l'on peut appeler n. Cette ligne doit être divisée en plusieurs partitions que je nomme K. Par exemple, si j'ai 2K partitions, j'aurais donc 3n parties, ces n parties ne sont pas obligatoirement égales, du moment que l'on ai 13 parties au total.

Je voudrais exécuter un programme qui lise K et n, suivant le nombre de partitions que je veux mettre.
Merci beaucoup.

[OG]

Tiens je croyais que le fortran n'existait plus ?
Ok je sors

O.G.

[MB]

Je ne comprends pas trop l'explication de ton problème. Les notations ne me semblent pas très limpides.

[alvarogill]

Y a t'il quelqu'un qui sait programmer en FORTRAN?
Je pense que le problême est assez comprehensif.
Je veux tout simplement diviser une ligne en k partition

[Arnaud]

Ton problème n'est pas assez précis, en tout cas pas pour moi.
Si tu veux "splitter" ta ligne, il faut une condition claire, ou un délimiteur facile à reconnaitre ( une virgule ? un trait-d'union ? le nombre de caractères dans la ligne ? ).

Tu peux regarder et adapter le code de cette page :

http://www.developpez.net/forums/d25056 ... ion-split/

[alvarogill]

En fait, je suis en train de travailler sur un programme qui calcule les effets aerodynamiques d'une aile d'un avion.
Je ne vais pas mettre tout le programme car il est assez long, mais je vous envoie la partie qui traite de sa géométrie.

Code: Tout sélectionner

   subroutine geo(b,c,S,AR,ib,jb,dx,dy,dgap,alfa)
   dimension bb(13),alf(5),sn(5),cs(5),sno(5),cso(5)
   dimension qf(5,16,3),qc(4,13,3),ds(4,13)
   dimension vortic(50,13),qw(50,14,3),a1(4,13),alamda(2)
   dimension vort1(50,13),qw1(50,14,3)
   common vortic,qw,vort1,qw1,qf,a1
   common it,alf,sno,cso,bb,qc,ds,alamda,dxw
!
   pay      = 3.141592654
! ib = no. of chordwise boxes; jb = no. of spanwise boxes. [b](Traduction: ib=Nombre de carré au long de la corde;jb=nombre de carré au long de l'envergure.[/b]
   ib1      = ib+1                [b](ib=4 et jb=13) donc 13*4 = 52 carrés au total qui vont constituer mon aile[/b]
   jb1      = jb+1
   do 2 i=1,ib1
     sn(i)   = sin(alf(i))
     cs(i) = cos(alf(i))
2   continue
   ctg1   = tan(pay/2.0-alamda(1))
   ctg2   = tan(pay/2.0-alamda(2))
   ctip   = c+b*(ctg2-ctg1)                    [b](c c'  est la corde de l'avion; b c'est l'envergure,alamada (1)correspond a l'angle au bord d'attaque et  alamda (2)[/b]
   S      = b*(c+ctip)/2.0             [b] corresponde à l'angle de fuite de l'aile).[/b]
   AR      = 2.0*b*b/S                     
!
! Wing fixed vortices location (qf(i,j,(x,y,z))...)
!
   bj      = 0.0                   
                                                                            [b](dy=b/jb)
                                                                      (bb(j)=dy; donc pour moi bj correspond a l'abcisse de l'envergure, par exemple a l'origine bj=0 et
                                                                            a l'extremité de l'aile bj=6.[/b]


   do 3 j=1,jb1                                                             [b](il fait un cycle qui va du premier segment jusqu'au 14 eme segment puisque l'on a 13 carré.[/b]
     if(j.GT.1) bj   = bj+bb(j-1)                                     [b]  VOILA: en fait c'est ici que je voudrais que vous m'aidiez.
                                                                                    Je veux diviser cette envergure en N parties pour pouvoir effectuer après plusieurs
                                                                                      angle du bord d'attaque et de fuite de l'aile. [/b]
     z1   = 0.0
! dc1 = leading edge x; dc2 = trailing edge x.
     dc1   = bj*ctg1
     dc2   = bj*ctg2
     dx1   = (c+dc2-dc1)/ib
   
     do 1 i=1,ib
       qf(i,j,1)   = dc1+dx1*(i-0.75)
       qf(i,j,2)   = bj
       qf(i,j,3)   = z1-0.25*dx1*sn(i)
       z1         = z1-dx1*sn(i)
   
   
1     continue

J'ai envoyé le premier message avec un shema plus simple pour que l'on puisse mieux comprendre ma question.
Je vous remercie.

[alvarogill]

Pardon

c c' est la corde de l' aile et b c'est l'envergure de l'aile de l'avion

[Arnaud]

En fait avec le code je comprends encore moins ta question...( je ne connais pas la syntaxe fortran )
On ne doit pas avoir la même définition de "ligne".
Pour moi une "ligne" est une chaine de caractères, or il semblerait que dans ton programme, tu aies un problème dans une boucle while, sur un nombre...( bj en l'occurence ).

[alvarogill]

En fait, avec ligne, je voulais dire segment. Désolé, pour le malentendu.
Est-ce que même comme ça mathématiquement, serait-il possible de résoudre mon problême.
Peut-être qu'à partir de la j'arrive a trouver une lumière au fond du tunnel
Donc, j'ai une droite que je dois diviser en K parties. Mais cette droite et déjà divisé par n parties.
Dans mon exemple j'ai mis n=13 parties. Sachant que ces K parties doivent coincider parfaitement avec un des n parties.
Je ne sais pas si j' ai été plus clair... j'ai pensé utiliser une succession.

Grand merci

[alvarogill]

Guten tag Arnaud

Connaissez-vous quelqu'un qui sache programmer sur FORTRAN, sinon c'est pas grave.
Le problème c'est que, en ce moment je suis en train de préparer une thèse et il y a seulement peu de temps que je m'habitue petit à petit a ce langage (archaique, diront certains, ).

Vielen dank.

[Arnaud]

Salut,

Demande peut-être sur developpez.net, il y a un sous-forum spécialisé dans fortran.