Skip to content

GitLab

Commit ed7dff09 authored by ymeur's avatar ymeur
Browse files

Update of the XML files. 

YM
parent 9de4289f
Pipeline #157233 failed with stages
in 4 seconds
Hide whitespace changes
Inline Side-by-side
Showing with 247 additions and 85 deletions
TEST_CASE/LAM_ERA5/XML/nudging_dynamico.xml
View file @ ed7dff09
......@@ -5,50 +5,57 @@
<context id="icosagcm" >
<file_definition>
<file id="nudging_out" enabled="false" output_freq="1h" sync_freq="1ts" output_level="10" >
<field_group freq_offset="0ts" scale_factor="_reset_" add_offset="_reset_">
<field field_ref="ulon_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ulat_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="temp_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ps_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="p_nudged" operation="instant" scale_factor="_reset_" add_offset="_reset_"/>
</field_group>
</file>
<file id="nudging_test" enabled="true" output_freq="12ts" sync_freq="3h" output_level="10" >
<field id="coeff_i" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="ps_nudged_out" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="temp_nudged_out" domain_ref="i" axis_ref="lev" operation="instant" enabled="true" prec="8" />
<file_group id="nudging_files" enabled="false"> <!-- enabled will be set from model in nudging case -->
<file id="nudging_out" enabled="false" output_freq="1h" sync_freq="1ts" output_level="10" >
<field_group freq_offset="0ts" scale_factor="_reset_" add_offset="_reset_">
<field field_ref="ulon_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ulat_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="temp_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ps_nudged" scale_factor="_reset_" add_offset="_reset_"/>
</field_group>
</file>
<file id="nudging_test" enabled="false" output_freq="1ts" sync_freq="1ts" output_level="10" >
<field id="coeff_i" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="ps_nudged_out" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="temp_nudged_out" domain_ref="i" axis_ref="lev" operation="instant" enabled="true" prec="8" />
</file>
<file id="nudging_in" name="ERA5" enabled="true" output_freq="1h" mode="read" record_offset="1" time_counter_name="time">
<field_group operation="instant" freq_offset="1ts" >
<field id="ps_nudged_read" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
<file id="nudging_in0" name="ERA5" enabled="true" output_freq="1000y" mode="read" record_offset="0" time_counter_name="time">
<field_group operation="once" freq_offset="1ts" >
<field id="ps_nudged_read0" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read0" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read0" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read0" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read0" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read0" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read0" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read0" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
<file id="nudging_in" name="ERA5_062000" enabled="true" output_freq="1h" mode="read" record_offset="1" time_counter_name="time">
<field_group operation="instant" freq_offset="1ts" >
<field id="ps_nudged_read" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
<file id="nudging_in0" name="ERA5_062000" enabled="true" output_freq="1000y" mode="read" record_offset="0" time_counter_name="time">
<field_group operation="once" freq_offset="1ts" >
<field id="ps_nudged_read0" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read0" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read0" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read0" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read0" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read0" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read0" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read0" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
</file_group>
</file_definition>
<field_definition>
<field id="temp_nudged_horiz_interp" field_ref="temp_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
......@@ -59,11 +66,11 @@
<field id="H2O-l_nudged_horiz_interp" field_ref="H2O-l_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-s_nudged_horiz_interp" field_ref="H2O-s_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field_group id="fields_nudged" read_access="true">
<field_group id="nudging_fields" read_access="false"> <!-- read_access will be set from model -->
<field id="ps_nudged" field_ref="ps_nudged_read" domain_ref="from_nudged"/>
<!--
<field id="p_nudged_horiz_interp" field_ref="p_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged" />
-->
<field id="p_nudged_horiz_interp" field_ref="p_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged" read_access="false" /> <!-- read_access will be set from model -->
<field id="temp_nudged" field_ref="temp_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulon_nudged" field_ref="ulon_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulat_nudged" field_ref="ulat_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
......@@ -72,8 +79,9 @@
<field id="H2O-l_nudged" field_ref="H2O-l_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="H2O-s_nudged" field_ref="H2O-s_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="q_nudged" field_ref="q_nudged_horiz_interp" grid_ref="grid_q_from_nudged"/>
</field_group>
<field id="p_nudged" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
</field_group>
<field id="p_nudged" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
<field id="temp_nudged_horiz_interp0" field_ref="temp_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
......@@ -84,11 +92,11 @@
<field id="H2O-l_nudged_horiz_interp0" field_ref="H2O-l_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-s_nudged_horiz_interp0" field_ref="H2O-s_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field_group id="fields_nudged0" read_access="true">
<field_group id="nudging_fields0" read_access="false"> <!-- read_access will be set from model -->
<field id="ps_nudged0" field_ref="ps_nudged_read0" domain_ref="from_nudged"/>
<!--
<field id="p_nudged_horiz_interp0" field_ref="p_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged" />
-->
<field id="p_nudged_horiz_interp0" field_ref="p_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged" read_access="false"/> <!-- read_access will be set from model -->
<field id="temp_nudged0" field_ref="temp_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="ulon_nudged0" field_ref="ulon_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="ulat_nudged0" field_ref="ulat_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
......@@ -96,17 +104,15 @@
<field id="H2O-g_nudged0" field_ref="H2O-g_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="H2O-l_nudged0" field_ref="H2O-l_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="H2O-s_nudged0" field_ref="H2O-s_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
</field_group>
<field id="p_nudged0" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
</field_group>
<field id="p_nudged0" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
</field_definition>
<domain_definition>
<domain id="domain_nudged" type="rectilinear">
......@@ -122,35 +128,21 @@
<axis_definition>
<axis id="axis_nudged" axis_type="Z" convert_from_factor="100" /> <!-- ERA5 files are in milibar => pascal -->
<!--
<axis id="from_nudged" axis_ref="axis_nudged" />
-->
<!--
<axis id="from_nudged" axis_ref="lev">
<axis id="nudged_vert_interp_src_dest" axis_ref="lev">
<interpolate_axis order="1" coordinate_src="p_nudged_horiz_interp" coordinate_dst="p_nudged" />
</axis>
-->
<axis id="from_nudged" axis_ref="lev">
<axis id="nudged_vert_interp_dest" axis_ref="lev">
<interpolate_axis order="1" coordinate_dst="p_nudged" />
</axis>
</axis>
<!--
<axis id="from_nudged0" axis_ref="axis_nudged" />
-->
<!--
<axis id="from_nudged0" axis_ref="lev">
<interpolate_axis order="1" coordinate_src="p_nudged_horiz_interp0" coordinate_dst="p_nudged0" />
</axis>
-->
<axis id="from_nudged0" axis_ref="lev">
<interpolate_axis order="1" coordinate_dst="p_nudged0" />
</axis>
<axis id="from_nudged" axis_ref="axis_nudged" /> <!-- axis_ref will be set from model -->
<axis id="from_nudged0" axis_ref="axis_nudged" /> <!-- axis_ref will be set from model -->
</axis_definition>
<grid_definition>
......
xml/iodef.xml
View file @ ed7dff09
......@@ -22,6 +22,7 @@
<context id="icosagcm_input" src="./context_input_dynamico.xml" />
<context id="icosagcm" src="./context_dynamico.xml" />
<context id="icosagcm" src="./file_def_dynamico.xml" />
<context id="icosagcm" src="./nudging_dynamico.xml" />
<!--context id="icosagcm" src="./tracers_dynamico.xml" /-->
<!-- For running with LMDZ, following line is added by dynamicolmdz.driver while using libIGCM: -->
<!-- <context id="LMDZ" src="./context_lmdz.xml" /> -->
......
xml/iodef.xml.py
View file @ ed7dff09
......@@ -3,7 +3,6 @@ import argparse
parser = argparse.ArgumentParser(description='Prints iodef.xml')
parser.add_argument('--enable_LMDZ', action='store_true', help='Includes context_lmdz.xml' )
parser.add_argument('--enable_tracers', action='store_true', help='Includes tracers_dynamico.xml' )
parser.add_argument('--enable_nudging', action='store_true', help='Includes nuding_dynamico.xml' )
args = parser.parse_args()
if args.enable_LMDZ is None:
......@@ -19,11 +18,6 @@ if args.enable_tracers:
else:
tracers = ' <!--context id="icosagcm" src="./tracers_dynamico.xml" /--> '
if args.enable_nudging:
nudging = ' <context id="icosagcm" src="./nudging_dynamico.xml" />'
else:
nudging = ''
iodef="""
<?xml version="1.0"?>
<simulation>
......@@ -49,12 +43,12 @@ iodef="""
<context id="icosagcm_input" src="./context_input_dynamico.xml" />
<context id="icosagcm" src="./context_dynamico.xml" />
<context id="icosagcm" src="./file_def_dynamico.xml" />
<context id="icosagcm" src="./nudging_dynamico.xml" />
{tracers}
{LMDZ}
{nudging}
</simulation>
"""
print(iodef.format(tracers=tracers, LMDZ=LMDZ, nudging=nudging))
print(iodef.format(tracers=tracers, LMDZ=LMDZ))
xml/nudging_dynamico.xml.py → xml/nudging_dynamico.xml
View file @ ed7dff09
import argparse
parser = argparse.ArgumentParser(description='Prints nudging_dynamico.xml')
parser.add_argument('--input_mesh', help='Mesh type of input file',
choices=['rectilinear','unstructured'], default='rectilinear')
args = parser.parse_args()
nudging="""
<!-- =========================================================================================================== -->
<!-- nudging_dynamico.xml -->
<!-- Definition of nudging input/output -->
<!-- nudging_dynamico.xml -->
<!-- Definition of nudging input/output -->
<!-- =========================================================================================================== -->
<context id="icosagcm" >
<file_definition>
<file id="nudging_out" enabled="false" output_freq="3h" sync_freq="3h" output_level="10" >
<field field_ref="ulon"/>
<field field_ref="ulat"/>
<field field_ref="p" />
<field field_ref="temp" />
<field field_ref="ps" />
<field field_ref="q" />
</file>
<file_group id="nudging_files" enabled="false"> <!-- enabled will be set from model in nudging case -->
<file id="nudging_out" enabled="false" output_freq="1h" sync_freq="1ts" output_level="10" >
<field_group freq_offset="0ts" scale_factor="_reset_" add_offset="_reset_">
<field field_ref="ulon_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ulat_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="temp_nudged" scale_factor="_reset_" add_offset="_reset_"/>
<field field_ref="ps_nudged" scale_factor="_reset_" add_offset="_reset_"/>
</field_group>
</file>
<file id="nudging_test" enabled="false" output_freq="1ts" sync_freq="1ts" output_level="10" >
<field id="coeff_i" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="ps_nudged_out" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="temp_nudged_out" domain_ref="i" axis_ref="lev" operation="instant" enabled="true" prec="8" />
</file>
<file id="nudging_test" enabled="true" output_freq="12ts" sync_freq="3h" output_level="10" >
<field id="coeff_i" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="ps_nudged_out" domain_ref="i" operation="instant" enabled="true" prec="8"/>
<field id="temp_nudged_out" domain_ref="i" axis_ref="lev" operation="instant" enabled="true" prec="8" />
<file id="nudging_in" name="ERA5" enabled="true" output_freq="1h" mode="read" record_offset="1" time_counter_name="time">
<field_group operation="instant" freq_offset="1ts" >
<field id="ps_nudged_read" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
<file id="nudging_in0" name="ERA5" enabled="true" output_freq="1000y" mode="read" record_offset="0" time_counter_name="time">
<field_group operation="once" freq_offset="1ts" >
<field id="ps_nudged_read0" name="sp" domain_ref="domain_nudged" />
<field id="ulon_nudged_read0" name="u" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read0" name="v" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read0" name="t" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-g_nudged_read0" name="q" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-l_nudged_read0" name="clwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="H2O-s_nudged_read0" name="ciwc" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read0" name="q" grid_ref="grid_q_nudged" enabled="false"/>
</field_group>
</file>
<file id="nudging_in" name="nudging_in" enabled="true" output_freq="3h" mode="read" record_offset="1">
<field_group operation="instant" freq_offset="1ts" >
<field id="ps_nudged_read" name="ps" domain_ref="domain_nudged" />
<field id="ulon_nudged_read" name="U" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read" name="V" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="p_nudged_read" name="P" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read" name="T" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read" name="q" grid_ref="grid_q_nudged" />
</field_group>
</file>
<file id="nudging_in0" name="nudging_in" enabled="true" output_freq="1000y" mode="read" record_offset="0">
<field_group operation="once" freq_offset="1ts" >
<field id="ps_nudged_read0" name="ps" domain_ref="domain_nudged" />
<field id="ulon_nudged_read0" name="U" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="ulat_nudged_read0" name="V" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="p_nudged_read0" name="P" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged_read0" name="T" domain_ref="domain_nudged" axis_ref="axis_nudged" />
<field id="q_nudged_read0" name="q" grid_ref="grid_q_nudged" />
</field_group>
</file>
</file_group>
</file_definition>
<field_definition>
<field id="temp_nudged_horiz_interp" field_ref="temp_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="ulon_nudged_horiz_interp" field_ref="ulon_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="ulat_nudged_horiz_interp" field_ref="ulat_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="q_nudged_horiz_interp" field_ref="q_nudged_read" grid_ref="grid_q_from_nudged_horiz_interp"/>
<field id="H2O-g_nudged_horiz_interp" field_ref="H2O-g_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-l_nudged_horiz_interp" field_ref="H2O-l_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-s_nudged_horiz_interp" field_ref="H2O-s_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field_group id="fields_nudged" read_access="true">
<field_group id="nudging_fields" read_access="false"> <!-- read_access will be set from model -->
<field id="ps_nudged" field_ref="ps_nudged_read" domain_ref="from_nudged"/>
<field id="p_nudged_horiz_interp" field_ref="p_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged" />
<field id="temp_nudged" field_ref="temp_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulon_nudged" field_ref="ulon_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulat_nudged" field_ref="ulat_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="p_nudged_horiz_interp" field_ref="p_nudged_read" domain_ref="from_nudged" axis_ref="axis_nudged" read_access="false" /> <!-- read_access will be set from model -->
<field id="temp_nudged" field_ref="temp_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulon_nudged" field_ref="ulon_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulat_nudged" field_ref="ulat_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="ulat_nudged" field_ref="ulat_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="H2O-g_nudged" field_ref="H2O-g_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="H2O-l_nudged" field_ref="H2O-l_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="H2O-s_nudged" field_ref="H2O-s_nudged_horiz_interp" domain_ref="from_nudged" axis_ref="from_nudged"/>
<field id="q_nudged" field_ref="q_nudged_horiz_interp" grid_ref="grid_q_from_nudged"/>
</field_group>
<field id="p_nudged" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
</field_group>
<field id="p_nudged" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
<field id="temp_nudged_horiz_interp0" field_ref="temp_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="ulon_nudged_horiz_interp0" field_ref="ulon_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="ulat_nudged_horiz_interp0" field_ref="ulat_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="q_nudged_horiz_interp0" field_ref="q_nudged_read0" grid_ref="grid_q_from_nudged_horiz_interp"/>
<field id="H2O-g_nudged_horiz_interp0" field_ref="H2O-g_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-l_nudged_horiz_interp0" field_ref="H2O-l_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field id="H2O-s_nudged_horiz_interp0" field_ref="H2O-s_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged"/>
<field_group id="fields_nudged0" read_access="true">
<field_group id="nudging_fields0" read_access="false"> <!-- read_access will be set from model -->
<field id="ps_nudged0" field_ref="ps_nudged_read0" domain_ref="from_nudged"/>
<field id="p_nudged_horiz_interp0" field_ref="p_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged" />
<field id="p_nudged_horiz_interp0" field_ref="p_nudged_read0" domain_ref="from_nudged" axis_ref="axis_nudged" read_access="false"/> <!-- read_access will be set from model -->
<field id="temp_nudged0" field_ref="temp_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="ulon_nudged0" field_ref="ulon_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="ulat_nudged0" field_ref="ulat_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="q_nudged0" field_ref="q_nudged_horiz_interp0" grid_ref="grid_q_from_nudged0"/>
</field_group>
<field id="p_nudged0" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
<field id="H2O-g_nudged0" field_ref="H2O-g_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="H2O-l_nudged0" field_ref="H2O-l_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
<field id="H2O-s_nudged0" field_ref="H2O-s_nudged_horiz_interp0" domain_ref="from_nudged" axis_ref="from_nudged0"/>
</field_group>
<field id="p_nudged0" domain_ref="i" axis_ref="lev" /> <!--sent from model-->
</field_definition>
<domain_definition>
<domain id="domain_nudged" type="{input_mesh}">
<domain id="domain_nudged" type="rectilinear">
<generate_rectilinear_domain/>
</domain>
......@@ -103,16 +126,23 @@ nudging="""
</domain_definition>
<axis_definition>
<axis id="axis_nudged"/>
<axis id="from_nudged" axis_ref="lev">
<interpolate_axis order="1" coordinate_src="p_nudged_horiz_interp" coordinate_dst="p_nudged" extrapolate="true"/>
</axis>
<axis id="from_nudged0" axis_ref="lev">
<interpolate_axis order="1" coordinate_src="p_nudged_horiz_interp0" coordinate_dst="p_nudged0" extrapolate="true"/>
<axis id="axis_nudged" axis_type="Z" convert_from_factor="100" /> <!-- ERA5 files are in milibar => pascal -->
<axis id="nudged_vert_interp_src_dest" axis_ref="lev">
<interpolate_axis order="1" coordinate_src="p_nudged_horiz_interp" coordinate_dst="p_nudged" />
</axis>
<axis id="nudged_vert_interp_dest" axis_ref="lev">
<interpolate_axis order="1" coordinate_dst="p_nudged" />
</axis>
<axis id="from_nudged" axis_ref="axis_nudged" /> <!-- axis_ref will be set from model -->
<axis id="from_nudged0" axis_ref="axis_nudged" /> <!-- axis_ref will be set from model -->
</axis_definition>
<grid_definition>
......@@ -143,6 +173,3 @@ nudging="""
</grid_definition>
</context>
"""
print(nudging.format(input_mesh=args.input_mesh))
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment