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

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
7 results
Show changes
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 1585 additions and 0 deletions
MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/clean_mesonh_xyz 0 → 100644
View file @ 6bb35ca8
set -x
rm -f PGD* CPL* AZF* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex *.dat *.bin
MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/run_mesonh 0 → 100755
View file @ 6bb35ca8
#!/bin/sh
#SBATCH -J azf
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 16 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run1_16jan.eo%j #
#SBATCH -e run1_16jan.eo%j #
#SBATCH -t 02:00:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
# Nom de la machine
hostname
ln -sf ../006_preal/PGD* .
ln -sf ../006_preal/CPL* .
ln -sf ../004_arp2lfi/CPL* .
. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 8"
set -x
set -e
ls -lrt
time ${MPIRUN} MESONH${XYZ}
cd ../009_python
sbatch run_python
MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/009_python/plot_AZF2M.py 0 → 100644
View file @ 6bb35ca8
#!/usr/bin/env python3
"""
@author: Quentin Rodier
Creation : 07/01/2021
Last modifications
"""
import matplotlib as mpl
mpl.use('Agg')
from read_MNHfile import read_netcdf
from Panel_Plot import PanelPlot
from misc_functions import mean_operator, convert_date
import cartopy.crs as ccrs
import numpy as np
import os
import cartopy.io.shapereader as shpreader
import matplotlib.patches as mpatches
os.system('rm -f tempgraph*')
#
# User's parameter / Namelist
#
path=""
LnameFiles = ['AZF02.1.CEN4T.001.nc', 'AZF02.1.CEN4T.002.nc', 'AZF02.1.CEN4T.003.nc',
'AZF02.1.CEN4T.004.nc', 'AZF02.1.CEN4T.005.nc', 'AZF02.1.CEN4T.007.nc',
'AZF02.2.CEN4T.001.nc', 'AZF02.2.CEN4T.002.nc', 'AZF02.2.CEN4T.003.nc',
'AZF02.2.CEN4T.004.nc', 'AZF02.2.CEN4T.005.nc', 'AZF02.2.CEN4T.007.nc',
'AZF02.1.CEN4T.000.nc']
LG_AVION='/Flyers/Aircrafts/AVION/'
LG_AVIONT='/Flyers/Aircrafts/AVION/Point/'
LG_AVIONZT='/Flyers/Aircrafts/AVION/Vertical_profile/'
Dvar_input = {
'f1':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f2':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f3':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f4':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f5':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f6':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f7':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level','LONOR','LATOR','LAT','LON'],
'f8':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f9':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f10':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f11':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f12':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
'f13':[(LG_AVION,'time_flyer'),(LG_AVIONT,'ZS'), (LG_AVIONT,'P'), (LG_AVIONT,'LON'),(LG_AVIONT,'MER_WIND'),
(LG_AVIONT,'ZON_WIND'),(LG_AVIONT,'W'), (LG_AVIONT,'Th'), (LG_AVIONT,'Rv'),(LG_AVIONT,'Tke'),
(LG_AVIONT,'H_FLUX'),(LG_AVIONT,'LE_FLUX'), (LG_AVIONT,'Tke_Diss'), (LG_AVIONT,'Tsrad')]
}
# Read the variables in the files
Dvar = {}
Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
################################################################
######### PANEL 1
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 1 SV 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
Lplot = [Dvar['f1']['SVPP001'], Dvar['f2']['SVPP001'], Dvar['f3']['SVPP001'],
Dvar['f4']['SVPP001'], Dvar['f5']['SVPP001'], Dvar['f6']['SVPP001']]
lon = [Dvar['f1']['longitude']]*len(Lplot)
lat = [Dvar['f1']['latitude']]*len(Lplot)
Ltitle = ['SVPP001']*len(Lplot)
Lcbarlabel = ['kg/kg']*len(Lplot)
Lxlab = ['longitude']*len(Lplot)
Lylab = ['latitude']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.15E-5]*len(Lplot)
Lstep = [0.05E-6]*len(Lplot)
Lstepticks = [0.2E-6]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
Lprojection = [ccrs.PlateCarree()]*len(Lplot)
LaddWhite = [True]*len(Lplot)
Llevel = [0]*len(Lplot)
Ltime = [Dvar['f1']['date'], Dvar['f2']['date'], Dvar['f3']['date'], Dvar['f4']['date'], Dvar['f5']['date'], Dvar['f6']['date']]
Lcbformatlabel=[True]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
Lplot1 = [ Dvar['f1']['UT'], Dvar['f2']['UT'], Dvar['f3']['UT'], Dvar['f4']['UT'], Dvar['f5']['UT'], Dvar['f6']['UT']]
Lplot2 = [ Dvar['f1']['VT'], Dvar['f2']['VT'], Dvar['f3']['VT'], Dvar['f4']['VT'], Dvar['f5']['VT'], Dvar['f6']['VT']]
Ltitle = ['wind vectors at K=2']*len(Lplot)
Llegendval = [7.5]*len(Lplot)
Lcbarlabel = ['(m/s)']*len(Lplot1)
Larrowstep = [2]*len(Lplot1)
Lwidth = [0.002]*len(Lplot1)
Lcolor = ['black']*len(Lplot1)
Lscale = [100]*len(Lplot1)
fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
# Departements francais
departements_shp='departements-20180101.shp'
adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
# Add departements to each axes + scatter point of emission source
loncar, latcar = [1.439,1.5], [43.567, 43.9]
label=['AZF1','AZF2']
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
# Add a Rectangle displaying the domain of the model 2
for i in range(len(Lplot)):
fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
fig.tight_layout()
Panel.save_graph(1,fig)
################################################################
######### PANEL 2
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 1 SV 002', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
Lplot = [Dvar['f1']['SVPP002'], Dvar['f2']['SVPP002'], Dvar['f3']['SVPP002'],
Dvar['f4']['SVPP002'], Dvar['f5']['SVPP002'], Dvar['f6']['SVPP002']]
lon = [Dvar['f1']['longitude']]*len(Lplot)
lat = [Dvar['f1']['latitude']]*len(Lplot)
Ltitle = ['SVPP001']*len(Lplot)
Lcbarlabel = ['kg/kg']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.15E-5]*len(Lplot)
Lstep = [0.05E-6]*len(Lplot)
Lstepticks = [0.2E-6]*len(Lplot)
Lcbformatlabel=[True]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
# Add departements to each axes + scatter point of emission source
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
# Add a Rectangle displaying the domain of the model 2
for i in range(len(Lplot)):
fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
fig.tight_layout()
Panel.save_graph(2,fig)
################################################################
######### PANEL 3
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 1 ATC 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
Lplot = [Dvar['f1']['ATC001'], Dvar['f2']['ATC001'], Dvar['f3']['ATC001'],
Dvar['f4']['ATC001'], Dvar['f5']['ATC001'], Dvar['f6']['ATC001']]
Ltitle = ['ATC001']*len(Lplot)
Lcbarlabel = ['$m^{-3}$']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.15E-6]*len(Lplot)
Lstep = [0.05E-7]*len(Lplot)
Lstepticks = [0.2E-7]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
# Add departements to each axes + scatter point of emission source
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
# Add a Rectangle displaying the domain of the model 2
for i in range(len(Lplot)):
fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
fig.tight_layout()
Panel.save_graph(3,fig)
################################################################
######### PANEL 4
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 1 ATC 002', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
Lplot = [Dvar['f1']['ATC002'], Dvar['f2']['ATC002'], Dvar['f3']['ATC002'],
Dvar['f4']['ATC002'], Dvar['f5']['ATC002'], Dvar['f6']['ATC002']]
Ltitle = ['ATC002']*len(Lplot)
Lcbarlabel = ['$m^{-3}$']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.15E-6]*len(Lplot)
Lstep = [0.05E-7]*len(Lplot)
Lstepticks = [0.2E-7]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
# Add departements to each axes + scatter point of emission source
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
# Add a Rectangle displaying the domain of the model 2
for i in range(len(Lplot)):
fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
fig.tight_layout()
Panel.save_graph(4,fig)
################################################################
######### PANEL 5 : Domaine fils
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 2 SV 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
Lplot = [Dvar['f7']['SVPP001'], Dvar['f8']['SVPP001'], Dvar['f9']['SVPP001'],
Dvar['f10']['SVPP001'], Dvar['f11']['SVPP001'], Dvar['f12']['SVPP001']]
lon = [Dvar['f7']['longitude']]*len(Lplot)
lat = [Dvar['f7']['latitude']]*len(Lplot)
Ltitle = ['SVPP001']*len(Lplot)
Lcbarlabel = ['kg/kg']*len(Lplot)
Lxlab = ['longitude']*len(Lplot)
Lylab = ['latitude']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.15E-5]*len(Lplot)
Lstep = [0.05E-6]*len(Lplot)
Lstepticks = [0.2E-6]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
Lprojection = [ccrs.PlateCarree()]*len(Lplot)
LaddWhite = [True]*len(Lplot)
Llevel = [0]*len(Lplot)
Ltime = [Dvar['f7']['date'], Dvar['f8']['date'], Dvar['f9']['date'], Dvar['f10']['date'], Dvar['f11']['date'], Dvar['f12']['date']]
Lcbformatlabel=[True]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
Lplot1 = [ Dvar['f7']['UT'], Dvar['f8']['UT'], Dvar['f9']['UT'], Dvar['f10']['UT'], Dvar['f11']['UT'], Dvar['f12']['UT']]
Lplot2 = [ Dvar['f7']['VT'], Dvar['f8']['VT'], Dvar['f9']['VT'], Dvar['f10']['VT'], Dvar['f11']['VT'], Dvar['f12']['VT']]
Ltitle = ['wind vectors at K=2']*len(Lplot)
Llegendval = [7.5]*len(Lplot)
Lcbarlabel = ['(m/s)']*len(Lplot1)
Larrowstep = [4]*len(Lplot1)
Lwidth = [0.002]*len(Lplot1)
Lcolor = ['black']*len(Lplot1)
Lscale = [75]*len(Lplot1)
fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
# Departements francais
departements_shp='departements-20180101.shp'
adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
# Add departements to each axes + scatter point of emission source
loncar, latcar = [1.439,1.5], [43.567, 43.9]
label=['AZF1','AZF2']
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
fig.tight_layout()
Panel.save_graph(5,fig)
################################################################
######### PANEL 6
###############################################################
Panel = PanelPlot(2,3, [25,14],'Domaine 2 ATC 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
Lplot = [Dvar['f7']['ATC001'], Dvar['f8']['ATC001'], Dvar['f9']['ATC001'],
Dvar['f10']['ATC001'], Dvar['f11']['ATC001'], Dvar['f12']['ATC001']]
Ltitle = ['ATC001']*len(Lplot)
Lcbarlabel = ['$m^{-3}$']*len(Lplot)
Lminval = [0]*len(Lplot)
Lmaxval = [0.6E-6]*len(Lplot)
Lstep = [0.01E-6]*len(Lplot)
Lstepticks = [0.1E-6]*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel,
Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
# Add departements to each axes + scatter point of emission source
for i in range(len(Lplot)):
fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
fig.axes[i*2].scatter(loncar,latcar)
for lab, txt in enumerate(label):
fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
fig.tight_layout()
Panel.save_graph(6,fig)
################################################################
######### PANEL 7
###############################################################
Panel = PanelPlot(8,2, [14,20],'Time series from Aircraft', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
Lplot = [ Dvar['f13'][(LG_AVIONT,'ZS')]]
Ltime = [Dvar['f13'][(LG_AVION,'time_flyer')]/3600.0]
Ltitle = ['Orography']
Lxlab = ['Time (h)']
Lylab = ['ZS (m)']
Lylim = [(0, 350)]
Lxlim = [(9.0, 9.2)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim)
Lplot = [ Dvar['f13'][(LG_AVIONT,'P')]]
Ltitle = ['Pressure']
Lylab = ['P (Pa)']
Lylim = [(0, 95000)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'LON')]]
Ltitle = ['Longitude']
Lylab = ['Longitude']
Lylim = [(0, 2.5)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'ZON_WIND')]]
Ltitle = ['Zonal wind']
Lylab = ['u (m/s)']
Lylim = [(-1, 11)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'MER_WIND')]]
Ltitle = ['Meridional wind']
Lylab = ['v (m/s)']
Lylim = [(-3, 3)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'W')]]
Ltitle = ['Vertical velocity']
Lylab = ['w (m/s)']
Lylim = [(-0.1, 0.1)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'Th')]]
Ltitle = ['Potential Temperature']
Lylab = [r'$\theta$ (K)']
Lylim = [(290, 305)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'Rv')]]
Ltitle = ['Water vapor mixing ratio']
Lylab = ['Rv (kg/kg))']
Lylim = [(0, 0.01)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'Tke')]]
Ltitle = ['Turbulent Kinetic Energy']
Lylab = ['TKE ($m^2s^{-2}$)']
Lylim = [(0, 0.1)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'Tke_Diss')]]
Ltitle = ['Turbulent Kinetic Energy Dissipation']
Lylab = ['TKE Diss ($m^2s^{-2}$']
Lylim = [(0, 1000)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'H_FLUX')]]
Ltitle = ['Sensible Heat Flux H']
Lylab = ['H ($W/m^2$)']
Lylim = [(-0.7, 0.)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'LE_FLUX')]]
Ltitle = ['Latent Heat Flux LE']
Lylab = ['LE ($W/m^2$)']
Lylim = [(0, 2.0)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
Lplot = [ Dvar['f13'][(LG_AVIONT,'Tsrad')]]
Ltitle = ['Radiative surface temperature']
Lylab = ['Tsrad (K))']
Lylim = [(250, 1000)]
fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
fig.tight_layout()
Panel.save_graph(7,fig)
MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/009_python/run_python 0 → 100755
View file @ 6bb35ca8
#!/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 02: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/* .
ln -sf ~rodierq/SAVE/OUTILS/PYTHON/departements-20180101.* .
ln -sf ../007_run/AZF02.*.CEN4T.*.nc .
module purge
module load python/3.7.6
python3 plot_AZF2M.py
convert *.png AZF_2M.pdf
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/001_prep_pgd1/PRE_PGD1.nam 0 → 100644
View file @ 6bb35ca8
&NAM_CONFIO LCDF4=.TRUE.,
LLFIOUT=.FALSE.,
LLFIREAD=.FALSE. /
&NAM_PGDFILE CPGDFILE='BIOMAIDO_CH_8km', NHALO=0 /
&NAM_COVER YCOVER='ECOCLIMAP_v2.0',
YCOVERFILETYPE='DIRECT' /
&NAM_ZS YZS='gtopo30',
YZSFILETYPE='DIRECT' /
&NAM_ZSFILTER NZSFILTER=0 /
&NAM_PGD_GRID CGRID='CONF PROJ' /
&NAM_CONF_PROJ XLAT0=-20.93,
XLON0=55.370,
XRPK=0.,
XBETA=0. /
&NAM_CONF_PROJ_GRID XLATCEN=-20.93,
XLONCEN=55.370,
NIMAX=60,
NJMAX=60,
XDX=8000.,
XDY=8000. /
&NAM_ISBA CPHOTO='NIT',
YCLAY='CLAY_HWSD_MOY_v2',
YCLAYFILETYPE='DIRECT',
YSAND='SAND_HWSD_MOY_v2',
YSANDFILETYPE='DIRECT',
CISBA ="3-L",
NPATCH=12 /
&NAM_CH_EMISSIONS CCH_EMIS = 'AGGR', CCH_BIOEMIS = 'MEGA', CCH_DMSEMIS= 'DMSD' /
&NAM_DMS_PGD NDMS_NBR = 12,
CDMS_NAME(1) = 'DMS_JAN',
CDMS_FILE(1) = 'data_DMS_january.dat',
CDMS_FILETYPE(1) = 'ASCLLV',
CDMS_AREA(1) = 'SEA',
CDMS_NAME(2) = 'DMS_FEV',
CDMS_FILE(2) = 'data_DMS_february.dat',
CDMS_FILETYPE(2) = 'ASCLLV',
CDMS_AREA(2) = 'SEA',
CDMS_NAME(3) = 'DMS_MAR',
CDMS_FILE(3) = 'data_DMS_march.dat',
CDMS_FILETYPE(3) = 'ASCLLV',
CDMS_AREA(3) = 'SEA',
CDMS_NAME(4) = 'DMS_APR',
CDMS_FILE(4) = 'data_DMS_april.dat',
CDMS_FILETYPE(4) = 'ASCLLV',
CDMS_AREA(4) = 'SEA',
CDMS_NAME(5) = 'DMS_MAY',
CDMS_FILE(5) = 'data_DMS_may.dat',
CDMS_FILETYPE(5) = 'ASCLLV',
CDMS_AREA(5) = 'SEA',
CDMS_NAME(6) = 'DMS_JUN',
CDMS_FILE(6) = 'data_DMS_june.dat',
CDMS_FILETYPE(6) = 'ASCLLV',
CDMS_AREA(6) = 'SEA',
CDMS_NAME(7) = 'DMS_JUL',
CDMS_FILE(7) = 'data_DMS_july.dat',
CDMS_FILETYPE(7) = 'ASCLLV',
CDMS_AREA(7) = 'SEA',
CDMS_NAME(8) = 'DMS_AUG',
CDMS_FILE(8) = 'data_DMS_august.dat',
CDMS_FILETYPE(8) = 'ASCLLV',
CDMS_AREA(9) = 'SEA',
CDMS_NAME(9) = 'DMS_SEP',
CDMS_FILE(9) = 'data_DMS_september.dat',
CDMS_FILETYPE(9) = 'ASCLLV',
CDMS_AREA(10) = 'SEA',
CDMS_NAME(10) = 'DMS_OCT',
CDMS_FILE(10) = 'data_DMS_october.dat',
CDMS_FILETYPE(10) = 'ASCLLV',
CDMS_AREA(11) = 'SEA',
CDMS_NAME(11) = 'DMS_NOV',
CDMS_FILE(11) = 'data_DMS_november.dat',
CDMS_FILETYPE(11) = 'ASCLLV',
CDMS_AREA(11) = 'SEA',
CDMS_NAME(12) = 'DMS_DEC',
CDMS_FILE(12) = 'data_DMS_december.dat',
CDMS_FILETYPE(12) = 'ASCLLV',
CDMS_AREA(12) = 'SEA' /
&NAM_CH_EMIS_PGD NEMIS_PGD_NBR = 38,
CEMIS_PGD_NAME(1) = 'COAN',
NEMIS_PGD_TIME(1) = 0,
CEMIS_PGD_FILE(1) = 'MACCity_CO_juil2014.asc',
CEMIS_PGD_FILETYPE(1) = 'ASCLLV',
CEMIS_PGD_AREA(1) = 'LAN',
CEMIS_PGD_NAME(2) = 'NOXAN',
NEMIS_PGD_TIME(2) = 0,
CEMIS_PGD_FILE(2) = 'MACCity_NOX_juil2014.asc',
CEMIS_PGD_FILETYPE(2) = 'ASCLLV',
CEMIS_PGD_AREA(2) = 'LAN',
CEMIS_PGD_NAME(3) = 'C2H6AN',
NEMIS_PGD_TIME(3) = 0,
CEMIS_PGD_FILE(3) = 'MACCity_C2H6_juil2014.asc',
CEMIS_PGD_FILETYPE(3) = 'ASCLLV',
CEMIS_PGD_AREA(3) = 'LAN',
CEMIS_PGD_NAME(4) = 'C2H4AN',
NEMIS_PGD_TIME(4) = 0,
CEMIS_PGD_FILE(4) = 'MACCity_C2H4_juil2014.asc',
CEMIS_PGD_FILETYPE(4) = 'ASCLLV',
CEMIS_PGD_AREA(4) = 'LAN',
CEMIS_PGD_NAME(5) = 'C3H8AN',
NEMIS_PGD_TIME(5) = 0,
CEMIS_PGD_FILE(5) = 'MACCity_C3H8_juil2014.asc',
CEMIS_PGD_FILETYPE(5) = 'ASCLLV',
CEMIS_PGD_AREA(5) = 'LAN',
CEMIS_PGD_NAME(6) = 'C3H6AN',
NEMIS_PGD_TIME(6) = 0,
CEMIS_PGD_FILE(6) = 'MACCity_C3H6_juil2014.asc',
CEMIS_PGD_FILETYPE(6) = 'ASCLLV',
CEMIS_PGD_AREA(6) = 'LAN',
CEMIS_PGD_NAME(7) = 'BALKAN',
NEMIS_PGD_TIME(7) = 0,
CEMIS_PGD_FILE(7) = 'MACCity_BIGALK_juil2014.asc',
CEMIS_PGD_FILETYPE(7) = 'ASCLLV',
CEMIS_PGD_AREA(7) = 'LAN',
CEMIS_PGD_NAME(8) = 'BENEAN',
NEMIS_PGD_TIME(8) = 0,
CEMIS_PGD_FILE(8) = 'MACCity_BIGENE_juil2014.asc',
CEMIS_PGD_FILETYPE(8) = 'ASCLLV',
CEMIS_PGD_AREA(8) = 'LAN',
CEMIS_PGD_NAME(9) = 'CH2OAN',
NEMIS_PGD_TIME(9) = 0,
CEMIS_PGD_FILE(9) = 'MACCity_CH2O_juil2014.asc',
CEMIS_PGD_FILETYPE(9) = 'ASCLLV',
CEMIS_PGD_AREA(9) = 'LAN',
CEMIS_PGD_NAME(10) = 'BALDAN',
NEMIS_PGD_TIME(10) = 0,
CEMIS_PGD_FILE(10) = 'MACCity_BIGALD_juil2014.asc',
CEMIS_PGD_FILETYPE(10) = 'ASCLLV',
CEMIS_PGD_AREA(10) = 'LAN',
CEMIS_PGD_NAME(11) = 'MTOLAN',
NEMIS_PGD_TIME(11) = 0,
CEMIS_PGD_FILE(11) = 'MACCity_CH3OH_juil2014.asc',
CEMIS_PGD_FILETYPE(11) = 'ASCLLV',
CEMIS_PGD_AREA(11) = 'LAN',
CEMIS_PGD_NAME(12) = 'BOLAN',
NEMIS_PGD_TIME(12) = 0,
CEMIS_PGD_FILE(12) = 'MACCity_BIGOL_juil2014.asc',
CEMIS_PGD_FILETYPE(12) = 'ASCLLV',
CEMIS_PGD_AREA(12) = 'LAN',
CEMIS_PGD_NAME(13) = 'ACETAN',
NEMIS_PGD_TIME(13) = 0,
CEMIS_PGD_FILE(13) = 'MACCity_ACETONE_juil2014.asc',
CEMIS_PGD_FILETYPE(13) = 'ASCLLV',
CEMIS_PGD_AREA(13) = 'LAN',
CEMIS_PGD_NAME(14) = 'BONEAN',
NEMIS_PGD_TIME(14) = 0,
CEMIS_PGD_FILE(14) = 'MACCity_BIGONE_juil2014.asc',
CEMIS_PGD_FILETYPE(14) = 'ASCLLV',
CEMIS_PGD_AREA(14) = 'LAN',
CEMIS_PGD_NAME(15) = 'AROAN',
NEMIS_PGD_TIME(15) = 0,
CEMIS_PGD_FILE(15) = 'MACCity_ARO_juil2014.asc',
CEMIS_PGD_FILETYPE(15) = 'ASCLLV',
CEMIS_PGD_AREA(15) = 'LAN',
CEMIS_PGD_NAME(16) = 'SO2AN',
NEMIS_PGD_TIME(16) = 0,
CEMIS_PGD_FILE(16) = 'MACCity_SO2_juil2014.asc',
CEMIS_PGD_FILETYPE(16) = 'ASCLLV',
CEMIS_PGD_AREA(16) = 'LAN',
NEMIS_PGD_TIME(17) = 0,
CEMIS_PGD_NAME(17) = 'NH3AN',
CEMIS_PGD_FILE(17) = 'MACCity_NH3_juil2014.asc',
CEMIS_PGD_FILETYPE(17) = 'ASCLLV',
CEMIS_PGD_AREA(17) = 'LAN',
CEMIS_PGD_NAME(18) = 'BCAN',
NEMIS_PGD_TIME(18) = 0,
CEMIS_PGD_FILE(18) = 'MACCity_BC_juil2014.asc',
CEMIS_PGD_FILETYPE(18) = 'ASCLLV',
CEMIS_PGD_AREA(18) = 'LAN',
CEMIS_PGD_NAME(19) = 'OCAN',
NEMIS_PGD_TIME(19) = 0,
CEMIS_PGD_FILE(19) = 'MACCity_OC_juil2014.asc',
CEMIS_PGD_FILETYPE(19) = 'ASCLLV',
CEMIS_PGD_AREA(19) = 'LAN',
CEMIS_PGD_NAME(20) = 'COANSEA',
NEMIS_PGD_TIME(20) = 0,
CEMIS_PGD_FILE(20) = 'MACCity_CO_juil2014.asc',
CEMIS_PGD_FILETYPE(20) = 'ASCLLV',
CEMIS_PGD_AREA(20) = 'SEA',
CEMIS_PGD_NAME(21) = 'NOXANSEA',
NEMIS_PGD_TIME(21) = 0,
CEMIS_PGD_FILE(21) = 'MACCity_NOX_juil2014.asc',
CEMIS_PGD_FILETYPE(21) = 'ASCLLV',
CEMIS_PGD_AREA(21) = 'SEA',
CEMIS_PGD_NAME(22) = 'C2H6ANSEA',
NEMIS_PGD_TIME(22) = 0,
CEMIS_PGD_FILE(22) = 'MACCity_C2H6_juil2014.asc',
CEMIS_PGD_FILETYPE(22) = 'ASCLLV',
CEMIS_PGD_AREA(22) = 'SEA',
CEMIS_PGD_NAME(23) = 'C2H4ANSEA',
NEMIS_PGD_TIME(23) = 0,
CEMIS_PGD_FILE(23) = 'MACCity_C2H4_juil2014.asc',
CEMIS_PGD_FILETYPE(23) = 'ASCLLV',
CEMIS_PGD_AREA(23) = 'SEA',
CEMIS_PGD_NAME(24) = 'C3H8ANSEA',
NEMIS_PGD_TIME(24) = 0,
CEMIS_PGD_FILE(24) = 'MACCity_C3H8_juil2014.asc',
CEMIS_PGD_FILETYPE(24) = 'ASCLLV',
CEMIS_PGD_AREA(24) = 'SEA',
CEMIS_PGD_NAME(25) = 'C3H6ANSEA',
NEMIS_PGD_TIME(25) = 0,
CEMIS_PGD_FILE(25) = 'MACCity_C3H6_juil2014.asc',
CEMIS_PGD_FILETYPE(25) = 'ASCLLV',
CEMIS_PGD_AREA(25) = 'SEA',
CEMIS_PGD_NAME(26) = 'BALKANSEA',
NEMIS_PGD_TIME(26) = 0,
CEMIS_PGD_FILE(26) = 'MACCity_BIGALK_juil2014.asc',
CEMIS_PGD_FILETYPE(26) = 'ASCLLV',
CEMIS_PGD_AREA(26) = 'SEA',
CEMIS_PGD_NAME(27) = 'BENEANSEA',
NEMIS_PGD_TIME(27) = 0,
CEMIS_PGD_FILE(27) = 'MACCity_BIGENE_juil2014.asc',
CEMIS_PGD_FILETYPE(27) = 'ASCLLV',
CEMIS_PGD_AREA(27) = 'SEA',
CEMIS_PGD_NAME(28) = 'CH2OANSEA',
NEMIS_PGD_TIME(28) = 0,
CEMIS_PGD_FILE(28) = 'MACCity_CH2O_juil2014.asc',
CEMIS_PGD_FILETYPE(28) = 'ASCLLV',
CEMIS_PGD_AREA(28) = 'SEA',
CEMIS_PGD_NAME(29) = 'BALDANSEA',
NEMIS_PGD_TIME(29) = 0,
CEMIS_PGD_FILE(29) = 'MACCity_BIGALD_juil2014.asc',
CEMIS_PGD_FILETYPE(29) = 'ASCLLV',
CEMIS_PGD_AREA(29) = 'SEA',
CEMIS_PGD_NAME(30) = 'MTOLANSEA',
NEMIS_PGD_TIME(30) = 0,
CEMIS_PGD_FILE(30) = 'MACCity_CH3OH_juil2014.asc',
CEMIS_PGD_FILETYPE(30) = 'ASCLLV',
CEMIS_PGD_AREA(30) = 'SEA',
CEMIS_PGD_NAME(31) = 'BOLANSEA',
NEMIS_PGD_TIME(31) = 0,
CEMIS_PGD_FILE(31) = 'MACCity_BIGOL_juil2014.asc',
CEMIS_PGD_FILETYPE(31) = 'ASCLLV',
CEMIS_PGD_AREA(31) = 'SEA',
CEMIS_PGD_NAME(32) = 'ACETANSEA',
NEMIS_PGD_TIME(32) = 0,
CEMIS_PGD_FILE(32) = 'MACCity_ACETONE_juil2014.asc',
CEMIS_PGD_FILETYPE(32) = 'ASCLLV',
CEMIS_PGD_AREA(32) = 'SEA',
CEMIS_PGD_NAME(33) = 'BONEANSEA',
NEMIS_PGD_TIME(33) = 0,
CEMIS_PGD_FILE(33) = 'MACCity_BIGONE_juil2014.asc',
CEMIS_PGD_FILETYPE(33) = 'ASCLLV',
CEMIS_PGD_AREA(33) = 'SEA',
CEMIS_PGD_NAME(34) = 'AROANSEA',
NEMIS_PGD_TIME(34) = 0,
CEMIS_PGD_FILE(34) = 'MACCity_ARO_juil2014.asc',
CEMIS_PGD_FILETYPE(34) = 'ASCLLV',
CEMIS_PGD_AREA(34) = 'SEA',
CEMIS_PGD_NAME(35) = 'SO2ANSEA',
NEMIS_PGD_TIME(35) = 0,
CEMIS_PGD_FILE(35) = 'MACCity_SO2_juil2014.asc',
CEMIS_PGD_FILETYPE(35) = 'ASCLLV',
CEMIS_PGD_AREA(35) = 'SEA',
CEMIS_PGD_NAME(36) = 'NH3ANSEA',
NEMIS_PGD_TIME(36) = 0,
CEMIS_PGD_FILE(36) = 'MACCity_NH3_juil2014.asc',
CEMIS_PGD_FILETYPE(36) = 'ASCLLV',
CEMIS_PGD_AREA(36) = 'SEA',
CEMIS_PGD_NAME(37) = 'BCANSEA',
NEMIS_PGD_TIME(37) = 0,
CEMIS_PGD_FILE(37) = 'MACCity_BC_juil2014.asc',
CEMIS_PGD_FILETYPE(37) = 'ASCLLV',
CEMIS_PGD_AREA(37) = 'SEA',
CEMIS_PGD_NAME(38) = 'OCANSEA',
NEMIS_PGD_TIME(38) = 0,
CEMIS_PGD_FILE(38) = 'MACCity_OC_juil2014.asc',
CEMIS_PGD_FILETYPE(38) = 'ASCLLV',
CEMIS_PGD_AREA(38) = 'SEA' /
&NAM_MEGAN_PGD NMEGAN_NBR = 15,
CMEGAN_NAME(1) = 'EFISOP',
CMEGAN_FILE(1) = 'isop100.asc',
CMEGAN_FILETYPE(1) = 'ASCLLV',
CMEGAN_AREA(1) = 'LAN',
CMEGAN_NAME(2) = 'EFAPIN',
CMEGAN_FILE(2) = 'apin100.asc',
CMEGAN_FILETYPE(2) = 'ASCLLV',
CMEGAN_AREA(2) = 'LAN',
CMEGAN_NAME(3) = 'EFBPIN',
CMEGAN_FILE(3) = 'bpin100.asc',
CMEGAN_FILETYPE(3) = 'ASCLLV',
CMEGAN_AREA(3) = 'LAN',
CMEGAN_NAME(4) = 'EFCARE',
CMEGAN_FILE(4) = 'bcar100.asc',
CMEGAN_FILETYPE(4) = 'ASCLLV',
CMEGAN_AREA(4) = 'LAN',
CMEGAN_NAME(5) = 'EFLIMO',
CMEGAN_FILE(5) = 'limo100.asc',
CMEGAN_FILETYPE(5) = 'ASCLLV',
CMEGAN_AREA(5) = 'LAN',
CMEGAN_NAME(6) = 'EFMYRC',
CMEGAN_FILE(6) = 'myrc100.asc',
CMEGAN_FILETYPE(6) = 'ASCLLV',
CMEGAN_AREA(6) = 'LAN',
CMEGAN_NAME(7) = 'EFOCIM',
CMEGAN_FILE(7) = 'ocim100.asc',
CMEGAN_FILETYPE(7) = 'ASCLLV',
CMEGAN_AREA(7) = 'LAN',
CMEGAN_NAME(8) = 'EFSABI',
CMEGAN_FILE(8) = 'sabi100.asc',
CMEGAN_FILETYPE(8) = 'ASCLLV',
CMEGAN_AREA(8) = 'LAN',
CMEGAN_NAME(10) = 'LAI',
CMEGAN_FILE(10) = 'lai100.asc',
CMEGAN_FILETYPE(10) = 'ASCLLV',
CMEGAN_AREA(10) = 'LAN'
CMEGAN_NAME(11) = 'PFT1',
CMEGAN_FILE(11) = 'pft1.asc',
CMEGAN_FILETYPE(11) = 'ASCLLV',
CMEGAN_AREA(11) = 'LAN'
CMEGAN_NAME(12) = 'PFT4',
CMEGAN_FILE(12) = 'pft4.asc',
CMEGAN_FILETYPE(12) = 'ASCLLV',
CMEGAN_AREA(12) = 'LAN'
CMEGAN_NAME(13) = 'PFT9',
CMEGAN_FILE(13) = 'pft9.asc',
CMEGAN_FILETYPE(13) = 'ASCLLV',
CMEGAN_AREA(13) = 'LAN'
CMEGAN_NAME(14) = 'PFT13',
CMEGAN_FILE(14) = 'pft13.asc',
CMEGAN_FILETYPE(14) = 'ASCLLV',
CMEGAN_AREA(14) = 'LAN'
CMEGAN_NAME(15) = 'PFT14',
CMEGAN_FILE(15) = 'pft14.asc',
CMEGAN_FILETYPE(15) = 'ASCLLV',
CMEGAN_AREA(15) = 'LAN'
CMEGAN_NAME(9) = 'PFTru',
CMEGAN_FILE(9) = 'pftru.asc',
CMEGAN_FILETYPE(9) = 'ASCLLV',
CMEGAN_AREA(9) = 'LAN' /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/001_prep_pgd1/run_pgd.sh 0 → 100755
View file @ 6bb35ca8
#!/bin/sh
#SBATCH -J pgd_megan1
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 8 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 00:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
# Nom de la machine
hostname
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
ln -sf ~/SAVE/mesonh/PGD/* .
ln -sf ~/SAVE/EMISSIONS_MEGAN/DATA/DATA_EMIS/* .
ln -sf ~/SAVE/CHIMIE_FILES/EMISSIONS/* .
ln -sf ~/SAVE/EMISSIONS_MEGAN/* .
ln -sf ~/SAVE/EMISSIONS_DMS/* .
export MPIRUN="Mpirun -np 8"
set -x
set -e
time ${MPIRUN} PREP_PGD${XYZ}
cd ../002_prep_pgd2
sbatch run_pgd.sh
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/002_prep_pgd2/PRE_PGD1.nam 0 → 100644
View file @ 6bb35ca8
&NAM_CONFIO LCDF4=.TRUE.,
LLFIOUT=.FALSE.,
LLFIREAD=.FALSE. /
&NAM_INIFILE_CONF_PROJ IXOR=23,
IYOR=23,
IXSIZE=15,
IYSIZE=15,
IDXRATIO=4,
IDYRATIO=4 /
&NAM_PGD_GRID YINIFILE='BIOMAIDO_CH_8km',
YINIFILETYPE='MESONH' /
&NAM_PGDFILE CPGDFILE='BIOMAIDO_CH_2km', NHALO= 0 /
&NAM_CONF_PROJ_GRID /
&NAM_COVER YCOVER='ECOCLIMAP_v2.0',
YCOVERFILETYPE='DIRECT' /
&NAM_ZS YZS='gtopo30',
YZSFILETYPE='DIRECT' /
&NAM_ISBA CPHOTO='NIT',
YCLAY='CLAY_HWSD_MOY_v2',
YCLAYFILETYPE='DIRECT',
YSAND='SAND_HWSD_MOY_v2',
YSANDFILETYPE='DIRECT',
CISBA ="3-L",
NPATCH=12 /
&NAM_CH_EMISSIONS CCH_EMIS = 'AGGR', CCH_BIOEMIS = 'MEGA', CCH_DMSEMIS= 'DMSD' /
&NAM_DMS_PGD NDMS_NBR = 12,
CDMS_NAME(1) = 'DMS_JAN',
CDMS_FILE(1) = 'data_DMS_january.dat',
CDMS_FILETYPE(1) = 'ASCLLV',
CDMS_AREA(1) = 'SEA',
CDMS_NAME(2) = 'DMS_FEV',
CDMS_FILE(2) = 'data_DMS_february.dat',
CDMS_FILETYPE(2) = 'ASCLLV',
CDMS_AREA(2) = 'SEA',
CDMS_NAME(3) = 'DMS_MAR',
CDMS_FILE(3) = 'data_DMS_march.dat',
CDMS_FILETYPE(3) = 'ASCLLV',
CDMS_AREA(3) = 'SEA',
CDMS_NAME(4) = 'DMS_APR',
CDMS_FILE(4) = 'data_DMS_april.dat',
CDMS_FILETYPE(4) = 'ASCLLV',
CDMS_AREA(4) = 'SEA',
CDMS_NAME(5) = 'DMS_MAY',
CDMS_FILE(5) = 'data_DMS_may.dat',
CDMS_FILETYPE(5) = 'ASCLLV',
CDMS_AREA(5) = 'SEA',
CDMS_NAME(6) = 'DMS_JUN',
CDMS_FILE(6) = 'data_DMS_june.dat',
CDMS_FILETYPE(6) = 'ASCLLV',
CDMS_AREA(6) = 'SEA',
CDMS_NAME(7) = 'DMS_JUL',
CDMS_FILE(7) = 'data_DMS_july.dat',
CDMS_FILETYPE(7) = 'ASCLLV',
CDMS_AREA(7) = 'SEA',
CDMS_NAME(8) = 'DMS_AUG',
CDMS_FILE(8) = 'data_DMS_august.dat',
CDMS_FILETYPE(8) = 'ASCLLV',
CDMS_AREA(9) = 'SEA',
CDMS_NAME(9) = 'DMS_SEP',
CDMS_FILE(9) = 'data_DMS_september.dat',
CDMS_FILETYPE(9) = 'ASCLLV',
CDMS_AREA(10) = 'SEA',
CDMS_NAME(10) = 'DMS_OCT',
CDMS_FILE(10) = 'data_DMS_october.dat',
CDMS_FILETYPE(10) = 'ASCLLV',
CDMS_AREA(11) = 'SEA',
CDMS_NAME(11) = 'DMS_NOV',
CDMS_FILE(11) = 'data_DMS_november.dat',
CDMS_FILETYPE(11) = 'ASCLLV',
CDMS_AREA(11) = 'SEA',
CDMS_NAME(12) = 'DMS_DEC',
CDMS_FILE(12) = 'data_DMS_december.dat',
CDMS_FILETYPE(12) = 'ASCLLV',
CDMS_AREA(12) = 'SEA' /
&NAM_CH_EMIS_PGD NEMIS_PGD_NBR = 38,
CEMIS_PGD_NAME(1) = 'COAN',
NEMIS_PGD_TIME(1) = 0,
CEMIS_PGD_FILE(1) = 'MACCity_CO_juil2014.asc',
CEMIS_PGD_FILETYPE(1) = 'ASCLLV',
CEMIS_PGD_AREA(1) = 'LAN',
CEMIS_PGD_NAME(2) = 'NOXAN',
NEMIS_PGD_TIME(2) = 0,
CEMIS_PGD_FILE(2) = 'MACCity_NOX_juil2014.asc',
CEMIS_PGD_FILETYPE(2) = 'ASCLLV',
CEMIS_PGD_AREA(2) = 'LAN',
CEMIS_PGD_NAME(3) = 'C2H6AN',
NEMIS_PGD_TIME(3) = 0,
CEMIS_PGD_FILE(3) = 'MACCity_C2H6_juil2014.asc',
CEMIS_PGD_FILETYPE(3) = 'ASCLLV',
CEMIS_PGD_AREA(3) = 'LAN',
CEMIS_PGD_NAME(4) = 'C2H4AN',
NEMIS_PGD_TIME(4) = 0,
CEMIS_PGD_FILE(4) = 'MACCity_C2H4_juil2014.asc',
CEMIS_PGD_FILETYPE(4) = 'ASCLLV',
CEMIS_PGD_AREA(4) = 'LAN',
CEMIS_PGD_NAME(5) = 'C3H8AN',
NEMIS_PGD_TIME(5) = 0,
CEMIS_PGD_FILE(5) = 'MACCity_C3H8_juil2014.asc',
CEMIS_PGD_FILETYPE(5) = 'ASCLLV',
CEMIS_PGD_AREA(5) = 'LAN',
CEMIS_PGD_NAME(6) = 'C3H6AN',
NEMIS_PGD_TIME(6) = 0,
CEMIS_PGD_FILE(6) = 'MACCity_C3H6_juil2014.asc',
CEMIS_PGD_FILETYPE(6) = 'ASCLLV',
CEMIS_PGD_AREA(6) = 'LAN',
CEMIS_PGD_NAME(7) = 'BALKAN',
NEMIS_PGD_TIME(7) = 0,
CEMIS_PGD_FILE(7) = 'MACCity_BIGALK_juil2014.asc',
CEMIS_PGD_FILETYPE(7) = 'ASCLLV',
CEMIS_PGD_AREA(7) = 'LAN',
CEMIS_PGD_NAME(8) = 'BENEAN',
NEMIS_PGD_TIME(8) = 0,
CEMIS_PGD_FILE(8) = 'MACCity_BIGENE_juil2014.asc',
CEMIS_PGD_FILETYPE(8) = 'ASCLLV',
CEMIS_PGD_AREA(8) = 'LAN',
CEMIS_PGD_NAME(9) = 'CH2OAN',
NEMIS_PGD_TIME(9) = 0,
CEMIS_PGD_FILE(9) = 'MACCity_CH2O_juil2014.asc',
CEMIS_PGD_FILETYPE(9) = 'ASCLLV',
CEMIS_PGD_AREA(9) = 'LAN',
CEMIS_PGD_NAME(10) = 'BALDAN',
NEMIS_PGD_TIME(10) = 0,
CEMIS_PGD_FILE(10) = 'MACCity_BIGALD_juil2014.asc',
CEMIS_PGD_FILETYPE(10) = 'ASCLLV',
CEMIS_PGD_AREA(10) = 'LAN',
CEMIS_PGD_NAME(11) = 'MTOLAN',
NEMIS_PGD_TIME(11) = 0,
CEMIS_PGD_FILE(11) = 'MACCity_CH3OH_juil2014.asc',
CEMIS_PGD_FILETYPE(11) = 'ASCLLV',
CEMIS_PGD_AREA(11) = 'LAN',
CEMIS_PGD_NAME(12) = 'BOLAN',
NEMIS_PGD_TIME(12) = 0,
CEMIS_PGD_FILE(12) = 'MACCity_BIGOL_juil2014.asc',
CEMIS_PGD_FILETYPE(12) = 'ASCLLV',
CEMIS_PGD_AREA(12) = 'LAN',
CEMIS_PGD_NAME(13) = 'ACETAN',
NEMIS_PGD_TIME(13) = 0,
CEMIS_PGD_FILE(13) = 'MACCity_ACETONE_juil2014.asc',
CEMIS_PGD_FILETYPE(13) = 'ASCLLV',
CEMIS_PGD_AREA(13) = 'LAN',
CEMIS_PGD_NAME(14) = 'BONEAN',
NEMIS_PGD_TIME(14) = 0,
CEMIS_PGD_FILE(14) = 'MACCity_BIGONE_juil2014.asc',
CEMIS_PGD_FILETYPE(14) = 'ASCLLV',
CEMIS_PGD_AREA(14) = 'LAN',
CEMIS_PGD_NAME(15) = 'AROAN',
NEMIS_PGD_TIME(15) = 0,
CEMIS_PGD_FILE(15) = 'MACCity_ARO_juil2014.asc',
CEMIS_PGD_FILETYPE(15) = 'ASCLLV',
CEMIS_PGD_AREA(15) = 'LAN',
CEMIS_PGD_NAME(16) = 'SO2AN',
NEMIS_PGD_TIME(16) = 0,
CEMIS_PGD_FILE(16) = 'MACCity_SO2_juil2014.asc',
CEMIS_PGD_FILETYPE(16) = 'ASCLLV',
CEMIS_PGD_AREA(16) = 'LAN',
NEMIS_PGD_TIME(17) = 0,
CEMIS_PGD_NAME(17) = 'NH3AN',
CEMIS_PGD_FILE(17) = 'MACCity_NH3_juil2014.asc',
CEMIS_PGD_FILETYPE(17) = 'ASCLLV',
CEMIS_PGD_AREA(17) = 'LAN',
CEMIS_PGD_NAME(18) = 'BCAN',
NEMIS_PGD_TIME(18) = 0,
CEMIS_PGD_FILE(18) = 'MACCity_BC_juil2014.asc',
CEMIS_PGD_FILETYPE(18) = 'ASCLLV',
CEMIS_PGD_AREA(18) = 'LAN',
CEMIS_PGD_NAME(19) = 'OCAN',
NEMIS_PGD_TIME(19) = 0,
CEMIS_PGD_FILE(19) = 'MACCity_OC_juil2014.asc',
CEMIS_PGD_FILETYPE(19) = 'ASCLLV',
CEMIS_PGD_AREA(19) = 'LAN',
CEMIS_PGD_NAME(20) = 'COANSEA',
NEMIS_PGD_TIME(20) = 0,
CEMIS_PGD_FILE(20) = 'MACCity_CO_juil2014.asc',
CEMIS_PGD_FILETYPE(20) = 'ASCLLV',
CEMIS_PGD_AREA(20) = 'SEA',
CEMIS_PGD_NAME(21) = 'NOXANSEA',
NEMIS_PGD_TIME(21) = 0,
CEMIS_PGD_FILE(21) = 'MACCity_NOX_juil2014.asc',
CEMIS_PGD_FILETYPE(21) = 'ASCLLV',
CEMIS_PGD_AREA(21) = 'SEA',
CEMIS_PGD_NAME(22) = 'C2H6ANSEA',
NEMIS_PGD_TIME(22) = 0,
CEMIS_PGD_FILE(22) = 'MACCity_C2H6_juil2014.asc',
CEMIS_PGD_FILETYPE(22) = 'ASCLLV',
CEMIS_PGD_AREA(22) = 'SEA',
CEMIS_PGD_NAME(23) = 'C2H4ANSEA',
NEMIS_PGD_TIME(23) = 0,
CEMIS_PGD_FILE(23) = 'MACCity_C2H4_juil2014.asc',
CEMIS_PGD_FILETYPE(23) = 'ASCLLV',
CEMIS_PGD_AREA(23) = 'SEA',
CEMIS_PGD_NAME(24) = 'C3H8ANSEA',
NEMIS_PGD_TIME(24) = 0,
CEMIS_PGD_FILE(24) = 'MACCity_C3H8_juil2014.asc',
CEMIS_PGD_FILETYPE(24) = 'ASCLLV',
CEMIS_PGD_AREA(24) = 'SEA',
CEMIS_PGD_NAME(25) = 'C3H6ANSEA',
NEMIS_PGD_TIME(25) = 0,
CEMIS_PGD_FILE(25) = 'MACCity_C3H6_juil2014.asc',
CEMIS_PGD_FILETYPE(25) = 'ASCLLV',
CEMIS_PGD_AREA(25) = 'SEA',
CEMIS_PGD_NAME(26) = 'BALKANSEA',
NEMIS_PGD_TIME(26) = 0,
CEMIS_PGD_FILE(26) = 'MACCity_BIGALK_juil2014.asc',
CEMIS_PGD_FILETYPE(26) = 'ASCLLV',
CEMIS_PGD_AREA(26) = 'SEA',
CEMIS_PGD_NAME(27) = 'BENEANSEA',
NEMIS_PGD_TIME(27) = 0,
CEMIS_PGD_FILE(27) = 'MACCity_BIGENE_juil2014.asc',
CEMIS_PGD_FILETYPE(27) = 'ASCLLV',
CEMIS_PGD_AREA(27) = 'SEA',
CEMIS_PGD_NAME(28) = 'CH2OANSEA',
NEMIS_PGD_TIME(28) = 0,
CEMIS_PGD_FILE(28) = 'MACCity_CH2O_juil2014.asc',
CEMIS_PGD_FILETYPE(28) = 'ASCLLV',
CEMIS_PGD_AREA(28) = 'SEA',
CEMIS_PGD_NAME(29) = 'BALDANSEA',
NEMIS_PGD_TIME(29) = 0,
CEMIS_PGD_FILE(29) = 'MACCity_BIGALD_juil2014.asc',
CEMIS_PGD_FILETYPE(29) = 'ASCLLV',
CEMIS_PGD_AREA(29) = 'SEA',
CEMIS_PGD_NAME(30) = 'MTOLANSEA',
NEMIS_PGD_TIME(30) = 0,
CEMIS_PGD_FILE(30) = 'MACCity_CH3OH_juil2014.asc',
CEMIS_PGD_FILETYPE(30) = 'ASCLLV',
CEMIS_PGD_AREA(30) = 'SEA',
CEMIS_PGD_NAME(31) = 'BOLANSEA',
NEMIS_PGD_TIME(31) = 0,
CEMIS_PGD_FILE(31) = 'MACCity_BIGOL_juil2014.asc',
CEMIS_PGD_FILETYPE(31) = 'ASCLLV',
CEMIS_PGD_AREA(31) = 'SEA',
CEMIS_PGD_NAME(32) = 'ACETANSEA',
NEMIS_PGD_TIME(32) = 0,
CEMIS_PGD_FILE(32) = 'MACCity_ACETONE_juil2014.asc',
CEMIS_PGD_FILETYPE(32) = 'ASCLLV',
CEMIS_PGD_AREA(32) = 'SEA',
CEMIS_PGD_NAME(33) = 'BONEANSEA',
NEMIS_PGD_TIME(33) = 0,
CEMIS_PGD_FILE(33) = 'MACCity_BIGONE_juil2014.asc',
CEMIS_PGD_FILETYPE(33) = 'ASCLLV',
CEMIS_PGD_AREA(33) = 'SEA',
CEMIS_PGD_NAME(34) = 'AROANSEA',
NEMIS_PGD_TIME(34) = 0,
CEMIS_PGD_FILE(34) = 'MACCity_ARO_juil2014.asc',
CEMIS_PGD_FILETYPE(34) = 'ASCLLV',
CEMIS_PGD_AREA(34) = 'SEA',
CEMIS_PGD_NAME(35) = 'SO2ANSEA',
NEMIS_PGD_TIME(35) = 0,
CEMIS_PGD_FILE(35) = 'MACCity_SO2_juil2014.asc',
CEMIS_PGD_FILETYPE(35) = 'ASCLLV',
CEMIS_PGD_AREA(35) = 'SEA',
CEMIS_PGD_NAME(36) = 'NH3ANSEA',
NEMIS_PGD_TIME(36) = 0,
CEMIS_PGD_FILE(36) = 'MACCity_NH3_juil2014.asc',
CEMIS_PGD_FILETYPE(36) = 'ASCLLV',
CEMIS_PGD_AREA(36) = 'SEA',
CEMIS_PGD_NAME(37) = 'BCANSEA',
NEMIS_PGD_TIME(37) = 0,
CEMIS_PGD_FILE(37) = 'MACCity_BC_juil2014.asc',
CEMIS_PGD_FILETYPE(37) = 'ASCLLV',
CEMIS_PGD_AREA(37) = 'SEA',
CEMIS_PGD_NAME(38) = 'OCANSEA',
NEMIS_PGD_TIME(38) = 0,
CEMIS_PGD_FILE(38) = 'MACCity_OC_juil2014.asc',
CEMIS_PGD_FILETYPE(38) = 'ASCLLV',
CEMIS_PGD_AREA(38) = 'SEA' /
&NAM_MEGAN_PGD NMEGAN_NBR = 15,
CMEGAN_NAME(1) = 'EFISOP',
CMEGAN_FILE(1) = 'isop100.asc',
CMEGAN_FILETYPE(1) = 'ASCLLV',
CMEGAN_AREA(1) = 'LAN',
CMEGAN_NAME(2) = 'EFAPIN',
CMEGAN_FILE(2) = 'apin100.asc',
CMEGAN_FILETYPE(2) = 'ASCLLV',
CMEGAN_AREA(2) = 'LAN',
CMEGAN_NAME(3) = 'EFBPIN',
CMEGAN_FILE(3) = 'bpin100.asc',
CMEGAN_FILETYPE(3) = 'ASCLLV',
CMEGAN_AREA(3) = 'LAN',
CMEGAN_NAME(4) = 'EFCARE',
CMEGAN_FILE(4) = 'bcar100.asc',
CMEGAN_FILETYPE(4) = 'ASCLLV',
CMEGAN_AREA(4) = 'LAN',
CMEGAN_NAME(5) = 'EFLIMO',
CMEGAN_FILE(5) = 'limo100.asc',
CMEGAN_FILETYPE(5) = 'ASCLLV',
CMEGAN_AREA(5) = 'LAN',
CMEGAN_NAME(6) = 'EFMYRC',
CMEGAN_FILE(6) = 'myrc100.asc',
CMEGAN_FILETYPE(6) = 'ASCLLV',
CMEGAN_AREA(6) = 'LAN',
CMEGAN_NAME(7) = 'EFOCIM',
CMEGAN_FILE(7) = 'ocim100.asc',
CMEGAN_FILETYPE(7) = 'ASCLLV',
CMEGAN_AREA(7) = 'LAN',
CMEGAN_NAME(8) = 'EFSABI',
CMEGAN_FILE(8) = 'sabi100.asc',
CMEGAN_FILETYPE(8) = 'ASCLLV',
CMEGAN_AREA(8) = 'LAN',
CMEGAN_NAME(10) = 'LAI',
CMEGAN_FILE(10) = 'lai100.asc',
CMEGAN_FILETYPE(10) = 'ASCLLV',
CMEGAN_AREA(10) = 'LAN'
CMEGAN_NAME(11) = 'PFT1',
CMEGAN_FILE(11) = 'pft1.asc',
CMEGAN_FILETYPE(11) = 'ASCLLV',
CMEGAN_AREA(11) = 'LAN'
CMEGAN_NAME(12) = 'PFT4',
CMEGAN_FILE(12) = 'pft4.asc',
CMEGAN_FILETYPE(12) = 'ASCLLV',
CMEGAN_AREA(12) = 'LAN'
CMEGAN_NAME(13) = 'PFT9',
CMEGAN_FILE(13) = 'pft9.asc',
CMEGAN_FILETYPE(13) = 'ASCLLV',
CMEGAN_AREA(13) = 'LAN'
CMEGAN_NAME(14) = 'PFT13',
CMEGAN_FILE(14) = 'pft13.asc',
CMEGAN_FILETYPE(14) = 'ASCLLV',
CMEGAN_AREA(14) = 'LAN'
CMEGAN_NAME(15) = 'PFT14',
CMEGAN_FILE(15) = 'pft14.asc',
CMEGAN_FILETYPE(15) = 'ASCLLV',
CMEGAN_AREA(15) = 'LAN'
CMEGAN_NAME(9) = 'PFTru',
CMEGAN_FILE(9) = 'pftru.asc',
CMEGAN_FILETYPE(9) = 'ASCLLV',
CMEGAN_AREA(9) = 'LAN' /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/002_prep_pgd2/run_pgd.sh 0 → 100755
View file @ 6bb35ca8
#!/bin/sh
#SBATCH -J pgd_megan1
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 8 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 00:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
# Nom de la machine
hostname
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
ln -sf ../001_prep_pgd1/BIOM* .
ln -sf ~/SAVE/mesonh/PGD/* .
ln -sf ~/SAVE/EMISSIONS_MEGAN/DATA/DATA_EMIS/* .
ln -sf ~/SAVE/CHIMIE_FILES/EMISSIONS/* .
ln -sf ~/SAVE/EMISSIONS_MEGAN/* .
ln -sf ~/SAVE/EMISSIONS_DMS/* .
export MPIRUN="Mpirun -np 8"
set -x
set -e
time ${MPIRUN} PREP_PGD${XYZ}
cd ../003_prep_nest
sbatch run_prep_nest.sh
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/003_prep_nest/PRE_NEST_PGD1.nam 0 → 100644
View file @ 6bb35ca8
&NAM_CONFIO LCDF4=.TRUE.,
LLFIOUT=.FALSE.,
LLFIREAD=.FALSE. /
&NAM_CONF NHALO=0 /
&NAM_PGD1 YPGD1 = 'BIOMAIDO_CH_8km' /
&NAM_PGD2 YPGD2 = 'BIOMAIDO_CH_2km', IDAD = 1 /
&NAM_NEST_PGD YNEST = 'ne' /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/003_prep_nest/run_prep_nest.sh 0 → 100755
View file @ 6bb35ca8
#!/bin/sh
#SBATCH -J pgd_megan1
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 1 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 00:10:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
# Nom de la machine
hostname
ln -sf ../001_prep_pgd1/BIOM* .
ln -sf ../002_prep_pgd2/BIOM* .
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 1"
set -x
set -e
time ${MPIRUN} PREP_NEST_PGD${XYZ}
cd ../004_prep_real_case_d1
sbatch run_prep_real_case.sh
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/CAM1.nam 0 → 100644
View file @ 6bb35ca8
CAM2MESONH
transfer CAMS/ReLACS variables (default values) included aerosols
37 # NUMBER OF OPTIONAL GRIB VARIABLES
(A6,1X,I1,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2)
O3 1 1.00 go3 48.00
CO 1 1.00 co 28.00
H2O2 1 1.00 h2o2 34.00
NO 1 1.00 no 30.00
NO2 1 1.00 no2 46.00
HNO3 1 1.00 hno3 63.00
HNO4 1 1.00 ho2no2 79.00
N2O5 1 1.00 n2o5 108.00
NO3 1 1.00 no3 62.00
HO2 1 1.00 ho2 33.00
SO2 1 1.00 so2 64.10
DMS 1.1.00 dms 62.10
NH3 1 1.00 nh3 17.00
MO2 1 1.00 ch3o2 47.00
ETH 1 1.00 c2h6 30.02
ALKA 4 0.39 c3h8 44.03 0.33 ch3oh 31.01 0.96 c2h5oh 46.02 1.03 par 12.00
ALKE 3 0.88 c2h4 28.00 1.19 c3h6 42.03 1.20 ole 24.00
BIO 2 1.00 c5h8 68.10 0.99 c10h16 136.00
KET 1 0.28 ch3coch3 58.00
CARBO 2 0.96 ch3cocho 72.10 1.07 ispd 70.00
PAN 1 1.00 pan 121.00
ONIT 1 1.00 onit 77.00
HCHO 1 1.00 hcho 30.00
ALD 1 1.00 ald2 24.00
ORA1 1 1.00 hcooh 46.01
OP1 1 1.00 ch3ooh 48.00
OP2 1 1.00 rooh 47.00
SO4I 1 0.20 aermr11 1.00
SO4J 1 0.80 aermr11 1.00
NH3I 1 0.20 nh4 1.00
NH3J 1 0.80 nh4 1.00
NO3I 1 0.20 no3_a 1.00
NO3J 1 0.80 no3_a 1.00
BCI 2 0.20 aermr10 1.00 0.20 aermr09 1.00
OCI 2 0.20 aermr08 1.00 0.20 aermr07 1.00
BCJ 2 0.80 aermr10 1.00 0.80 aermr09 1.00
OCJ 2 0.80 aermr08 1.00 0.80 aermr07 1.00
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/CAM1_relacs.nam 0 → 100644
View file @ 6bb35ca8
CAM2MESONH
transfer CAMS/ReLACS variables (default values) included aerosols
37 # NUMBER OF OPTIONAL GRIB VARIABLES
(A6,1X,I1,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2)
O3 1 1.00 go3 48.00
CO 1 1.00 co 28.00
H2O2 1 1.00 h2o2 34.00
NO 1 1.00 no 30.00
NO2 1 1.00 no2 46.00
HNO3 1 1.00 hno3 63.00
HNO4 1 1.00 ho2no2 79.00
N2O5 1 1.00 n2o5 108.00
NO3 1 1.00 no3 62.00
HO2 1 1.00 ho2 33.00
SO2 1 1.00 so2 64.10
DMS 1.1.00 dms 62.10
NH3 1 1.00 nh3 17.00
MO2 1 1.00 ch3o2 47.00
ETH 1 1.00 c2h6 30.02
ALKA 4 0.39 c3h8 44.03 0.33 ch3oh 31.01 0.96 c2h5oh 46.02 1.03 par 12.00
ALKE 3 0.88 c2h4 28.00 1.19 c3h6 42.03 1.20 ole 24.00
BIO 2 1.00 c5h8 68.10 0.99 c10h16 136.00
KET 1 0.28 ch3coch3 58.00
CARBO 2 0.96 ch3cocho 72.10 1.07 ispd 70.00
PAN 1 1.00 pan 121.00
ONIT 1 1.00 onit 77.00
HCHO 1 1.00 hcho 30.00
ALD 1 1.00 ald2 24.00
ORA1 1 1.00 hcooh 46.01
OP1 1 1.00 ch3ooh 48.00
OP2 1 1.00 rooh 47.00
SO4I 1 0.20 aermr11 1.00
SO4J 1 0.80 aermr11 1.00
NH3I 1 0.20 nh4 1.00
NH3J 1 0.80 nh4 1.00
NO3I 1 0.20 no3_a 1.00
NO3J 1 0.80 no3_a 1.00
BCI 2 0.20 aermr10 1.00 0.20 aermr09 1.00
OCI 2 0.20 aermr08 1.00 0.20 aermr07 1.00
BCJ 2 0.80 aermr10 1.00 0.80 aermr09 1.00
OCJ 2 0.80 aermr08 1.00 0.80 aermr07 1.00
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/PRE_REAL1.nam 0 → 100644
View file @ 6bb35ca8
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20190328.06.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_06.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_VER_GRID NKMAX=45,
LTHINSHELL=.TRUE.,
YZGRID_TYPE='FUNCTN',
ZDZGRD=60.,
ZDZTOP=1000.,
ZZMAX_STRGRD=500.,
ZSTRGRD=7.,
ZSTRTOP=10. /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/PRE_REAL1.nam_06 0 → 100644
View file @ 6bb35ca8
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20190328.06.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_06.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_VER_GRID NKMAX=45,
LTHINSHELL=.TRUE.,
YZGRID_TYPE='FUNCTN',
ZDZGRD=60.,
ZDZTOP=1000.,
ZZMAX_STRGRD=500.,
ZSTRGRD=7.,
ZSTRTOP=10. /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/run_prep_real_case.sh 0 → 100755
View file @ 6bb35ca8
#!/bin/sh
#SBATCH -J prep_megan
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 1 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 01:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
ln -sf ../003_prep_nest/BIOM* .
ln -sf ~/SAVE/GRIB_KTEST/CAMS_20190328_* .
ln -sf ~/SAVE/GRIB_KTEST/ecmwf.OD.20190328.* .
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 1"
set -x
set -e
cp PRE_REAL1.nam_06 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cp PRE_REAL1.nam_12 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cd ../005_prep_real_case_d2
sbatch run_prep_real_case.sh

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