changes

CLM集合卡曼滤波数据同化算法

1, Water.f90

          if(tss(i) > tfrz)then
             smpmax = -1.5e5        !-1.e10
                   s_node = max(wliq(i)/(dzmm(i)*eff_porosity(i)),0.01)
                                --------------------------------- previously use "vol_liq"
             smp_node = max(smpmax, -phi0*s_node**(-bsw))
             rresis(i)=(1.-smp_node/smpmax)/(1.+phi0/smpmax)
             rootr(i) = rootfr(i)*rresis(i)
             roota    = roota + rootr(i)