Replacing only negative Z values in GeoTIFF file, grdmath?

My code most likely could be written in one single command using grdmath LE instead of grdclip.