ONE OF THE MOST IMPORTANT SCIENTIFIC CONTRIBUTION OF THE FERMI SATELLITE HAS BEEN TO SHOW THAT THE GAMMA-RAY EMISSION FROM SEVERAL SUPERNOVA REMNANTS (SNRS) IS OF HADRONIC ORIGIN THEREBY SUPPORTING THE IDEA THAT SNRS ARE THE SOURCES OF GALACTIC COSMIC RAYS (CRS). ALMOST ALL GEV--BRIGHT SNRS EXPAND INTO PARTIALLY-NEUTRAL MEDIA AND/OR ARE ASSOCIATED WITH DENSE MOLECULAR CLOUDS (MCS). IN THE PRESENCE OF NEUTRAL HYDROGEN THE SHOCK DYNAMICS THE SPECTRUM OF ACCELERATED PARTICLES AND THE CR--MAGNETIC FIELD COUPLING ARE RADICALLY DIFFERENT THAN IN FULLY-IONIZED PLASMAS. HOWEVER THE EFFECTS OF NEUTRALS ARE USUALLY NEGLECTED IN CALCULATIONS OF CR ACCELERATION AND PROPAGATION. WE WILL USE SELF-CONSISTENT KINETIC SIMULATIONS OF COLLISIONLESS SHOCKS PROPAGATING IN PARTIALLY-IONIZED PLASMAS TO ASSESS THE IMPACT OF PHENOMENA SUCH AS CHARGE EXCHANGE IONIZATION AND ION--NEUTRAL DAMPING ON SHOCK DYNAMICS CR ACCELERATION AND MAGNETIC FIELD AMPLIFICATION WHICH ARE NEEDED TO CHARACTERIZE THE HADRONIC GAMMA-RAY EMISSION OBSERVED IN SNR/MC SYSTEMS.