I need to model a volcanic blast wave. This mean that the domain cannot be prismatic, as the bottom corresponds to ground coordinates.
Is there any way to do that? Which is the best way?
An easy task: lets hace an square (-1<x>1; -1<y<1) with just 3 ponits at each coordinate direction. I need something like: for x=0 make transform y=-1 to -0.5, y=0 to y=0.3 and y=1 to y=1. How to do that from makeGrid.m? Or even from a new script, lets say how the wave numbers become affected?