_RL cosz(sNx,sNy)

      call ASTRO ( 20040321,nhms, xlats,xlons, im2*jm2, cosz,ra )

      call ASTRO ( nymd,nhms, xlats,xlons, im2*jm2, cosz,ra )

       radswt(i,j,bi,bj) = S0*(1.0/ra**2)*cosz(i,j)

      _RL getcon, cosz(irun), alat(irun), alon(irun), ra

      _RL  cosz(m,ndim),fcld(m,ndim,np),taucld(m,ndim,np,2)

           fm=cosz(i,j)/dm

     .             albirdr,albirdf,albuvdr,albuvdf,cosz,

       if( cosz(i).gt.0.0 ) then

     *                 taual,rsirbm,rsirdf,rsuvbm,rsuvdf,cosz,

     *     rsuvbm(m,ndim),rsuvdf(m,ndim),cosz(m,ndim),co2

         csm(i,j)=35./sqrt(1224.*cosz(i,j)*cosz(i,j)+1.)

      call CLDSCALE(m,n,ndim,np,cosz,fcld,taucld,ict,icb,

      subroutine CLDSCALE (m,n,ndim,np,cosz,fcld,taucld,ict,icb,

      _RL lats(istrip),lons(istrip),cosz(istrip),icest(istrip)

       call ASTRO ( 20040321,nhms,lats,lons,istrip,cosz,ra )

       call ASTRO ( nymd,nhms,lats,lons,istrip,cosz,ra )

       cosz(i) = max(cosz(i),0.0001 _d 0)

     I   tmpnlay, sh(1,nlay), cd, cosz, pardr, pardf,

        call ASTRO(nymd,nhms,lats,lons,im2*jm2,cosz,radius)

        call GETALB(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp,

      _RL lats(sNx,sNy), lons(sNx,sNy), cosz(sNx,sNy)

      subroutine GETALB(sec,month,day,cosz,snodep,fraci,fracg,im,jm,

      _RL cosz(im,jm),fraci(im,jm),fracg(im,jm)

       ALBOC(I,J) = A0 + (A1 + (A2 +  A3*cosz(I,J))*cosz(I,J))*cosz(I,J)

      call GRD2MSC(cosz,im,jm,igrd,zenith,nchp,nchpland)

        call ASTRO(20040321,nhms,lats,lons,im2*jm2,cosz,radius)

      _RL cosz(sNx,sNy)

      _RL    cosz(irun), alat(irun), alon(irun), ra

       _RL  cosz, frac, fluxi

          cosz = ( sin(delta)*sin(lat)+              !average zenith angle

          IF (cosz.LE.5. _d -3) cosz= 5. _d -3

          fluxi = solar*(1. _d 0-albedo)*cosz*frac*parfrac