From allocate_and_insert_fields a radiation specific internal subroutine has been added that allocates and inserts radiation specific scalar fields  from templates into states(1). This took some inspiration from the  sediments model code above it in that all the fields are actually  templates in the flml and and names change as needed.
From allocate_and_insert_fields a radiation specific internal subroutine
has been added that allocates and inserts radiation specific scalar fields 
from templates into states(1). This took some inspiration from the 
sediments model code above it in that all the fields are actually 
templates in the flml and and names change as needed.

This includes the solution, diagnostic and coupling to fluids fields.

The coupling fields will be calculated from fluids fields via the already 
there diagnositc algorithms. This is useful as it breaks a direct link 
between radiation and fluids and allows for the possibility of 
the radiation and fluids to have different function spaces (or even meshes) 
as they could have very different resolution requirements.

This will insert the fields as needed for each 
neutral particle object initiated in the flml (INT is an integer): 
 - NeutralParticleKeff//object_name
 - NeutralParticleFluxGroupINT//object_name (for all energy groups)
 - NeutralParticleFluxOldGroupINT//object_name (for all energy groups) 
 - NeutralParticlePrescribedSourceINT//object_name (for all source)
 - RadRegionFission//object_name
 - RadRegionPower//object_name
 - RadMaterialVolumeFractionINT//object_name
 - RadMaterialTemperature//object_name
 - RadMaterialFission//object_name
 - RadMaterialPower//object_name
 - RadPorosity//object_name
 - RadPorousTemperature//object_name
 - RadPorousFission//object_name
 - RadPorousPower//object_name
 - RadTotalFission//object_name
 - RadTotalPower//object_name
 - DelayedNeutronPrecursorGroupINTRegion//object_name (for all delayed groups)
 - DelayedNeutronPrecursorGroupINTPorous//object_name (for all delayed groups)
 - DelayedNeutronPrecursorGroupINTMaterial//object_name (for all delayed groups)

All these fields are associated with a particular named function space  
specific to radiation by default.

Then the allocate of auxillary scalar fields cycles such as to do 
nothing for fields starting with 'NeutralParticle' as we dont want 
Old or Iterated fields appearing for radiation that will depend 
on fluids options


FluidityProject/multifluids,FluidityProject/multifluids,FluidityProject/multifluids,FluidityProject/multifluids