Init_atmosphere has trouble reading grid.nc when creating static.nc

[munez]

Hi,
I am trying to create the static.nc files from the original mesh files with the init_atmosphere of MPAS-8.2.2. The program keeps giving "CRITICAL ERROR: Could not open input file 'x1.40962.grid.nc' to read mesh fields" error even though the grid file is present in the same working directory. The streams.init_atmosphere and the namelist has been attached, along with the output log files.
Thank you!

 

[mgduda]

Everything in your namelist and streams files looks reasonable to me, so perhaps there's something peculiar about the x1.40962.grid.nc file. Could you post the output of the following commands, run from within your working directory?

ls -lL x1.40962.grid.nc

(the first option to ls is a little-L)

ncdump -k x1.40962.grid.nc

ncdump -h x1.40962.grid.nc

Thanks!

 

[munez]

The outputs are as below:

(mpas) mm@ip 120km % ls -lL x1.40962.grid.nc
-rw-r--r--  1 mm  staff  58001928 Jun  3  2013 x1.40962.grid.nc
(mpas) mm@ip 120km % ncdump -k x1.40962.grid.nc
64-bit offset
(mpas) mm@ip 120km % ncdump -h x1.40962.grid.nc
netcdf x1.40962.grid {
dimensions:
    nCells = 40962 ;
    nEdges = 122880 ;
    nVertices = 81920 ;
    maxEdges = 10 ;
    maxEdges2 = 20 ;
    TWO = 2 ;
    vertexDegree = 3 ;
    nVertLevels = 1 ;
    Time = UNLIMITED ; // (0 currently)
variables:
    double latCell(nCells) ;
    double lonCell(nCells) ;
    double meshDensity(nCells) ;
    double xCell(nCells) ;
    double yCell(nCells) ;
    double zCell(nCells) ;
    int indexToCellID(nCells) ;
    double latEdge(nEdges) ;
    double lonEdge(nEdges) ;
    double xEdge(nEdges) ;
    double yEdge(nEdges) ;
    double zEdge(nEdges) ;
    int indexToEdgeID(nEdges) ;
    double latVertex(nVertices) ;
    double lonVertex(nVertices) ;
    double xVertex(nVertices) ;
    double yVertex(nVertices) ;
    double zVertex(nVertices) ;
    int indexToVertexID(nVertices) ;
    int cellsOnEdge(nEdges, TWO) ;
    int nEdgesOnCell(nCells) ;
    int nEdgesOnEdge(nEdges) ;
    int edgesOnCell(nCells, maxEdges) ;
    int edgesOnEdge(nEdges, maxEdges2) ;
    double weightsOnEdge(nEdges, maxEdges2) ;
    double dvEdge(nEdges) ;
    double dv1Edge(nEdges) ;
    double dv2Edge(nEdges) ;
    double dcEdge(nEdges) ;
    double angleEdge(nEdges) ;
    double areaCell(nCells) ;
    double areaTriangle(nVertices) ;
    int cellsOnCell(nCells, maxEdges) ;
    int verticesOnCell(nCells, maxEdges) ;
    int verticesOnEdge(nEdges, TWO) ;
    int edgesOnVertex(nVertices, vertexDegree) ;
    int cellsOnVertex(nVertices, vertexDegree) ;
    double kiteAreasOnVertex(nVertices, vertexDegree) ;

// global attributes:
        :on_a_sphere = "YES             " ;
        :sphere_radius = 1. ;
        :np = 40962 ;
        :n_scvt_iterations = 1000 ;
        :eps = 1.e-10 ;
        :Convergence = "L2" ;
}

 

[mgduda]

It looks like the x1.40962.grid.nc file that you have may be quite old (or at least, not produced using our current tools). Although I can't point to a specific reason why that mesh shouldn't work, it might be worth trying to download a version of the x1.40962 mesh file from the MPAS-A mesh download page (MPAS-Atmosphere mesh downloads).

If you're still seeing errors when using the x1.40962.grid.nc file from the download page, please don't hesitate to follow up here and we'll be glad to take another look at the issue!