Skip to content

GitLab

Commit 71b9bed8 authored by ymeur's avatar ymeur
Browse files

LAM ERA5 configuration example 

YM
parent ac3dec79
Pipeline #156109 failed with stages
in 4 minutes and 53 seconds
Hide whitespace changes
Inline Side-by-side
Showing with 793 additions and 0 deletions
TEST_CASE/LAM_ERA5/CREATE_ETAT0/context_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/context_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/context_input_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/context_input_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/earth_const.def 0 → 100644
View file @ 71b9bed8
#planetary radius : real (default=37122E6)
radius=6.371229E6
# gravity : real (default=9.80616)
g=9.80665
#planetary rotation rate (default=7.292E-5)
omega=7.292E-5
# kappa=Rd/cpp (default=0.2857143)
kappa=0.2857143
#kappa=0.285773692052
# Heat capacity at constant pressure real (default=1004.70885)
#cpp=1004.5
cpp= 1004.70885783307
# Reference pressure : real (default=101325.)
preff=101325.
# Scale factor for small planet experiment : real (default=1.)
scale_factor=1.
# Precision of NetCDF output
netcdf_prec = double
# solarlong0 = -9999.999 for normal run
# solarlong0 = 1000 for aquaplanet run
#solarlong0 = _AUTO_ : DEFAULT = -9999.999
TEST_CASE/LAM_ERA5/CREATE_ETAT0/field_def_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/field_def_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/file_def_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/file_def_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/iodef.xml 0 → 120000
View file @ 71b9bed8
../XML/iodef.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/nudging_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/nudging_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/CREATE_ETAT0/run.def 0 → 100644
View file @ 71b9bed8
# $Id$
#
INCLUDEDEF=run_dynamico.def
INCLUDEDEF=earth_const.def
TEST_CASE/LAM_ERA5/CREATE_ETAT0/run_dynamico.def 0 → 100644
View file @ 71b9bed8
mpi_threading_mode=funneled
# ------------------------------- Mesh ---------------------------------
metric_type = icosa_area
area_center_lon = 3
area_center_lat = 47
area_radius=800
#---------------- guided ----------------
#guided_type = none
guided_type = nudging
itau_nudging = 120
nudging_zone = area
# Number of subdivision on a main triangle (nbp) : integer (default=40)
nbp=40
# nbp 20 40 80 160
# T-edge length (km) 500 250 120 60
# sub splitting of main rhombus : integer (default=1)
nsplit_i=2
nsplit_j=2
omp_level_size=1
# read_metric : read grid description from start.nc file
read_metric= n
# optim_it : mesh optimisation : number of iteration : integer (default=0)
optim_it= 0
# ---------------------------------- Time---------------------------------
# timestep : real (default=480s)
dt = 30.
# advection called every itau_adv time steps : integer (default 2)
# standard : umax=100m/s vs c=340m/s (ratio 1:3)
# in JW06 umax=35m/s vs c=340m/s (ratio 1:10)
itau_adv=1
# run length : real (default=dt*itaumax)
run_length= 900
# etat0_start_iteration_reset : set y to reset iteration from start.nc
etat0_start_iteration_reset= y
# output field period : integer (default none)
write_period=86400
itau_check_conserv=100
# ---------------------------------- Misc --------------------------------
# number of tracer (nqtot) : integer (default 1)
nqtot=4
# pression value where output is interpolated : real (default=0, no output)
out_pression_level=85000
start_file_name=start
restart_file_name=restart
# Use of oasis or not (needed to initialize and finalize coupling mode)
using_oasis= n
# ---------------------------------- Dissipation ---------------------------------
# dissipation time graddiv : real (default=5000)
tau_graddiv = 3600
# number of iteration for graddiv : integer (default=1)
nitergdiv = 1
# dissipation time nxgradrot (default=5000)
tau_gradrot = 21600
# number of iteration for nxgradrot : integer (default=1)
nitergrot=2
# dissipation time divgrad (theta) (default=5000)
tau_divgrad=10800
# number of iteration for divgrad : integer (default=1)
niterdivgrad=2
# -------------------------------- Physics -------------------------------------
# kind of physics : string : none, dcmip (default=none)
#physics=phys_external
physics=none
# itau_physics: dt_phys=itau_physics*dt here dt=450s => dt_phys=900s
itau_physics=30
# Following parameters are changed by options set in dynamico.card
create_etat0_limit= y
etat0= database
#etat0 = start_file
iflag_phys= 1
# heldsz_p0 is only read if etat0=held_suarez (aquaplanet simulation)
heldsz_p0=101080
TEST_CASE/LAM_ERA5/ERA_062000.nc_header 0 → 100644
View file @ 71b9bed8
netcdf ERA5_062000 {
dimensions:
time = 720 ;
level = 37 ;
latitude = 69 ;
longitude = 89 ;
variables:
short ciwc(time, level, latitude, longitude) ;
ciwc:scale_factor = 1.32719542301207e-08 ;
ciwc:add_offset = 0.000434868852304135 ;
ciwc:_FillValue = -32767s ;
ciwc:missing_value = -32767s ;
ciwc:units = "kg kg**-1" ;
ciwc:long_name = "Specific cloud ice water content" ;
short clwc(time, level, latitude, longitude) ;
clwc:scale_factor = 2.21635814610459e-08 ;
clwc:add_offset = 0.000726211910152629 ;
clwc:_FillValue = -32767s ;
clwc:missing_value = -32767s ;
clwc:units = "kg kg**-1" ;
clwc:long_name = "Specific cloud liquid water content" ;
float latitude(latitude) ;
latitude:units = "degrees_north" ;
latitude:long_name = "latitude" ;
int level(level) ;
level:units = "millibars" ;
level:long_name = "pressure_level" ;
float longitude(longitude) ;
longitude:units = "degrees_east" ;
longitude:long_name = "longitude" ;
short q(time, level, latitude, longitude) ;
q:scale_factor = 2.64000697934743e-07 ;
q:add_offset = 0.00864792981254272 ;
q:_FillValue = -32767s ;
q:missing_value = -32767s ;
q:units = "kg kg**-1" ;
q:long_name = "Specific humidity" ;
q:standard_name = "specific_humidity" ;
short r(time, level, latitude, longitude) ;
r:scale_factor = 0.00229625616917222 ;
r:add_offset = 67.7087207798439 ;
r:_FillValue = -32767s ;
r:missing_value = -32767s ;
r:units = "%" ;
r:long_name = "Relative humidity" ;
r:standard_name = "relative_humidity" ;
short sp(time, latitude, longitude) ;
sp:scale_factor = 0.42606110661804 ;
sp:add_offset = 89552.8182194467 ;
sp:_FillValue = -32767s ;
sp:missing_value = -32767s ;
sp:units = "Pa" ;
sp:long_name = "Surface pressure" ;
sp:standard_name = "surface_air_pressure" ;
short t(time, level, latitude, longitude) ;
t:scale_factor = 0.00168413112723563 ;
t:add_offset = 257.49478629137 ;
t:_FillValue = -32767s ;
t:missing_value = -32767s ;
t:units = "K" ;
t:long_name = "Temperature" ;
t:standard_name = "air_temperature" ;
short t2m(time, latitude, longitude) ;
t2m:scale_factor = 0.000684527350485633 ;
t2m:add_offset = 290.570580587887 ;
t2m:_FillValue = -32767s ;
t2m:missing_value = -32767s ;
t2m:units = "K" ;
t2m:long_name = "2 metre temperature" ;
int time(time) ;
time:units = "hours since 1900-01-01 00:00:00.0" ;
time:long_name = "time" ;
time:calendar = "gregorian" ;
short u(time, level, latitude, longitude) ;
u:scale_factor = 0.00179951866151881 ;
u:add_offset = 11.6168523158646 ;
u:_FillValue = -32767s ;
u:missing_value = -32767s ;
u:units = "m s**-1" ;
u:long_name = "U component of wind" ;
u:standard_name = "eastward_wind" ;
short v(time, level, latitude, longitude) ;
v:scale_factor = 0.00187790934275837 ;
v:add_offset = -2.20490318806982 ;
v:_FillValue = -32767s ;
v:missing_value = -32767s ;
v:units = "m s**-1" ;
v:long_name = "V component of wind" ;
v:standard_name = "northward_wind" ;
short z(time, latitude, longitude) ;
z:scale_factor = 0.374993085544687 ;
z:add_offset = 11500.6079136135 ;
z:_FillValue = -32767s ;
z:missing_value = -32767s ;
z:units = "m**2 s**-2" ;
z:long_name = "Geopotential" ;
z:standard_name = "geopotential" ;
// global attributes:
:Conventions = "CF-1.6" ;
:history = "Fri Jan 7 11:03:02 2022: ncks --fix_rec_dmn time ERA5_062000.nc ERA5_fix_062000.nc\nFri Dec 17 12:37:36 2021: ncks -4 --mk_rec_dmn time ERA5_062000.nc ERA5_062000.nc4\nFri Dec 17 12:35:51 2021: ncks -A ERA5_QR_062000.nc ERA5_062000.nc\n2021-12-15 11:27:27 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data9/adaptor.mars.internal-1639567586.8756938-8273-3-26f80432-6499-4aae-a0ef-9673aaf54ec8.nc /cache/tmp/26f80432-6499-4aae-a0ef-9673aaf54ec8-adaptor.mars.internal-1639565273.4298384-8273-3-tmp.grib" ;
:history_of_appended_files = "Fri Dec 17 12:35:51 2021: Appended file ERA5_QR_062000.nc had following \"history\" attribute:\n2021-12-15 11:27:27 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data9/adaptor.mars.internal-1639567586.8756938-8273-3-26f80432-6499-4aae-a0ef-9673aaf54ec8.nc /cache/tmp/26f80432-6499-4aae-a0ef-9673aaf54ec8-adaptor.mars.internal-1639565273.4298384-8273-3-tmp.grib\nFri Dec 17 12:34:56 2021: Appended file ERA5_QR_062000.nc had following \"history\" attribute:\n2021-12-15 11:27:27 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data9/adaptor.mars.internal-1639567586.8756938-8273-3-26f80432-6499-4aae-a0ef-9673aaf54ec8.nc /cache/tmp/26f80432-6499-4aae-a0ef-9673aaf54ec8-adaptor.mars.internal-1639565273.4298384-8273-3-tmp.grib\nFri Dec 17 11:11:14 2021: Appended file ERA5_Q_062000.nc had following \"history\" attribute:\n2021-12-14 18:20:41 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data3/adaptor.mars.internal-1639505797.8218899-16873-6-11757d01-baee-4b75-9710-8acd06e29d85.nc /cache/tmp/11757d01-baee-4b75-9710-8acd06e29d85-adaptor.mars.internal-1639500782.495831-16873-7-tmp.grib\nFri Dec 17 11:11:00 2021: Appended file ERA5_PHIS_062000.nc had following \"history\" attribute:\n2021-12-14 18:16:35 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data9/adaptor.mars.internal-1639505793.7545817-20389-2-e322aeb3-8b88-4a81-a2f3-975ddd9ad6fb.nc /cache/tmp/e322aeb3-8b88-4a81-a2f3-975ddd9ad6fb-adaptor.mars.internal-1639505734.8771975-20389-3-tmp.grib\nFri Dec 17 11:10:52 2021: Appended file ERA5_TS_062000.nc had following \"history\" attribute:\n2021-12-14 17:42:36 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data6/adaptor.mars.internal-1639503752.642435-18115-10-f0920f01-e4f3-4a0c-b9cb-9562589ddd4f.nc /cache/tmp/f0920f01-e4f3-4a0c-b9cb-9562589ddd4f-adaptor.mars.internal-1639503641.6824787-18115-3-tmp.grib\nFri Dec 17 11:10:17 2021: Appended file ERA5_PS_062000.nc had following \"history\" attribute:\n2021-12-14 16:58:32 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data7/adaptor.mars.internal-1639501110.8745177-19594-13-52d5faeb-d077-41bb-8a28-59b04f8aeac7.nc /cache/tmp/52d5faeb-d077-41bb-8a28-59b04f8aeac7-adaptor.mars.internal-1639501054.0349972-19594-13-tmp.grib\nFri Dec 17 11:09:36 2021: Appended file ERA5_UVT_062000.nc had following \"history\" attribute:\n2021-12-13 12:00:05 GMT by grib_to_netcdf-2.23.0: /opt/ecmwf/mars-client/bin/grib_to_netcdf -S param -o /cache/data2/adaptor.mars.internal-1639396572.7511187-9755-14-b351a23c-bccc-4f09-aed7-ccb19f5c59cc.nc /cache/tmp/b351a23c-bccc-4f09-aed7-ccb19f5c59cc-adaptor.mars.internal-1639393619.0348225-9755-15-tmp.grib\n" ;
:NCO = "netCDF Operators version 4.9.1 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)" ;
}
TEST_CASE/LAM_ERA5/RUN/context_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/context_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/context_input_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/context_input_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/earth_const.def 0 → 100644
View file @ 71b9bed8
#planetary radius : real (default=37122E6)
radius=6.371229E6
# gravity : real (default=9.80616)
g=9.80665
#planetary rotation rate (default=7.292E-5)
omega=7.292E-5
# kappa=Rd/cpp (default=0.2857143)
kappa=0.2857143
#kappa=0.285773692052
# Heat capacity at constant pressure real (default=1004.70885)
#cpp=1004.5
cpp= 1004.70885783307
# Reference pressure : real (default=101325.)
preff=101325.
# Scale factor for small planet experiment : real (default=1.)
scale_factor=1.
# Precision of NetCDF output
netcdf_prec = double
# solarlong0 = -9999.999 for normal run
# solarlong0 = 1000 for aquaplanet run
#solarlong0 = _AUTO_ : DEFAULT = -9999.999
TEST_CASE/LAM_ERA5/RUN/field_def_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/field_def_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/file_def_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/file_def_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/iodef.xml 0 → 120000
View file @ 71b9bed8
../XML/iodef.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/nudging_dynamico.xml 0 → 120000
View file @ 71b9bed8
../XML/nudging_dynamico.xml
\ No newline at end of file
TEST_CASE/LAM_ERA5/RUN/run.def 0 → 100644
View file @ 71b9bed8
# $Id$
#
INCLUDEDEF=run_dynamico.def
INCLUDEDEF=earth_const.def
TEST_CASE/LAM_ERA5/RUN/run_dynamico.def 0 → 100644
View file @ 71b9bed8
mpi_threading_mode=funneled
# ------------------------------- Mesh ---------------------------------
metric_type = icosa_area
area_center_lon = 3
area_center_lat = 47
area_radius=800
#---------------- guided ----------------
#guided_type = none
guided_type = nudging
itau_nudging = 120
nudging_zone = area
# Number of subdivision on a main triangle (nbp) : integer (default=40)
nbp=40
# nbp 20 40 80 160
# T-edge length (km) 500 250 120 60
# sub splitting of main rhombus : integer (default=1)
nsplit_i=4
nsplit_j=4
omp_level_size=1
# read_metric : read grid description from start.nc file
read_metric= n
# optim_it : mesh optimisation : number of iteration : integer (default=0)
optim_it= 0
# ---------------------------------- Time---------------------------------
# timestep : real (default=480s)
dt = 30.
# advection called every itau_adv time steps : integer (default 2)
# standard : umax=100m/s vs c=340m/s (ratio 1:3)
# in JW06 umax=35m/s vs c=340m/s (ratio 1:10)
itau_adv=1
# run length : real (default=dt*itaumax)
#run_length=864000
run_length=2592000
# etat0_start_iteration_reset : set y to reset iteration from start.nc
etat0_start_iteration_reset= y
# output field period : integer (default none)
write_period=10800
itau_check_conserv=100
# ---------------------------------- Misc --------------------------------
# number of tracer (nqtot) : integer (default 1)
nqtot=4
# pression value where output is interpolated : real (default=0, no output)
out_pression_level=85000
start_file_name=start
restart_file_name=restart
# Use of oasis or not (needed to initialize and finalize coupling mode)
using_oasis= n
# ---------------------------------- Dissipation ---------------------------------
# dissipation time graddiv : real (default=5000)
tau_graddiv = 360
# number of iteration for graddiv : integer (default=1)
nitergdiv = 1
# dissipation time nxgradrot (default=5000)
tau_gradrot = 2160
# number of iteration for nxgradrot : integer (default=1)
nitergrot=2
# dissipation time divgrad (theta) (default=5000)
tau_divgrad=1080
# number of iteration for divgrad : integer (default=1)
niterdivgrad=2
# -------------------------------- Physics -------------------------------------
# kind of physics : string : none, dcmip (default=none)
#physics=phys_external
physics=none
# itau_physics: dt_phys=itau_physics*dt here dt=450s => dt_phys=900s
itau_physics=15
# Following parameters are changed by options set in dynamico.card
#create_etat0_limit= y
#etat0= database
etat0 = start_file
iflag_phys= 1
# heldsz_p0 is only read if etat0=held_suarez (aquaplanet simulation)
heldsz_p0=101080
TEST_CASE/LAM_ERA5/XML/context_dynamico.xml 0 → 100644
View file @ 71b9bed8
<!-- =========================================================================================================== -->
<!-- DYNAMICO context -->
<!-- context_dynamico.xml : Configuration file for DYNAMICO for reading and writing files using XIOS -->
<!-- -->
<!-- =========================================================================================================== -->
<context id="icosagcm" >
<calendar type="d360" start_date="2000-01-01 00:00:00" time_origin="2000-01-01 00:00:00"/>
<!-- =========================================================================================================== -->
<!-- Definition of all existing variables in field_def_dynamico.xml -->
<!-- =========================================================================================================== -->
<field_definition src="./field_def_dynamico.xml" />
<!-- =========================================================================================================== -->
<!-- Definition of restart files and other input files -->
<!-- =========================================================================================================== -->
<file_definition type="one_file" par_access="collective">
<file id="output" name="output" enabled=".FALSE." output_freq="2h" sync_freq="2h" output_level="10" >
<field_group group_ref="standard_output"/>
</file>
<!-- =========================================================================================================== -->
<!-- Writing restart file -->
<!-- =========================================================================================================== -->
<file id="restart" enabled="true">
<field_group id="group_restart" operation="once">
<field id="it_restart" name="iteration" grid_ref="scalar" />
<field id="q_restart" name="q" grid_ref="grid_q" prec="8"/>
<field_group domain_ref="i" prec="8">
<field id="xcell_restart" name="xcell" />
<field id="ycell_restart" name="ycell" />
<field id="zcell_restart" name="zcell" />
<field id="ps_restart" name="ps" />
<field id="phis_restart" name="phis" />
<field_group axis_ref="lev" >
<field id="theta_rhodz_restart" name="theta_rhodz" />
<field id="u_restart" name="u" domain_ref="u"/>
<field id="ulon_restart" name="ulon" />
<field id="ulat_restart" name="ulat" />
<field id="p_restart" name="p" />
<field_group id="tracers_restart">
<!-- tracers will be add atomatically here-->
</field_group>
</field_group>
<field_group axis_ref="levp1" >
<field id="geopot_restart" name="geopot" />
<field id="W_restart" name="W" />
</field_group>
</field_group>
</field_group>
</file>
<!-- =========================================================================================================== -->
<!-- Writing initial state file (when not previously read from start file) -->
<!-- =========================================================================================================== -->
<file id="start0" mode="write" enabled="false">
<field_group id="group_start0" operation="once">
<field id="it_start0" name="iteration" grid_ref="scalar" />
<field id="q_start0" name="q" grid_ref="grid_q" prec="8"/>
<field_group domain_ref="i" prec="8">
<field id="xcell_start0" name="xcell" />
<field id="ycell_start0" name="ycell" />
<field id="zcell_start0" name="zcell" />
<field id="ps_start0" name="ps" />
<field id="phis_start0" name="phis" />
<field_group axis_ref="lev" >
<field id="theta_rhodz_start0" name="theta_rhodz" />
<field id="u_start0" name="u" domain_ref="u"/>
<field id="ulon_start0" name="ulon" />
<field id="ulat_start0" name="ulat" />
<field id="p_start0" name="p" />
<field_group id="tracers_start0">
<!-- tracers will be add atomatically here-->
</field_group>
</field_group>
</field_group>
</field_group>
</file>
<!-- =========================================================================================================== -->
<!-- reading start file -->
<!-- =========================================================================================================== -->
<file id="start" enabled="false" mode="read">
<field_group id="group_start" operation="once">
<field id="it_start" name="iteration" grid_ref="scalar" />
<field_group domain_ref="i" prec="8">
<field id="xcell_start" name="xcell" />
<field id="ycell_start" name="ycell" />
<field id="zcell_start" name="zcell" />
<field id="ps_start" name="ps" />
<field id="phis_start" name="phis" />
<field_group axis_ref="lev" >
<field id="theta_rhodz_start" name="theta_rhodz" />
<field id="u_start" name="u" domain_ref="u"/>
<field id="ulon_start" name="ulon" />
<field id="ulat_start" name="ulat" />
<field_group id="tracers_start">
<!-- tracers will be add atomatically here-->
</field_group>
</field_group>
<field_group axis_ref="levp1" enabled="false">
<field id="geopot_start" name="geopot" />
<field id="W_start" name="W" />
</field_group>
</field_group>
<field id="q_start" name="q" grid_ref="grid_q" prec="8"/>
</field_group>
</file>
<!-- =========================================================================================================== -->
<!-- reading start file and perform interpolation (restarting from a different resolution) -->
<!-- =========================================================================================================== -->
<file id="start_interp" enabled="false" mode="read">
<field_group operation="once">
<field id="it_start_interp" name="iteration" grid_ref="scalar" />
<field_group domain_ref="i_start_interp_read" prec="8">
<field id="ps_start_interp_read" name="ps" />
<field id="phis_start_interp_read" name="phis" />
<field_group axis_ref="lev_start_interp_read" >
<field id="theta_rhodz_start_interp_read" name="theta_rhodz" />
<field id="ulon_start_interp_read" name="ulon" />
<field id="ulat_start_interp_read" name="ulat" />
<field id="p_start_interp_read" name="p" />
<field_group id="tracers_start_interp_read">
<!-- tracers will be add atomatically here-->
</field_group>
</field_group>
</field_group>
</field_group>
</file>
</file_definition>
<field_definition>
<field_group id="start_interp">
<field_group domain_ref="i_start_interp" >
<field id="ps_start_interp" field_ref="ps_start_interp_read" />
<field id="phis_start_interp" field_ref="phis_start_interp_read" />
<field id="p_start_interp_horiz" field_ref="p_start_interp_read" axis_ref="lev_start_interp_read"/>
<field_group axis_ref="lev_start_interp" >
<field id="theta_rhodz_start_interp" field_ref="theta_rhodz_start_interp_read" />
<field id="ulon_start_interp" field_ref="ulon_start_interp_read" />
<field id="ulat_start_interp" field_ref="ulat_start_interp_read" />
<field_group id="tracers_start_interp">
<!-- tracers will be add atomatically here-->
</field_group>
</field_group>
</field_group>
</field_group>
<field id="p_start_interp" domain_ref="i" axis_ref="lev" />
</field_definition>
<domain_definition>
<domain id="i_start_interp_read" type="unstructured" >
<generate_rectilinear_domain/>
</domain>