[MITgcm-support] how to represent advective forcing

[William Boos]

I am using the MIT GCM in 2D (latitude-height), and would like to 
represent a dynamical forcing (such as a zonally propagating wave) 
passing through the model.  My thought is to use a vertical velocity 
forcing in the advection tendencies for temperature, humidity, and 
momentum.   Is there some standard way to do this, or would I need to 
modify the advection code (e.g. specify my applied vertical velocity in 
subroutines THERMODYNAMICS and CALC_COMMON_FACTORS)?

Or is there a better method for representing a full-field dynamical forcing?

Regards,
Bill Boos