Move freeze curve code to registered Julia package

Additional changes:

• Types no longer have Param values by default. They must be added by the user or by the parameterize method which automatically wraps all numerical quantities with Param.
• Prognostic on-grid variables now have flux terms j[x]defined automatically.
• Refactored heat conduction and non linear diffusion implementations.