README               May 27, 2021

                     by Hiromasa Yoshimura
                        Meteorological Research Institute, Japan Meteorological Agency


Output data in shallow water model experiments in Yoshimura (2021).
(Preliminary version.)


[Paper]

  Yoshimura, H.: Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow water model, 2021.


[Directory]

  (N  : Truncation wavenumber.)
  (IxJ: Number of longitudinal and latitudinal grid points.)
  (DFS: Double Fourier Series)
  (SH : Spherical Harmonics)
      
  Williamson2/        Williamson test case 2   (Output time interval is 1 hour)

    N63_128x64/
      dfsold_grid0/       Old DFS model using Grid [0].  N=63. IxJ=128x64
      dfs_grid0/          New DFS model using Grid [0].  N=63. IxJ=128x64
      dfs_grid1/          New DFS model using Grid [1].  N=63. IxJ=129x64
      dfs_grid-1/         New DFS model using Grid [-1]. N=63. IxJ=127x64
      sh/                 SH model.                      N=62. IxJ=128x64

  Williamson5/        Williamson test case 5   (Output time interval is 24 hour)

    N42_128x64/
      dfs_grid0/          New DFS model using Grid [0].  N=42. IxJ=128x64.
      dfs_grid1/          New DFS model using Grid [1].  N=42. IxJ=129x64.
      dfs_grid-1/         New DFS model using Grid [-1]. N=42. IxJ=127x64.
      sh/                 SH model.                      N=42. IxJ=128x64.

    N63_128x64/
      dfsold_grid0/       Old DFS model using Grid [0].  N=63. IxJ=128x64.
      dfs_grid0/          New DFS model using Grid [0].  N=63. IxJ=128x64.
      dfs_grid1/          New DFS model using Grid [1].  N=63. IxJ=129x64.
      dfs_grid-1/         New DFS model using Grid [-1]. N=63. IxJ=127x64.
      sh/                 SH model.                      N=62. IxJ=128x64.

    N639_1920x960/
      dfs_grid0/          New DFS model using Grid [0].  N=639. IxJ=1920x960.
      dfs_grid1/          New DFS model using Grid [1].  N=639. IxJ=1921x960.
      dfs_grid-1/         New DFS model using Grid [-1]. N=639. IxJ=1919x960.
      sh/                 SH model.                      N=639. IxJ=1920x960.

    N959_1920x960/
      dfsold_grid0/       Old DFS model using Grid [0].  N=959. IxJ=1920x960.
      dfs_grid0/          New DFS model using Grid [0].  N=959. IxJ=1920x960.
      dfs_grid1/          New DFS model using Grid [1].  N=959. IxJ=1921x960.
      dfs_grid-1/         New DFS model using Grid [-1]. N=959. IxJ=1919x960.
      sh/                 SH model.                      N=958. IxJ=1920x960.

  Galewsky/           Galewsky test case   (Output time interval is 24 hour)

    N10239_30720x15360/
      dfs_grid0/          New DFS model using Grid [0].  N=10239. IxJ=30720x15360.
      dfs_grid1/          New DFS model using Grid [1].  N=10239. IxJ=30721x15360.
      dfs_grid-1/         New DFS model using Grid [-1]. N=10239. IxJ=30719x15360.
      sh/                 SH model.                      N=10239. IxJ=30720x15360.
    
  truncation_error/   Check of the error due to meridional wavenumber truncation


[Data]

  *.nc : NetCDF file
  *.ncctl: Grads contol file for NetCDF file (text file).

  data.nc, data.ncctl: Simulated lon-lat data.
    Variables: h (height), u (zonal wind), v (meridional wind), vor (vorticity), div (divergence).

  vor.nc, vor.ncctl: Simulated lon-lat vorticity data after 144-hour integration.
                     (For Galewsky test case.)
    Variables: vor (vorticity)

  norm.nc, norm.ncctl: Error of predicted height (for Williamson test case 2).
    Variables: norm_L1 (L1 norm), norm_L2 (L2 norm), norm_Lmax (Lmax norm).

  spectrum.nc, spectrum.ncctl: Kinetic Energy spectrum.
    Variables: KE_VOR (Kinetic Energy (Vorticity)), KE_DIV (Kinetic Energy (Divergence)).

  In 'truncation_error' directory
    orszag.nc, orszag.ncctl : Orszag's method to calculate expansion coefficients
    cheong.nc, cheong.ncctl : Cheong's method to calculate expansion coefficients
    new.nc   , new.ncctl    : New method to calculate expansion coefficients
      lon: Wavenumber
      lat: Latitude
      t=1: Original values
      t=2: Truncated values


  - Output 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/
    http://cola.gmu.edu/grads/
    

[Web site]

  https://climate.mri-jma.go.jp/pre/Yoshimura_DFS_SW_Testcase_2021/
  ( Final version will be put at
    https://climate.mri-jma.go.jp/pub/archives/Yoshimura_DFS_SW_Testcase_2021/ )