!
! This is a subroutine to save the model statement of cur
      subroutine save_statement()
      include 'p1werm.inc'
      include 's1phys.inc'
      include 's1agg.inc'
      include 's1dbh.inc'
      include 's1dbc.inc'
      include 's1sgeo.inc'
      include 'h1hydro.inc'
      include 'm1subr.inc'
      include 's1surf.inc'
      include 's1psd.inc'
      include 's1layr.inc' 
      include 'h1scs.inc'
      include 'h1db1.inc'
      include 'command.inc'
      include 'w1clig.inc'
      include 'h1temp.inc' 
      include 'h1balance.inc'

      include 'airpact/prevdayMan.inc'
      include 'soil/prevday.inc'	  
      include 'hydro/htheta.inc'
      include 'airpact/prevdaySoil.inc'
      include 'airpact/prevdayHydro.inc'
      include 'airpact/prevdayClim.inc'
      include 'airpact/emit.inc'

! local variable
      integer isr
      integer lay
! update soil hydric and soil texture info 
! after call propsaxt()
      do isr =1, nsubr 

      do lay=1, nslay(isr)

!    layer thickness
      aszlyt0(lay,isr) = aszlyt(lay,isr)
! Surface physical properties      
      asfsan0(lay,isr) = asfsan(lay,isr)
      asfcla0(lay,isr) = asfcla(lay,isr)
      asfvfs0(lay,isr) = asfvfs(lay,isr)
      asfsil0(lay,isr) = asfsil(lay,isr)
      asvroc0(lay,isr) = asvroc(lay,isr)
      asfvcs0(lay,isr) = asfvcs(lay,isr)
      asfcs0(lay,isr) = asfcs(lay, isr)
      asfms0(lay,isr) = asfms(lay,isr)
      asffs0(lay, isr) = asffs(lay, isr)
      asdwblk0(lay,isr) = asdwblk(lay,isr)
     
! soil chemical properties
      asfom0(lay,isr) = asfom(lay,isr) !soil organic matter content
      aseags0(lay,isr) = aseags(lay,isr)
      as0ph0(lay, isr) = as0ph(lay,isr)
      asfcce0(lay,isr) = asfcce(lay,isr)
      asfcec0(lay, isr) = asfcec(lay,isr)
      asfcle0(lay,isr) = asfcle(lay,isr)
! aggregate properties
      aslagm0(lay,isr) = aslagm(lay,isr)
      aslagn0(lay,isr) = aslagn(lay,isr)
      aslagx0(lay,isr) = aslagx(lay,isr)
      as0ags0(lay,isr) = as0ags(lay,isr)
      asdagd0(lay, isr) = asdagd(lay,isr)
      aseags0(lay,isr) = aseags(lay,isr)
! soil hydrologic properties
      asdblk0(lay, isr) = asdblk(lay,isr) !already defined
!      ahrwc0(lay,isr) = ahrwc(lay,isr)
      ah0cb0(lay, isr) = ah0cb(lay,isr) 
      aheaep0(lay, isr) = aheaep(lay,isr) 
      ahrsk0(lay,isr) = ahrsk(lay,isr)

      ahrwcs0(lay,isr) = ahrwcs(lay,isr)
      ahrwcw0(lay,isr) = ahrwcw(lay,isr)
      ahrwcf0(lay,isr) = ahrwcf(lay,isr)
      ahrwca0(lay,isr) = ahrwca(lay,isr)
   
      asdsblk0(lay,isr) = asdsblk(lay,isr)
      asdpart0(lay,isr) = asdpart(lay,isr)
      ahfredsat0(lay,isr) = ahfredsat(lay,isr)
      
! mean daily soil temperature
      ahtsav0(lay,isr) = ahtsav(lay,isr)
! max daily soil temperature
      ahtsmx0(lay, isr) = ahtsmx(lay,isr)
! min daily soil temperature
      ahtsmn0(lay,isr) = ahtsmn(lay,isr)
 
! fraction of soil water in layer which is frozen
      ahfice0(lay,isr) = ahfice(lay,isr)  

      end do
! snow depth ?
      ahzsnd0(isr) = ahzsnd(isr)
!water infiltration depth (mm) 
      ahzwid0(isr) = ahzwid(isr)
! accumlated surface evaporation since last complete rewetting
      ahzeasurf0(isr) = ahzeasurf(isr)
      ahzsno0(isr) = ahzsno(isr)
      ahtsno0(isr) = ahtsno(isr)
      ahfsnfrz0(isr) = ahfsnfrz(isr) 
      
      presswc0(isr) = presswc(isr)
      pressnow0(isr) = pressnow(isr)
!     presday(isr) = initday(isr) 	 

! hydrologic balance variables	
      cumprecip0(isr) = cumprecip(isr)
      cumrunoff0(isr) = cumrunoff(isr)
      cumevap0(isr) = cumevap0(isr)
      cumtrans0(isr) = cumtrans(isr)
      cumdrain0(isr) = cumdrain(isr)
      hprevrotation0(isr) = hprevrotation(isr)

! irrigation variables
      ahzirr0(isr) = ahzirr(isr)
      ahratirr0(isr) = ahratirr(isr)
      ahdurirr0(isr) = ahdurirr(isr) 
      ahlocirr0(isr) = ahlocirr(isr)
      am0monirr0(isr) = am0monirr(isr)  
      ahmadirr0(isr) = ahmadirr(isr)
      ahminirr0(isr) = ahminirr(isr)
      ahndayirr0(isr) = ahndayirr(isr)
      ahmintirr0(isr) = ahmintirr(isr)
      
! daily relative humidity
      awrrh0 = awrrh
      ahzrun0(isr) = ahzrun(isr)
      ahzsmt0(isr) = ahzsmt(isr)
      ah0cng0(isr) = ah0cng(isr)
      ah0cnp0(isr) = ah0cnp(isr)

!biomass
      
! save Bulk Density and Praticle Density with texture based calculation
!  asdpart0(1,isr) = asdpart(1,isr)	       
! soil surface variable
! crust properties
      asfcr0(isr) = asfcr(isr)
      aszcr0(isr) = aszcr(isr)
      asflos0(isr) = asflos(isr)
      asmlos0(isr) = asmlos(isr)
      asdcr0(isr) = asdcr(isr)
      asecr0(isr) = asecr(isr)
! surface roughness properties
      aslrr0(isr) = aslrr(isr)
      aslrro0(isr) = aslrro(isr)
      asargo0(isr) = asargo(isr)
      aszrgh0(isr) = aszrgh(isr)
      asxrgs0(isr) = asxrgs(isr)
      asxrgw0(isr) = asxrgw(isr)
      as0rrk0(isr) = as0rrk(isr)
      aslrrc0(isr) = aslrrc(isr)
    
      asxdks0(isr) = asxdks(isr)

      asxdkh0(isr) = asxdkh(isr)
      asf10an0(isr) = asf10an(isr)
      asf10en0(isr) = asf10en(isr)
      asf10bk0(isr) = asf10bk(isr)	 	  
        
      asfald0(isr) = asfald(isr)
      asfalw0(isr) = asfalw(isr)
      amrslp0(isr) = amrslp(isr)
      SFCov0(isr) = SFCov(isr)
      bedrock_depth0(isr) = bedrock_depth(isr)
      restrict_depth0(isr) = restrict_depth(isr)
      
! save soil emit as previous day emit
      aegtp(isr) = aegt(isr)
      aegtssp(isr) = aegtss(isr) 
      aegt10p(isr) = aegt10(isr)  
      end do 
      end