1 |       double precision function profblack(abmin,abmax,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
  17 |       real a, b, xx,lum0
  18 | 
  19 | 
  20 | 
  21 | !!!!!!!!!! valeurs pour parabole entre 0 et 1
  22 |       a = -4
  23 |       b= 1
  24 |       xx =0
  25 | 
  26 | !!!!!!!!!! valeurs pour parabole entre b et b/10
  27 | c      lum0 = 5.d+0
  28 | cc      a= - lum0 /4.d+0
  29 | c      a = -(lum0 / 10.d+0 * 4.d+0) 
  30 | c      b=lum0/10.d+0  + lum0 
  31 | c      xx=0
  32 | !!!!!!!!!!!
  33 |       profblack = a* (abscisse-xx-rfictif)**2 +b 
  34 | 
  35 | !!!!!!!!!!nvelle valeur richard
  36 | !!      profblack = 0.5d+0*(abscisse+1-rfictif)/3.1416
  37 | !!      if (abscisse .eq. rfictif) then
  38 | !!         profblack = 0.d+0
  39 | !!      endif
  40 | 
  41 | 
  42 | !!! pour richard: test fig a et b 
  43 | !!! fichier rayon 0 0.1 0.2 0.3 ...... 1.0 1.0      
  44 | cccccccccc      xx = abscisse-rfictif
  45 | ccccccccc      profblack = xx * (1-xx)*4
  46 | 
  47 |       return
  48 | 
  49 |       end
  50 |