Télécharger la plaquette Cast3M


fr en

Télécharger crbary.eso

Retour à la liste

Numérotation des lignes : 

crbary

  1. C CRBARY    SOURCE    KICH      22/11/21    21:15:03     11505          
  2.       subroutine crbary(ith)
  3. C
  4. C
  5.       implicit integer(i-n)
  6.       implicit real*8(a-h,o-z)
  7.  
  8. -INC  CCASSIS
  9.  
  10.       PARAMETER(MASDIM=64)
  11.       common/CCONNE/iwrk3,ipmodl1,xmultl,icle1,ihgsel
  12.       common/CCONNE/jconl,ihg1,ihg2,ihg,nbthr,ixlong1
  13.       common/CCONN1/d,pt1(3),pt2(3)
  14.       common/CCONN1/xmn(masdim),ymn(masdim),zmn(masdim)
  15.       common/CCONN1/xmx(masdim),ymx(masdim),zmx(masdim)
  16.       common/CCONN1/hmxt(masdim),xlg2m(masdim)
  17.  
  18.       segment hgsele
  19.          real*8 xmult,ymult,zmult
  20.          real*8 hmaxt,xlong2,tmax,tmin,xlong2m
  21.          integer nels,nbpb,ipass
  22.          integer nbzt,indt,khug
  23.       endsegment
  24.       segment hg
  25.         integer ielh(nbpb,2)
  26.         real*8 hmax(nbpb)
  27.         real*8 xll(nbpb)
  28.         integer inoa(nbpb,nsous+1)
  29.       endsegment
  30.       segment conl
  31.          integer iconl(nbpb)
  32.       endsegment
  33.       segment hg1
  34.         real*8 hcoor(3*nbpb)
  35.       endsegment
  36.       pointeur hg2.hg1
  37.       hgsele = ihgsel
  38.       icle = icle1
  39.       ixlong = ixlong1
  40.       hg = ihg
  41.       ipmodl = ipmodl1
  42.       nbthr1 = nbthr
  43.       if (ith.gt.nbthr1) goto 9999
  44.       ires = mod(nbpb,nbthr1)
  45.       if (ires.eq.0) then
  46.          ilon = nbpb / nbthr1
  47.          ideb = (ith - 1) * ilon + 1
  48.       else
  49.          if (ith.le.ires) then
  50.             ilon = nbpb / nbthr1 + 1
  51.             ideb = (ith - 1) * ilon + 1
  52.          else
  53.             ilon = nbpb / nbthr1
  54.             ideb = (ires * (ilon + 1)) + (ith - ires - 1) * ilon + 1
  55.          endif
  56.       endif
  57.       ifin = ideb + ilon - 1
  58.       if (ifin.ge.ideb) then
  59.           call crbar1(iwrk3,ipmodl,ideb,ifin,xmultl,icle,d,pt1,pt2
  60.      +      ,jconl,ihg1,ihg2,xmn,ymn,zmn,xmx,ymx,zmx
  61.      +      ,hmxt,ihg,ith,ixlong,xlong2,xlg2m)
  62.       endif
  63.  9999 continue
  64.       return
  65.       end
  66.  
  67.  
  68.  
  69.  

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