diff --git a/Analysis/plot_eddy_contours.py b/Analysis/plot_eddy_contours.py
index e99e586ca84ba5e26d69d726a14ef7f45cc10f84..a34c56d3bdc96e5193cd93a52636f15849419021 100755
--- a/Analysis/plot_eddy_contours.py
+++ b/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:
diff --git a/shp_tr_create.f90 b/shp_tr_create.f90
index f52af7b2a41025f352e5865d04c262333ab18733..b23df255141ae4db5ea5456929d7f510e2c788fe 100644
--- a/shp_tr_create.f90
+++ b/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)