noah  1.6.0
 All Data Structures Files Functions
set_soilveg_snippet.f90
Go to the documentation of this file.
1 
4 
7 
8 module set_soilveg_snippet_mod
9 
10  implicit none
11 
12  private
13 
14  public set_soilveg
15 
16 contains
17 
26 subroutine set_soilveg(isot,ivet, maxsmc, bb, satpsi, iret)
27  implicit none
28 
29  integer, intent(in) :: isot,ivet
30  real, dimension(30), intent(out) :: maxsmc, bb, satpsi
31  integer, intent(out) :: iret
32 
33 ! set vegetation-dependent params (May 2021, UFS uses ivet=1)
34 ! Draper, not needed for now, but might need SNUPX
35 ! for SWE-> SCF calculation later
36 !
37 ! if(ivet.eq.1)then
38 
39  !defined_veg=20
40 ! might want this later
41  ! SNUPX =(/0.080, 0.080, 0.080, 0.080, 0.080, 0.020,
42  !* 0.020, 0.060, 0.040, 0.020, 0.010, 0.020,
43  !* 0.020, 0.020, 0.013, 0.013, 0.010, 0.020,
44  !& 0.020, 0.020, 0.000, 0.000, 0.000, 0.000,
45  !& 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
46 
47 ! endif
48 
49 ! set soil-dependent params (May 2021, UFS uses isot=1)
50 
51  if (isot .eq. 1) then
52 
53 ! using stasgo table
54  bb =(/4.05, 4.26, 4.74, 5.33, 5.33, 5.25, &
55  6.77, 8.72, 8.17, 10.73, 10.39, 11.55,&
56  5.25, 4.26, 4.05, 4.26, 11.55, 4.05, &
57  4.05, 0.00, 0.00, 0.00, 0.00, 0.00, &
58  0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
59 ! Draper, these are provided for reference only, and
60 ! may be useful for later SMC updates
61 ! DRYSMC=(/0.010, 0.025, 0.010, 0.010, 0.010, 0.010,
62 ! & 0.010, 0.010, 0.010, 0.010, 0.010, 0.010,
63 ! & 0.010, 0.010, 0.010, 0.010, 0.010, 0.010,
64 ! & 0.010, 0.000, 0.000, 0.000, 0.000, 0.000,
65 ! & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
66 
67  maxsmc=(/0.395, 0.421, 0.434, 0.476, 0.476, 0.439, &
68  0.404, 0.464, 0.465, 0.406, 0.468, 0.457, &
69  0.464, 0.421, 0.200, 0.421, 0.457, 0.200, &
70  0.395, 0.000, 0.000, 0.000, 0.000, 0.000, &
71  0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
72 
73  satpsi=(/0.035, 0.0363, 0.1413, 0.7586, 0.7586, 0.3548, &
74  0.1349, 0.6166, 0.2630, 0.0977, 0.3236, 0.4677,&
75  0.3548, 0.0363, 0.0350, 0.0363, 0.4677, 0.0350,&
76  0.0350, 0.00, 0.00, 0.00, 0.00, 0.00, &
77  0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
78 
79  !defined_soil=19
80  else
81  print *, 'set_soilveg_snippet not coded for soil type ', isot
82  iret = -1
83  return
84  endif
85 
86  iret = 0
87 
88 end subroutine set_soilveg
89 
90 end module set_soilveg_snippet_mod
set_soilveg_snippet_mod
Below was extracted from namelist_soilveg.f and set_soilveg.f (couldn't get above to compile for doxy...
Definition: set_soilveg_snippet.f90:8
set_soilveg_snippet_mod::set_soilveg
subroutine, public set_soilveg(isot, ivet, maxsmc, bb, satpsi, iret)
This subroutine initializes soil and vegetation parameters needed in global_cycle/land_increment.f90.
Definition: set_soilveg_snippet.f90:26