No subject

I guess a different route would be to allow thin walls.,AFAIK, thin walls can be handled by the model.
You'd have to implement that in ini_masks_etc.F after Ro_surf/R_low and hFacW/S are already known, but before hFacW/S are exchanged in line 287, i.e. if Ro_surf(i,:)<=R_low(i-1,:) =>  hFacW(i:,:)=0. It's worth a try, isn't it?
It should lead to kSurfW/S=Nr+1 and then the corresponding rStarFacW/S=0 and there is no issue, right?

Martin
</pre>
            </blockquote>
          </blockquote>
        </blockquote>
      </blockquote>
    </blockquote>
  </body>
</html>

--------------060203020505050407060509--