Fix initialization of reference pressure for physics when using MPAS dynamical core

[kuanchihwang]

Tag name (required for release branches):

None

Originator(s):

PeterHjortLauritzen, nusbaume, kuanchihwang

Description (include the issue title, and the keyword ['closes', 'fixes', 'resolves'] followed by the issue number):

Presently, MPAS dynamical core is not initializing reference pressure for physics in the same way as others. This PR backports the fix from ESCOMP/CAM#1169.

To confirm the fix, observe log entries similar to the following in atm.log.<job-id>.<date>-<time>. The reference pressure at surface now starts from 1000 hPa. Previously, it starts from 1013.25 hPa.

dyn_debug_print (0): Reference layer information:
dyn_debug_print (0): ----- | -------------- | --------------
dyn_debug_print (0): Index |     Height (m) | Pressure (hPa)
dyn_debug_print (0): ----- |   44999.999819 |       1.431394
dyn_debug_print (0):     1 |   43953.595504 |       1.655836
dyn_debug_print (0): ----- |   42907.191189 |       1.880278
dyn_debug_print (0):     2 |   41877.531046 |       2.177235
dyn_debug_print (0): ----- |   40847.870904 |       2.474192
dyn_debug_print (0):     3 |   39835.232205 |       2.867591
dyn_debug_print (0): ----- |   38822.593507 |       3.260989

... (SNIPPED) ...

dyn_debug_print (0): ----- |    1292.109873 |     856.014954
dyn_debug_print (0):    30 |     997.820649 |     887.686997
dyn_debug_print (0): ----- |     703.531425 |     919.359041
dyn_debug_print (0):    31 |     476.271636 |     945.093654
dyn_debug_print (0): ----- |     249.011847 |     970.828267
dyn_debug_print (0):    32 |     124.505924 |     985.414133
dyn_debug_print (0): ----- |       0.000000 |    1000.000000

Closes #315.

Describe any changes made to build system:

None

Describe any changes made to the namelist:

None

List any changes to the defaults for the input datasets (e.g. boundary datasets):

None

List all files eliminated and why:

None

List all files added and what they do:

None

List all existing files that have been modified, and describe the changes:

• M src/dynamics/mpas/dyn_comp.F90
  • Sort statements
• M src/dynamics/mpas/dyn_grid.F90
  • Properly call std_atm_pres to initialize reference pressure for physics
  • Wire up history support
  • More explicit memory management
  • Update code comments
• M src/utils/std_atm_profile.F90
  • Allow custom surface pressure in std_atm_pres
  • Remove unused variable in std_atm_pres

[nusbaume]

Looks good to me, thanks @kuanchihwang!

[peverwhee]

thanks @kuanchihwang !