Merge pull request swat-model#171 from tugraskan/change_snow

Snow change (#149)

Author: odav

src/actions.f90

@@ -4,7 +4,7 @@ subroutine actions (ob_cur, ob_num, idtbl)
       use time_module
       use aquifer_module
       use hru_module, only : hru, cn2, fertno3, fertnh3, fertorgn, fertorgp, fertsolp,   &
-        ihru, ipl, isol,  phubase, sol_sumno3, sol_sumsolp, qtile
+        ihru, ipl, isol,  phubase, sol_sumno3, sol_sumsolp, qtile, snodb
       use soil_module
       use plant_module
       use plant_data_module
@@ -91,6 +91,7 @@ subroutine actions (ob_cur, ob_num, idtbl)
       real :: stor_m3 = 0.
       character(len=1) :: action = ""      !         |
       character(len=40) :: lu_prev = ""    !         |
+      character(len=40) :: snow_prev = ""  !         |
 
       do iac = 1, d_tbl%acts
         action = "n"
@@ -944,6 +945,22 @@ subroutine actions (ob_cur, ob_num, idtbl)
               end do
               !pcom(j)%dtbl(idtbl)%num_actions(iac) = pcom(j)%dtbl(idtbl)%num_actions(iac) + 1
             !end if
+            
+          !snow change - change snow parameters for HRU
+          case ("snow_change")
+            j = d_tbl%act(iac)%ob_num
+            if (j == 0) j = ob_cur
+            d_tbl%snow_chg_mx(iac) = d_tbl%snow_chg_mx(iac) + 1
+            ilu = d_tbl%act_typ(iac)
+            hru(j)%dbs%snow = ilu
+            
+            snow_prev = hru(j)%dbsc%snow
+            hru(j)%dbsc%snow = d_tbl%act(iac)%file_pointer
+            hru(j)%sno = snodb(ilu)
+            !! write to lu_change output file
+            write (3612,*) j, time%yrc, time%mo, time%day_mo,  " SNOW_CHANGE ",        &
+                    snow_prev, hru(j)%dbsc%snow, "   0   0"
+                            
           !land use change - contouring
           case ("p_factor")
             j = d_tbl%act(iac)%ob_num

src/conditional_module.f90

@@ -34,6 +34,7 @@ module conditional_module
         character(len=25), dimension(:,:), allocatable :: alt           ! condition alternatives
         type (actions_var), dimension(:), allocatable :: act            ! actions
         integer, dimension(:), allocatable :: lu_chg_mx                 ! max times lu change can occur
+        integer, dimension(:), allocatable :: snow_chg_mx               ! max times snow change can occur
         character(len=1), dimension(:,:), allocatable :: act_outcomes   ! action outcomes ("y" to perform action; "n" to not perform)
         character(len=1), dimension(:), allocatable :: act_hit          ! "y" if all condition alternatives (rules) are met; "n" if not
         integer, dimension(:), allocatable :: act_typ                   ! pointer to action type (ie plant, fert type, tillage implement, release type, etc)

src/dtbl_lum_read.f90

@@ -12,7 +12,7 @@ subroutine dtbl_lum_read
       use plant_data_module
       use constituent_mass_module
       use hydrograph_module, only : sp_ob
-      use hru_module, only : hru
+      use hru_module, only : hru, snodb
 
       implicit none
 
@@ -59,6 +59,8 @@ subroutine dtbl_lum_read
             allocate (dtbl_lum(i)%act_hit(dtbl_lum(i)%alts))
             allocate (dtbl_lum(i)%act_typ(dtbl_lum(i)%acts), source = 0)
             allocate (dtbl_lum(i)%act_app(dtbl_lum(i)%acts), source = 0)
+            allocate (dtbl_lum(i)%lu_chg_mx(dtbl_lum(i)%acts), source = 0)
+            allocate (dtbl_lum(i)%snow_chg_mx(dtbl_lum(i)%acts), source = 0)
             allocate (dtbl_lum(i)%act_outcomes(dtbl_lum(i)%acts,dtbl_lum(i)%alts))
 
             !read conditions and condition alternatives
@@ -230,6 +232,22 @@ subroutine dtbl_lum_read
                       exit
                     end if
                   end do
+                  
+                case ("lu_change")
+                  do idb = 1, db_mx%landuse
+                    if (dtbl_lum(i)%act(iac)%file_pointer == lum(idb)%name) then
+                      dtbl_lum(i)%act_typ(iac) = idb
+                      exit
+                    end if
+                  end do
+                  
+                case ("snow_change")
+                  do idb = 1, db_mx%sno
+                    if (dtbl_lum(i)%act(iac)%file_pointer == snodb(idb)%name) then
+                      dtbl_lum(i)%act_typ(iac) = idb
+                      exit
+                    end if
+                  end do
                 end select
 
                 !xwalk conditions and actions for days since last action

src/dtbl_scen_read.f90

@@ -8,6 +8,7 @@ subroutine dtbl_scen_read
       use fertilizer_data_module
       use input_file_module
       use conditional_module
+      use hru_module, only : snodb
 
       implicit none
 
@@ -52,6 +53,7 @@ subroutine dtbl_scen_read
             allocate (dtbl_scen(i)%act_typ(dtbl_scen(i)%acts), source = 0)
             allocate (dtbl_scen(i)%act_app(dtbl_scen(i)%acts), source = 0)
             allocate (dtbl_scen(i)%lu_chg_mx(dtbl_scen(i)%acts), source = 0)
+            allocate (dtbl_scen(i)%snow_chg_mx(dtbl_scen(i)%acts), source = 0)
 
             allocate (dtbl_scen(i)%act_outcomes(dtbl_scen(i)%acts,dtbl_scen(i)%alts))
 
@@ -83,6 +85,13 @@ subroutine dtbl_scen_read
                       exit
                     end if
                   end do
+                case ("snow_change")
+                  do ilum = 1, db_mx%sno
+                    if (dtbl_scen(i)%act(iac)%file_pointer == snodb(ilum)%name) then
+                      dtbl_scen(i)%act_typ(iac) = ilum
+                      exit
+                    end if
+                  end do
                 end select
 
             end do