code: def test_3d_ellipsoidal_cs_depth():
      crs = CRS(
          {
              "$schema": "https://proj.org/schemas/v0.2/projjson.schema.json",
              "type": "GeographicCRS",
              "name": "WGS 84 (geographic 3D)",
              "datum": {
                  "type": "GeodeticReferenceFrame",
                  "name": "World Geodetic System 1984 ensemble",
                  "ellipsoid": {
                      "name": "WGS 84",
                      "semi_major_axis": 6378137,
                      "inverse_flattening": 298.257223563,
                  },
              },
              "coordinate_system": {
                  "subtype": "ellipsoidal",
                  "axis": [
                      {
                          "name": "Geodetic latitude",
                          "abbreviation": "Lat",
                          "direction": "north",
                          "unit": {
                              "type": "AngularUnit",
                              "name": "degree minute second hemisphere",
                              "conversion_factor": 0.0174532925199433,
                          },
                      },
                      {
                          "name": "Geodetic longitude",
                          "abbreviation": "Long",
                          "direction": "east",
                          "unit": {
                              "type": "AngularUnit",
                              "name": "degree minute second hemisphere",
                              "conversion_factor": 0.0174532925199433,
                          },
                      },
                      {
                          "name": "Ellipsoidal depth",
                          "abbreviation": "d",
                          "direction": "down",
                          "unit": "metre",
                      },
                  ],
              },
              "area": "World",
              "bbox": {
                  "south_latitude": -90,
                  "west_longitude": -180,
                  "north_latitude": 90,
                  "east_longitude": 180,
              },
          }
      )
      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",
          },
          {
              "standard_name": "height_above_reference_ellipsoid",
              "long_name": "Ellipsoidal depth",
              "units": "metre",
              "positive": "down",
              "axis": "Z",
          },
      ]
  
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 46 for calling a global function with example usage: test_3d_ellipsoidal_cs_depth() without return types