ReqCycTorque(probdata)

function [model avail] = ReqCycTorque(probdata)

if isfield(probdata,'pmodel2')
                %create a list with available models
                modellist=fieldnames(probdata.pmodel2);
                comparelist={'ZA','MAe','MAn','MAp'};
                for k=1:length(comparelist);
                    if sum(strcmp(modellist,comparelist(k)))~=0
                        avail(k)=1;
                    else
                        avail(k)=0;
                    end
                end

                %define the models
                if avail(1)==1
                model.ZA=probdata.pmodel2.ZA;
                else
                    modelZA=0;
                end
                if avail(2)==1
                model.MAe=probdata.pmodel2.MAe;
                else
                    model.MAe=0;
                end
                if avail(3)==1
                model.MAn=probdata.pmodel2.MAn;
                else
                    model.MAn=0;
                end
                if avail(4)==1
                model.MAp=probdata.pmodel2.MAp;
                else
                    model.MAp=0;
                end
else
    avail=[0 0 0 0];
    model.ZA=0;
    model.MAe=0;
    model.MAn=0;
    model.MAp=0;
end