This dataset contains the output of the Mars Global Climate Model simulations. This model has a horizontal resolution of 64 x 32 grid points (longitude x latitude) and a vertical resolution of 32 sigma levels and 20 ground levels. The file "dataset_Kobayashi_2026.tar.gz" includes atmospheric, surface, and subsurface variables in ASCII format. The file "plev.asc" includes vertical grids of pressure levels. The file "glev.asc" includes vertical grids of ground levels. The following provides an explanation of variables and their resolutions in the parentheses (nx, ny, nz, nt, nd). where nx, ny, and nz are the numbers of longitudinal, latitudinal, and vertical grids; nt is the number of data points per sol; and nd is the number of sols. The variable starting "colwv2-1h" represents the total water vapor column abundance in m (64, 32, 1, 24, 10). The variable starting "PS-1h" represents the surface pressure in hPa (64, 32, 1, 24, 10). The variable starting "omg-1h-p" represents the vertical velocity in hPa/s (64, 32, 70, 24, 10). The variable starting "QL2-1h-p" represents the mass mixing ratio of water ice clouds in kg/kg (64, 32, 70, 24, 10). The variable starting "RHOAIR-1h-p" represents the air density in kg/m3 (64, 32, 70, 24, 10). The variable starting "RHOWTR-1h-p" represents the water vapor density in kg/m3 (64, 32, 70, 24, 10). The variable starting "U-1h-p" represents the zonal wind in m/s (64, 32, 70, 24, 10). The variable starting "STF-p" represents the mass stream function in kg/s (64, 32, 70, 1, 49). The variable starting "U-12h" represents the zonal wind at local noon in m/s (64, 32, 1, 1, 10). The variable starting "V-12h" represents the meridional wind at local noon in m/s (64, 32, 1, 1, 10). The variable starting "wflux2s-12h" without "1yr" represents the water vapor flux at local noon in kg/m2/s (64, 32, 1, 1, 10). The variable starting "wflux2s-1h" represents the water vapor flux in kg/m2/s (64, 32, 1, 24, 10). The variable starting "wflux2s-8h...1yr" represents the water vapor flux at local 8 a.m. in kg/m2/s (64, 32, 1, 1, 669). The variable starting "wflux2s-12h...1yr" represents the water vapor flux at local noon in kg/m2/s (64, 32, 1, 1, 669). The variable starting "wflux2s...1yr" represents the water vapor flux in kg/m2/s (64, 32, 1, 1, 669). The variable starting "wga2c" represents the total adsorbed water column abundance in kg/m2 (64, 32, 20, 1, 1). The variable starting "wgv2c" represents the total water vapor column abundance in kg/m2 (64, 32, 20, 1, 1). The variable starting "wgi2" without "1h" represents the water ice mass per unit volume of the regolith in kg/m3 (64, 32, 20, 1, 1). The variable starting "TG-1h" represents the ground temperature in K (64, 32, 20, 24, 10). The variable starting "wga2-1h" represents the adsorbed water mass per unit volume of the regolith in kg/m2 (64, 32, 20, 24, 10). The variable starting "wgv2-1h" represents the water vapor mass per unit volume of the regolith in kg/m2 (64, 32, 20, 24, 10). The variable starting "wgi2-1h" represents the pore ice mass per unit volume of the regolith in kg/m2 (64, 32, 20, 24, 10). The data containing the term "inactive" represents the case of the inactive regolith. The data containing the term "homo" represents the case of the active, homogeneous regolith. The data containing the term "inhomo" represents the case of the active, inhomogeneous regolith. The data containing the term "nocloud" represents the case of the inactive regolith without cloud formation. The data containing the term "Ls150" represents the 10-sols simulation from Ls = 150°. The data containing the term "Ls180" represents the 10-sols simulation from Ls = 180°. The data containing the term "Ls210" represents the 10-sols simulation from Ls = 210°. The data containing the term "Ls186" represents the 10-sols simulation from Ls = 186°. The data containing the term "1yr" represents the 1-year simulation starting from Ls = 0°. For the variables with nz = 70, the vertical grid has been converted to pressure levels, as shown in "plev.asc". For the variables with nz = 20, the vertical grid is defined by "glev.asc".