| 1 | \chapter{Advanced simulations}\label{advance} |
|---|
| 2 | |
|---|
| 3 | \vk |
|---|
| 4 | In this chapter, useful advice to perform more sophisticated simulations is provided to advanced users of the model. |
|---|
| 5 | |
|---|
| 6 | \mk |
|---|
| 7 | \section{Running nested simulations}\label{nests} |
|---|
| 8 | |
|---|
| 9 | \paragraph{Parameter files} In case you run simulations with \ttt{max\_dom} nested domains, you have to set \ttt{max\_dom} parameters wherever there is a ``," in the \ttt{namelist.input} template in chapter~\ref{zeparam}. Below we reproduce an example of the resulting syntax of the \ttt{\&time\_control}, \ttt{\&domains} and \ttt{\&bdy\_control} categories in \ttt{namelist.input}. We recommend running hydrostatic nested simulations\footnote{Non-hydrostatic nested simulations are sometimes unstable at boundaries; this should be fixed in future versions of the model.} by setting \ttt{non\_hydrostatic = F} in \ttt{\&dynamics} in \ttt{namelist.input}. If you run a simulation with, say, $3$ domains, please ensure that you defined three files \ttt{callphys.def}, \ttt{callphys\_d2.def} and \ttt{callphys\_d3.def} (one per nest). Usually all settings in these files are similar except \ttt{iradia}. |
|---|
| 10 | |
|---|
| 11 | \scriptsize |
|---|
| 12 | \codesource{OMG_namelist.input} |
|---|
| 13 | \normalsize |
|---|
| 14 | |
|---|
| 15 | \paragraph{Preprocessing steps} The additional settings for nests in \ttt{namelist.input} must correspond to the ones in the \ttt{namelist.wps} file which define the domain settings for \ttt{geogrid.exe}. A typical \ttt{namelist.wps} file for nested simulations is given below. An automated generation of \ttt{namelist.wps} from \ttt{namelist.input} is provided in the \ttt{runmeso} script as for single-domain simulations. Defining several domains yield one output per domain: e.g. for three domains \ttt{geogrid.exe} yields \ttt{geo\_em.d01.nc}, \ttt{geo\_em.d02.nc}, \ttt{geo\_em.d03.nc}\ldots, \ttt{real.exe} yields \ttt{wrfinput\_d02}, \ttt{wrfinput\_d03}, \ldots |
|---|
| 16 | |
|---|
| 17 | \scriptsize |
|---|
| 18 | \codesource{namelist.wps_NEST} |
|---|
| 19 | \normalsize |
|---|
| 20 | |
|---|
| 21 | \paragraph{Useful remarks} The model presently supports 3 nests, but more nests can be included by adaptating \ttt{runmeso} and the following files: |
|---|
| 22 | \begin{verbatim} |
|---|
| 23 | $LMDMOD/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis3.inc |
|---|
| 24 | $LMDMOD/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq3.inc |
|---|
| 25 | $LMDMOD/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/duplicate3 |
|---|
| 26 | $LMDMOD/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/generate3 |
|---|
| 27 | $LMDMOD/LMD_MM_MARS/SRC/WRFV2/mars_lmd/makegcm* ## search for 'nest' |
|---|
| 28 | \end{verbatim} |
|---|
| 29 | |
|---|
| 30 | grid points - 1 est divisé par nproc et ratio dx |
|---|
| 31 | exemple: 4 domaines et ratio 3 : il faut multiple de 12 (+1) |
|---|
| 32 | |
|---|
| 33 | attention seulement 4 processeurs |
|---|
| 34 | |
|---|
| 35 | nests >> commencer par faire un test nest par nest en en rajoutant un n-1 à chaque succès au niveau n |
|---|
| 36 | -- une instabilité dans le nest 2 peut provenir du nest 1. stabiliser le nest 1 tout d'abord. |
|---|
| 37 | |
|---|
| 38 | resolution à virgule au moins précis à deux chiffres entre WPS et WRF |
|---|
| 39 | |
|---|
| 40 | \mk |
|---|
| 41 | \section{Running simulations with tracers} |
|---|
| 42 | |
|---|
| 43 | mars: number of corresponding tracers |
|---|
| 44 | |
|---|
| 45 | |
|---|
| 46 | |
|---|
| 47 | \mk |
|---|
| 48 | \section{Running simulations with the new physics} |
|---|
| 49 | different callphys.def |
|---|
| 50 | the step datafile.h is not needed anymore ! use callphys.def. |
|---|
| 51 | traceur.def |
|---|
| 52 | run.def |
|---|
| 53 | different callphys.def |
|---|
| 54 | makemeso -p |
|---|
| 55 | |
|---|
| 56 | \mk |
|---|
| 57 | \section{Running Large-Eddy Simulations} |
|---|
| 58 | |
|---|
| 59 | \mk |
|---|
| 60 | \section{Controlling which fields to output} |
|---|
| 61 | %\section{geogrid.tbl} |
|---|
| 62 | |
|---|
| 63 | \scriptsize |
|---|
| 64 | \codesource{Registry.EM.extract} |
|---|
| 65 | \normalsize |
|---|
| 66 | |
|---|
| 67 | ne pas déclarer un 2D en 3D et cioce versa |
|---|
| 68 | FAIRE registry.bash UNEFOISFINI |
|---|
| 69 | |
|---|
| 70 | \mk |
|---|
| 71 | \section{Interpolating outputs on altitude and pressure levels}\label{postproc} |
|---|
| 72 | \ttt{api} |
|---|
| 73 | |
|---|
| 74 | \mk |
|---|
| 75 | \section{Generating maps for winds and meteorological fields}\label{plots} |
|---|
| 76 | |
|---|
| 77 | \ttt{python} |
|---|
| 78 | \ttt{idl} ? |
|---|
| 79 | |
|---|
