diff options
-rw-r--r-- | interface.ccl | 8 | ||||
-rw-r--r-- | schedule.ccl | 1 | ||||
-rw-r--r-- | src/setup.c | 8 |
3 files changed, 17 insertions, 0 deletions
diff --git a/interface.ccl b/interface.ccl index 363489f..03af59d 100644 --- a/interface.ccl +++ b/interface.ccl @@ -69,3 +69,11 @@ CCTK_REAL sf_coordinate_descriptors[nsurfaces] TYPE=scalar sf_origin_theta sf_origin_phi sf_delta_theta sf_delta_phi } "Surface coordinate descriptors" + + + +CCTK_INT sf_shape_descriptors[nsurfaces] TYPE=scalar +{ + sf_ntheta sf_nphi + sf_nghoststheta sf_nghostsphi +} "Surface shape descriptors" diff --git a/schedule.ccl b/schedule.ccl index 0e28d1d..b7feabc 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -5,6 +5,7 @@ STORAGE: sf_valid STORAGE: sf_info STORAGE: sf_radius sf_origin STORAGE: sf_coordinate_descriptors +STORAGE: sf_shape_descriptors diff --git a/src/setup.c b/src/setup.c index 982bbe9..f30d1ce 100644 --- a/src/setup.c +++ b/src/setup.c @@ -50,6 +50,14 @@ void SphericalSurfaceInfo_Setup (CCTK_ARGUMENTS) + /* copy parameters into grid functions */ + sf_ntheta[n] = ntheta[n]; + sf_nphi[n] = nphi[n]; + sf_nghoststheta[n] = nghoststheta[n]; + sf_nghostsphi[n] = nghostsphi[n]; + + + /* coordinates in the theta direction */ /* avoid_sf_origin_theta = 1 */ if (symmetric_z[n]) { |