commented out and added thoughts to a bunch of places where there seems to be a hardwired east-west periodic boundary condition. I think the original intent of the code was for global modeling and hence the fact that it is hardwired,  but I am not sure if another boundary condition is needed to replace that one. That is a question to yet be answered.