Preparing for two collections of shapefiles. We create the
procedure `write_nml` so we will be able to call it twice. But it
needs to access the namelist. So the namelist must be at the module level.