Problem with bathymetry map and contoured data

I am trying to map ocean data over the Pacific Basin with an ETOPO2 file for the coastlines (which is grey land with a black outline). The coastlines and land load fine and then I add the data I want to map.... Where the data field is nan, the black outline of the coast remains, but the fill in grey color goes away. How do I fix this?

Here is my piece of code that graphs the land and data...

[cc,hh] = contourf(flipud(rot90(LON)),flipud(rot90(LAT)),flipud(rot90(Zf)),[1 1], 'k');
set(hh, 'fill', 'on', 'LineWidth',2);
set(get(hh, 'children'), 'facecolor', [128 128 128]/255);
contourf(lon2_2d,lat_2d,sst_2d,40,'linestyle','none');