Add functionality to output GRF variogram data to file

Author: oddvarlia

src/aps/algorithms/APSModel.py

@@ -1140,6 +1140,33 @@ def getAllGaussFieldNamesUsed(self) -> List[str]:
                     gfAllZones.append(gf)
         return copy.copy(gfAllZones)
 
+    def getAllGaussFieldVariogramsUsed(self, job_name: str) -> str:
+        content = job_name + '\n'
+        for key, zoneModel in self.__zoneModelTable.items():
+            (zone_number, region_number) = key
+
+            gfNames = zoneModel.used_gaussian_field_names
+            for gfname in gfNames:
+                if self.__debug_level >= Debug.VERY_VERBOSE:
+                    print(
+                        f'--- Get variogram data for zone:{zone_number} region:{region_number} for {gfname}'
+                    )
+                content += (
+                    str(zone_number) + ' ' + str(region_number) + '  ' + gfname + '  '
+                )
+                variogram_type = zoneModel.getVariogramType(gfname)
+                main_range = zoneModel.getMainRange(gfname)
+                perp_range = zoneModel.getPerpRange(gfname)
+                main_range = zoneModel.getMainRange(gfname)
+                vert_range = zoneModel.getVertRange(gfname)
+                azimuth = zoneModel.getAzimuthAngle(gfname)
+                content += str(variogram_type) + '  '
+                content += str(main_range) + '  '
+                content += str(perp_range) + '  '
+                content += str(vert_range) + '  '
+                content += str(azimuth) + '\n'
+        return content
+
     @property
     def zone_parameter(self) -> str:
         return self.getZoneParamName()
@@ -1462,6 +1489,7 @@ def write_model(
         attributes_file_name: Optional[FilePath] = None,
         param_file_name: Optional[FilePath] = None,
         probability_distribution_file_name: Optional[FilePath] = None,
+        variogram_output_file_name: Optional[FilePath] = None,
         current_job_name: Optional[str] = None,
         debug_level: Debug = Debug.OFF,
     ) -> None:
@@ -1473,6 +1501,16 @@ def write_model(
         top = ET.Element('APSModel', {'version': self.__aps_model_version})
         root_updated = self.XMLAddElement(top, fmu_attributes)
         write_string_to_file(model_file_name, root_updated, debug_level=debug_level)
+        if variogram_output_file_name is not None:
+            if current_job_name is None:
+                current_job_name = 'apsgui_job_name'
+            content = self.getAllGaussFieldVariogramsUsed(current_job_name)
+            write_string_to_file(
+                variogram_output_file_name,
+                content,
+                debug_level=debug_level,
+            )
+
         if len(fmu_attributes) == 0:
             return
         if attributes_file_name is not None:

src/aps/rms_jobs/export_fmu_config_files.py

@@ -41,7 +41,7 @@ def run(project, **kwargs):
         prob_name = job_name + aps_config.ert_config_probs_extension()
         ert_field_name = job_name + aps_config.ert_config_fields_extension()
         aps_fmuconfig_name = job_name + aps_config.fmu_master_config_extension()
-
+        variogram_file_name = job_name + '_variograms.txt'
         global_variables_file_path = aps_config.global_variables_file()
 
         # APS model file
@@ -72,6 +72,10 @@ def run(project, **kwargs):
             aps_config.ert_model_dir_absolute() + '/' + ert_field_name
         )
 
+        variogram_output_file_name = (
+            aps_config.ert_model_dir_absolute() + '/' + variogram_file_name
+        )
+
         # Check that default directory exist for export of APS model file
         export_file_full_path = Path(model_file_name).absolute()
         export_file_dir = export_file_full_path.parent
@@ -108,6 +112,10 @@ def run(project, **kwargs):
                 print(
                     f'-- Use in ERT distribution config:       {probability_distribution_file_name}'
                 )
+            if variogram_output_file_name:
+                print(
+                    f'-- Use optionally to specify localization ranges in ERT:       {variogram_output_file_name}'
+                )
             if aps_model.fmu_mode == 'FIELDS':
                 print(
                     f'-- Use in ERT main config:               {ert_field_keyword_file_name}'
@@ -119,6 +127,7 @@ def run(project, **kwargs):
             attributes_file_name=param_file_name_alternative,
             param_file_name=param_file_name,
             probability_distribution_file_name=probability_distribution_file_name,
+            variogram_output_file_name=variogram_output_file_name,
             current_job_name=job_name,
             debug_level=debug_level,
         )