This repository holds the Zhang McFarlane convective scheme (Zhang and McFarlane, 1995) in CAM 5.0 (Neale et al. 1995), modified to "protect" convection from entraining dry air in selected locations. Each pocket is ~ 2.5 x 2.5 degrees in size. Inside this pocket, convection only entrains saturated air.
The code works by asking if the current location is inside a protected location. If yes, then it implements protections in two separate places within the code by:
i) Transforming the entraining CAPE closure (Neale et al. 2012) into an undilute CAPE in the subroutine 'closure'.
ii) Make the updrafts and downdrafts interact with saturated air in the subroutine 'cldprp'.
The code uses Newton's method to invert the plume entropy to deduce temperature, which sometimes fails to converge. In case of failure, the code switches to Halley's method (which has not yet failed).
Search for !FA within the code to examine all the modifications.
To run these codes in your version of CAM/CESM, place:
i) user_nl_cam in your case directory (CASE_DIR)
ii) namelist_definition.xml, zm_conv.F90 and zm_conv_intr.F90 in CASE_DIR/SourceMods/src.cam, and build.