1 |       double precision function proftemp(abscisse, rfictif)
   2 | 
   3 | c*******************************************************************
   4 | c     cette function renvoie la valeur de la luminance 
   5 | c     pour une abscisse donne dans un volume.
   6 | c     Elle utilise un profil sous maille
   7 | c
   8 | c     in:   abscisse: min max et position
   9 | c           pivots du profil (si interpolation)
  10 | c
  11 | c     out:  la luminance noire = profblack
  12 | c*******************************************************************
  13 | c     profblack = (abscisse**2-abmin**2)/(abmax**2-abmin**2)
  14 | c*******************************************************************
  15 |       implicit none
  16 |       double precision abmin, abmax, abscisse, rfictif,a,c
  17 |       double precision tmin,tmax,xcentre,x2,xx,x1
  18 | 
  19 | 
  20 | 
  21 | !!! pour richard: test fig a et b 
  22 | !!! fichier rayon 0 0.1 0.2 0.3 ...... 1.0 1.0      
  23 |       xx = abscisse-rfictif
  24 | !!!      proftemp = xx * (1-xx)*4
  25 | !!!profils en parabole
  26 |       tmin=500
  27 |       tmax=2500
  28 |       xcentre=0.1
  29 |       x2=0.2
  30 |       a=(tmin-tmax)/(xcentre*(xcentre-x2))
  31 |       c=tmax
  32 |       proftemp = a*xx*(xx-x2)+c
  33 | 
  34 | !!!profils en pointe
  35 | !      tmin=500
  36 | !      tmax=2500
  37 | !      x1=0.45
  38 | !      xcentre=0.5
  39 | !      x2=0.55
  40 | !      if ((xx.ge.x1).and.(xx.le.xcentre)) then
  41 | !         a = (tmax-tmin)/(xcentre-x1)
  42 | !         c = tmax - a*xcentre
  43 | !         proftemp= a*xx+c 
  44 | !      elseif ((xx.gt.xcentre).and.(xx.le.x2)) then
  45 | !         a = (tmin-tmax)/(x2-xcentre)
  46 | !         c=  tmax -a*xcentre
  47 | !         proftemp=a*xx+c
  48 | !      else
  49 | !         proftemp= tmin
  50 | !      endif
  51 |       
  52 |       !print *, xx,proftemp
  53 | 
  54 | 
  55 | 
  56 |       return
  57 | 
  58 |       end
  59 |