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


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

c  
c  Wiscombe, W.J. and and S.G. Warren, 1980: "A model for the spectral
c  albedo of snow. I: pure snow."  J. Atmospheric Sciences, 37, 2712-2733.
c  (based on there Mie/delta-eddington model with Reff=125um )

      data (sr(i),i=1,110) /
     a.970,.970,.970,.970,.970,.970,.970,.970,.970,.970,.970,
     a.971,.972,.973,.974,.975,.975,.975,.975,.976,.977,.978,
     a.979,.980,.980,.980,.980,.980,.980,.980,.980,.981,.982,
     a.983,.984,.985,.985,.985,.985,.986,.987,.988,.989,.990,
     a.990,.990,.990,.990,.990,.990,.990,.990,.990,.990,.990,
     a.988,.987,.986,.985,.983,.982,.981,.980,.978,.977,.976,
     a.975,.973,.972,.971,.970,.968,.967,.966,.965,.965,.965,
     a.965,.965,.963,.962,.961,.960,.953,.952,.951,.950,.950,
     a.950,.950,.950,.948,.945,.943,.940,.938,.935,.933,.930,
     a.928,.925,.923,.920,.918,.915,.913,.910,.908,.907,.906/
      data (sr(i),i=111,220) /
     a.905,.903,.902,.901,.900,.898,.895,.893,.890,.885,.880,
     a.875,.870,.865,.860,.855,.850,.843,.842,.841,.840,.835,
     a.830,.825,.820,.810,.800,.790,.780,.771,.762,.753,.745,
     a.739,.733,.727,.720,.718,.715,.713,.710,.715,.720,.725,
     a.730,.733,.735,.738,.740,.742,.743,.744,.745,.745,.745,
     a.745,.745,.743,.742,.741,.740,.732,.724,.716,.710,.692,
     a.674,.656,.640,.630,.620,.610,.600,.584,.568,.552,.535,
     a.524,.513,.502,.490,.490,.490,.490,.490,.490,.490,.490,
     a.490,.492,.493,.494,.495,.495,.495,.495,.495,.495,.495,
     a.495,.495,.493,.492,.491,.490,.488,.485,.483,.480,.470/
      data (sr(i),i=221,330) /
     a.460,.450,.440,.424,.408,.392,.375,.320,.264,.198,.150,
     a.134,.118,.101,.085,.074,.063,.051,.040,.036,.032,.028,
     a.025,.025,.025,.025,.025,.028,.030,.033,.035,.039,.043,
     a.047,.050,.053,.055,.058,.060,.063,.065,.068,.070,.073,
     a.075,.078,.080,.083,.085,.088,.090,.093,.095,.098,.100,
     a.104,.108,.112,.115,.119,.123,.126,.130,.135,.140,.145,
     a.150,.153,.155,.158,.160,.163,.165,.168,.170,.173,.175,
     a.178,.180,.184,.188,.192,.195,.204,.213,.222,.230,.238,
     a.246,.254,.260,.254,.248,.242,.235,.206,.177,.248,.120,
     a.104,.088,.072,.055,.049,.043,.037,.030,.038,.035,.033/
      data (sr(i),i=331,mxwv) /
     a.030,.028,.025,.023,.020,.023,.022,.021,.020,.018,.017,
     a.016,.015,.013,.012,.011,.010,.010,.010,.010,.010,.010,
     a.010,.010,.010,.012,.013,.014,.015,.017,.018,.019,.020,
     a.022,.023,.024,.025,.025,.025,.025,.025,.027,.030,.032,
     a.035,.040,.045,.050,.055,.059,.063,.067,.070,.074,.078,
     a.082,.085,.089,.093,.097,.100,.102,.103,.104,.105,.109,
     a.113,.117,.120,.124,.128,.132,.135,.130,.125,.120,.115,
     a.108,.100,.092,.085,.081,.077,.073,.070,.068,.065,.063,
     a.062,.058,.055,.053,.050,.048,.047,.046,.045,.042,.040,
     a.038,.035,.033,.032,.031,.030,.030,.030,.030,.030,.028,
     a.027,.026,.025,.025,.025,.025,.025,.025,.025,.025,.025,
     a.025,.025,.025,.025,.025,.025,.025,.025,.025,.025,.025,
     a.025,.025,.025,.025,.025,.023,.022,.021,.020,.018,.017,
     a.016,.015,.013,.012,.011,.010,.008,.007,.006,.005,.003,
     a.002,.001,.000,264*0.0/

      do 1 i=1,mxwv
      r(i)=sr(i)
1     continue
      return
      end