nutrient_allocate Subroutine

  subroutine nutrient_allocate(mb, mcrdb, meap, ml)
    integer, intent(in) :: mb, mcrdb, meap, ml
    allocate(ano3(mb, meap, ml))
    allocate(anora(mb, meap, ml))
    allocate(anors(mb, meap, ml))
    allocate(bn1(mcrdb))
    allocate(bn3(mcrdb))
    allocate(bnu1(mcrdb))
    allocate(bnu2(mcrdb))
    allocate(bp1(mcrdb))
    allocate(bp2(mcrdb))
    allocate(bn2(mcrdb))
    allocate(bp3(mcrdb))
    allocate(bpu1(mcrdb))
    allocate(bpu2(mcrdb))
    allocate(dflow(mb, meap, 20))
    allocate(plab(mb, meap, ml))
    allocate(pma(mb, meap, ml))
    allocate(pms(mb, meap, ml))
    allocate(porg(mb, meap, ml))
    ano3 = 0.
    anora = 0.
    anors = 0.
    bn1 = 0.
    bn3 = 0.
    bnu1 = 0.
    bnu2 = 0.
    bp1 = 0.
    bp3 = 0.
    bpu1 = 0.
    bpu2 = 0.
    bn2 = 0.
    bp2 = 0.
    dflow = 0.
    plab = 0.
    pma = 0.
    pms = 0.
    porg = 0.
  end subroutine nutrient_allocate