<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hello,</div><div><br></div><div>I am trying to overlay two different grids, which are both on different projections to see the differences between the grids.</div><div><br></div><div>One grid (rectilinear, lat-lon) has attributes:</div><div><div>dimensions:</div><div>        lat = 258 ;</div><div>        lon = 600 ;</div><div>variables:</div><div>        double mask(lat, lon) ;</div><div>                mask:_FillValue = 1.00000002004088e+20 ;</div><div>                mask:long_name = "Land Mask" ;</div><div>                mask:missing_value = 1.e+20f ;</div><div>                mask:standard_name = "land_binary_mask" ;</div><div>                mask:units = "1" ;</div><div>        double lat(lat) ;</div><div>                lat:long_name = "latitude" ;</div><div>                lat:units = "degrees_north" ;</div><div>                lat:standard_name = "latitude" ;</div><div>        double lon(lon) ;</div><div>                lon:long_name = "longitude" ;</div><div>                lon:units = "degrees_east" ;</div><div>                lon:standard_name = "longitude" ;</div></div><div><br></div><div>The other grid (cylindrical, rotated mercator) has attributes:</div><div><div>dimensions:</div><div>        y = 278 ;</div><div>        x = 334 ;</div><div>variables:</div><div>        float mask(y, x) ;</div><div>                mask:_FillValue = 1.e+20f ;</div><div>                mask:coordinates = "lat lon" ;</div><div>                mask:grid_mapping = "rotated_mercator" ;</div><div>                mask:long_name = "Land Area Fraction" ;</div><div>                mask:missing_value = 1.e+20f ;</div><div>                mask:standard_name = "land_area_fraction" ;</div><div>                mask:units = "%" ;</div><div>        double lat(y, x) ;</div><div>                lat:_Coordinate_AxisType = "Lat" ;</div><div>                lat:long_name = "Latitude on Cross Points" ;</div><div>                lat:standard_name = "latitude" ;</div><div>                lat:units = "degrees_north" ;</div><div>                lat:grid_mapping = "rcm_map" ;</div><div>        double lon(y, x) ;</div><div>                lon:_Coordinate_AxisType = "Lon" ;</div><div>                lon:long_name = "Longitude on Cross Points" ;</div><div>                lon:standard_name = "longitude" ;</div><div>                lon:units = "degrees_east" ;</div><div>                lon:grid_mapping = "rcm_map" ;</div><div>        char oblique_mercator ;</div><div>                oblique_mercator:grid_mapping_name = "oblique_mercator" ;</div><div>                oblique_mercator:latitude_of_projection_origin = 46. ;</div><div>                oblique_mercator:longitude_of_projection_origin = 263. ;</div><div>                oblique_mercator:_CoordinateTransformType = "Projection" ;</div><div>                oblique_mercator:_CoordinateAxisTypes = "GeoX GeoY" ;</div><div>                oblique_mercator:azimuth_of_central_line = 90. ;</div><div>                oblique_mercator:scale_factor_at_projection_origin = 1. ;</div><div>                oblique_mercator:false_easting = 0. ;</div><div>                oblique_mercator:false_northing = 0. ;</div><div>        char rotated_mercator ;</div><div>                rotated_mercator:grid_mapping_name = "rotated_mercator" ;</div><div>                rotated_mercator:azimuth_of_central_line = 90. ;</div><div>                rotated_mercator:latitude_of_projection_origin = 46.5 ;</div><div>                rotated_mercator:longitude_of_projection_origin = 263. ;</div><div>                rotated_mercator:scale_factor_at_projection_origin = 1. ;</div><div>                rotated_mercator:false_easting = 0. ;</div><div>                rotated_mercator:false_northing = 0. ;</div><div>        double x(x) ;</div><div>                x:_CoordinateAxisType = "GeoX" ;</div><div>                x:axis = "X" ;</div><div>                x:long_name = "x-coordinate in Cartesian system" ;</div><div>                x:standard_name = "projection_x_coordinate" ;</div><div>                x:units = "m" ;</div><div>        double y(y) ;</div><div>                y:_CoordinateAxisType = "GeoY" ;</div><div>                y:axis = "Y" ;</div><div>                y:long_name = "y-coordinate in Cartesian system" ;</div><div>                y:standard_name = "projection_y_coordinate" ;</div><div>                y:units = "m" ;</div><div><br></div></div><div>I have been following something similar to <a href="https://www.ncl.ucar.edu/Applications/Scripts/datagrid_3.ncl">https://www.ncl.ucar.edu/Applications/Scripts/datagrid_3.ncl</a> from <a href="https://www.ncl.ucar.edu/Applications/datagrid.shtml">NCL Graphics: Drawing locations of data values via markers or lines (ucar.edu)</a>, but it seems to only plot the lat-lon grid and not also the rotated mercator grid.</div><div><br></div><div>Is this possible to plot two different projections of a grid as an overlay?</div><div><br></div><div>Thanks,</div><div><br></div><div>BJF</div></div></div></div></div></div>