corrected the error with CPO initial_volume

[srivatsvyas]

1. Made the change discussed to line 194 regarding initial volume fraction. This is reflected in the CPO file output.
2. Added an if else for pole deflection in line 99. Created a variable dt to store time of simulation in line 76. If dt = 0, then create a normal rotation matrix. If dt != 0, constrained the pole deflection to 0.1PIthree rather than 2PIthree. Did not test this as had problems with the output and piezometer addtion.
3. Added both the piezometer and paleowattmeter in line-718.(just substituted the values in the equation for the paleowattmeter directly so will look messy, sorry :(! ). After making the change to the initial volume fraction, this part does not cause a floating point exception anymore but would be better to check from your side (I tried the paleopiezometer, dont feel comfortable deriving all the terms for the paleowattmeter yet, its there but am still a bit unsure for the equation works so didnt test it now).
4. TO-DOs left: the threshold condition for discriminating dislocation creep and diffusion creep.

[MFraters]

Thanks!

1. So is the issue with the initial grain sizes being zero fixed?
2. hmm, interesting approach. Did you check whether limits the rotation to 10 degrees? If it does, we can leave it for now, but it probably should become an input parameter for the function later.
3. I actually think it is good to split up the equation into multiple parts and then combine them. please add a comment explaining the equation you are solving. Also, I doubt you need a long double for this. a normal double should be fine.
4. I can take a look at that.

[MFraters]

dt is usally used for a time difference, you can just call the variable time.