1. Bug fixed in background diffusion. 2.Uvs lake setup updated

setup/Uvs1979-2017_driver.dat

@@ -73,6 +73,8 @@ N_LatentFlux    -1
 N_Ustar         -1
 N_surfrad       -1
 N_cloud         -1
+N_NetRad        -1
+N_SurfTemp      -1
 #
 #-----------------------------------------------------------------------------------------
 #                              TIME INTEGRATION PARAMETERS
@@ -100,6 +102,8 @@ spinup_period 0.
 # 31536000 is year
 dt          20.
 call_Flake  0
+cp_period   0.
+control_point 0
 #
 #----------------------------------------------------------------------------------------
 #                                   PHYSICAL PARAMETERS

setup/Uvs1979-2017_setup.dat

@@ -149,6 +149,17 @@ zero_model 0
 thermokarst_meth_prod 0.
 soil_meth_prod 1.
 outflpar 0 #2 #0
+nManning 0.
+horvisc  0.
+pgrad    0.
+backdiff0 0.
+backdiff 0
+botfric  1
+#
+VmaxCH4aeroboxid 0.
+khsCH4           0.
+khsO2            0.
+r0methprod       0.
 #
 sensflux0 100.
 momflux0  1.e-15
@@ -163,14 +174,15 @@ deadvol 0.
 #----------------------------------------------------------------------------------------
 #
 T_profile 8
-0.    15.23 0. 0.0012  0.080   2.94E-1
-1.    15.33 0. 0.0012  0.080   2.93E-1
-3.    15.29 0. 0.0013  0.080   2.63E-1
-5.    14.92 0. 0.0016  0.082   2.63E-1
-7.    11.88 0. 0.0010  0.144   2.45E-1
-9.    9.92  0. 0.0023  0.171   1.0E-1
-11.   9.12  0. 0.1076  0.175   0.0E-0
-12.   8.96  0. 0.1513  0.173   0.0E-0
+0.    15.23 0. 0.0012  0.080   2.94E-1 0.
+1.    15.33 0. 0.0012  0.080   2.93E-1 0.
+3.    15.29 0. 0.0013  0.080   2.63E-1 0.
+5.    14.92 0. 0.0016  0.082   2.63E-1 0.
+7.    11.88 0. 0.0010  0.144   2.45E-1 0.
+9.    9.92  0. 0.0023  0.171   1.0E-1  0.
+11.   9.12  0. 0.1076  0.175   0.0E-0  0.
+12.   8.96  0. 0.1513  0.173   0.0E-0  0.
+T_soilprofile -1
 #.
 #----------------------------------------------------------------------------------------
 #                               BOUNDARY CONDITIONS: TRIBUTARIES AND EFFLUENTS
@@ -415,6 +427,7 @@ accum_begin 2003060100
 accum_end   2004060100
 rtemp 1
 -999. -999. -999.
+zserout -999.
 #
 ngrid_out 12
 0.2
@@ -442,6 +455,7 @@ ngridsoil_out 11
 8.
 9.
 10.
+ngridice_out -1
 #
 #----------------------------------------------------------------------------------------
 # NOTE: VARIABLES, DENOTED BY ASTERISK *, ARE USED ONLY IN STANDALONE RUNS OF THE MODEL

source/model/phys_func.f90

@@ -1975,7 +1975,7 @@ END FUNCTION FP_THETA
     else
       x = max(x, xmin)
       y = area*1.E-6 ! converting from m**2 to km**2
-      wst%lamw_back(i) = C1*(y)**C2 * x**C3
+      wst%lamw_back(i) = z*(y)**C2 * x**C3
       wst%lamw_back(i) = wst%lamw_back(i) * 1.E-4 !converting from cm**2/s to m**2/s
       wst%lamw_back(i) = wst%lamw_back(i) * cw_m_row0
       wst%lamw_back(i) = wst%lamw_back(i) * 5.E-1 !Calibration multiplyer