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
Showing
with 1811 additions and 0 deletions
#!/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)
#!/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
&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' /
#!/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
&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' /
#!/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
&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' /
#!/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
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
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
&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 /
&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 /
&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.12.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_12.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_12'/
&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 /
#!/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
&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_D1_20190328_06.spa',
HATMFILETYPE='MESONH',
HPGDFILE='BIOMAIDO_CH_2km.nestne',
CINIFILE='ECMWF_D2_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&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 /
&NAM_LUNIT2_SPA CINIFILE = "ECMWF_D1_20190328_06",
CINIFILEPGD = "BIOMAIDO_CH_8km.nestne",
YDOMAIN = "BIOMAIDO_CH_2km.nestne",
YSPANBR = '' /
&NAM_CONFIO LCDF4 = T, NIO_ABORT_LEVEL=1, NGEN_ABORT_LEVEL=1 /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
#!/bin/sh
#SBATCH -J prep_d2_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 00: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 ../004_prep_real_case_d1/ECMWF* .
. ~/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} SPAWNING${XYZ}
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cd ../006_mesonh_283TO3003
sbatch run_mesonh.sh
&NAM_CONFIO LCDF4 = .TRUE.,
LLFIREAD = .FALSE.,
LLFIOUT = .FALSE. /
&NAM_DIAG_FILE YINIFILE(1) = 'A2830.1.MEGA2.001',
YINIFILEPGD(1) = 'BIOMAIDO_CH_8km.nestne',
YSUFFIX = 'diag' /
&NAM_DIAG LTPZH = .TRUE,
LMSLP = .TRUE.,
LTHW = .TRUE.,
CISO='PREVTK',
LCHEMDIAG = .TRUE. /
&NAM_DIAG_SURFn N2M=2,
LSURF_BUDGET=.TRUE.,
LCOEF=.TRUE. /
&NAM_DIAG_SURF_ATMn LFRAC=.TRUE.,
LDIAG_GRID=.TRUE. /
&NAM_CONF CCONF = "START", NVERB=1, NMODEL = 2,
CEXP = "A2830", CSEG = "MEGAD" , NHALO=3,
CSPLIT="BSPLITTING" /
&NAM_CONFIO LCDF4 = .TRUE., LLFIOUT = .FALSE., NIO_ABORT_LEVEL=1, NGEN_ABORT_LEVEL=1 /
&NAM_DYN XSEGLEN = 21600.,
LCORIO = T,
LNUMDIFU = F,
XALKTOP = 0.005,
XALZBOT = 15000. /
&NAM_SERIESn NFREQSERIES = 30 /
&NAM_LUNITn CINIFILE = "ECMWF_D1_20190328_06",
CINIFILEPGD="BIOMAIDO_CH_8km.nestne",
CCPLFILE(1) = "ECMWF_D1_20190328_12" /
&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 5, XWAY(2) = 2. /
&NAM_DYNn XTSTEP = 30.,
CPRESOPT = "ZRESI",
NITR = 12,
LHORELAX_UVWTH = F, LHORELAX_RV = F, LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01 /
&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",NWENO_ORDER=5,CTEMP_SCHEME="RK53",
CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
&NAM_PARAMn CCLOUD = "LIMA", CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "NONE", CSCONV="EDKF" /
&NAM_PARAM_RADn XDTRAD = 600., XDTRAD_CLONLY = 600., NRAD_COLNBR = 400 /
&NAM_PARAM_KAFRn XDTCONV = 360., NICE = 1, LCHTRANS = F,
LREFRESH_ALL = T, LDOWN = T /
&NAM_PARAM_MFSHALLn CMF_UPDRAFT = "EDKF", CMF_CLOUD = "DIRE" /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSIG_CONV = F, LTURB_FLX=F /
&NAM_CH_MNHCn LUSECHEM = .TRUE.,
LCH_CONV_LINOX = .FALSE.,
LCH_INIT_FIELD = .FALSE.,
LCH_CONV_SCAV = .TRUE.,
CCHEM_INPUT_FILE = "ReLACS_MACCity_aero.nam",
NCH_SUBSTEPS = 2,
LCH_TUV_ONLINE = .FALSE.,
CCH_TUV_LOOKUP = "PHOTO.TUV5",
CCH_TUV_CLOUDS = "CHAN",
XCH_TUV_TUPDATE = 300.,
XCH_TUV_DOBNEW = 250.,
CCH_VEC_METHOD = 'MAX'
NCH_VEC_LENGTH = 50 /
&NAM_CH_SOLVERn CSOLVER = "EXQSSA",
XSLOW = 90.,
XFAST = 0.1,
XDTMIN = 0.1,
XDTMAX = 45.,
XDTFIRST = 10. /
&NAM_BACKUP XBAK_TIME(1,1)=120. ,
XBAK_TIME(1,2)=7200. ,
XBAK_TIME(1,3)=10800. ,
XBAK_TIME(1,4)=14400. ,
XBAK_TIME(1,5)=18000. ,
XBAK_TIME(1,6)=21600. ,
XBAK_TIME(1,7)=25200. ,
XBAK_TIME(1,8)=28800. ,
XBAK_TIME(1,9)=32400. ,
XBAK_TIME(1,10)=36000. ,
XBAK_TIME(1,11)=39600. ,
XBAK_TIME(1,12)=43200. /
&NAM_PARAM_LIMA LACTI = T,
HINI_CCN='AER',
HTYPE_CCN(1) = 'C',
HTYPE_CCN(2) = 'C',
HTYPE_CCN(3) = 'C',
NMOD_CCN = 3, ! Nombre de modes de CCN
LACTIT=T,
LSEDC=T,
LSCAV=T, ! Scavenging
LAERO_MASS=F,
LCCN_HOM=F, ! Concentration en CCN homogène selon z
CCCN_MODES='CAMS_JPP', ! Choix des modes de CCN
XCCN_CONC(1) = 300., ! Concentration du mode 1 de CCN (/cm3)
XCCN_CONC(2) = 140., ! Concentration du mode 2 de CCN (/cm3)
XCCN_CONC(3) = 50., ! Concentration du mode 3 de CCN (/cm3)
XCCN_CONC(4) = 0.01, ! Concentration du mode 4 de CCN (/cm3)
LNUCL = T,
LSEDI = T,
LHHONI = F,
NMOD_IFN = 2, ! Nombre de modes d'IFN
CPRISTINE_ICE_LIMA = 'PLAT',
CHEVRIMED_ICE_LIMA = 'GRAU',
NIND_SPECIE = 1,
LMEYERS = F,
NMOD_IMM = 1, ! Nombre de modes de CCN agissant par IMM
LIFN_HOM = T, ! Concentration en IFN homogène selon z
CIFN_SPECIES = '', !
CINT_MIXING = 'MACC_JPP', ! DM1, DM2, BC, O, MACC, default (Phillips 2008)
XIFN_CONC(1) = 10., ! Concentration du mode 1 d'IFN (/L)
XIFN_CONC(2) = 10. / ! Concentration du mode 2 d'IFN (/L)
&NAM_CH_ORILAM LORILAM = .TRUE.,
LVARSIGI = .FALSE.,
LVARSIGJ = .FALSE.,
LRGFIX = .FALSE.,
LSEDIMAERO = .TRUE.,
LDEPOS_AER(1) = .TRUE.,
LDEPOS_AER(2) = .TRUE.,
XN0IMIN = 5.,
XN0JMIN = 0.001,
XINIRADIUSI= 0.015,
XINIRADIUSJ= 0.055,
XINISIGI = 1.400,
XINISIGJ = 1.600,
CNUCLEATION = "MAATTANEN_NEUTRAL",
LCONDENSATION = .TRUE.,
LCOAGULATION = .TRUE.,
LMODE_MERGING = .TRUE.,
XRADIUS_NUCL = 15E-3,
XSIGMA_NUCL = 1.2 /
&NAM_DIAG_SURFn N2M = 0, LSURF_BUDGET = F /
&NAM_ISBAn CSCOND="NP89", CC1DRY="DEF",
CSOILFRZ="DEF", CDIFSFCOND="DEF",
CSNOWRES="DEF" /
&NAM_SSOn CROUGH="Z01D" /
&NAM_DIAG_ISBAn /
&NAM_DUST LDUST = .FALSE.,
LRGFIX_DST = .FALSE.,
LVARSIG = .FALSE.,
LSEDIMDUST = .TRUE.,
LDEPOS_DST(1) = .TRUE.,
LDEPOS_DST(2) = .TRUE.,
LDEPOS_DST(3) = .TRUE.,
NMODE_DST = 3 /
&NAM_SALT LSALT = .TRUE.,
LRGFIX_SLT = .FALSE.,
LVARSIG_SLT = .FALSE.,
LSEDIMSALT = .TRUE.,
LDEPOS_SLT(1) = .TRUE.,
LDEPOS_SLT(2) = .TRUE.,
CRGUNITS = 'NUMB',
NMODE_SLT = 8 /
&NAM_SURF_SLT CEMISPARAM_SLT='OvB21b' /
&NAM_SURF_DST CEMISPARAM_DST='AMMA', CVERMOD='CMDVER' /
&NAM_SEAFLUXn CSEA_ALB="TA96", LPRECIP = T /
&NAM_CH_SURFn LCH_SURF_EMIS=T, LCH_DMSEMIS=T /
&NAM_CH_CONTROLn CCHEM_SURF_FILE= "ReLACS_MACCity_aero.nam" /
&NAM_CH_ISBAn LCH_BIO_FLUX=T, CCH_DRY_DEP="WES89", CPARAMBVOC="MEGAN",
XDAILYPAR=300., XDAILYTEMP=290. /
&NAM_CH_SEAFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_WATFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_TEBn CCH_DRY_DEP="WES89" /
&NAM_CHS_ORILAM LCH_AERO_FLUX = .TRUE.,
LCO2PM=.FALSE.,
XEMISRADIUSI = 0.015,
XEMISRADIUSJ = 0.055,
XEMISSIGI = 1.40,
XEMISSIGJ = 1.60,
CRGUNIT = "NUMB" /
&NAM_BLANK /
&NAM_NEBn LSIGMAS=T, LSUBG_COND=F /
&NAM_PARAM_ICEn CSUBG_AUCV_RC='NONE'/
&NAM_LUNITn CINIFILE = "ECMWF_D2_20190328_06",
CINIFILEPGD="BIOMAIDO_CH_2km.nestne" /
&NAM_DYNn CPRESOPT = "ZRESI",
NITR = 12,
LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083,LHORELAX_RH=F /
&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",NWENO_ORDER=5,CTEMP_SCHEME="RK53",
CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
&NAM_PARAMn CCLOUD = "LIMA", CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "NONE", CSCONV= "EDKF" /
&NAM_PARAM_MFSHALLn CMF_UPDRAFT='EDKF',CMF_CLOUD='DIRE', LMIXUV=T /
&NAM_PARAM_RADn XDTRAD = 600., XDTRAD_CLONLY = 600., NRAD_COLNBR = 400 /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LTURB_FLX=F /
&NAM_CH_MNHCn LUSECHEM = .TRUE.,
LCH_CONV_LINOX = .FALSE.,
LCH_INIT_FIELD = .FALSE.,
LCH_CONV_SCAV = .TRUE.,
CCHEM_INPUT_FILE = "ReLACS_MACCity_aero.nam",
NCH_SUBSTEPS = 2,
LCH_TUV_ONLINE = .FALSE.,
CCH_TUV_LOOKUP = "PHOTO.TUV5",
CCH_TUV_CLOUDS = "CHAN",
XCH_TUV_TUPDATE = 300.,
XCH_TUV_DOBNEW = 250.,
CCH_VEC_METHOD = 'MAX'
NCH_VEC_LENGTH = 50 /
&NAM_CH_SOLVERn CSOLVER = "EXQSSA",
XSLOW = 90.,
XFAST = 0.1,
XDTMIN = 0.1,
XDTMAX = 45.,
XDTFIRST = 10. /
&NAM_DIAG_SURFn N2M = 0, LSURF_BUDGET = F /
&NAM_ISBAn CSCOND="NP89",
CC1DRY="DEF", CSOILFRZ="DEF", CDIFSFCOND="DEF", CSNOWRES="DEF" /
&NAM_SSOn CROUGH="Z01D" /
&NAM_DIAG_ISBAn /
&NAM_SEAFLUXn CSEA_ALB="TA96", LPRECIP = T /
&NAM_CH_SURFn LCH_SURF_EMIS=T /
&NAM_CH_CONTROLn CCHEM_SURF_FILE= "ReLACS_MACCity_aero.nam" /
&NAM_CH_ISBAn LCH_BIO_FLUX=T, CCH_DRY_DEP="WES89", CPARAMBVOC="MEGAN",
XDAILYPAR=150., XDAILYTEMP=290. /
&NAM_CH_SEAFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_WATFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_TEBn CCH_DRY_DEP="WES89" /
&NAM_NEBn LSIGMAS=T , LSUBG_COND = F /
&NAM_PARAM_ICEn CSUBG_AUCV_RC='NONE'/
\ No newline at end of file