      subroutine sand(r)
      real r,sr
      integer l,i
      dimension sr(1501),r(1501)
c
c     sand average reflectance
c     warning : values of dry sand ground reflectance are given
c     between 0.4 and 2.2 microns. outside this interval the
c     values are set to 0.
c
 
      data (sr(l),l=1,128)/  58*0.,
     a .00000, .04550, .09100, .09100, .09100, .09100, .09100,
     a .09100, .09100, .09100, .09100, .09100, .09100, .09300,
     a .09500, .09500, .09500, .09500, .09500, .09500, .09500,
     a .09500, .09500, .09500, .09500, .09500, .09500, .09600,
     a .09700, .09700, .09700, .09850, .10000, .10150, .10300,
     a .10300, .10300, .10500, .10700, .10700, .10700, .10700,
     a .10700, .10700, .10700, .10850, .11000, .11200, .11400,
     a .11550, .11700, .11900, .12100, .12100, .12100, .12300,
     a .12500, .12500, .12500, .12500, .12500, .12650, .12800,
     a .12700, .12600, .12850, .13100, .13250, .13400, .13400/
      data (sr(l),l=129,198)/  
     a .13400, .13400, .13400, .13400, .13400, .13400, .13400,
     a .13400, .13400, .13750, .14100, .14300, .14500, .14700,
     a .14900, .14900, .14900, .15150, .15400, .15550, .15700,
     a .15850, .16000, .16150, .16300, .16300, .16300, .16450,
     a .16600, .16750, .16900, .17100, .17300, .17500, .17700,
     a .17900, .18100, .18300, .18500, .18700, .18900, .19200,
     a .19500, .19700, .19900, .20150, .20400, .20400, .20400,
     a .20600, .20800, .21100, .21400, .21700, .22000, .22200,
     a .22400, .22600, .22800, .23050, .23300, .23400, .23500,
     a .23700, .23900, .24050, .24200, .24400, .24600, .24600/
      data (sr(l),l=199,268)/  
     a .24600, .24700, .24800, .24950, .25100, .25250, .25400,
     a .25500, .25600, .25800, .26000, .26150, .26300, .26300,
     a .26300, .26300, .26300, .26450, .26600, .26800, .27000,
     a .27150, .27300, .27450, .27600, .27750, .27900, .27900,
     a .27900, .28050, .28200, .28300, .28400, .28400, .28400,
     a .28400, .28400, .28500, .28600, .28600, .28600, .28600,
     a .28600, .28800, .29000, .29100, .29200, .29200, .29200,
     a .29200, .29200, .29200, .29200, .29200, .29200, .29200,
     a .29200, .29200, .29200, .29350, .29500, .29500, .29500,
     a .29650, .29800, .29800, .29800, .29950, .30100, .30100/
      data (sr(l),l=269,338)/  
     a .30100, .30350, .30600, .30600, .30600, .30800, .31000,
     a .31000, .31000, .31200, .31400, .31400, .31400, .31400,
     a .31400, .31500, .31600, .31600, .31600, .31600, .31600,
     a .31600, .31600, .31600, .31600, .31750, .31900, .32000,
     a .32100, .32100, .32100, .32100, .32100, .32100, .32100,
     a .32250, .32400, .32400, .32400, .32400, .32400, .32400,
     a .32400, .32300, .32200, .32200, .32200, .32200, .32200,
     a .32100, .32000, .32000, .32000, .31850, .31700, .31700,
     a .31700, .31700, .31700, .31850, .32000, .32000, .32000,
     a .32000, .32000, .32150, .32300, .32300, .32300, .32500/
      data (sr(l),l=339,408)/
     a .32700, .32700, .32700, .32700, .32700, .32800, .32900,
     a .32900, .32900, .33050, .33200, .33200, .33200, .33200,
     a .33200, .33200, .33200, .33250, .33300, .33400, .33500,
     a .33650, .33800, .33950, .34100, .34300, .34500, .34650,
     a .34800, .35050, .35300, .35300, .35300, .35300, .35300,
     a .35300, .35300, .35400, .35500, .35500, .35500, .35500,
     a .35500, .35500, .35500, .35400, .35300, .35300, .35300,
     a .35300, .35300, .35300, .35300, .35300, .35300, .35300,
     a .35300, .35450, .35600, .35600, .35600, .35600, .35600,
     a .35600, .35600, .35800, .36000, .36000, .36000, .36000/
      data (sr(l),l=409,464)/  
     a .36000, .36150, .36300, .36300, .36300, .36300, .36300,
     a .36500, .36700, .36700, .36700, .36700, .36700, .36850,
     a .37000, .37000, .37000, .37000, .37000, .37000, .37000,
     a .37000, .37000, .37000, .37000, .36850, .36700, .36700,
     a .36700, .36550, .36400, .36300, .36200, .36200, .36200,
     a .36200, .36200, .36050, .35900, .35900, .35900, .35900,
     a .35900, .35900, .35900, .35900, .35900, .36100, .36300,
     a .36300, .36300, .36300, .36300, .36300, .36300, .36300/
        data (sr(l),l=465,534)/
     a .36300, .36400, .36500, .36700, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .37050,
     a .37200, .37200, .37200, .37200, .37200, .37200, .37200,
     a .37350, .37500, .37500, .37500, .37500, .37500, .37500,
     a .37500, .37500, .37500, .37500, .37500, .37500, .37500,
     a .37500, .37500, .37500, .37500, .37500, .37500, .37500,
     a .37500, .37500, .37500, .37500, .37500, .37500, .37500,
     a .37500, .37500, .37600, .37700, .37700, .37700, .37700,
     a .37700, .37700, .37700, .37700, .37700, .37700, .37700,
     a .37700, .37700, .37700, .37700, .37700, .37700, .37700/
        data (sr(l),l=535,604)/
     a .37700, .37700, .37700, .37700, .37700, .37800, .37900,
     a .37900, .37900, .37900, .37900, .37900, .37900, .38050,
     a .38200, .38200, .38200, .38200, .38200, .38200, .38200,
     a .38350, .38500, .38500, .38500, .38500, .38500, .38500,
     a .38500, .38500, .38500, .38500, .38500, .38650, .38800,
     a .38800, .38800, .38800, .38800, .38800, .38800, .38800,
     a .38800, .38800, .38800, .38800, .38800, .38950, .39100,
     a .39100, .39100, .39100, .39100, .39100, .39100, .39100,
     a .39100, .39100, .39100, .39100, .39100, .39100, .39100,
     a .39100, .39100, .39100, .39100, .39100, .39100, .39250/
        data (sr(l),l=605,674)/
     a .39400, .39400, .39400, .39400, .39400, .39400, .39400,
     a .39400, .39400, .39400, .39400, .39400, .39400, .39500,
     a .39600, .39600, .39600, .39600, .39600, .39600, .39600,
     a .39600, .39600, .39600, .39600, .39450, .39300, .39300,
     a .39300, .39300, .39300, .39300, .39300, .39100, .38900,
     a .38750, .38600, .38600, .38600, .38400, .38200, .38200,
     a .38200, .38000, .37800, .37800, .37800, .37600, .37400,
     a .37150, .36900, .36900, .36900, .36900, .36900, .37000,
     a .37100, .37100, .37100, .37100, .37100, .37100, .37100,
     a .37100, .37100, .37100, .37100, .37100, .37100, .37100/
        data (sr(l),l=675,744)/
     a .37100, .37100, .37100, .37100, .37100, .37250, .37400,
     a .37400, .37400, .37550, .37700, .37700, .37700, .37700,
     a .37700, .37700, .37700, .37800, .37900, .37900, .37900,
     a .37900, .37900, .37900, .37900, .38100, .38300, .38400,
     a .38500, .38650, .38800, .38800, .38800, .38800, .38800,
     a .39000, .39200, .39200, .39200, .39200, .39200, .39350,
     a .39500, .39500, .39500, .39500, .39500, .39500, .39500,
     a .39500, .39500, .39400, .39300, .39300, .39300, .39300,
     a .39300, .39050, .38800, .38800, .38800, .38800, .38800,
     a .38800, .38800, .38800, .38800, .38650, .38500, .38500/
        data (sr(l),l=745,1501)/
     a .38500, .38500, .38500, .38300, .38100, .38100, .38100,
     a .38100, .38100, .38100, .38100, .38100, .38100, .38100,
     a .38100, .38100, .38100, .38100, .38100, .37750, .37400,
     a .37400, .37400, .37400, .37400, .37400, .37400, .37400,
     a .37400, .37400, .37400, .37400, .37400, .37400, .37400,
     a .37300, .37200, .37050, .36900, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .36900,
     a .36900, .36900, .36900, .36900, .36900, .36900, .36900,
     a 680*0./
      do 1 i=1,1501
      r(i)=sr(i)
    1 continue
      return
      end

