To obtain data for more isobaric levels in MPAS model

[met-sree]

Hi everyone,
I am new to the MPAS model. I have installed the model and successfully run a test case. In the diagnostic output file (as specified in stream_list.atmosphere.diagnostics), I obtained data at certain isobaric levels: 50, 100, 200, 250, 500, 700, 850, and 925 hPa.

However, I would like to have data at additional isobaric levels in between these. Is it possible to output more isobaric levels? If so, could you please provide suggestions or share a detailed manual or documentation that explains how to configure this?

Thank you,
Sree

 

[Ming Chen]

Sree,

It is totally feasible to output data on more isobaric levels. Please take a look at the codes below:

(1) core_atmosphere/diagnostics/Registry_isobaric.xml

You can add the variables on specific pressure levels to this Registry file. For example, suppose you need temperature at 300 hPa, then you can add the following code to the Registry file:

<var name="temperature_300hPa" type="real" dimensions="nCells Time" units="K"
description="Temperature vertically interpolated to 300 hPa"/>

2) core_atmosphere/diagnostics/isobaric_diagnostics.F

You need to modify this code to obtain temperature_300hPa. A good example to follow is the variable "temperature_250hPa". Just take a look how "temperature_250hPa" is calculated, and the same method can be used to calculate "temperature_300hPa"

Please try and let me know if you have any issues.

 

[met-sree]

Dear Ming Chen,
Thank you for your suggestion. I followed the procedure for adding 1000 hPa, based on the 925 hPa information. I modified both the registry file (Registry_isobaric.xml) and the Fortran code (isobaric_diagnostics.F).

I cleaned the previously installed model and recompiled it with the modifications mentioned above. However, the stream_list.atmosphere.diagnostics file still looks the same as the default one, and there is no information for 1000 hPa."

For your reference, I have attached all the files at the link below. Please take a look and provide any suggestions

Thank you,
Sree

error1 - Google Drive

drive.google.com