getting periodic_channel_dx40km_LxLy1040km to work
changes to script_Cartesian_IC