1 |       program mprofblack
   2 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   3 | !Trace dans le fichier prof_black.out le profil
   4 | !sous maille utilise pour la luminance et qui
   5 | ! a ete impose dans: profblack.f
   6 | 
   7 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   8 |       implicit none
   9 |       real*8 rfictif,abscisse,abmax,abmin,y
  10 |       real *8 profblack
  11 |       real *8 blanu, blat, dblae
  12 | 
  13 | 
  14 | 
  15 | 
  16 | 
  17 | 
  18 | 
  19 | 
  20 | 
  21 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  22 |       rfictif = 0.d+0
  23 |       abmax = 0.d+0
  24 |       abmin = 0.d+0
  25 |       open (24, file='prof_black.out', status='unknown')
  26 |       
  27 | 
  28 | !      do abscisse=-0.5,0.5,0.01
  29 | !      y = profblack(abmin,abmax,abscisse, rfictif)
  30 | !      write (24,*) abscisse, y
  31 | !      enddo
  32 | 
  33 | !!      do blat= 1000, 0, -1
  34 | !!      call black_dt(3250.d+0, blat, dblae)
  35 | !!      write (24,*) blat, dblae
  36 | !!      enddo 
  37 | 
  38 | !       call black_dt(557500.d+0,5.7447417d+0,dblae)
  39 | !       write (*,*) '557500.   ','5.7447417', dblae
  40 |       do blanu= 100000, 500000, 2000
  41 |       call black_dt(blanu, 5.d+0, dblae)
  42 | !      write (24,*) blat, dblae
  43 |       enddo       
  44 | 
  45 |       close(24)
  46 |       end