script for total establishment plots