c $Header: /weru/cvs/weps/weps.src/asd/test/tstasd.for,v 1.1.1.1 1999-03-12 17:05:17 wagner Exp $
c
c     program tstasd.for

*$noereference
      include 'p1werm.inc'
      include 'm1subr.inc'
      include 's1layr.inc'
      include 's1agg.inc'
      include 'manage/asd.inc'
*$reference

      integer sr,lay,i
      real    massf(msieve+1,mnsz)

      nsubr = 2
      nslay(1) = 2
      nslay(2) = 2


      call asdini()
c     -------------------------------------------------------------------------
c
      do 20 sr=1, nsubr
         do 10 lay=1, nslay(sr)
            aslagm(lay,sr) = 10*sr+5*lay
            as0ags(lay,sr) = 4 + lay
            aslagn(lay,sr) = 0.01 * sr
            aslagx(lay,sr) = 100.0 + sr
10       continue
20    continue

      print*,'initial data - lognormal case:', logcas
      do 300 sr=1, nsubr
         print*,'subregion',sr
         do 200 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
200      continue
300    continue

      do 500 sr=1, nsubr

         call asd2m(aslagn(1,sr), aslagx(1,sr), aslagm(1,sr),
     &   as0ags(1,sr), nslay, massf)

         print*,'----1------subregion',sr,' after asd2m --------------'
         do 210 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
210      continue

         call m2asd(massf, nslay,
     &   aslagn(1,sr), aslagx(1,sr), aslagm(1,sr), as0ags(1,sr))

         print*,'subregion',sr,' after m2asd'
         do 220 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
220      continue

         call asd2m(aslagn(1,sr), aslagx(1,sr), aslagm(1,sr),
     &   as0ags(1,sr), nslay, massf)

         print*,'----2------subregion',sr,' after asd2m --------------'
         do 211 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
211      continue

         call m2asd(massf, nslay,
     &   aslagn(1,sr), aslagx(1,sr), aslagm(1,sr), as0ags(1,sr))

         print*,'subregion',sr,' after m2asd'
         do 221 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
221      continue

         call asd2m(aslagn(1,sr), aslagx(1,sr), aslagm(1,sr),
     &   as0ags(1,sr), nslay, massf)

         print*,'----3------subregion',sr,' after asd2m --------------'
         do 212 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
212      continue

         call m2asd(massf, nslay,
     &   aslagn(1,sr), aslagx(1,sr), aslagm(1,sr), as0ags(1,sr))

         print*,'subregion',sr,' after m2asd'
         do 222 lay=1, nslay(sr)
            print*, aslagm(lay,sr), as0ags(lay,sr),
     &              aslagn(lay,sr), aslagx(lay,sr)
            print*, (massf(i,lay),i=1,nsieve+1)
222      continue

500   continue
      stop
      end
c
c$Log: not supported by cvs2svn $
c Revision 1.2  1995/09/13  15:50:02  wagner
c Necessary changes made to allow FORTRAN src files (*.for) to use the
c extended FORTRAN include statement rather than the MICROSOFT $INCLUDE
c directive as previously used.  This is required to allow use of other
c FORTRAN compilers.
c
c Changes have been made to the prologue.mk, epilogue.mk, and the Unix
c master startup.mk files as well as the src files.
c
c Revision 1.1.1.1  1995/01/18  04:19:57  wagner
c Initial checkin
c
c Revision 1.3  1992/10/10  21:44:14  wagner
c Changed names appropriate for submodel name change
c from TILLAGE to MANAGEMENT.
c
c Revision 1.2  1992/04/16  21:41:37  wagner
c Uses common memory now.
c