Quentin 14/08/2024: adapt test cases with modif of BACKUP management +...

Quentin 14/08/2024: adapt test cases with modif of BACKUP management + cleaning (profiles, LFI old files) + add new plots for recent test cases

MY_RUN/INTEGRATION_CASES/HPC/16JAN/PYTHON/plot_16JAN.py

@@ -17,7 +17,7 @@ os.system('rm -f tempgraph*')
 #  User's parameter / Namelist
 #
 #
-path="../RUN/NETCDF/"
+path="../RUN/"
 LnameFiles = ['16JA1.1.WENO5.004dg.nc', '16JA2.2.WENO5.002dg.nc']
 
 Dvar_input = {

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_1

@@ -15,7 +15,10 @@
   LCLD_COV=T, 
   NCAPE=1,
   LRADAR=T,
-  LDIAG(:)=.FALSE. /
+  LDIAG(:)=.FALSE.,
+  LISOAL=T, XISOAL(1)=9000,
+  LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
+
 &NAM_DIAG_FILE      YINIFILE(1) = "16JA1.1.WENO5.004" ,
                     YINIFILEPGD(1) = "16JAN98_36km.neste1",
                     YSUFFIX='dg' /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_2

 &NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
-&NAM_DIAG CISO='', LVAR_RS=F, LVAR_LS=F, LVAR_PR=T, LTOTAL_PR=T /
+&NAM_DIAG CISO='', LVAR_RS=F, LVAR_LS=F, LVAR_PR=T, LTOTAL_PR=T,
+  LISOAL=T, XISOAL(1)=9000,
+  LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
 &NAM_BLANK /
 &NAM_DIAG_FILE  YINIFILE(1) = "16JA2.2.WENO5.002",
                 YINIFILEPGD(1) = "16JAN98_9km.neste1",

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/PYTHON/plot_ARMCU_LES_DEAR.py

@@ -17,7 +17,7 @@ os.system('rm -f tempgraph*')
 #
 #  User's parameter / Namelist
 #
-path="../DEAR/NETCDF/"
+path="../DEAR/"
 
 LnameFiles = ['ARM__.1.CEN4T.000.nc' ]
 

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/PYTHON/plot_ARMCU_LES_HM21.py

@@ -17,7 +17,7 @@ os.system('rm -f tempgraph*')
 #
 #  User's parameter / Namelist
 #
-path="../HM21/NETCDF/"
+path="../HM21/"
 
 LnameFiles = ['ARM__.1.CEN4T.000.nc' ]
 

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/EXSEG1.nam

@@ -30,7 +30,7 @@
 &NAM_DYN XSEGLEN = 3600, XASSELIN = 0.2, LCORIO = T, LNUMDIFU = T,
          XALKTOP = 0.003, XALZBOT = 15000. /
 &NAM_NESTING NDAD(2) = 1,NDTRATIO(2) =5 ,XWAY(2) = 2.  /
-&NAM_BACKUP XBAK_TIME(1,1)=0., XBAK_TIME(1,2)=60.,XBAK_TIME(1,3)=600.,XBAK_TIME(1,4)=900.,XBAK_TIME(1,5)=1800.,XBAK_TIME(1,6)=2700.,XBAK_TIME(1,7)=3600. /
+&NAM_BACKUP XBAK_TIME(1,1)=10., XBAK_TIME(1,2)=60.,XBAK_TIME(1,3)=600.,XBAK_TIME(1,4)=900.,XBAK_TIME(1,5)=1800.,XBAK_TIME(1,6)=2700.,XBAK_TIME(1,7)=3600. /
 &NAM_PASPOL
 LPASPOL   =  T  ,
 NRELEASE     =   2        ,

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/001_prep_pgd1/run_pgd.sh

@@ -14,8 +14,6 @@ set -e
 # Nom de la machine
 hostname
 
-
-#. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PIERRE
 . ~/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/* .

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/002_prep_pgd2/run_pgd.sh

@@ -14,8 +14,6 @@ set -e
 # Nom de la machine
 hostname
 
-
-#. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PIERRE
 . ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
 
 ln -sf ../001_prep_pgd1/BIOM* .

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/003_prep_nest/run_prep_nest.sh

@@ -17,13 +17,11 @@ 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
-#. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PIERRE
 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/run_prep_real_case.sh

@@ -11,15 +11,11 @@ ulimit -c 0
 ulimit -s unlimited
 # Arrete du job des la premiere erreur
 set -e
-# Nom de la machine
-hostname
-
 
 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-DEBUG
-#. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PIERRE
+. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
 export MPIRUN="Mpirun -np 1"
 
 set -x

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/005_prep_real_case_d2/run_prep_real_case.sh

@@ -11,13 +11,10 @@ ulimit -c 0
 ulimit -s unlimited
 # Arrete du job des la premiere erreur
 set -e
-# Nom de la machine
-hostname
-
 
 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-RELACS
+. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
 
 export MPIRUN="Mpirun -np 1"
 

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/run_mesonh.sh

 #!/bin/sh
-#SBATCH -J run_megan1
+#SBATCH -J run_biomaido
 #SBATCH -N 1            # nodes number (=NBP)
 #SBATCH -n 64           # CPUs number (on all nodes) (=NBP*TPN)
 #SBATCH -o run_megan.eo%j   #
@@ -11,11 +11,8 @@ 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-RELACS
+. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
 
 ln -sf ../003_prep_nest/BIOM* .
 ln -sf ../004_prep_real_case_d1/ECMWF* .

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/PYTHON/plot_BIOMAIDO_DMS.py

+#!/usr/bin/env python3
+"""
+@author: Quentin Rodier
+Creation : 07/01/2021
+
+Last modifications
+"""
+import matplotlib as mpl
+mpl.use('Agg')
+import cartopy.crs as ccrs
+try:
+	import MNHPy
+	from MNHPy.read_MNHfile import read_netcdf
+	from MNHPy.Panel_Plot import PanelPlot
+except:
+	from read_MNHfile import read_netcdf
+	from Panel_Plot import PanelPlot
+import os
+
+Lpath=["../006_mesonh_283TO3003/"]
+LnameFiles=['A2830.1.MEGAD.006.nc']
+
+variables=['CO','O3','NO2','NO','BIOP','HCHO','BIO','latitude','longitude',
+           'ETH', 'CARBO','PHO','SO2','SULF','DMS','CH4','XO2','PAN',
+           'DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R',
+           'SLTM01','SLTM02','SLTM03','SLTM04','SLTM05','SLTM06','SLTM07','SLTM08']
+Dvar_input = {'f1': variables}
+
+#  Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=Lpath, removeHALO=True, get_data_only=True)
+
+os.system('rm -f tempgraph*')
+
+f=Dvar['f1']
+nb_subplot=4
+
+lon = [Dvar['f1']['longitude']]*nb_subplot
+lat = [Dvar['f1']['latitude']]*nb_subplot
+Lxlab = ['longitude']*nb_subplot
+Lylab = ['latitude']*nb_subplot
+Lprojection = [ccrs.PlateCarree()]*nb_subplot
+Llvl = [0]*nb_subplot
+
+################################################################
+#########          PANEL 1
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['BIOP'],f['BIO'],f['CO'],f['HCHO']]
+Ltitle = ['BIOP','BIO','CO','HCHO']
+Lminval=[0,      0,    30,    0.001 ]
+Lmaxval=[15.0,  1000, 110,    1  ]
+Lstep=[0.25,      20,   2,    0.04 ]
+Lstepticks=[1,   100,  10,   0.1 ]
+Lfacconv=[1E12,1E12,1E9,1E9    ]
+Lcbarlabel = ['ppt','ppt','ppb','ppb']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20], 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(1,fig)
+
+################################################################
+#########          PANEL 2
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['O3'],f['NO2'],f['NO'],f['ETH']]
+Ltitle = ['O3','NO2','NO','ETH']
+Lminval=[15,      0,    30,    125 ]
+Lmaxval=[25.0,  1000, 1000,    175  ]
+Lstep=[0.1,      20,   20,    0.25 ]
+Lstepticks=[1,   100,  100,   5 ]
+Lfacconv=[1E9,1E12,1E12,1E12    ]
+Lcbarlabel = ['ppb','ppt','ppt','ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20], 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(2,fig)
+
+################################################################
+#########          PANEL 3
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['CARBO'],f['PHO'],f['SO2'],f['SULF']]
+Ltitle = ['CARBO','PHO','SO2','SULF']
+Lminval=[0,      0,    50,    0 ]
+Lmaxval=[400,  0.7, 1000,    1  ]
+Lstep=[5,      0.01,   10,    0.01 ]
+Lstepticks=[100,   0.1,  100,   0.1 ]
+Lfacconv=[1E12,1E12,1E12,1E12    ]
+Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(3,fig)
+
+
+################################################################
+#########          PANEL 4
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['DMS'],f['CH4'],f['XO2'],f['PAN']]
+Ltitle = ['DMS','CH4','XO2','PAN']
+Lminval=[0,      0,    0.1,    0 ]
+Lmaxval=[300,  0.1, 0.6,       500  ]
+Lstep=[2,      0.001,   0.01,  5 ]
+Lstepticks=[50,   0.01,  0.1,   100 ]
+Lfacconv=[1E12,1E12,1E12,1E12    ]
+Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(4,fig)
+
+################################################################
+#########          PANEL 5
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['DEAERM31C'],f['DEAERM32C'],f['DEAERM31R'],f['DEAERM32R']]
+Ltitle = ['DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R']
+Lminval=[0,      0,    0,    0 ]
+Lmaxval=[10,  40, 1,    1  ]
+Lstep=[0.005,     ]*4
+Lstepticks=[0.5  ]*4
+Lfacconv=[1E12,1E12,1E12,1E12    ]
+Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(5,fig)
+
+################################################################
+#########          PANEL 6
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['SLTM01'],f['SLTM02'],f['SLTM03'],f['SLTM04']]
+Ltitle = ['SLTM01','SLTM02','SLTM03','SLTM04']
+Lminval=[0]*4
+Lmaxval=[3]*4
+Lstep=[0.01      ]*4
+Lstepticks=[1  ]*4
+Lfacconv=[1E18,1E18,1E18,1E18    ]
+Lcbarlabel = ['1E-6 ppt','1E-6 ppt','1E-6 ppt','1E-6 ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(6,fig)
+
+################################################################
+#########          PANEL 7
+###############################################################
+Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
+
+Lplot=[f['SLTM05'],f['SLTM06'],f['SLTM07'],f['SLTM08']]
+Ltitle = ['SLTM05','SLTM06','SLTM07','SLTM08']
+Lminval=[0 ]*4
+Lmaxval=[1  ]*4
+Lstep=[0.01  ]*4
+Lstepticks=[0.1  ]*4
+Lfacconv=[1E18,1E21,1E21,1E21    ]
+Lcbarlabel = ['1E-6 ppt','1E-9 ppt','1E-9 ppt','1E-9 ppt']*nb_subplot
+LaddWhite_cm=[True]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
+
+Panel.save_graph(7,fig)
+

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/PYTHON/run_python

@@ -11,12 +11,12 @@ 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-PIERRE
+. ~rodierq/DEV_57/MNH-git_open_source-lfs-0a57991/conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-0-ECRAD140-MPIAUTO-O2
 
 ln -sf ${SRC_MESONH}/src/LIB/Python/* .
 
 module purge
 module load python/3.7.6
 
-python3 plot_MEGAN_REUNION.py
-convert *.png MEGAN_REUNION.pdf
+python3 plot_BIOMAIDO_DMS.py
+convert *.png BIOMAIDO_DMS.pdf

MY_RUN/INTEGRATION_CASES/HPC/BLAZE_OROG/02_mesonh/run_mesonh

@@ -19,12 +19,15 @@ export MPIRUN="Mpirun -np 128"
 
 set -x
 set -e
+
+# Prepare the Fuel map (input of the MESONH step)
 #module purge
 #module load python/3.7.6
-
 #python3 fuelmap.py
 #module purge
 
+ln -sf ~rodierq/SAVE/GRIB_KTEST/Fuel* .
+
 . ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
 
 ln -sf ../01_prep_ideal_case/PGDFireTest.* .

MY_RUN/INTEGRATION_CASES/HPC/BLAZE_OROG/03_python/plot_BLAZE.py

@@ -43,10 +43,10 @@ Ltitle = ['u', 'v','w','théta',]
 Lcbarlabel = ['m/s','m/s','m/s','K']
 Lxlab = ['x (m)']*len(Lplot)
 Lylab = ['y (m)']*len(Lplot)
-Lminval = [-1.5, -3.5, -1.5, 286]
-Lmaxval = [3.5,3.5,1.5,300]
-Lstep = [0.1,0.1,0.1,0.25]
-Lstepticks = [0.5,0.5,0.5,1]
+Lminval = [-1.5, -3.5, -1.5, 299]
+Lmaxval = [3.5,3.5,1.5,300.5]
+Lstep = [0.1,0.1,0.1,0.05]
+Lstepticks = [0.5,0.5,0.5,0.25]
 Lfacconv = [1]*len(Lplot)
 Lcolormap = ['gist_rainbow_r','RdBu_r','RdBu_r','gist_rainbow_r']*len(Lplot)
 Llvl = [1]*len(Lplot)
@@ -66,8 +66,8 @@ Panel1.save_graph(1,fig1)
 Panel1 = PanelPlot(2,2, [25,15],'K=28', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.01)
 
 Llvl = [28]*len(Lplot)
-Lminval = [-2, -2., -1.5, 285]
-Lmaxval = [4,2.,4,288]
+Lminval = [-2, -2., -1.5, 299.5]
+Lmaxval = [4,2.,4,302.5]
 Lstep = [0.1,0.1,0.1,0.05]
 
 fig2 = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 

MY_RUN/INTEGRATION_CASES/HPC/CHARMEX/NO_MEGAN/plot_chem_06.ncl

@@ -8,7 +8,7 @@ begin
 ; open file and read in data
 ;=================================================;
 
-  a = addfile("../RUN/OUTFILES/CHARM.1.SEG02.003.nc", "r")
+  a = addfile("OUTFILES/CHARM.1.SEG02.003.nc", "r")
 ;=================================================;
 ; Get informations on variable sizes
 ; dims are dims-2 to remove non-physical values

MY_RUN/INTEGRATION_CASES/HPC/CHARMEX/NO_MEGAN/plot_chem_12.ncl

@@ -8,7 +8,7 @@ begin
 ; open file and read in data
 ;=================================================;
 
-  a = addfile("../RUN/OUTFILES/CHARM.1.SEG02.009.nc", "r")
+  a = addfile("OUTFILES/CHARM.1.SEG02.009.nc", "r")
 ;=================================================;
 ; Get informations on variable sizes
 ; dims are dims-2 to remove non-physical values