Skip to content
Snippets Groups Projects
Commit 193d25bb authored by lkugler's avatar lkugler
Browse files

fix syntax

parent aa9e2632
No related branches found
No related tags found
No related merge requests found
...@@ -41,20 +41,14 @@ def obskind_read(): ...@@ -41,20 +41,14 @@ def obskind_read():
##################### #####################
# Global variables # Global variables
obs_kind_nrs = obskind_read() # DART internal indices
# position on earth for RTTOV ray geometry # position on earth for RTTOV ray geometry
lat0 = 45. lat0 = 45.
lon0 = 0. 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_az = "180.0"
sat_zen = "45.0" sat_zen = "45.0"
obs_kind_nrs = obskind_read()
def degr_to_rad(degr): def degr_to_rad(degr):
"""Convert to DART convention (radians) """Convert to DART convention (radians)
......
...@@ -18,8 +18,6 @@ lon0_center = 0. ...@@ -18,8 +18,6 @@ lon0_center = 0.
radius_earth_meters = 6.371*1E6 radius_earth_meters = 6.371*1E6
def square_array_from_domaincenter(n_obs, distance_between_obs_km): def square_array_from_domaincenter(n_obs, distance_between_obs_km):
""" """
Create equally spaced grid for satellite observations every 4 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): ...@@ -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 dx_4km_in_degree = distance_between_obs_meters/m_per_degree_x
lon = lon0_center + ix*dx_4km_in_degree lon = lon0_center + ix*dx_4km_in_degree
coords.append((lat, lon)) coords.append((lat, lon))
return coords
def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True): def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True):
"""Observations spread evenly over domain """Observations spread evenly over domain
...@@ -91,4 +89,4 @@ def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True): ...@@ -91,4 +89,4 @@ def evenly_on_grid(n_obs, omit_covloc_radius_on_boundary=True):
for j in range(n_obs_x): for j in range(n_obs_x):
coords.append((lats[skip_gridpoints+i*gridpoints_between_obs, skip_gridpoints+j*gridpoints_between_obs], 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])) lons[skip_gridpoints+i*gridpoints_between_obs, skip_gridpoints+j*gridpoints_between_obs]))
return coords
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment