c=============================================================================


      subroutine   sand (r) 2
      parameter (mxwv=751)
      dimension sr(mxwv),r(mxwv)

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
c     Staetter, R., and M. Schroeder, 1978: "Spectral characteristics of
c     natural surfaces", Proceeding of the tenth Int. Conf. on Earth
c     Obs. from Space, 6-11 March 1978, (ESA-SP, 134)


      data (sr(i),i=1,140) / 30*0.,
     a.091,.091,.091,.091,.091,.091,.095,.095,.095,.095,.095,
     a.095,.095,.097,.097,.100,.103,.103,.107,.107,.107,.107,
     a.110,.114,.117,.121,.121,.125,.125,.125,.128,.126,.131,
     a.134,.134,.134,.134,.134,.134,.141,.145,.149,.149,.154,
     a.157,.160,.163,.163,.166,.169,.173,.177,.181,.185,.189,
     a.195,.199,.204,.204,.208,.214,.220,.224,.228,.233,.235,
     a.239,.242,.246,.246,.248,.251,.254,.256,.260,.263,.263,
     a.263,.266,.270,.273,.276,.279,.279,.282,.284,.284,.284,
     a.286,.286,.286,.290,.292,.292,.292,.292,.292,.292,.292,
     a.295,.295,.298,.298,.301,.301,.306,.306,.310,.310,.314/
      data (sr(i),i=141,240) /
     a.314,.314,.316,.316,.316,.316,.316,.319,.321,.321,
     a.321,.321,.324,.324,.324,.324,.322,.322,.322,.320,
     a.320,.317,.317,.317,.320,.320,.320,.323,.323,.327,
     a.327,.327,.329,.329,.332,.332,.332,.332,.333,.335,
     a.338,.341,.345,.348,.353,.353,.353,.353,.355,.355,
     a.355,.355,.353,.353,.353,.353,.353,.353,.356,.356,
     a.356,.356,.360,.360,.360,.363,.363,.363,.367,.367,
     a.367,.370,.370,.370,.370,.370,.370,.367,.367,.364,
     a.362,.362,.362,.359,.359,.359,.359,.359,.363,.363,
     a.363,.363,.363,.365,.369,.369,.369,.369,.369,.372/
      data (sr(i),i=241,340) /
     a.372,.372,.372,.375,.375,.375,.375,.375,.375,.375,
     a.375,.375,.375,.375,.375,.375,.375,.375,.377,.377,
     a.377,.377,.377,.377,.377,.377,.377,.377,.377,.377,
     a.379,.379,.379,.379,.382,.382,.382,.382,.385,.385,
     a.385,.385,.385,.385,.388,.388,.388,.388,.388,.388,
     a.388,.391,.391,.391,.391,.391,.391,.391,.391,.391,
     a.391,.391,.394,.394,.394,.394,.394,.394,.394,.396,
     a.396,.396,.396,.396,.396,.393,.393,.393,.393,.389,
     a.386,.386,.382,.382,.378,.378,.374,.369,.369,.369,
     a.371,.371,.371,.371,.371,.371,.371,.371,.371,.371/
      data (sr(i),i=341,mxwv) /
     a.374,.374,.377,.377,.377,.377,.379,.379,.379,.379,
     a.383,.385,.388,.388,.388,.392,.392,.392,.395,.395,
     a.395,.395,.395,.393,.393,.393,.388,.388,.388,.388,
     a.388,.385,.385,.385,.381,.381,.381,.381,.381,.381,
     a.381,.381,.374,.374,.374,.374,.374,.374,.374,.374,
     a.372,.369,.369,.369,.369,.369,.369,.369,.369,.369,
     a.369,.369,.369,.369,.369,.369,.369,.369,.369,.369,
     a.369,340*0./
      do 1 i=1,mxwv
      r(i)=sr(i)
    1 continue
      return
      end