1/nlayers;
0362 temp   = 0.5*tlayer;
0363 <span class="keyword">for</span> i=1:nlayers
0364     Msiglay(i,:)= temp*Msiglay(i,:);
0365     temp        = temp + tlayer;
0366 <span class="keyword">end</span>
0367 
0368 Msiglev = -1*ones(nlevels,nObcs);
0369 temp    = 0.0;
0370 <span class="keyword">for</span> i=1:nlevels
0371     Msiglev(i,:)= temp*Msiglev(i,:);
0372     temp        = temp + tlayer;
0373 <span class="keyword">end</span>
0374 <span class="comment">%}</span>
0375 netcdf.putVar(nc,siglay_varid,Ms