R: Reading a netCDF file as a raster
<p>I'm trying to read into R a netCDF file as a raster. The netCDF file describes mean annual temperature in the ocean as a function of longitude, latitude, and depth. I'm concerned with the surface ocean (i.e., the first level in the netCDF file) and <code>t_an</code> is the name of the variable in the netCDF file. Therefore, I use the following code:</p>

<pre><code>MyRast &lt;- raster("Temperature.nc", level = 1, varname = "t_an")
</code></pre>

<p>This gives me the following warning:</p>

<pre><code>Warning message:
In .getCRSfromGridMap4(atts) : cannot process these parts of the CRS: epsg_code=EPSG:4326
</code></pre>

<p>As you can see, the netCDF file has CRS EPSG 4326 (or WGS 84), yet the raster that is created has the following CRS:</p>

<pre><code>+proj=longlat +lon_0=0 +a=6378137 +rf=298.257232666016
</code></pre>

<p>Any ideas how I read in the file netCDF with the correct CRS?</p>