Zonal Refinement

This feature allows a reduction of the number of points in x/y by a factor of 2 at specified k indices referred to as interfaces. Using this method, simulations can reach altitudes of ~500 km.

New Input Options

k_interface1:Index for interface refinement k_interface2:Index for interface refinement ...: i_atmos_var:1-Use variable atmospheric composition (background.gas), 0-do not

For high altitude simulations, it is important to allow variable composition. Namely gas constant (R_gas) and specific heat ratio (gam_gas) as functions of altitude in the file background.gas. A sample is included in (cgcam/examples/4Dwind).

Constraints

Interfaces must be above the fade region for terrain cases. The number of points in x/y must be divisible at each interface. Due to boundary extrapolations, the number of points per processor is required to be more than 1. Therefore for large cases, it is recommended that the user specify n_pencil_y as the processor limit is approached. A new program tst_decomp -n [num_procs] is useful to determine the number of processors and decomposition.

Output

Each zone collects at writes its data to a separate file. Therefore to visualize x-z or y-z planes, multiple files are loaded. Additional zones can be added on restart by modifying the first header file.