Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found
Select Git revision

Target

Select target project
  • mesonh/mesonh-code
  • quentin.rodier/mesonh-code-fork
  • benoit.vie/mesonh-code
  • joris.pianezze/mesonh-code
  • 8qq4g5s7/mesonh-code
  • jean.baptiste.filippi/meso-nh-fire-code
  • fdl68d9p/mesonh-code-sophia
  • jean.wurtz/mesonh-code
8 results
Select Git revision
Show changes
Commits on Source (586)
486 additional commits have been omitted to prevent performance issues.
Hide whitespace changes
Inline Side-by-side
Showing
with 456 additions and 170 deletions
.gitignore
View file @ af717e63
conf/profile_mesonh conf/profile_mesonh
conf/profile_mesonh-* conf/profile_mesonh-*
exe/* exe/*
MY_RUN/INTEGRATION_CASES/**/*.des
MY_RUN/INTEGRATION_CASES/**/*.dir
MY_RUN/INTEGRATION_CASES/**/*.hdr
MY_RUN/INTEGRATION_CASES/**/*.nc
MY_RUN/INTEGRATION_CASES/**/ecmwf.OD.????????.??
MY_RUN/INTEGRATION_CASES/**/file_for_xtransfer
MY_RUN/INTEGRATION_CASES/**/OUTPUT_LISTING*
MY_RUN/INTEGRATION_CASES/**/PRESSURE
MY_RUN/KTEST/*/*/*.des MY_RUN/KTEST/*/*/*.des
MY_RUN/KTEST/*/*/*.lfi MY_RUN/KTEST/*/*/*.lfi
MY_RUN/KTEST/*/*/*.nc MY_RUN/KTEST/*/*/*.nc
MY_RUN/KTEST/*/*/dir.??????:??:??
MY_RUN/KTEST/*/*/dir_save
MY_RUN/KTEST/*/*/dirconv
MY_RUN/KTEST/*/*/dirextract
MY_RUN/KTEST/*/*/ecmwf.OD.????????.?? MY_RUN/KTEST/*/*/ecmwf.OD.????????.??
MY_RUN/KTEST/*/*/file_for_xtransfer MY_RUN/KTEST/*/*/file_for_xtransfer
MY_RUN/KTEST/*/*/gmeta
MY_RUN/KTEST/*/*/gtopo30.* MY_RUN/KTEST/*/*/gtopo30.*
MY_RUN/KTEST/*/*/CLAY_HWSD_MOY.* MY_RUN/KTEST/*/*/CLAY_HWSD_MOY.*
MY_RUN/KTEST/*/*/DATA?1 MY_RUN/KTEST/*/*/DATA?1
...@@ -22,15 +25,14 @@ MY_RUN/KTEST/*/*/OUTPUT_LISTING* ...@@ -22,15 +25,14 @@ MY_RUN/KTEST/*/*/OUTPUT_LISTING*
MY_RUN/KTEST/*/*/PRESSURE MY_RUN/KTEST/*/*/PRESSURE
MY_RUN/KTEST/*/*/REMAP*FFT* MY_RUN/KTEST/*/*/REMAP*FFT*
MY_RUN/KTEST/*/*/SAND_HWSD_MOY.* MY_RUN/KTEST/*/*/SAND_HWSD_MOY.*
MY_RUN/KTEST/*/*/output_save*
MY_RUN/KTEST/*/*/pipe_name MY_RUN/KTEST/*/*/pipe_name
MY_RUN/KTEST/*/*/gshhs(?).rim MY_RUN/KTEST/*/*/__pycache__
MY_RUN/KTEST/*/*/gshhs(?).zip MY_RUN/KTEST/*/*/misc_functions.py
MY_RUN/KTEST/*/*/rangs(?).cat MY_RUN/KTEST/*/*/Panel_Plot.py
MY_RUN/KTEST/*/*/rangs(?).cel MY_RUN/KTEST/*/*/read_MNHfile.py
MY_RUN/KTEST/*/*/rangs(?).zip MY_RUN/KTEST/*/*/tutorial_readme
MY_RUN/KTEST/*/*/visu*.png MY_RUN/KTEST/*/*/*.png
MY_RUN/KTEST/*/*/zsection*.png MY_RUN/KTEST/*/*/*.pdf
MY_RUN/KTEST/007_16janvier/012_spectre/spectra_16JAN.1.12B18.001_* MY_RUN/KTEST/007_16janvier/012_spectre/spectra_16JAN.1.12B18.001_*
MY_RUN/KTEST/009_ICARTT/001_pgd1/*.asc MY_RUN/KTEST/009_ICARTT/001_pgd1/*.asc
MY_RUN/KTEST/009_ICARTT/002_arp2lfi/ecmwf.OD.20040810.18-V2 MY_RUN/KTEST/009_ICARTT/002_arp2lfi/ecmwf.OD.20040810.18-V2
...@@ -57,3 +59,9 @@ src/LIB/oasis3-mct* ...@@ -57,3 +59,9 @@ src/LIB/oasis3-mct*
!src/LIB/oasis3-mct*.tar.gz !src/LIB/oasis3-mct*.tar.gz
src/LIB/toy* src/LIB/toy*
!src/LIB/toy*.tar.gz !src/LIB/toy*.tar.gz
vscode_mod
.DS_Store
*.vts
*.vtu
*.mp4
A-INSTALL
View file @ af717e63
# #
# Version of PACKAGE MESONH "Open distribution" # Version of PACKAGE MESONH "Open distribution"
# PACK-MNH-V5-6-2 # PACK-MNH-V5-7-1
# DATE : 28/11/2023 # DATE : 04/09/2024
# VERSION : MESONH MASDEV5_6 + BUG-2 # VERSION : MESONH MASDEV5_7 + BUG-1
# #
# MAP # MAP
# #
...@@ -28,7 +28,7 @@ ...@@ -28,7 +28,7 @@
# e) cleaning previous compiled version # e) cleaning previous compiled version
# #
# #
# NEW on MNH-56X : For conda package for graphic output in test-case examples , read # NEW since MNH-56X : For conda package for graphic output in test-case examples , read
# the 'README_MNH_CONDA' file in this (root) directory # the 'README_MNH_CONDA' file in this (root) directory
# #
# #
...@@ -84,14 +84,14 @@ ...@@ -84,14 +84,14 @@
# #
# or directly # or directly
# #
# http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_MESONH/MNH-V5-6-2.tar.gz # http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_MESONH/MNH-V5-7-1.tar.gz
# #
# Then untar the file "MNH-V5-6-2.tar.gz" where you want to. # Then untar the file "MNH-V5-7-1.tar.gz" where you want to.
# For example, in your home directory: # For example, in your home directory:
# #
cd ~ cd ~
tar xvfz MNH-V5-6-2.tar.gz tar xvfz MNH-V5-7-1.tar.gz
# #
# Process now to the chapter to configure the MesoNH package. # Process now to the chapter to configure the MesoNH package.
...@@ -134,7 +134,7 @@ git lfs install ...@@ -134,7 +134,7 @@ git lfs install
# only access) by following the next link: # only access) by following the next link:
# #
http://mesonh.aero.obs-mip.fr/mesonh56/GitSources?action=AttachFile&do=get&target=anongitmesonh.key http://mesonh.aero.obs-mip.fr/mesonh57/GitSources?action=AttachFile&do=get&target=anongitmesonh.key
# #
# and save the file in your $HOME/.ssh/ directory. # and save the file in your $HOME/.ssh/ directory.
...@@ -165,68 +165,68 @@ git config --global http.sslverify false ...@@ -165,68 +165,68 @@ git config --global http.sslverify false
# This is necessary to disable the certificate checks because a self-signed # This is necessary to disable the certificate checks because a self-signed
# certificate was used for the LFS server. # certificate was used for the LFS server.
# #
# c) Cloning the Meso-NH Source repository on the developpement branch MNH-56-branch # c) Cloning the Meso-NH Source repository on the developpement branch MNH-57-branch
# ---------------------------------------------------------------------------------- # ----------------------------------------------------------------------------------
# #
# Finally you can clone the Meso-NH Git repository with the following command: # Finally you can clone the Meso-NH Git repository with the following command:
# #
git lfs clone anongit@anongit_mesonh:/gitrepos/MNH-git_open_source-lfs.git -b MNH-56-branch MNH-V5-6-2 git clone anongit@anongit_mesonh:/gitrepos/MNH-git_open_source-lfs.git -b MNH-57-branch MNH-V5-7-1
# #
# that will create the MNH-V5-6-2 directory containing a clone (copy) of the # that will create the MNH-V5-7-1 directory containing a clone (copy) of the
# Meso-NH package on the remote developpement branch MNH-56-branch # Meso-NH package on the remote developpement branch MNH-57-branch
# #
# #
# d) Checking out a given version of MESONH # d) Checking out a given version of MESONH
# ----------------------------------------- # -----------------------------------------
# #
# Once the repository is cloned, it's better for you to checkout your own branch # Once the repository is cloned, it's better for you to checkout your own branch
# (by default, you are on HEAD of the MNH-56-branch development branch ). # (by default, you are on HEAD of the MNH-57-branch development branch ).
# #
# To create your local branch corresponding to the V5-6-2 version, type: # To create your local branch corresponding to the V5-7-1 version, type:
# #
cd MNH-V5-6-2 cd MNH-V5-7-1
git checkout -b MYB-MNH-V5-6-2 PACK-MNH-V5-6-2 git checkout -b MYB-MNH-V5-7-1 PACK-MNH-V5-7-1
# #
# MYB-MNH-V5-6-2 is the name of the local branch you created # MYB-MNH-V5-7-1 is the name of the local branch you created
# and # and
# PACK-MNH-V5-6-2 is the remote/origin tag on which it is based. # PACK-MNH-V5-7-1 is the remote/origin tag on which it is based.
# #
# The advantage of this way of downloading the package is that in the future # The advantage of this way of downloading the package is that in the future
# you could check/update quickly differences with the new version of the # you could check/update quickly differences with the new version of the
# package without having to download entirely the full package. # package without having to download entirely the full package.
# #
# Suppose that a new version, for example "PACK-MNH-V5-6-3", is announced. # Suppose that a new version, for example "PACK-MNH-V5-7-1", is announced.
# #
# To see the differences with your working copy, do: # To see the differences with your working copy, do:
# #
git fetch git fetch
git diff HEAD PACK-MNH-V5-6-3 git diff HEAD PACK-MNH-V5-7-1
# #
# To go to the new version, you can, for example, create a new local branch: # To go to the new version, you can, for example, create a new local branch:
# #
git checkout -b MYB-MNH-V5-6-3 PACK-MNH-V5-6-3 git checkout -b MYB-MNH-V5-7-1 PACK-MNH-V5-7-1
# #
# At any time, you can also check for "uptodate" changes in the Git branch # At any time, you can also check for "uptodate" changes in the Git branch
# dedicated to the MNH56 version before the official release of the "bugN+1" # dedicated to the MNH57 version before the official release of the "bugN+1"
# bugfix version. # bugfix version.
# #
git fetch git fetch
git diff HEAD MNH-56-branch git diff HEAD MNH-57-branch
# #
# And, test this development (not yet official) version by going to this branch: # And, test this development (not yet official) version by going to this branch:
# #
git checkout --track origin/MNH-56-branch git checkout --track origin/MNH-57-branch
# #
# e) Cloning the Meso-NH Documentation repository # e) Cloning the Meso-NH Documentation repository
# #
...@@ -257,7 +257,7 @@ git clone anongit@anongit_mesonh:/gitrepos/MNH-DOC.git ...@@ -257,7 +257,7 @@ git clone anongit@anongit_mesonh:/gitrepos/MNH-DOC.git
# use the "./configure" script like this # use the "./configure" script like this
# #
cd ~/MNH-V5-6-2/src cd ~/MNH-V5-7-1/src
./configure ./configure
. ../conf/profile_mesonh . ../conf/profile_mesonh
...@@ -308,7 +308,7 @@ export OPTLEVEL=O2 # Compile in O2, 4 times faster then DEBUG, but less ...@@ -308,7 +308,7 @@ export OPTLEVEL=O2 # Compile in O2, 4 times faster then DEBUG, but less
# and then source/load the new generate file # and then source/load the new generate file
. ../conf/profile_mesonh.LXifort.MNH-V5-6-2.MPIAUTO.O2 . ../conf/profile_mesonh.LXifort.MNH-V5-7-1.MPIAUTO.O2
# #
# REM: # REM:
...@@ -333,7 +333,7 @@ export OPTLEVEL=O2 # Compile in O2, 4 times faster then DEBUG, but less ...@@ -333,7 +333,7 @@ export OPTLEVEL=O2 # Compile in O2, 4 times faster then DEBUG, but less
# go to the directory "src" # go to the directory "src"
# #
cd ~/MNH-V5-6-2/src cd ~/MNH-V5-7-1/src
# #
# if you have not already configured your MESONH environment # if you have not already configured your MESONH environment
...@@ -562,7 +562,7 @@ make examples ...@@ -562,7 +562,7 @@ make examples
# #
cd $WORKDIR cd $WORKDIR
cd MNH-V5-6-2/src cd MNH-V5-7-1/src
./configure ./configure
...@@ -576,8 +576,8 @@ cd MNH-V5-6-2/src ...@@ -576,8 +576,8 @@ cd MNH-V5-6-2/src
# - On JEAN-ZAY ( HPE ) the compilation is in interactive : # - On JEAN-ZAY ( HPE ) the compilation is in interactive :
cd MNH-V5-6-2/src cd MNH-V5-7-1/src
. ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIINTEL-O2 . ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIINTEL-O2
make -j16 |& tee error$XYZ make -j16 |& tee error$XYZ
make installmaster make installmaster
...@@ -612,12 +612,12 @@ sbatch job_make_examples_BullX_occigen ...@@ -612,12 +612,12 @@ sbatch job_make_examples_BullX_occigen
# - ssh irene-fr : for Intel SkyLake processors # - ssh irene-fr : for Intel SkyLake processors
# On Intel processors the MPI use is OPENMPI/4.1.4 # On Intel processors the MPI use is OPENMPI/4.1.4
# the configure will generate a # the configure will generate a
# profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2 # profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIAUTO-O2
# #
# - ssh irene-amd : for AMD , processors # - ssh irene-amd : for AMD , processors
# On AMD processors the MPI use is OPENMPI/4.1.4 # On AMD processors the MPI use is OPENMPI/4.1.4
# the configure will generate a # the configure will generate a
# profile_mesonh-LXifort-R8I4-MNH-V5-6-2-AMD-MPIAUTO-O2 # profile_mesonh-LXifort-R8I4-MNH-V5-7-1-AMD-MPIAUTO-O2
# #
# #
# At TGCC they use 'one login' for multi-project allocation . # At TGCC they use 'one login' for multi-project allocation .
...@@ -642,14 +642,14 @@ id -ng ...@@ -642,14 +642,14 @@ id -ng
echo $CCCHOME echo $CCCHOME
--> /ccc/.../home/genXXXX/{your_login} --> /ccc/.../home/genXXXX/{your_login}
# If all is OK , you could install the code in interactive as usally # If all is OK, you can install the code interactively as usual
# #
# - install the PACKAGE in your ${CCCHOME} ( default 20Go of quota ) # - install the PACKAGE in your ${CCCHOME} ( default 20Go of quota )
# for example on Intel frontal # for example on Intel frontal
cd ${CCCHOME} cd ${CCCHOME}
cd MNH-V5-6-2/src cd MNH-V5-7-1/src
./configure ./configure
. ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2 . ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIAUTO-O2
make -j16 |& tee error$XYZ make -j16 |& tee error$XYZ
make installmaster make installmaster
...@@ -689,14 +689,14 @@ sbatch job_make_examples_Atos_HPCF ...@@ -689,14 +689,14 @@ sbatch job_make_examples_Atos_HPCF
# - At Meteo-France DSI on belenos # - At Meteo-France DSI on belenos
# #
# to install the whole package on your "$HOME" directory # to install the whole package on your "$HOME" directory
# untar the file "MNH-V5-6-2.tar.gz" from its location : # untar the file "MNH-V5-7-1.tar.gz" from its location :
cd ~ cd ~
tar xvf $MESONH/MNH-V5-6-2.tar.gz tar xvf $MESONH/MNH-V5-7-1.tar.gz
# run the "./configure" command : # run the "./configure" command :
cd MNH-V5-6-2/src cd MNH-V5-7-1/src
./configure ./configure
# #
...@@ -795,7 +795,7 @@ scandollar ...@@ -795,7 +795,7 @@ scandollar
## OUTPUT :: ## OUTPUT ::
># read default config file :: ---> CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-6-2/conf/post/confdollar_aeropc_default ># read default config file :: ---> CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-7-1/conf/post/confdollar_aeropc_default
># >#
># read user config file :: ---> CONFIG=confdollar ># read user config file :: ---> CONFIG=confdollar
># >#
...@@ -817,7 +817,7 @@ scandollar 0* ...@@ -817,7 +817,7 @@ scandollar 0*
## OUTPUT :: ## OUTPUT ::
># >#
># read default config file :: ---> CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-6-2/conf/post/confdollar_aeropc_default ># read default config file :: ---> CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-7-1/conf/post/confdollar_aeropc_default
># >#
># read user config file :: ---> CONFIG=confdollar ># read user config file :: ---> CONFIG=confdollar
># >#
...@@ -891,22 +891,22 @@ cp -R 007_16janvier_scandollar /.../your_directory ...@@ -891,22 +891,22 @@ cp -R 007_16janvier_scandollar /.../your_directory
# #
# use this "profile_mesonh" : # use this "profile_mesonh" :
. /home/rech/mnh/rmnh007/DEV/MNH-V5-6-2/conf/profile_mesonh-SX8-MNH-V5-6-2-MPIAUTO-O4 . /home/rech/mnh/rmnh007/DEV/MNH-V5-7-1/conf/profile_mesonh-SX8-MNH-V5-7-1-MPIAUTO-O4
# And the examples are here ( link to my $WORKDIR in actually ) # And the examples are here ( link to my $WORKDIR in actually )
/home/rech/mnh/rmnh007/DEV/MNH-V5-6-2/MY_RUN/KTEST/007_16janvier_scandollar /home/rech/mnh/rmnh007/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
# #
# On vargas # On vargas
# --------- # ---------
# use this "profile_mesonh" : # use this "profile_mesonh" :
. /workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-6-2/conf/profile_mesonh-AIX64-MNH-V5-6-2-MPIAUTO-O2 . /workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-7-1/conf/profile_mesonh-AIX64-MNH-V5-7-1-MPIAUTO-O2
# and examples here : # and examples here :
/workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-6-2/MY_RUN/KTEST/007_16janvier_scandollar /workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
# #
# - At CINES on JADE : # - At CINES on JADE :
...@@ -914,11 +914,11 @@ cp -R 007_16janvier_scandollar /.../your_directory ...@@ -914,11 +914,11 @@ cp -R 007_16janvier_scandollar /.../your_directory
# #
# use # use
. /work/escobar/DEV/MNH-V5-6-2/conf/profile_mesonh-LXifort-MNH-V5-6-2-MPIICE-O2 . /work/escobar/DEV/MNH-V5-7-1/conf/profile_mesonh-LXifort-MNH-V5-7-1-MPIICE-O2
# and the exemples # and the exemples
/work/escobar/DEV/MNH-V5-6-2/MY_RUN/KTEST/007_16janvier_scandollar /work/escobar/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
# #
# - At ECMWF on cxa : # - At ECMWF on cxa :
...@@ -926,11 +926,11 @@ cp -R 007_16janvier_scandollar /.../your_directory ...@@ -926,11 +926,11 @@ cp -R 007_16janvier_scandollar /.../your_directory
# #
# use # use
. /c1a/ms_perm/au5/MNH-V5-6-2/conf/profile_mesonh-AIX64-MNH-V5-6-2-MPIAUTO-O2 . /c1a/ms_perm/au5/MNH-V5-7-1/conf/profile_mesonh-AIX64-MNH-V5-7-1-MPIAUTO-O2
# and the examples # and the examples
/c1a/ms_perm/au5/MNH-V5-6-2/MY_RUN/KTEST/007_16janvier_scandollar /c1a/ms_perm/au5/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
# #
...@@ -989,56 +989,42 @@ ulimit -s unlimited ...@@ -989,56 +989,42 @@ ulimit -s unlimited
# ------------------------------------------------------------- # -------------------------------------------------------------
# #
# If you want to use coupled(inline) run with FOREFIRE and MESONH # If you want to use coupled(inline) run with FOREFIRE and MESONH
# or use routines to do high frequency parallel file dumps
# you could compile the interfaced/coupling routine by activating # you could compile the interfaced/coupling routine by activating
# this variable before any compilation # this variable before any compilation.
# FOREFIRE available at https://github.com/forefireAPI/forefire
# follow instruction to compile with option -DMPI_COUPLING to run
# you have to set #FOREFIREHOME to the directory where you have installed
# the forefire package so it can find extras, such as web interface
# After compilation you must link :
# $MESONH/EXE/libForeFire.so with $FOREFIREHOME/lib/libforefireL.{so,dylib,dll}
# for the MESONH binary to find the FOREFIRE library.
export MNH_FOREFIRE=1.0 export MNH_FOREFIRE=1.0
# and then the configure
./configure
#
# after this compile as usually
#
make
etc ...
#
# The FOREFIRE API package himself must be compiled independently from MesoNH
# The git repository is here https://github.com/forefireAPI/firefront/tree/2014.01
# it could be cloned by
git clone -b 2014.01 https://github.com/forefireAPI/firefront.git
# it depend on netcdf and scons for is compilation
# the "libForeFIre.so" generate must by referenced in the LD_LIBRARY_PATH
# or move/linked to the 'exe' directory of MesoNH
#
# see d) for cleaning previously version if needed # see d) for cleaning previously version if needed
# b) MNH_RTTOV for optional radiative computation # b) MNH_RTTOV for optional radiative computation
# -------------------------------------- # --------------------------------------
# #
# The RTTOV package was not included into the open source version of Meso-NH # The RTTOV 13.2 package was not included into the open source version of Meso-NH
# because it needs a licence agrement. # because it needs a licence agrement.
# #
# ----------------------------------
# OPTION 1: Use version 13.2 of RTTOV
# -----------------------------------
#
# Run the 'configure' script preceded with the setting of the MNH_RTTOV variable: # Run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
# #
cd MNH.../src/ cd $SRC_MESONH/src/
export MNH_RTTOV=1 export MNH_RTTOV=1
export VER_RTTOV=13.2 export VER_RTTOV=13.2
# #
# Compile the HDF5 library
#
make cdf
#
# Download the RTTOV package rttov132.tar.xz by following the instructions given on https://nwpsaf.eu/site/software/rttov/ # Download the RTTOV package rttov132.tar.xz by following the instructions given on https://nwpsaf.eu/site/software/rttov/
# #
# Install the RTTOV package rttov132.tar.xz # Install the RTTOV package rttov132.tar.xz
cd MNH.../src/LIB cd $SRC_MESONH/src/LIB
mkdir RTTOV-13.2 mkdir RTTOV-13.2
cd RTTOV-13.2 cd RTTOV-13.2
tar xJf rttov132.tar.xz tar xJf rttov132.tar.xz
...@@ -1049,72 +1035,15 @@ HDF5_PREFIX = $(SRC_MESONH)/src/dir_obj${XYZ}/MASTER/NETCDF-${VERSION_CDFF} ...@@ -1049,72 +1035,15 @@ HDF5_PREFIX = $(SRC_MESONH)/src/dir_obj${XYZ}/MASTER/NETCDF-${VERSION_CDFF}
FFLAGS_HDF5 = -D_RTTOV_HDF $(FFLAG_MOD)$(HDF5_PREFIX)/include FFLAGS_HDF5 = -D_RTTOV_HDF $(FFLAG_MOD)$(HDF5_PREFIX)/include
LDFLAGS_HDF5 = -L$(HDF5_PREFIX)/lib64 -lhdf5hl_fortran -lhdf5_hl -lhdf5_fortran -lhdf5 -lsz -laec -lz -ldl LDFLAGS_HDF5 = -L$(HDF5_PREFIX)/lib64 -lhdf5hl_fortran -lhdf5_hl -lhdf5_fortran -lhdf5 -lsz -laec -lz -ldl
" "
cd src cd ../src
../build/Makefile.PL RTTOV_HDF=1 ../build/Makefile.PL RTTOV_HDF=1
make ARCH=ifort # Use Intel "ifort" compiler; other options: gfortran, NAG, pgf90, IBM make ARCH=ifort # Use Intel "ifort" compiler; other options: gfortran, NAG, pgf90, IBM
# #
# And then for the Meso-NH compilation, do # And then for the Meso-NH compilation, do
# #
cd MNH.../src/ cd $SRC_MESONH/src/
make make
# ----------------------------------
# OPTION 2: Use version 11.3 of RTTOV
# -----------------------------------
# Download the RTTOV package rttov113.tar.gz by following the instructions given on https://nwpsaf.eu/site/software/rttov/
#
# Install the RTTOV package rttov113.tar.gz
cd MNH.../src/LIB
mkdir RTTOV-11.3
cd RTTOV-11.3
tar xvfz rttov113.tar.gz
cd src
make ARCH=ifort # Use Intel "ifort" compiler; other options: gfortran, NAG, pgf90, IBM
#
# And then for the compilation, run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
#
cd MNH.../src/
export MNH_RTTOV=1
export VER_RTTOV=11.3
./configure
etc ...
# ----------------------------------
# OPTION 3: Use version 8.7 of RTTOV
# ----------------------------------
# For already(old) licencied MesoNH users (MNH-4-X version with research licence see here: http://mesonh.aero.obs-mip.fr/mesonh410/UserInformation)
#
# the package could be reloaded in this way
#
# - With cvs access
#
cd MNH.../src/LIB
cvs up -rPACK-MNH-V4-10-3 -d -P RTTOV
#
# - With WEB access (with WEB login/pass as usually) the RTTOV package could also be retrieve in tarball with wget like this:
#
cd MNH.../src/LIB
wget --http-user=USER --http-password=PASS 'http://mesonh.aero.obs-mip.fr/cgi-bin/mesonh_interne/viewcvs.cgi/MNH-VX-Y-Z/src/LIB/RTTOV/?view=tar&pathrev=PACK-MNH-V4-10-3' -O RTTOV.tar.gz
tar xvfz RTTOV.tar.gz
# /!\ don't forget the ''. Otherwise, this will not work!
#
# And then for the compilation, run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
#
cd MNH.../src/
export MNH_RTTOV=1
export VER_RTTOV=8.7
./configure
etc ...
# c) MNH_ECRAD for optional compilation of new ECRAD radiative library from ECMWF # c) MNH_ECRAD for optional compilation of new ECRAD radiative library from ECMWF
# -------------------------------------- # --------------------------------------
# #
......
LIBTOOLS/tools/lfi2cdf/src/lfi2cdf.f90
View file @ af717e63
!MNH_LIC Copyright 1994-2023 CNRS, Meteo-France and Universite Paul Sabatier !MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1. !MNH_LIC for details. version 1.
...@@ -59,8 +59,8 @@ program LFI2CDF ...@@ -59,8 +59,8 @@ program LFI2CDF
ALLOCATE(TIMEZ) !Used by IO_WRITE_FIELD ALLOCATE(TIMEZ) !Used by IO_WRITE_FIELD
NIO_VERB = NVERB_WARNING NIO_VERB = NVERB_INFO
NGEN_VERB = NVERB_WARNING NGEN_VERB = NVERB_INFO
NIO_ABORT_LEVEL = NVERB_FATAL NIO_ABORT_LEVEL = NVERB_FATAL
NGEN_ABORT_LEVEL = NVERB_FATAL NGEN_ABORT_LEVEL = NVERB_FATAL
LVERB_OUTLST = .FALSE. LVERB_OUTLST = .FALSE.
......
LIBTOOLS/tools/lfi2cdf/src/mode_util.f90
View file @ af717e63
!MNH_LIC Copyright 1994-2023 CNRS, Meteo-France and Universite Paul Sabatier !MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1. !MNH_LIC for details. version 1.
...@@ -69,6 +69,7 @@ CONTAINS ...@@ -69,6 +69,7 @@ CONTAINS
SUBROUTINE parse_infiles(infiles, outfiles, KNFILES_OUT, nbvar_infile, nbvar_tbr, nbvar_calc, nbvar_tbw, & SUBROUTINE parse_infiles(infiles, outfiles, KNFILES_OUT, nbvar_infile, nbvar_tbr, nbvar_calc, nbvar_tbw, &
tpreclist, options, runmode) tpreclist, options, runmode)
USE MODD_DIM_n, ONLY: NIMAX_ll, NJMAX_ll, NKMAX USE MODD_DIM_n, ONLY: NIMAX_ll, NJMAX_ll, NKMAX
use modd_io, only: nio_verb
USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT, NGRIDUNKNOWN USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT, NGRIDUNKNOWN
use mode_io_tools_nc4, only: IO_Dimids_guess_nc4 use mode_io_tools_nc4, only: IO_Dimids_guess_nc4
...@@ -126,10 +127,11 @@ CONTAINS ...@@ -126,10 +127,11 @@ CONTAINS
kcdf_id = INFILES(1)%TFILE%NNCID kcdf_id = INFILES(1)%TFILE%NNCID
END IF END IF
PRINT *,'MESONH 3D, 2D articles DIMENSIONS used :' WRITE( cmnhmsg(1), '( A )' ) 'MESONH 3D, 2D articles DIMENSIONS used :'
PRINT *,'DIMX =',NIMAX_ll+2*JPHEXT WRITE( cmnhmsg(2), '( "DIMX = ", I0 )' ) NIMAX_ll + 2 * JPHEXT
PRINT *,'DIMY =',NJMAX_ll+2*JPHEXT WRITE( cmnhmsg(3), '( "DIMY = ", I0 )' ) NJMAX_ll + 2 * JPHEXT
PRINT *,'DIMZ =',NKMAX +2*JPVEXT WRITE( cmnhmsg(4), '( "DIMZ = ", I0 )' ) NKMAX + 2 * JPVEXT
call Print_msg( NVERB_INFO, 'IO', 'parse_infiles' )
! Phase 1 : build articles list to convert. ! Phase 1 : build articles list to convert.
! !
...@@ -270,10 +272,14 @@ CONTAINS ...@@ -270,10 +272,14 @@ CONTAINS
END DO END DO
maxvar = nbvar_tbr+nbvar_calc maxvar = nbvar_tbr+nbvar_calc
DO ji=1,nbvar_tbr+nbvar_calc
print *,ji,'name=',tpreclist(ji)%name,' calc=',tpreclist(ji)%calc,' tbw=',tpreclist(ji)%tbw,& if ( nio_verb >= NVERB_INFO ) then
' tbr=',tpreclist(ji)%tbr,' found=',tpreclist(ji)%found do ji = 1, nbvar_tbr + nbvar_calc
END DO write ( cmnhmsg(1), '( I0, " name=", A, "calc=", L1, " tbw=", L1, " tbr=", L1, " found=", L1 )' ) &
ji, tpreclist(ji)%name, tpreclist(ji)%calc, tpreclist(ji)%tbw, tpreclist(ji)%tbr, tpreclist(ji)%found
call Print_msg( NVERB_INFO, 'IO', 'parse_infiles' )
end do
end if
ELSE ELSE
! Entire file is converted ! Entire file is converted
...@@ -977,7 +983,8 @@ END DO ...@@ -977,7 +983,8 @@ END DO
USE MODD_CONFZ, ONLY: NB_PROCIO_R USE MODD_CONFZ, ONLY: NB_PROCIO_R
USE MODD_DIM_n, ONLY: NIMAX_ll, NJMAX_ll, NKMAX USE MODD_DIM_n, ONLY: NIMAX_ll, NJMAX_ll, NKMAX
USE MODD_GRID, ONLY: XBETA, XRPK, XLAT0, XLON0, XLATORI, XLONORI USE MODD_GRID, ONLY: XBETA, XRPK, XLAT0, XLON0, XLATORI, XLONORI
USE MODD_GRID_n, ONLY: LSLEVE, XXHAT, XYHAT, XZHAT USE MODD_GRID_n, ONLY: LSLEVE, XXHAT, XXHATM, XYHAT, XYHATM, XZHAT, XZHATM, &
XXHAT_ll, XYHAT_ll, XXHATM_ll, XYHATM_ll
USE MODD_IO, ONLY: LIOCDF4 USE MODD_IO, ONLY: LIOCDF4
USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT
USE MODD_PARAMETERS_ll, ONLY: JPHEXT_ll=>JPHEXT, JPVEXT_ll=>JPVEXT USE MODD_PARAMETERS_ll, ONLY: JPHEXT_ll=>JPHEXT, JPVEXT_ll=>JPVEXT
...@@ -985,6 +992,7 @@ END DO ...@@ -985,6 +992,7 @@ END DO
USE MODE_IO_FILE, ONLY: IO_FILE_OPEN, IO_FILE_CLOSE USE MODE_IO_FILE, ONLY: IO_FILE_OPEN, IO_FILE_CLOSE
USE MODE_IO_MANAGE_STRUCT, ONLY: IO_FILE_ADD2LIST USE MODE_IO_MANAGE_STRUCT, ONLY: IO_FILE_ADD2LIST
USE MODE_SET_GRID, ONLY: INTERP_VERGRID_TO_MASSPOINTS
TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: infiles TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: infiles
TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: outfiles TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: outfiles
...@@ -997,6 +1005,7 @@ END DO ...@@ -997,6 +1005,7 @@ END DO
character(len=:), allocatable :: yunits character(len=:), allocatable :: yunits
INTEGER :: idx, IRESP2 INTEGER :: idx, IRESP2
integer :: iiu, iju, iku
integer :: inb_procio_r_save integer :: inb_procio_r_save
INTEGER(KIND=CDFINT) :: ioldmode INTEGER(KIND=CDFINT) :: ioldmode
INTEGER(KIND=CDFINT) :: istatus INTEGER(KIND=CDFINT) :: istatus
...@@ -1019,7 +1028,11 @@ END DO ...@@ -1019,7 +1028,11 @@ END DO
CALL IO_FILE_ADD2LIST(INFILES(1)%TFILE,HINFILE,'MNH','READ',HFORMAT='NETCDF4') CALL IO_FILE_ADD2LIST(INFILES(1)%TFILE,HINFILE,'MNH','READ',HFORMAT='NETCDF4')
CALL IO_FILE_OPEN(INFILES(1)%TFILE) CALL IO_FILE_OPEN(INFILES(1)%TFILE)
nbvar_infile = INFILES(1)%TFILE%NNCNAR istatus = NF90_INQUIRE( infiles(1)%tfile%nncid, nvariables=nbvar_infile )
if (istatus /= NF90_NOERR) then
call Print_msg( NVERB_FATAL, 'IO', 'OPEN_FILES', 'NF90_INQUIRE for ' &
// trim(infiles(1)%tfile%cname) // '.nc: ' // NF90_STRERROR(istatus) )
end if
!Open fallback file if provided !Open fallback file if provided
if ( options( OPTFALLBACK )%set ) then if ( options( OPTFALLBACK )%set ) then
...@@ -1095,18 +1108,33 @@ END DO ...@@ -1095,18 +1108,33 @@ END DO
if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'STORAGE_TYPE',CSTORAGE_TYPE,IRESP2) if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'STORAGE_TYPE',CSTORAGE_TYPE,IRESP2)
if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'STORAGE_TYPE not found') if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'STORAGE_TYPE not found')
! !
ALLOCATE(XXHAT(NIMAX_ll+2*JPHEXT)) iiu = nimax_ll + 2 * JPHEXT
ALLOCATE( XXHAT(iiu), XXHATM(iiu) )
CALL IO_Field_read(INFILES(1)%TFILE,'XHAT',XXHAT,IRESP2) CALL IO_Field_read(INFILES(1)%TFILE,'XHAT',XXHAT,IRESP2)
!If not found in main file, try the fallback one !If not found in main file, try the fallback one
if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'XHAT',XXHAT,IRESP2) if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'XHAT',XXHAT,IRESP2)
if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'XHAT not found') if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'XHAT not found')
ALLOCATE(XYHAT(NJMAX_ll+2*JPHEXT)) iju = njmax_ll + 2 * JPHEXT
ALLOCATE( XYHAT(iju), XYHATM(iju) )
CALL IO_Field_read(INFILES(1)%TFILE,'YHAT',XYHAT,IRESP2) CALL IO_Field_read(INFILES(1)%TFILE,'YHAT',XYHAT,IRESP2)
!If not found in main file, try the fallback one !If not found in main file, try the fallback one
if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'YHAT',XYHAT,IRESP2) if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'YHAT',XYHAT,IRESP2)
if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'YHAT not found') if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'YHAT not found')
! Interpolations of positions to mass points
XXHATM(1:IIU-1) = 0.5 * XXHAT(1:IIU-1) + 0.5 * XXHAT(2:IIU)
XXHATM( IIU) = 1.5 * XXHAT( IIU) - 0.5 * XXHAT(IIU-1)
XYHATM(1:IJU-1) = 0.5 * XYHAT(1:IJU-1) + 0.5 * XYHAT(2:IJU)
XYHATM( IJU) = 1.5 * XYHAT( IJU) - 0.5 * XYHAT(IJU-1)
! Set global domain boundaries (hypothesis: only 1 process)
XXHAT_ll => XXHAT
XYHAT_ll => XYHAT
XXHATM_ll => XXHATM
XYHATM_ll => XYHATM
CALL IO_Field_read(INFILES(1)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2) CALL IO_Field_read(INFILES(1)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2)
!If not found in main file, try the fallback one !If not found in main file, try the fallback one
if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2) if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2)
...@@ -1146,12 +1174,16 @@ END DO ...@@ -1146,12 +1174,16 @@ END DO
! !
IF (TRIM(CPROGRAM_ORIG)/='PGD' .AND. TRIM(CPROGRAM_ORIG)/='NESPGD' .AND. TRIM(CPROGRAM_ORIG)/='ZOOMPG' & IF (TRIM(CPROGRAM_ORIG)/='PGD' .AND. TRIM(CPROGRAM_ORIG)/='NESPGD' .AND. TRIM(CPROGRAM_ORIG)/='ZOOMPG' &
.AND. .NOT.(TRIM(CPROGRAM_ORIG)=='REAL' .AND. CSTORAGE_TYPE=='SU') ) THEN !condition to detect PREP_SURFEX .AND. .NOT.(TRIM(CPROGRAM_ORIG)=='REAL' .AND. CSTORAGE_TYPE=='SU') ) THEN !condition to detect PREP_SURFEX
ALLOCATE(XZHAT(NKMAX+2*JPVEXT)) iku = nkmax + 2 * JPVEXT
ALLOCATE( XZHAT(iku), XZHATM(iku) )
CALL IO_Field_read(INFILES(1)%TFILE,'ZHAT',XZHAT,IRESP2) CALL IO_Field_read(INFILES(1)%TFILE,'ZHAT',XZHAT,IRESP2)
!If not found in main file, try the fallback one !If not found in main file, try the fallback one
if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'ZHAT',XZHAT,IRESP2) if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'ZHAT',XZHAT,IRESP2)
if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'ZHAT not found') if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'ZHAT not found')
! Interpolations of vertical positions to mass points
CALL INTERP_VERGRID_TO_MASSPOINTS( XZHAT, XZHATM )
ALLOCATE(LSLEVE) ALLOCATE(LSLEVE)
CALL IO_Field_read(INFILES(1)%TFILE,'SLEVE',LSLEVE,IRESP2) CALL IO_Field_read(INFILES(1)%TFILE,'SLEVE',LSLEVE,IRESP2)
!If not found in main file, try the fallback one !If not found in main file, try the fallback one
...@@ -1264,7 +1296,7 @@ END DO ...@@ -1264,7 +1296,7 @@ END DO
CALL IO_FILE_OPEN(outfiles(idx)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG) CALL IO_FILE_OPEN(outfiles(idx)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG)
END IF END IF
PRINT *,'--> Converted to file: ', TRIM(houtfile) call Print_msg( NVERB_INFO, 'IO', 'parse_infiles', '--> Converted to file: ' // trim(houtfile) )
END SUBROUTINE OPEN_FILES END SUBROUTINE OPEN_FILES
......
MY_RUN/BENCH/003_MEDATL_FULLHD/003_mesonh/run_mesonh_xyz_SEGN.ihm
View file @ af717e63
#MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier #MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
#MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence #MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
#MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt #MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
#MNH_LIC for details. version 1. #MNH_LIC for details. version 1.
...@@ -54,7 +54,6 @@ cat > EXSEG1.nam << EOF ...@@ -54,7 +54,6 @@ cat > EXSEG1.nam << EOF
&NAM_DYNn &NAM_DYNn
XTSTEP = ${XTSTEP} , XTSTEP = ${XTSTEP} ,
CPRESOPT = "CRESI", CPRESOPT = "CRESI",
NITR = 12,
LITRADJ = T, LITRADJ = T,
XRELAX = 1., XRELAX = 1.,
LHORELAX_UVWTH = T, LHORELAX_UVWTH = T,
......
MY_RUN/INTEGRATION_CASES/HPC/16JAN/PYTHON/plot_16JAN.py 0 → 100644
View file @ af717e63
#!/usr/bin/env python3
"""
@author: Quentin Rodier
Creation : 07/01/2021
Last modifications
"""
import matplotlib as mpl
mpl.use('Agg')
import cartopy.crs as ccrs
from read_MNHfile import read_netcdf
from Panel_Plot import PanelPlot
import os
os.system('rm -f tempgraph*')
#
# User's parameter / Namelist
#
#
path="../RUN/"
LnameFiles = ['16JA1.1.WENO5.004dg.nc', '16JA2.2.WENO5.002dg.nc']
Dvar_input = {
'f1':['MRV700HPA','THT850HPA','UT850HPA','VT850HPA','UT700HPA','VT700HPA', 'ALT_PRESSURE','ALT_U','ALT_V', 'ZS', 'latitude', 'longitude'],
'f2':['MRV700HPA','THT850HPA','UT850HPA','VT850HPA','UT700HPA','VT700HPA', 'ALT_PRESSURE', 'ZS', 'ALT_U','ALT_V','latitude', 'longitude']
}
# Read the variables in the files
Dvar = {}
Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
################################################################
######### PANEL 1
###############################################################
Panel1 = PanelPlot(2,2, [15,10],'16 JAN domaine 1 ', minmaxpad=1.05, minmaxTextSize=7.7, colorbaraspect=14)
Lplot = [ Dvar['f1']['ZS'],Dvar['f1']['THT850HPA'], Dvar['f1']['MRV700HPA'],Dvar['f1']['ALT_PRESSURE']]
lon = [Dvar['f1']['longitude']]*len(Lplot)
lat = [Dvar['f1']['latitude']]*len(Lplot)
Ltitle = ['Orography', 'Potential Temperature at 850hPa', 'Water vapor mixing at 700hPa','Pressure at z = 9000m']
Lcbarlabel = ['m','K', 'g/kg', 'hPa']
Lxlab = ['longitude']*len(Lplot)
Lylab = ['latitude']*len(Lplot)
Lminval = [0, 280, 0.5, 270]
Lmaxval = [1500, 293, 4.2, 303]
Lstep = [10, 0.25, 0.1, 1]
Lstepticks = [500, 2, 0.5, 5]
Lfacconv = [1.0, 1.0, 1.0, 1./100.0]
Lcolormap = ['terrain', 'gist_rainbow_r', 'gist_rainbow_r', 'gist_rainbow_r']
Lprojection = [ccrs.PlateCarree()]*len(Lplot)
Llvl = [0, 0, 0, 0]
fig1 = Panel1.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv)
Lplot1 = [ Dvar['f1']['UT850HPA'], Dvar['f1']['UT700HPA'], Dvar['f1']['ALT_U']]
Lplot2 = [ Dvar['f1']['VT850HPA'], Dvar['f1']['VT700HPA'], Dvar['f1']['ALT_V']]
Ltitle = ['Wind at 850hPa', 'Wind at 700hPa', 'Wind at 9000m']
Lxlab = ['longitude']*len(Lplot1)
Lylab = ['latitude']*len(Lplot1)
Llegendval = [20,20,40]
Llegendlabel = ['(m/s)']*len(Lplot1)
Larrowstep = [2]*len(Lplot1)
Lwidth = [0.002]*len(Lplot1)
Lcolor = ['black']*len(Lplot1)
Lprojection = [ccrs.PlateCarree()]*len(Lplot1)
Llvl = [0]*len(Lplot1)
fig2 = Panel1.pvector(Lxx=lon, Lyy=lat, Lvar1=Lplot1, Lvar2=Lplot2, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab,
Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep, Lproj=Lprojection,
Lcolor=Lcolor, Llegendval=Llegendval, Llegendlabel=Llegendlabel, Lid_overlap=[2,4,6], ax=fig1.axes)
Panel1.save_graph(1,fig2)
################################################################
######### PANEL 2
###############################################################
Panel2 = PanelPlot(2,2, [15,10],'16 JAN domaine 2 ', minmaxpad=1.05, minmaxTextSize=8, colorbaraspect=14)
Lplot = [ Dvar['f2']['ZS'],Dvar['f2']['THT850HPA'], Dvar['f2']['MRV700HPA'],Dvar['f2']['ALT_PRESSURE']]
lon = [Dvar['f2']['longitude']]*len(Lplot)
lat = [Dvar['f2']['latitude']]*len(Lplot)
Ltitle = ['Orography', 'Potential Temperature at 850hPa', 'Water vapor mixing at 700hPa','Pressure at z = 9000m']
Lcbarlabel = ['m','K', 'g/kg', 'hPa']
Lxlab = ['longitude']*len(Lplot)
Lylab = ['latitude']*len(Lplot)
Lminval = [0, 285, 0.9, 280]
Lmaxval = [2600, 293, 4, 294]
Lstep = [10, 0.5, 0.1, 0.5]
Lstepticks = [500, 1, 0.2, 5]
Lfacconv = [1.0, 1.0, 1.0, 1./100.0]
Lcolormap = ['terrain', 'gist_rainbow_r', 'gist_rainbow_r', 'gist_rainbow_r']
Lprojection = [ccrs.PlateCarree()]*len(Lplot)
Llvl = [0]*len(Lplot)
fig1 = Panel2.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv)
Lplot1 = [ Dvar['f2']['UT850HPA'], Dvar['f2']['UT700HPA'], Dvar['f2']['ALT_U']]
Lplot2 = [ Dvar['f2']['VT850HPA'], Dvar['f2']['VT700HPA'], Dvar['f2']['ALT_V']]
Ltitle = ['Wind at 850hPa', 'Wind at 700hPa', 'Wind at 9000m']
Llegendval = [20,20,40]
Lxlab = ['longitude']*len(Lplot1)
Lylab = ['latitude']*len(Lplot1)
Llegendlabel = ['(m/s)']*len(Lplot1)
Larrowstep = [2]*len(Lplot1)
Lwidth = [0.002]*len(Lplot1)
Lcolor = ['black']*len(Lplot1)
Lprojection = [ccrs.PlateCarree()]*len(Lplot1)
Llvl = [0]*len(Lplot1)
fig2 = Panel2.pvector(Lxx=lon, Lyy=lat, Lvar1=Lplot1, Lvar2=Lplot2, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab,
Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep, Lproj=Lprojection,
Lcolor=Lcolor, Llegendval=Llegendval, Llegendlabel=Llegendlabel, Lid_overlap=[2,4,6], ax=fig1.axes)
Panel2.save_graph(2,fig2)
MY_RUN/INTEGRATION_CASES/HPC/16JAN/PYTHON/run_python 0 → 100755
View file @ af717e63
#!/bin/sh
#SBATCH -J tools
#SBATCH -N 1 # nodes number
#SBATCH -n 1 # CPUs number (on all nodes)
#SBATCH -o tools.eo%j #
#SBATCH -e tools.eo%j #
#SBATCH -t 01:00:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
ln -sf ${SRC_MESONH}/src/LIB/Python/* .
module purge
module load python/3.7.6
python3 plot_16JAN.py
convert *.png 16JAN.pdf
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_1 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_DIAG
CISO='TKPREV',
LVAR_RS=T,
LVAR_TURB=F,
NCONV_KF=1,
LVAR_MRW=T,
LVAR_MRSV=F,
LTRAJ=F,
LTPZH=T,
CBLTOP="RICHA",
LMSLP=T,
LAGEO=T,
LTHW=T,
LCLD_COV=T,
NCAPE=1,
LRADAR=T,
LDIAG(:)=.FALSE.,
LISOAL=T, XISOAL(1)=9000,
LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
&NAM_DIAG_FILE YINIFILE(1) = "16JA1.1.WENO5.004" ,
YINIFILEPGD(1) = "16JAN98_36km.neste1",
YSUFFIX='dg' /
&NAM_DIAG_SURFn N2M=2 LSURF_BUDGET=T /
&NAM_DIAG_ISBAn LPGD=F LSURF_EVAP_BUDGET=T /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_2 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_DIAG CISO='', LVAR_RS=F, LVAR_LS=F, LVAR_PR=T, LTOTAL_PR=T,
LISOAL=T, XISOAL(1)=9000,
LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
&NAM_BLANK /
&NAM_DIAG_FILE YINIFILE(1) = "16JA2.2.WENO5.002",
YINIFILEPGD(1) = "16JAN98_9km.neste1",
YSUFFIX = "dg" /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG1.nam_run1 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_LUNITn CINIFILE = "15JAN_12_MNH" ,
CINIFILEPGD = "16JAN98_36km.neste1",
CCPLFILE(1) = "15JAN_18_MNH",
CCPLFILE(2) = "16JAN_00_MNH" /
&NAM_CONFZ
NZ_VERB=5 , NZ_PROC=0 , NB_PROCIO_R=1 , NB_PROCIO_W=1
/
&NAM_DYNn XTSTEP = 120., CPRESOPT = "ZRESI", NITR = 8,
LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000.,XT4DIFTH = 5000. /
&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K" ,CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
&NAM_PARAMn CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "KAFR", CCLOUD= "KESS",CSCONV="KAFR" /
&NAM_PARAM_RADn XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR= 400 /
&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T, LDOWN = T /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM",
LTURB_DIAG = F, LTURB_FLX = F /
&NAM_CONF CCONF = "START", NMODEL = 1, NVERB = 5,
CEXP = "16JA1 ", CSEG = "WENO5",NHALO=3 /
&NAM_DYN XSEGLEN = 43200., LCORIO = T, LNUMDIFU = F,LNUMDIFTH = F,
XALKTOP = 0.001, XALZBOT = 15000. /
&NAM_BACKUP XBAK_TIME(1,1)=10800., XBAK_TIME(1,2)=21600.,XBAK_TIME(1,3)=32400.,XBAK_TIME(1,4)=43200./
&NAM_DIAG_SURFn /
&NAM_ISBA CALBEDO = "DRY" /
&NAM_ISBAn CSCOND = "NP89",
CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
&NAM_SSOn CROUGH="Z01D"/
&NAM_DIAG_ISBAn /
&NAM_SEAFLUXn CSEA_ALB="UNIF" /
&NAM_BLANK /
&NAM_SGH_ISBAn CRUNOFF = "WSAT"/
&NAM_NEBn LSUBG_COND= F /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG1.nam_run2 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_LUNITn CINIFILE = "16JA1.1.WENO5.004",
CINIFILEPGD = "16JAN98_36km.neste1",
CCPLFILE(1) = "16JAN_06_MNH" /
&NAM_DYNn XTSTEP = 120., CPRESOPT = "ZRESI", NITR = 8,
LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000. /
&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
&NAM_PARAMn CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "KAFR", CSCONV="KAFR",CCLOUD= "KESS" /
&NAM_PARAM_RADn XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR = 400 /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
LDOWN = T /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM" /
&NAM_CH_MNHCn LUSECHEM = F /
&NAM_CONF CCONF = "START", NVERB=5, NMODEL = 2,
CEXP = "16JA2", CSEG = "WENO5",NHALO=3 /
&NAM_DYN XSEGLEN = 21600., LCORIO = T, LNUMDIFU = F,
XALKTOP = 0.001, XALZBOT = 14500. /
&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 2, XWAY(2) = 2. /
&NAM_BACKUP XBAK_TIME(1,1) = 10800., XBAK_TIME(1,2) = 21600. ,
XBAK_TIME(2,1) = 10800., XBAK_TIME(2,2) = 21600. /
&NAM_DIAG_SURFn /
&NAM_ISBA CALBEDO = "DRY" /
&NAM_ISBAn CSCOND = "NP89",
CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
&NAM_SSOn CROUGH="Z01D"/
&NAM_DIAG_ISBAn /
&NAM_SGH_ISBAn
CRUNOFF = "WSAT" /
&NAM_SEAFLUXn CSEA_ALB="UNIF" /
&NAM_BLANK /
&NAM_NEBn LSUBG_COND = F /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG2.nam_run2 0 → 100644
View file @ af717e63
&NAM_LUNITn CINIFILE = "16JAN_00_MNH2",
CINIFILEPGD = "16JAN98_9km.neste1" /
&NAM_DYNn CPRESOPT = "ZRESI", NITR = 8, XRELAX = 1.,
LHORELAX_UVWTH = F, LHORELAX_RV = F, LHORELAX_RC= F, LHORELAX_RR= F,
LHORELAX_RS= F, LHORELAX_RI= F, LHORELAX_RG= F, LHORELAX_TKE= F,
NRIMX = 0, NRIMY = 0, LVE_RELAX = T, XT4DIFU = 1500. /
&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
&NAM_PARAMn CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW", CDCONV = "KAFR" ,
CSCONV="KAFR" /
&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
NRAD_COLNBR = 400 /
&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
LDOWN = T /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
&NAM_TURBn XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
LTURB_FLX = F, LTURB_DIAG = F /
&NAM_CH_MNHCn LUSECHEM = F /
&NAM_DIAG_SURFn /
&NAM_ISBA CALBEDO = "DRY" /
&NAM_ISBAn CSCOND = "NP89",
CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
&NAM_SSOn CROUGH="Z01D"/
&NAM_DIAG_ISBAn /
&NAM_SGH_ISBAn
CRUNOFF = "WSAT" /
&NAM_SEAFLUXn CSEA_ALB="UNIF" /
&NAM_BLANK /
&NAM_NEBn LSUBG_COND = F /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_NEST_PGD1.nam 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_PGD1 YPGD1 = '16JAN98_36km' /
&NAM_PGD2 YPGD2 = '16JAN98_9km', IDAD = 1 /
&NAM_NEST_PGD YNEST = 'e1' /
"
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_PGD1.nam_1 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_PGDFILE CPGDFILE='16JAN98_36km' /
&NAM_CONF_PROJ XLAT0=50., XLON0=-2.8, XRPK=0.58, XBETA=0. /
&NAM_CONF_PROJ_GRID XLATCEN=48., XLONCEN=-1.,
NIMAX=80, NJMAX=60, XDX=36000., XDY=36000. /
&NAM_COVER YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
&NAM_ZS YZS='gtopo30', YZSFILETYPE='DIRECT' /
&NAM_ISBA YCLAY='CLAY_HWSD_MOY', YCLAYFILETYPE='DIRECT' ,
YSAND='SAND_HWSD_MOY', YSANDFILETYPE='DIRECT' /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_PGD1.nam_2 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_PGDFILE CPGDFILE='16JAN98_9km' ,NHALO=20/
&NAM_PGD_GRID YINIFILE='16JAN98_36km' YINIFILETYPE='MESONH' /
&NAM_INIFILE_CONF_PROJ IXOR=40, IYOR=16, IXSIZE=18, IYSIZE=15, IDXRATIO=4, IDYRATIO=4 /
&NAM_COVER YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
&NAM_ZS YZS='gtopo30', YZSFILETYPE='DIRECT' /
&NAM_ISBA YCLAY='CLAY_HWSD_MOY', YCLAYFILETYPE='DIRECT' ,
YSAND='SAND_HWSD_MOY', YSANDFILETYPE='DIRECT' /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_1 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980115.12' , HATMFILETYPE='GRIBEX',
HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='15JAN_12_MNH' /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
ZDZGRD=60., ZDZTOP=700.,
ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_2 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980115.18' , HATMFILETYPE='GRIBEX',
HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='15JAN_18_MNH' /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
ZDZGRD=60., ZDZTOP=700.,
ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_3 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980116.00' , HATMFILETYPE='GRIBEX',
HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='16JAN_00_MNH' /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
ZDZGRD=60., ZDZTOP=700.,
ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_4 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=T, LLFIOUT=F, LLFIREAD=F /
&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980116.06' , HATMFILETYPE='GRIBEX',
HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='16JAN_06_MNH' /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
ZDZGRD=60., ZDZTOP=700.,
ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_dom2 0 → 100644
View file @ af717e63
&NAM_CONFIO LCDF4=F LLFIOUT=F LLFIREAD=F /
&NAM_FILE_NAMES HATMFILE ='16JA1.1.WENO5.004.spa04' , HATMFILETYPE='MESONH',
HPGDFILE ='16JAN98_9km.neste1' , CINIFILE='16JAN_00_MNH2' /
&NAM_REAL_CONF NVERB=10 /
&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
ZDZGRD=60., ZDZTOP=700.,
ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /

Documentation | Mentions légales | CGU | Accessibilité : non conforme