def calc_dh(DTheta, F0, rho, cp, dt):
    """Gets change in height for a change in time per phil's notes?

    Arguments:
    DTheta -- change in potential temperature, K
    F0 -- surface heat flux, w/m2
    rho -- density of air, kg/m3
    cp -- heat capacity of air
    dt -- change in time, s

    Returns:
    dh -- change in height, m  
    
    """     
    dh = 1.0*(dt)*((0.2*F0)/(rho*cp*DTheta))
    return dh