@feathern OK I have convinced myself that reference_type = 5 produces the same simulation as reference_type = 2, just rescaled. I added two input files to show how this works. One is the original "model_15" from our parameter survey and the other is the non-dimensional equivalent.
Note the different conventions in Rossby number.
I checked volume averaged output non-d numbers, differential rotation, meridional circulation, and radial energy flux balance and they all seem the same (after rescaling) up to < 1%.
I think the 1% difference comes from the heating being normalized slightly differently for ref type 5 (it doesn't get rescaled in Transport_Dependencies for ref type 5 to match the conductive flux loss out the boundaries). This makes the input Rayleigh numbers slightly different.
Anyway, I think we're good to go! And these input files should show people how to use the new reference type.
@feathern OK I have convinced myself that reference_type = 5 produces the same simulation as reference_type = 2, just rescaled. I added two input files to show how this works. One is the original "model_15" from our parameter survey and the other is the non-dimensional equivalent.
Note the different conventions in Rossby number.
I checked volume averaged output non-d numbers, differential rotation, meridional circulation, and radial energy flux balance and they all seem the same (after rescaling) up to < 1%.
I think the 1% difference comes from the heating being normalized slightly differently for ref type 5 (it doesn't get rescaled in Transport_Dependencies for ref type 5 to match the conductive flux loss out the boundaries). This makes the input Rayleigh numbers slightly different.
Anyway, I think we're good to go! And these input files should show people how to use the new reference type.