gridquantity3D.m
dbcb156c5ec93Dmpimode

gridquantity3D.m
View Options

			classdef gridquantity3D < h5quantity3D

			properties
			centered
			fileversion
			end

			methods
			function obj=gridquantity3D(filename, dataset, nr, nthet, nz, scale, centered)
			if nargin<5
			scale=1;
			end
			obj=obj@h5quantity3D(filename, dataset, nr, nthet, nz, scale);
			if nargin<7
			obj.centered=false;
			else
			obj.centered=centered;
			end
			try
			obj.fileversion = h5readatt(filename,'/data','file_version');
			catch
			obj.fileversion=1;
			end
			end

			function quantity=val(obj,indices)
			if strcmp(indices{1},':')
			r=1:obj.nr;
			else
			r=indices{1};
			end
			if strcmp(indices{2},':')
			thet=1:obj.nthet;
			else
			thet=indices{2};
			end
			if strcmp(indices{3},':')
			z=1:obj.nz;
			else
			z=indices{3};
			end
			if strcmp(indices{4},':')
			t=1:obj.nt;
			else
			t=indices{4};
			end
			if obj.centered
			quantity=zeros((obj.nz-1)*(obj.nr-1),obj.nthet,length(t));
			else
			quantity=zeros((obj.nz)*(obj.nr),obj.nthet,length(t));
			end
			for i=1:length(t)
			quantity(:,:,i) = h5read(obj.filename, obj.dataset,[1 1 t(i)],[Inf Inf 1]) ;
			end
			if obj.centered
			quantity=reshape(quantity,obj.nr-1,obj.nz-1,obj.nthet,[]);
			quantity=padarray(quantity,[1,1],0,'post');
			else
			quantity=reshape(quantity,obj.nr,obj.nz,obj.nthet,[]);
			end
			quantity=permute(quantity,[1,3,2,4]);
			quantity=quantity(r,thet,z,:)*obj.scale;
			end
			end
			end