A2 - Metric Scalars

Computing Metric Scalars

Get Metric

%% Alcubierre
gridSize = [1 20 20 20];
worldCenter = (gridSize+1)./2;
velocity = 0.5;
R = 5;
sigma = 0.5;
Metric = metricGet_AlcubierreComoving(gridSize,worldCenter,velocity,R,sigma);

Get Scalars

[expansionScalar, shearScalar, vorticityScalar] = getScalars(Metric);

Plot Scalars

surfq(expansionScalar(1,:,:,10),'EdgeColor','none')
title('Expansion Scalar')

surfq(shearScalar(1,:,:,10),'EdgeColor','none')
title('Shear Scalar')

surfq(vorticityScalar(1,:,:,10),'EdgeColor','none')
title('Vorticity Scalar')

When copying this code to your own script, add figure() in between each call to create new figures for each plot.