Inability to construct noisy version of a parametric gate

Usually I can over-ride some standard gate (e.g. X) with some Kraus operators using define_noisy_gate. The same doesn't seem to work with parametric gates such as RY. If I try to run the following:

from pyquil import get_qc, Program
from pyquil.gates import RY, MEASURE
import numpy as np


# bit flip kraus ops
prob = 0.1
kraus_ops = [np.sqrt(1 - prob) * np.array([[1, 0], [0, 1]]),
             np.sqrt(prob) * np.array([[0, 1], [1, 0]])]

p = Program(RY(np.pi/2, 0))
p.define_noisy_gate('RY', [0], kraus_ops)

ro = p.declare('ro', 'BIT', 1)
p += MEASURE(0, ro[0])
p.wrap_in_numshots_loop(1000)

qc = get_qc('1q-qvm')
results = qc.run(p)

I get the following error msg

QVMError: The value of QVM::GATE is #<CL-QUIL:PARAMETERIZED-GATE {100354A5C3}>, which is not of type CL-QUIL:STATIC-GATE.

The QVM returned the above error. This could be due to a bug in the server or a
bug in your code. If you suspect this to be a bug in pyQuil or Rigetti Forest,
then please describe the problem in a GitHub issue at:
    https://github.com/rigetti/pyquil/issues

Is there any way to construct a noisy version of a parametric gate?

Reproduce the failure with Program.define_noisy_gate, RY, and the 1q-qvm example in the issue. Trace how the parameterized RY instruction reaches the QVM and determine whether noisy construction is supported; done means the behavior is fixed or the supported limitation is clearly established and covered by a regression test.