A classical result by Kreweras (1965) allows one to compute the number of plane partitions of a given skew shape and bounded parts as certain determinants. We prove that these determinants expand as polynomials with nonnegative coefficients. This result can be reformulated in terms of order polynomials of cell posets of skew shapes, and explains important positivity phenomena about the Ehrhart polynomials of shard polytopes, matroids, and order polytopes. Among other applications, we generalize a positivity statement from Schubert calculus by Fomin and Kirillov (1997) from straight shapes to skew shapes. We show that all shard polytopes are Ehrhart positive and, stronger, that all fence posets, including the zig-zag poset, and all circular fence posets have order polynomials with nonnegative coefficients. We discuss a general method for proving positivity which reduces to showing positivity of the linear terms of the order polynomials. We propose positivity conjectures on other relevant classes of posets.