GWDO scheme error

[amandacarolsq]

Hello everyone,
I'm trying to run a limited-area (regional) simulation with MPAS, but the model fails during initialization with the following errors:

Setting nominalMinDc to 60000.0 based on namelist option config_len_disp
ERROR: *******************************************************************************
ERROR: The GWDO scheme requires valid var2d, con, oa{1,2,3,4}, and ol{1,2,3,4} fields,
ERROR: but these fields appear to be zero everywhere in the model input.
ERROR: Either set config_gwdo_scheme = 'off' in the &physics namelist, or generate
ERROR: the GWDO static fields with the init_atmosphere core.
ERROR: *******************************************************************************

Checking consistency of limited-area settings...
- config_apply_lbcs = T
- Maximum value in bdyMaskCell = 0
ERROR: config_apply_lbcs = true, but no boundary cells found in the bdyMaskCell field.
ERROR: Please ensure that config_apply_lbcs = false for global simulations.

Checking consistency of dynamics settings...
----- done checking dynamics settings -----

ERROR: Please correct issues with the model input fields and/or namelist.
CRITICAL ERROR: Core init failed for core atmosphere

Could anyone please help me with a solution to this issue?
Thank you!

 

[Ming Chen]

The error message indicates some fundamental issues in this case.

Would you please clarify step by step how you set up this regional mpas case?

It will be helpful to upload your namelist.atmosphere and namelist.init_atmopshere for me to take look.

Thanks.

 

[amandacarolsq]

Thank you for your reply.

I have attached both namelist.atmosphere and namelist.init_atmosphere as requested.

Please let me know if you need any further details about the case setup.

I appreciate your help.

 

[Ming Chen]

Hi,

Thank you for uploading the files. However, these files don't provide helpful information for me to figure out what is wrong.

As I said in my previous post, please clarify step by step how you set up this case:

(1) what is your mesh structure?

(2) how did you create static, and initial and lateral boundary conditions for this case?

(3) what is your input data that is used to drive the MPAS run?

 

[amandacarolsq]

Hi,

(1) what is your mesh structure?


Mesh structure:

Name: SaoPaulo
Type: ellipse
Point: -22.0, -48.0
Semi-major-axis: 450000
Semi-minor-axis: 350000
Orientation-angle: 0
(2) how did you create static, and initial and lateral boundary conditions for this case?
I used:

./create_region SaoPaulo.pts x1.163842.grid
gpmetis -minconn -contig -niter=200 SaoPaulo.graph.info 64

Then, I modified namelist.init_atmosphere and streams.init_atmosphere to generate the SaoPaulo.static.nc file.

The streams.init_atmosphere file was configured as follows:

<streams>
<immutable_stream name="input"
                  type="input"
                  filename_template="SaoPaulo.grid.nc"
                  input_interval="initial_only" />


<immutable_stream name="output"
                  type="output"
                  filename_template="SaoPaulo.static.nc"
                  packages="initial_conds"
                  output_interval="initial_only" />


<immutable_stream name="surface"
                  type="output"
                  filename_template="#filename_templates#"
                  filename_interval="none"
                  packages="sfc_update"
                  output_interval="86400" />


</streams>

After that, I ran namelist.init_atmosphere and streams.init_atmosphere again to generate the initial conditions and the lateral boundary conditions, modifying the namelists accordingly for each case.

(3) what is your input data that is used to drive the MPAS run?

I am using GFS data.