(<span class="string">'writing vertical schematization data: sigma layers, sigma levels'</span>);
0351 inc = 1./real(nlayers);
0352 siglev = 0:-inc:-1;
0353 <span class="keyword">for</span> i=1:nlevels
0354     Msiglev(i,1:nObcs) = siglev(i);
0355 <span class="keyword">end</span>;
0356 <span class="keyword">for</span> i=1:nlayers
0357     Msiglay(i,1:nObcs) = mean(siglev(i:i+1));
0358 <span class="keyword">end</span>;
0359 <span class="comment">%{</span>
0360 Msiglay = -1*ones(nlayers,nObcs);
0361 tlayer = 