README.txt           December 3, 2025

                     by Hiromasa Yoshimura
                        Meteorological Research Institute, Japan Meteorological Agency

This directory contains output data from experiments using a two-dimensional x–z nonhydrostatic model.  
The data are used in the following paper.

[Paper]

  Yoshimura, H.: A nonhydrostatic spectral atmospheric model using non-constant coefficient semi-implicit
                 and vertically conservative semi-Lagrangian schemes
                  (Geoscientific Model Development, 2026)

[Directory]

  [Mountain wave test case (St-MIP; Satomura et al., 2003)]

    A1/        Case A1.
    A2/        Case A2.
    A3/        Case A3.
    A4/        Case A4.
    A4_prev/   Case A4. Previous method is used for the boundary condition.
    D1/        Case D1.
    D2/        Case D2.

  [Extremely steep mountain test case (Yoshimura, 2026)]

    S1/         Case S1.
    S1_winit1/  Case S1. Different initial condition is used for w.

  [Mountain wave test cases (Schär et al., 2002)]

    schar1_0/             Case Schär1. dt=5s.
    schar1/               Case Schär1. dt=30s.
    schar1_2/             Case Schär1. dt=60s.

    schar1_0_not_SLTPC/   Case Schär1. dt=5s.  SLTPC is not used.
    schar1_not_SLTPC/     Case Schär1. dt=30s. SLTPC is not used.
    schar1_2_not_SLTPC/   Case Schär1. dt=60s. SLTPC is not used.

    schar2_0/             Case Schär2. dt=5s.
    schar2/               Case Schär2. dt=30s.
    schar2_2/             Case Schär2. dt=60s.

    schar2_0_not_SATEM/   Case Schär2. dt=5s.  SATEM is not used.
    schar2_not_SATEM/     Case Schär2. dt=30s. SATEM is not used.
    schar2_2_not_SATEM/   Case Schär2. dt=60s. SATEM is not used.

    schar2_0_not_Hermite/ Case Schär2. dt=5s.  Hermite interpolation for SLTPC is not used.
    schar2_not_Hermite/   Case Schär2. dt=30s. Hermite interpolation for SLTPC is not used.
    schar2_2_not_Hermite/ Case Schär2. dt=60s. Hermite interpolation for SLTPC is not used.

  [Density current test case (Straka et al., 1993)]

    dc_50/       dx=dz=50m.  dt=3s.
    dc_100/      dx=dz=100m. dt=6s.
    dc_200/      dx=dz=200m. dt=12s.
    dc_400/      dx=dz=400m. dt=25s.

  [Warm bubble test case, Gaussian profile (Robert, 1993)]

    wb_g/              Default settings (quadratic grid is used).
    wb_g_lineargrid/   Linear grid is used.
    wb_g_cubicgrid/    Cubic grid is used.
    wb_g_halfdt/       dt=2.5s is used instead of dt=5.0s.
    wb_g_0/            dt=2.5s and dx=dz=5.0m are used instead of dt=5.0s and dx=dz=10.0m.

  [Gravity wave test case (Skamarock and Klemp, 1994)]

    gw_dt6/      dt=6s.
    gw/          dt=30s.

    gw_dt6_prev/ dt=6s.  Previous method for vertical interpolation is used.
    gw_prev/     dt=30s. Previous method for vertical interpolation is used.


[Data]

  *.nc : NetCDF file
  *.ncctl: GrADS control file for NetCDF file (text file).

  [Model half levels]
   phihf.nc: Geopotential (m**2/s**2)
   w.nc    : Vertical wind w (m/s)

  [Model full levels]
   phi.nc   : Geopotential (m**2/s**2)
   theta.nc : Potential temperature perturbation (K)

    - The original REAL(4) binary files and GrADS control files were converted to NetCDF files by using CDO.
      > cdo -f nc4 -z zip_5 import_binary aaa.ctl aaa.nc
      https://code.mpimet.mpg.de/projects/cdo/

  *.py    : Python script for creating figures.
  *.gexec : GrADS script for creating figures.
  *.png   : Figures created by the scripts.
    
[Website]

  https://climate.mri-jma.go.jp/pub/archives/Yoshimura_Nonhydro_2D_XZ_Testcase/