river_allocate Subroutine

  subroutine river_allocate(mb, mch, mhyd)
    integer, intent(in) :: mb, mch, mhyd
    allocate(c1_qd(mb))
    allocate(c1_ssf(mb))
    allocate(c2_qd(mb))
    allocate(c2_ssf(mb))
    allocate(c3_qd(mb))
    allocate(c3_ssf(mb))
    allocate(c4_qd(mb))
    allocate(c4_ssf(mb))
    allocate(chd(mch))
    allocate(chk(2, mb))
    allocate(chl(2, mb))
    allocate(chn(mb))
    allocate(chnn(mch))
    allocate(chs(mb))
    allocate(chss(mch))
    allocate(chw(2, mb))
    allocate(ovn(mb))
    allocate(pet_day(mb)) ! for transmission losses
    allocate(phi(20, mb))
    allocate(qdilast(mhyd))
    allocate(qdinp(mhyd, 366))
    allocate(qdolast(mhyd))
    allocate(qdout(mhyd, 366))
    allocate(qsilast(mhyd))
    allocate(qsolast(mhyd))
    allocate(qssinp(mhyd, 366))
    allocate(qssout(mhyd, 366))
    allocate(roc2(mb))
    allocate(roc4(mb))
    allocate(rt(mb))
    allocate(sdtsav(mb))
    allocate(srch(20, mch))
    allocate(varoute(19, mhyd))
    allocate(xkm_qd(mb))
    allocate(xkm_ssf(mb))
    allocate(runsub_m3s(mb))

    runsub_m3s = 0.
    chd = 0.
    chk = 0.37
    chl = 0.
    chn = 0.075
    chnn = 0.075 * chnnc0
    chs = 0.
    chss = 0.
    chw = 0.
    ovn = 0.15
    pet_day = 0.
    phi = 0.
    qdilast = 0.
    qdinp = 0.
    qdolast = 0.
    qdout = 0.
    qsilast = 0.
    qsolast = 0.
    qssinp = 0.
    qssout = 0.
    roc2 = roc2_0
    roc4 = roc4_0
    rt = 0.
    sdtsav = 0.
    srch = 0.
    varoute = 0.
  end subroutine river_allocate