Error in compute_layer_mean: pressure should increase with index

[edougherty]

Hi there,

I am trying to run MPAS using a 15-3km elliptical mesh. When I get to the stage of running the model, I get this error:

ERROR: Error in compute_layer_mean: pressure should increase with index

My static file appears to look reasonable, but they are issues in my init.nc file:

/glade/derecho/scratch/doughert/MPAS/6_Sept_2013_PR_15-3km_ellip2/PR_old.init.nc

When I look at the field 'surface pressure', the values are all nans. I have checked my intermediate file format:

/glade/derecho/scratch/doughert/MPAS/6_Sept_2013_PR_15-3km_ellip2/FILE*

And everything looks OK. Does anyone know what might be going on? For reference, my namelists are located here:

/glade/derecho/scratch/doughert/MPAS/6_Sept_2013_PR_15-3km_ellip2

 

[Ming Chen]

What data did you ungrib to provide initial condition for MAPS run?

 

[edougherty]

I actually used pywinter(pywinter — pywinter 2.0.6 documentation) to create the intermediate files using high-resolution CESM data.

 

[Ming Chen]

I guess the CESM data you processed is netCDF format. In this case pywinter is a good tool to handle the CESM data. However, please double check the output, and make sure the data are vertically monotonic , i.e., the vertical levels must be from bottom-up or from top-down.