diff --git a/previews/PR1080/.documenter-siteinfo.json b/previews/PR1080/.documenter-siteinfo.json index 72a2b38835..06fd684834 100644 --- a/previews/PR1080/.documenter-siteinfo.json +++ b/previews/PR1080/.documenter-siteinfo.json @@ -1 +1 @@ -{"documenter":{"julia_version":"1.9.4","generation_timestamp":"2024-05-17T15:32:39","documenter_version":"1.3.0"}} \ No newline at end of file +{"documenter":{"julia_version":"1.9.4","generation_timestamp":"2024-05-17T15:39:17","documenter_version":"1.3.0"}} \ No newline at end of file diff --git a/previews/PR1080/Boundaries/BoundaryFilePreparation/index.html b/previews/PR1080/Boundaries/BoundaryFilePreparation/index.html index 7eac3fcf90..271c287d1f 100644 --- a/previews/PR1080/Boundaries/BoundaryFilePreparation/index.html +++ b/previews/PR1080/Boundaries/BoundaryFilePreparation/index.html @@ -51,4 +51,4 @@ # hh_offset is 0 ; DTG is SURFEX_INI| /nobackup/smhid9/sm_esbol/hm_home/ice_36h1_4/20110901_12/SURFXINI.lfi 000|2011090112 /nobackup/smhid9/sm_esbol/hm_home/ice_36h1_4/20110901_12/ELSCFHARMALBC000 /nobackup/smhid9/sm_esbol/hm_home/ice_36h1_4/g05a/archive/2011/09/01/12/fc20110901_12+000 scp smhi:/data/arkiv/field/f_archive/hirlam/G05_60lev/201109/G05_201109011200+000H00M -003|2011090115 /nobackup/smhid9/sm_esbol/hm_home/ice_36h1_4/20110901_12/ELSCFHARMALBC001 /nobackup/smhid9/sm_esbol/hm_home/ice_36h1_4/g05a/archive/2011/09/01/12/fc20110901_12+003 scp smhi:/data/arkiv/field/f_archive/hirlam/G05_60lev/201109/G05_201109011200+003H00M
In this example an scp from smhi will be executed if the expected file is not in BDDIR. There are a few environment variables that one can play with in sms/confi_exp.h that deals with the initial and boundary files
HOST_MODEL
: Tells the origin of your boundary data * ifs
: ecmwf data * hir
: hirlam data * ald
: Output from aladin physics, this also covers arpege data after fullpos processing. * ala
: Output from alaro physics * aro
: Output from arome physicsBDINT
: Interval of boundaries in hoursBDLIB
: Name of the forcing experiment. SetECMWF
to use MARS dataRCRa
to use RCRa data from ECFSBDDIR
: The path to the boundary file. In the default location BDDIR=$HM_DATA/${BDLIB}/archive/@YYYY@/@MM@/@DD@/@HH@
the file retrieved from e.g. MARS will be stored in a separate directory. On could also consider to configure this so that all the retrieved files are located in your working directory $WRK
. Locally this points to the directory where you have all your common boundary HIRLAM or ECMWF files.INT_BDFILE
: is the full path of the interpolated boundary files. The default setting is to let the boundary file be removed by directing it to $WRK
.INT_SINI_FILE
: The full path of the initial surfex file. There are a few optional environment variables that could be used that are not visible in config_exp.h
EXT_BDDIR
: External location of boundary data. If not set rules are depending on HOST_MODEL
EXT_ACCESS
: Method for accessing external data. If not set rules are depending on HOST_MODEL
BDCYCLE
: Assimilation cycle interval of forcing data, default is 6h.More about this can be bounds in the Boundary_strategy.pl
script.
The bdstrategy file is parsed by the script ExtractBD
.
scr/ExtractBD
Checks if data are on BDDIR
otherwise copy from EXT_BDDIR
. The operation performed can be different depending on HOST
and HOST_MODEL
. IFS data at ECMWF are extracted from MARS, RCR data are copied from ECFS.In case data should be retrieved from MARS there is also a stage step. When calling MARS with the stage command we ask MARS to make sure data are on disk. In HARMONIE we ask for all data for one day of r forecasts ( normally four cycles ) at the time.
The use of near real time aerosols require the presence of aerosol fields in the boundary files.
BDAERO
: Origin of the aerosol fieldsnone
: no aerosols (default configuration)cams
: aerosol from CAMS.A bdstrategycams file is generated. After the data is retrieved, the files are merge with the files from the HOSTMODEL to get the final boundary conditions files.
To be able to start the model we need the variables defining the model state.
Optional:
For the surface we need the different state variables for the different tiles. The scheme selected determines the variables.
Boundary files (coupling files) for HARMONIE are prepared in two different ways depending on the nesting procedure defined by HOST_MODEL
.
If you use data from HIRLAM or ECMWF gl_grib_api
will be called to generate boundaries. The generation can be summarized in the following steps:
HOST_MODEL
gl_grib_api
is called by the script scr/gl_bd
where we make different choices depending on PHYSICS
and HOST_MODEL
When starting a forecast there are options to whether e.g. cloud properties and TKE should be read from the initial/boundary file through NREQIN
and NCOUPLING
. At the moment these fields are read from the initial file but not coupled to. gl reads them if they are available in the input files and sets them to zero otherwise. For a Non-Hydrostatic run the non-hydrostatic pressure departure and the vertical divergence are demanded as an initial field. The pressure departure is by definition zero if you start from a non-hydrostatic mode and since the error done when disregarding the vertical divergence is small it is also set to zero in gl. There are also a choice in the forecast model to run with Q in gridpoint or in spectral space.
It's possible to use an input file without e.g. the uppermost levels. By setting LDEMAND_ALL_LEVELS=.FALSE.
the missing levels will be ignored. This is used at some institutes to reduce the amount of data transferred for the operational runs.
If you use data generated by HARMONIE you will use fullpos to generate boundaries and initial conditions. Here we will describe how it's implemented in HARMONIE but there are also good documentation on the gmapdoc site.
In HARMONIE it is done by the script scr/E927
. It contains the following steps:
Fetcht climate files. Fullpos needs a climate file and the geometry definition for both the input and output domains.
Set different moist variables in the namelists depending if your run AROME or ALADIN/ALARO.
Check if input data has Q in gridpoint or spectral space.
Demand NH variables if we run NH.
Determine the number of levels in the input file and extract the correct levels from the definition in scr/Vertical_level.pl
Run fullpos
E927 is also called from 4DVAR when the resolution is changed between the inner and outer loops.
For SURFEX we have to fill the different tiles with correct information from the input data. This is called the PREP step in the SURFEX context. scr/Prep_ini_surfex
creates an initial SURFEX file from an FA file if you run with SURFACE=surfex
.
Settings
This document was generated with Documenter.jl version 1.3.0 on Friday 17 May 2024. Using Julia version 1.9.4.