code: def test_ellipsoidal_cs__geodetic():
      crs = CRS.from_epsg(4326)
      crs.cs_to_cf() == [
          {
              "standard_name": "latitude",
              "long_name": "geodetic latitude coordinate",
              "units": "degrees_north",
              "axis": "Y",
          },
          {
              "standard_name": "longitude",
              "long_name": "geodetic longitude coordinate",
              "units": "degrees_east",
              "axis": "X",
          },
      ]
  
imports: import numpy
  import pytest
  from numpy.testing import assert_almost_equal
  from packaging import version
  from pyproj import CRS
  from pyproj.crs import ProjectedCRS
  from pyproj.crs._cf1x8 import _try_list_if_string
  from pyproj.crs.coordinate_operation import LambertAzimuthalEqualAreaConversion, LambertCylindricalEqualAreaConversion, MercatorAConversion, OrthographicConversion, PolarStereographicAConversion, PolarStereographicBConversion, SinusoidalConversion, StereographicConversion, VerticalPerspectiveConversion
  from pyproj.exceptions import CRSError
  from test.conftest import PROJ_GTE_901, PROJ_LOOSE_VERSION
label_desc: Use image node_id 45 for calling a global function with example usage: test_ellipsoidal_cs__geodetic() without return types