refactor write to a file
print number of populations in RMES input