Changeset 2285

Timestamp:

01/10/22 11:29:46 (3 years ago)

Author:

ymipsl

Message:

Fix for vertical interpolation :
fixed pressure level to dynamic pressure level (i.e. for sigma hybrid coordinate)

YM

File:

• XIOS/trunk/src/transformation/axis_algorithm_interpolate.cpp (modified) (3 diffs)

XIOS/trunk/src/transformation/axis_algorithm_interpolate.cpp

@@ -375 +375 @@
   
   
-  if (!has_src)
+  if (!has_src && !has_dst)
   {
     vecAxisValue.resize(1);
@@ -385 +385 @@
   else 
   {
-    CField* field = CField::get(coordinate_);
+    CField* field ;
+    if (has_src) field = CField::get(coordinate_);
+    else field = CField::get(coordinateDST_);
     CGrid* grid = field->grid;
 
@@ -453 +455 @@
           if (iteIndex != itIndex)
           {
-            vecAxisValue[indexMask](axisValueSize) = (*dataAuxInputs[0])(itIndex->second);
+            if (has_src)
+            {
+              vecAxisValue[indexMask](axisValueSize) = (*dataAuxInputs[0])(itIndex->second);
+            }
+            else
+            {
+              vecAxisValue[indexMask](axisValueSize) = axisSrc_-> value(axisValueSize) ;
+            }
             ++axisValueSize;
           }