code: def test_milisecond_temporal_crs__coordinate_system():
      crs = CRS(
          'TIMECRS["GPS milliseconds",'
          'TDATUM["GPS time origin",TIMEORIGIN[1980-01-01T00:00:00.0Z]],'
          "CS[TemporalCount,1],"
          'AXIS["(T)",future,TIMEUNIT["millisecond",0.001]]]'
      )
      assert crs.cs_to_cf() == [
          {
              "standard_name": "time",
              "long_name": "time",
              "calendar": "proleptic_gregorian",
              "units": "millisecond since 1980-01-01T00:00:00.0",
              "axis": "T",
          }
      ]
  
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 51 for calling a global function with example usage: test_milisecond_temporal_crs__coordinate_system() without return types