You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: Murat_modv3D.m
    path of file: ./repos/MuRAT/bin
    
    the code of the file until where you have to start completion: function [modvProp,modvP,modvEqSpace,modvPlot]...
                            =...
    Murat_modv3D(FPath,FLabel,modvOriginal,origin,mLat,mLon,nLat,nLon,nzc)
% function [modvP,modvPropagation,modvEqS,modvPlot]	=...
%     Murat_modv3D(FPath,FLabel,modvOriginal,origin,mLat,mLon,nLat,nLon)
% CREATES the velocity models when using a 3D as input
%
% Input parameters:
%    FPath,FLabel:          Paths to save test velocity figure

    
function [modvProp,modvP,modvEqSpace,modvPlot]...
                            =...
    Murat_modv3D(FPath,FLabel,modvOriginal,origin,mLat,mLon,nLat,nLon,nzc)
% function [modvP,modvPropagation,modvEqS,modvPlot]	=...
%     Murat_modv3D(FPath,FLabel,modvOriginal,origin,mLat,mLon,nLat,nLon)
% CREATES the velocity models when using a 3D as input
%
% Input parameters:
%    FPath,FLabel:          Paths to save test velocity figure
%    modvOriginal:          original 3D velocity model
%    origin:                origin of the spatial grid
%    mLat:                  min/max latitudes of earthquakes
%    mLon:                  min/max longitudes of earthquakes
%    nLat:                  number nodes latitude
%    nLon:                  number nodes longitude
%
% Output parameters:
%    modvP:                 propagation velocity model in standard format
%    modvPropagation:       propagation velocity model in matrix format
%    modvEqS:               equally-spaced inversion velocity model
%    modvPlot:              equally-spaced ploting velocity model (lat/lon)

%% Rescale original lat/lon modV in cartesian coords
% Transform nodes velocity model in meters with WGS84
wgs84                       =   wgs84Ellipsoid("m");
[d_o,az_o]                  =   distance(origin(1),origin(2),...
    modvOriginal(:,1),modvOriginal(:,2),wgs84);

% Transform min/max lat/lon earthquakes in meters with WGS84
[d_ll,az_ll]                =   distance(origin(1),origin(2),...
    mLat',mLon',wgs84);
mLatLon                     =   [d_ll.*sin(az_ll*2*pi/360)...
    d_ll.*cos(az_ll*2*pi/360)];

%% Original velocity model in meters - uneven spacing
modvOrC                     =   sortrows([d_o.*sin(az_o*2*pi/360)...
    d_o.*cos(az_o*2*pi/360) modvOriginal(:,3:4)],[1,2,-3]);
xOrC                        =   modvOrC(:,1);
yOrC                        =   modvOrC(:,2);
zOrC                        =   modvOrC(:,3);
vOrC                        =   modvOrC(:,4);


%% Velocity model for ray tracing - inversion
% For ray-tracing, you need to space evenly with the same number of
% x-y nodes as in input - here we extend