This directory contains the necessary code and input to calculate the
balanced Sawyer-Eliassen circulation response to a given imposed
(radiative) heat forcing, based on Pendergrass and Willoughby (2009,
MWR, available from https://doi.org/10.1175/2008MWR2657.1).

The code is kindly provided courtesy of Hugh E. Willoughby
(willough@fiu.edu). Minor modifications from its originalversion
have been made to 1) read an input heat forcing instead of
prescribe it analytically and 2) to reduce the maximum radius of
the model to 1000 km (from 1500 km).

See associated publication (Ruppert et al. 2020) for further detail.

Directory/file listing:

CODE: the model, solver code.

DATA: the input forcing necessary to run the model.

DATA/WRF_SND.txt: the input base-state sounding.

DATA/Test_0/: directory containing both the heat forcing file and the
resulting circulation response (including *png images).

DATA/Test_0/seq_force_lwradanom.txt: the prescribed radiative heat
forcing (see CODE/LQ_LIB.f90 for how to read in Fortran).

DATA/Test_0/SetUp.txt: prescribed settings for initial vortex, etc.

FCIRC.exe is the executable that calculates the balanced Sawyer-Eliassen
circulation response to the imposed (radiative) heat forcing.