diff --git a/dartwrf/create_obsseq.py b/dartwrf/create_obsseq.py
index 351be1a6840957c8715659a3ee7c032ceb321357..9581d664d2cf45c5e798c827291da5f485b38d7b 100755
--- a/dartwrf/create_obsseq.py
+++ b/dartwrf/create_obsseq.py
@@ -41,20 +41,14 @@ def obskind_read():
 
 #####################
 # Global variables
+obs_kind_nrs = obskind_read()  # DART internal indices
 
 # position on earth for RTTOV ray geometry
 lat0 = 45.
 lon0 = 0.
-
-# position on Earth for DART, domain center when coords_from_domaincenter=True
-lat0_center = 45.
-lon0_center = 0.
-
 sat_az = "180.0"
 sat_zen = "45.0"
 
-obs_kind_nrs = obskind_read()
-
 
 def degr_to_rad(degr):
     """Convert to DART convention (radians)
diff --git a/dartwrf/obs/calculate_obs_locations.py b/dartwrf/obs/calculate_obs_locations.py
index 9e30bde6677fc70882dff34471b080101099fdb0..d4dc65cfd8947b01d6f2468d1fa0e5ecf77e56d7 100755
--- a/dartwrf/obs/calculate_obs_locations.py
+++ b/dartwrf/obs/calculate_obs_locations.py
@@ -18,8 +18,6 @@ lon0_center = 0.
 
 radius_earth_meters = 6.371*1E6
 
-
-
 def square_array_from_domaincenter(n_obs, distance_between_obs_km):
     """
     Create equally spaced grid for satellite observations every 4 km
@@ -45,7 +43,7 @@ def square_array_from_domaincenter(n_obs, distance_between_obs_km):
             dx_4km_in_degree = distance_between_obs_meters/m_per_degree_x
             lon = lon0_center + ix*dx_4km_in_degree
             coords.append((lat, lon))
-
+    return coords
 
 def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True):
     """Observations spread evenly over domain
@@ -91,4 +89,4 @@ def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True):
         for j in range(n_obs_x):
             coords.append((lats[skip_gridpoints+i*gridpoints_between_obs, skip_gridpoints+j*gridpoints_between_obs],
                             lons[skip_gridpoints+i*gridpoints_between_obs, skip_gridpoints+j*gridpoints_between_obs]))
-
+    return coords
\ No newline at end of file