Missing comma and wrong parameter name while in discretization:
"discretization = MOLFiniteDifference([x=>dx, y=>dy], t, approx_order=order grid_type=center_align)"
Which should be :
"discretization = MOLFiniteDifference([x=>dx, y=>dy], t, approx_order=order, grid_align=center_align)"