fixed bug in interpolation procedure (spherical vector)

Author: mo-lormi

src/atlas/interpolation/method/sphericalvector/SphericalVector.cc

@@ -150,14 +150,14 @@ void SphericalVector::do_execute(const Field& sourceField, Field& targetField,
   ATLAS_TRACE("atlas::interpolation::method::SphericalVector::do_execute()");
 
   if (targetField.size() == 0) {
+    haloExchange(sourceField);
     return;
   }
 
   const auto fieldType = sourceField.metadata().getString("type", "");
   if (fieldType != "vector") {
     auto metadata = Metadata();
     Method::do_execute(sourceField, targetField, metadata);
-
     return;
   }
 
@@ -194,14 +194,14 @@ void SphericalVector::do_execute_adjoint(Field& sourceField,
       "atlas::interpolation::method::SphericalVector::do_execute_adjoint()");
 
   if (targetField.size() == 0) {
+    adjointHaloExchange(sourceField);
     return;
   }
 
   const auto fieldType = sourceField.metadata().getString("type", "");
   if (fieldType != "vector") {
     auto metadata = Metadata();
     Method::do_execute_adjoint(sourceField, targetField, metadata);
-
     return;
   }