nuserplot.m

用于计算脉冲激光照射下物体的温度分布

function userplot(Rnode,Znode,Tnode,ptype)
%   画出浓度分布随R、Z、T的变化曲线
%   输入：R轴节点序号Rnode,Z轴节点序号Znode,时间轴节点序号Tnode,画图类型ptype
%       1<=Rnode<=M1 ; 1<=Znode<=L1; 1<= Tnode<=td,
%       ptype=1:给定R、T 浓度随Z轴变化曲线,线性坐标
%       ptype=2:给定R、T 浓度随Z轴变化曲线,N对数坐标
%       ptype=3:给定Z、T 浓度随Z轴变化曲线,线性坐标
%       ptype=4:给定Z、R 浓度随T轴变化曲线,线性坐标

global rd zd R Z r z rn zn td tn L1 M1 TT;


switch ptype
    case 1,
        plot(z([1:end]),[N([L1*(Rnode-1)+1:L1*(Rnode-1)+L1],Tnode)]',':r'),xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3)'),
        title('the Zn Concentration along z axis'),
    case 2,
        axis([0 zd 1 10^8]),
        semilogy(z([1:end]),[N([L1*(Rnode-1)+1:L1*(Rnode-1)+L1],Tnode)]',':r'),
        xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3) in logarithm'),
        title('the Zn Concentration along z axis(in logarithm)'),