0001 function dist = sigma_geo(nlev,p_sigma) 0002 0003 kb = nlev; 0004 if(p_sigma ==1) 0005 for k=1:nlev 0006 dist(k) = -((k-1)/float(kb-1))^p_sigma; 0007 end; 0008 else 0009 for k=1:(kb+1)/2 0010 dist(k) = -((k-1)/float((kb+1)/2-1))^p_sigma/2; 0011 end; 0012 do k=(kb+1)/2+1:kb 0013 dist(k) = ((kb-k)/float((kb+1)/2-1))^p_sigma/2-1.0; 0014 end; 0015 end; 0016