Rename DBF field `date_index` to `days_1950`

Analysis/plot_eddy_contours.py

@@ -208,7 +208,7 @@ def open_shpt(shp_coll_dir):
         except shapefile.ShapefileException:
             print(f"Shapefile {layer} not found in {shp_coll_dir}.")
 
-    k1 = readers["extremum"].record(0).date_index
+    k1 = readers["extremum"].record(0).days_1950
     fname = path.join(shp_coll_dir, "ishape_last.txt")
     
     try:

shp_tr_create.f90

@@ -21,7 +21,7 @@ contains
          hshp%extremum)
     call dbf_add_field_03(hshp%extr_ssh, hshp%extremum, 'ssh', ftdouble, &
          nwidth = 13, ndecimals = 6)
-    call dbf_add_field_03(hshp%extr_date, hshp%extremum, 'date_index', &
+    call dbf_add_field_03(hshp%extr_date, hshp%extremum, 'days_1950', &
          ftinteger, nwidth = 5, ndecimals = 0)
     call dbf_add_field_03(hshp%extr_eddy_index, hshp%extremum, 'eddy_index', &
          ftinteger, nwidth = 5, ndecimals = 0)
@@ -40,7 +40,7 @@ contains
          'r_eq_area', ftdouble, nwidth = 10, ndecimals = 4)
     call dbf_add_field_03(hshp%out_ssh, hshp%outermost, 'ssh', ftdouble, &
          nwidth = 13, ndecimals = 6)
-    call dbf_add_field_03(hshp%out_date, hshp%outermost, 'date_index', &
+    call dbf_add_field_03(hshp%out_date, hshp%outermost, 'days_1950', &
          ftinteger, nwidth = 5, ndecimals = 0)
     call dbf_add_field_03(hshp%out_eddy_index, hshp%outermost, 'eddy_index', &
          ftinteger, nwidth = 5, ndecimals = 0)
@@ -53,7 +53,7 @@ contains
          'r_eq_area', ftdouble, nwidth = 10, ndecimals = 4)
     call dbf_add_field_03(hshp%max_speed_ssh, hshp%max_speed, 'ssh', &
          ftdouble, nwidth = 13, ndecimals = 6)
-    call dbf_add_field_03(hshp%max_speed_date, hshp%max_speed, 'date_index', &
+    call dbf_add_field_03(hshp%max_speed_date, hshp%max_speed, 'days_1950', &
          ftinteger, nwidth = 5, ndecimals = 0)
     call dbf_add_field_03(hshp%max_speed_eddy_index, hshp%max_speed, &
          'eddy_index', ftinteger, nwidth = 5, ndecimals = 0)