Problems with S+ and nested GLM models


I was wondering if someone might be able to help. My department has just changed statistic programs to S+ and I'm trying to run a multi-factorial GLM with one nested factor (design below). What was simple task in our old software has stumped everyone in my department. Does anyone know how we can run this test in S+? I'm wondering if there is a code line from R that we could insert into S+.

zone - fixed
Height- fixed
Location(zone)- random
zone x height- fixed
location(zone) x height- random

Any advice or help would be appreciated.