diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl b/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl
new file mode 100755
index 0000000000000000000000000000000000000000..12eeedc8e2098dd3b9af358058400b0dc719d700
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl
@@ -0,0 +1,2 @@
+#!/bin/bash
+rm -f ARM__.1.CEN4T.000* *.png
diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl b/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl
new file mode 100644
index 0000000000000000000000000000000000000000..77b79867f4133f56b2c3c5cc8c1743a285598eaa
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl
@@ -0,0 +1,512 @@
+;================================================;
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
+; ================================================;
+;================================================;
+begin
+;=================================================;
+; open file and read in data
+;=================================================;
+ fichier1 = addfile("ARM__.1.CEN4T.000.nc", "r")
+;==================================================;
+; Open the workstation
+;==================================================;
+ type = "png"
+ wks = gsn_open_wks(type,"visu_ARM")
+;=================================================;
+; lecture des différents champs
+;=================================================;
+kmax=100
+zhat = fichier1->ZHAT(:) ; ZHAT
+; Unstagger zhat (from grid 4 to 1)
+ nzh=new(kmax,double)
+ do k=0,kmax-2
+ nzh(k)=(zhat(k)+zhat(k+1))/2.
+ end do
+ nzh(kmax-1)=2*zhat(kmax-1)-zhat(kmax-2)
+TIME=new((/180/),double)
+TIME(0)=300
+do j=1,179
+ TIME(j)=TIME(j-1)+300
+end do
+
+group_LES_budgets = fichier1=>LES_budgets/Cartesian/Not_time_averaged/Not_normalized/cart
+
+vname="SBG_TKE"
+SBGTKE = group_LES_budgets->$vname$(:,:)
+SBGTKE@long_name="SBG_TKE"
+SBGTKE!1="nzh"
+SBGTKE!0="TIME"
+SBG_TKE=transpose(SBGTKE)
+
+vname="MEAN_TH"
+MEANTH = group_LES_budgets->$vname$(:,:) ; MEAN_TH
+MEANTH@long_name="MEAN_TH"
+MEANTH@units="K"
+MEANTH!1="nzh"
+MEANTH!0="TIME"
+MEAN_TH=transpose(MEANTH)
+
+vname="MEAN_U"
+MEANU = group_LES_budgets->$vname$(:,:)
+MEANU!1="nzh"
+MEANU!0="TIME"
+MEAN_U=transpose(MEANU)
+MEAN_U@long_name="MEAN_U"
+
+vname="MEAN_V"
+MEANV = group_LES_budgets->$vname$(:,:)
+MEANV!1="nzh"
+MEANV!0="TIME"
+MEAN_V=transpose(MEANV)
+MEAN_V@long_name="MEAN_V"
+
+vname="MEAN_RC"
+MEANRC = group_LES_budgets->$vname$(:,:)
+MEANRC!1="nzh"
+MEANRC!0="TIME"
+MEAN_RC=transpose(MEANRC)
+MEAN_RC@long_name="MEAN_RC"
+
+vname="MEAN_RR"
+MEANRR = group_LES_budgets->$vname$(:,:)
+MEANRR!1="nzh"
+MEANRR!0="TIME"
+MEAN_RR=transpose(MEANRR)
+MEAN_RR@long_name="MEAN_RR"
+
+vname="SBG_WTHL"
+SBGWTHL = group_LES_budgets->$vname$(:,:)
+SBGWTHL!1="nzh"
+SBGWTHL!0="TIME"
+SBG_WTHL=transpose(SBGWTHL)
+SBG_WTHL@long_name="SBG_WTHL"
+
+vname="SBG_WRT"
+SBGWRT = group_LES_budgets->$vname$(:,:)
+SBGWRT!1="nzh"
+SBGWRT!0="TIME"
+SBG_WRT=transpose(SBGWRT)
+SBG_WRT@long_name="SBG_WRT"
+
+
+vname="THLUP_MF"
+ THLUPMF= group_LES_budgets->$vname$(:,:)
+THLUPMF!1="nzh"
+THLUPMF!0="TIME"
+THLUP_MF=transpose(THLUPMF)
+THLUP_MF@long_name="THLUP_MF"
+
+vname="RTUP_MF"
+RTUPMF = group_LES_budgets->$vname$(:,:)
+RTUPMF!1="nzh"
+RTUPMF!0="TIME"
+RTUP_MF=transpose(RTUPMF)
+RTUP_MF@long_name="RTUP_MF"
+
+vname="RVUP_MF"
+RVUPMF = group_LES_budgets->$vname$(:,:)
+RVUPMF!1="nzh"
+RVUPMF!0="TIME"
+RVUP_MF=transpose(RVUPMF)
+RVUP_MF@long_name="RVUP_MF"
+
+vname="RCUP_MF"
+RCUPMF = group_LES_budgets->$vname$(:,:)
+RCUPMF!1="nzh"
+RCUPMF!0="TIME"
+RCUP_MF=transpose(RCUPMF)
+RCUP_MF@long_name="RCUP_MF"
+
+vname="RIUP_MF"
+RIUPMF = group_LES_budgets->$vname$(:,:)
+RIUPMF!1="nzh"
+RIUPMF!0="TIME"
+RIUP_MF=transpose(RIUPMF)
+RIUP_MF@long_name="RIUP_MF"
+
+vname="WUP_MF"
+WUPMF = group_LES_budgets->$vname$(:,:)
+WUPMF!1="nzh"
+WUPMF!0="TIME"
+WUP_MF=transpose(WUPMF)
+WUP_MF@long_name="WUP_MF"
+
+vname="MAFLX_MF"
+MAFLXMF = group_LES_budgets->$vname$(:,:)
+MAFLXMF!1="nzh"
+MAFLXMF!0="TIME"
+MAFLX_MF=transpose(MAFLXMF)
+MAFLX_MF@long_name="MAFLX_MF"
+
+vname="DETR_MF"
+DETRMF = group_LES_budgets->$vname$(:,:)
+DETRMF!1="nzh"
+DETRMF!0="TIME"
+DETR_MF=transpose(DETRMF)
+DETR_MF@long_name="DETR_MF"
+
+vname="ENTR_MF"
+ENTRMF = group_LES_budgets->$vname$(:,:)
+ENTRMF!1="nzh"
+ENTRMF!0="TIME"
+ENTR_MF=transpose(ENTRMF)
+ENTR_MF@long_name="ENTR_MF"
+
+vname="FRCUP_MF"
+FRCUPMF = group_LES_budgets->$vname$(:,:)
+FRCUPMF!1="nzh"
+FRCUPMF!0="TIME"
+FRCUP_MF=transpose(FRCUPMF)
+FRCUP_MF@long_name="FRCUP_MF"
+
+vname="THVUP_MF"
+THVUPMF = group_LES_budgets->$vname$(:,:)
+THVUPMF!1="nzh"
+THVUPMF!0="TIME"
+THVUP_MF=transpose(THVUPMF)
+THVUP_MF@long_name="THVUP_MF"
+
+vname="WTHL_MF"
+WTHLMF = group_LES_budgets->$vname$(:,:)
+WTHLMF!1="nzh"
+WTHLMF!0="TIME"
+WTHL_MF=transpose(WTHLMF)
+WTHL_MF@long_name="WTHL_MF"
+
+vname="WRT_MF"
+WRTMF = group_LES_budgets->$vname$(:,:)
+WRTMF!1="nzh"
+WRTMF!0="TIME"
+WRT_MF=transpose(WRTMF)
+WRT_MF@long_name="WRT_MF"
+
+
+vname="WTHV_MF"
+WTHVMF = group_LES_budgets->$vname$(:,:)
+WTHVMF!1="nzh"
+WTHVMF!0="TIME"
+WTHV_MF=transpose(WTHVMF)
+WTHV_MF@long_name="WTHV_MF"
+
+
+vname="WU_MF"
+WUMF = group_LES_budgets->$vname$(:,:)
+WUMF!1="nzh"
+WUMF!0="TIME"
+WU_MF=transpose(WUMF)
+WU_MF@long_name="WU_MF"
+
+
+vname="WV_MF"
+WVMF = group_LES_budgets->$vname$(:,:)
+WVMF!1="nzh"
+WVMF!0="TIME"
+WV_MF=transpose(WVMF)
+WV_MF@long_name="WV_MF"
+
+;=================================================;
+; Altitude des niveaux modèles
+;=================================================;
+; Unstagger zhat (from grid 4 to 1)
+ nzh=new(kmax,double)
+ do k=0,kmax-2
+ nzh(k)=(zhat(k)+zhat(k+1))/2.
+ end do
+ nzh(kmax-1)=2*zhat(kmax-1)-zhat(kmax-2)
+;=================================================;
+; Set some other basic resources
+;=================================================;
+ resmap = True
+ resmap@gsnFrame = False
+ resmap@gsnDraw = False
+ resmap@gsnMaximize = True
+ resmap@gsnPaperOrientation = "portrait"
+ resmap@gsnSpreadColors = True ; use full range of colormap
+ resmap@tiYAxisString =" "
+ resmap@cnFillOn = True ; turn on color fill
+ resmap@cnLinesOn = False ; turn off contour lines
+ ;resmap@tmXBLabelStride = 2 ; to reduce the number of labels on xaxis
+ ;resmap@lbLabelStride = 2. ; to reduce the number of labels on labelbar
+ resmap@sfYArray = nzh(0:76)
+ resmap@sfXArray = TIME
+ resmap@tmXBTickSpacingF = 10800.
+ resmap@tmXBMode = "Explicit"
+ resmap@tmXBValues = (/10800,21600,32400,43200,54000/)
+ resmap@tmXBLabels =(/3,6,9,12,15/)
+resmap@tiXAxisPosition="Left"
+resmap@tiXAxisFontHeightF=0.015
+;=================================================;
+; TRACE
+;=================================================;
+ gsn_define_colormap(wks,"rainbow") ; Choose colormap
+ colors = gsn_retrieve_colormap(wks) ; retrieve color map for editing
+ colors(2,:) = (/ 1, 1, 1 /) ; replace the first color with white color (les deux premières sont background /foreground donc c'est bien la 2 qu'il faut changer)
+ gsn_define_colormap(wks,colors) ; redefine colormap to workstation, color map now includes a gray
+
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.25
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=2.5
+
+plot = gsn_csm_contour(wks,SBG_TKE(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=298
+opts@cnMaxLevelValF=316
+
+plot = gsn_csm_contour(wks,MEAN_TH(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.5
+opts@cnMinLevelValF=3
+opts@cnMaxLevelValF=11
+
+plot = gsn_csm_contour(wks,MEAN_U(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.2
+opts@cnMinLevelValF=0.
+opts@cnMaxLevelValF=3.2
+
+plot = gsn_csm_contour(wks,MEAN_V(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00001
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=0.0002
+
+plot = gsn_csm_contour(wks,MEAN_RC(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.000001
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=0.00002
+
+plot = gsn_csm_contour(wks,MEAN_RR(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.02
+opts@cnMinLevelValF=-0.3
+opts@cnMaxLevelValF=0.12
+
+plot = gsn_csm_contour(wks,SBG_WTHL(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00005
+opts@cnMinLevelValF=-0.0001
+opts@cnMaxLevelValF=0.0004
+
+plot = gsn_csm_contour(wks,SBG_WRT(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=301
+opts@cnMaxLevelValF=318
+
+plot = gsn_csm_contour(wks,THLUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.001
+opts@cnMinLevelValF=0.0035
+opts@cnMaxLevelValF=0.017
+plot = gsn_csm_contour(wks,RTUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.001
+opts@cnMinLevelValF=0.001
+opts@cnMaxLevelValF=0.017
+
+plot = gsn_csm_contour(wks,RVUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00025
+opts@cnMinLevelValF=0.00025
+opts@cnMaxLevelValF=0.00425
+
+plot = gsn_csm_contour(wks,RCUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "AutomaticLevels"
+plot = gsn_csm_contour(wks,RIUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.25
+opts@cnMinLevelValF=0.250
+opts@cnMaxLevelValF=5.5
+
+plot = gsn_csm_contour(wks,WUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.02
+opts@cnMinLevelValF=0.02
+opts@cnMaxLevelValF=0.32
+
+plot = gsn_csm_contour(wks,MAFLX_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.025
+opts@cnMinLevelValF=0.025
+opts@cnMaxLevelValF=0.45
+
+plot = gsn_csm_contour(wks,DETR_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00025
+opts@cnMinLevelValF=0.00025
+opts@cnMaxLevelValF=0.005
+
+plot = gsn_csm_contour(wks,ENTR_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.01
+opts@cnMinLevelValF=0.01
+opts@cnMaxLevelValF=0.16
+
+plot = gsn_csm_contour(wks,FRCUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=303
+opts@cnMaxLevelValF=319
+
+plot = gsn_csm_contour(wks,THVUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.075
+opts@cnMaxLevelValF=0.055
+
+plot = gsn_csm_contour(wks,WTHL_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00001
+opts@cnMinLevelValF=0.00001
+opts@cnMaxLevelValF=0.00019
+
+plot = gsn_csm_contour(wks,WRT_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.016
+opts@cnMaxLevelValF=0.075
+plot = gsn_csm_contour(wks,WTHV_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.120
+opts@cnMaxLevelValF=-0.005
+
+plot = gsn_csm_contour(wks,WU_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.004
+opts@cnMinLevelValF=-0.024
+opts@cnMaxLevelValF=0.048
+plot = gsn_csm_contour(wks,WV_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+end
diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl b/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl
new file mode 100755
index 0000000000000000000000000000000000000000..7fa36570be699763a33bdad50638a95b35f01eb7
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl
@@ -0,0 +1,37 @@
+#!/bin/bash
+#MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+#MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+#MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+#MNH_LIC for details. version 1.
+set -x
+
+ln -sf ../002_mesonh/ARM__.1.CEN4T.000.* .
+
+rm -f visu_ARM.*.png
+
+ncl plot_arm.ncl
+${POSTRUN} display visu_ARM.000001.png
+${POSTRUN} display visu_ARM.000002.png
+${POSTRUN} display visu_ARM.000003.png
+${POSTRUN} display visu_ARM.000004.png
+${POSTRUN} display visu_ARM.000005.png
+${POSTRUN} display visu_ARM.000006.png
+${POSTRUN} display visu_ARM.000007.png
+${POSTRUN} display visu_ARM.000008.png
+${POSTRUN} display visu_ARM.000009.png
+${POSTRUN} display visu_ARM.000010.png
+${POSTRUN} display visu_ARM.000011.png
+${POSTRUN} display visu_ARM.000012.png
+${POSTRUN} display visu_ARM.000013.png
+${POSTRUN} display visu_ARM.000014.png
+${POSTRUN} display visu_ARM.000015.png
+${POSTRUN} display visu_ARM.000016.png
+${POSTRUN} display visu_ARM.000017.png
+${POSTRUN} display visu_ARM.000018.png
+${POSTRUN} display visu_ARM.000019.png
+${POSTRUN} display visu_ARM.000020.png
+${POSTRUN} display visu_ARM.000021.png
+${POSTRUN} display visu_ARM.000022.png
+${POSTRUN} display visu_ARM.000023.png
+${POSTRUN} display visu_ARM.000024.png
+exit 0
diff --git a/MY_RUN/KTEST/005_ARM/Makefile b/MY_RUN/KTEST/005_ARM/Makefile
index ea3b47d02c482ebfd287cf86e058ea1f156a6b8c..e9802cb77dfeb5f17cc60876cc4024efe8a6c41f 100644
--- a/MY_RUN/KTEST/005_ARM/Makefile
+++ b/MY_RUN/KTEST/005_ARM/Makefile
@@ -1,9 +1,15 @@
all:
cd 001_prep_ideal && run_prep_ideal_case_xyz
cd 002_mesonh && run_mesonh_xyz
+ ifneq "$(MNH_PYTHON)" "NO"
cd 003_python && run_python
+ endif
+ cd 004_ncl && run_ncl
clean:
cd 001_prep_ideal && clean_prep_ideal_case_xyz
cd 002_mesonh && clean_mesonh_xyz
+ ifneq "$(MNH_PYTHON)" "NO"
cd 003_python && clean_python
+ endif
+ cd 004_ncl && clean_ncl
diff --git a/MY_RUN/KTEST/014_LIMA/003_ncl/clean_ncl b/MY_RUN/KTEST/014_LIMA/003_ncl/clean_ncl
new file mode 100755
index 0000000000000000000000000000000000000000..ef83c16a5de17b683dd873050ef83dedfa0a005d
--- /dev/null
+++ b/MY_RUN/KTEST/014_LIMA/003_ncl/clean_ncl
@@ -0,0 +1,2 @@
+#!/bin/bash
+rm -f *.nc *.png
diff --git a/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl b/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
new file mode 100644
index 0000000000000000000000000000000000000000..753a1d2ec97a0c960d674b36dea14f611165c718
--- /dev/null
+++ b/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
@@ -0,0 +1,126 @@
+;================================================;
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
+; ================================================;
+;================================================;
+begin
+;=================================================;
+; open file and read in data
+;=================================================;
+ fichier1 = addfile("XPREF.1.SEG01.000.nc", "r")
+;==================================================;
+; Open the workstation
+;==================================================;
+ type = "png"
+ wks = gsn_open_wks(type,"visu_LIMA")
+;=================================================;
+; lecture des différents champs
+;=================================================;
+
+group_Budgets = fichier1=>Budgets
+group_RI = group_Budgets=>RI
+group_RS = group_Budgets=>RS
+group_RG = group_Budgets=>RG
+group_CICE = group_Budgets=>CICE
+group_CIFNFREE01 = group_Budgets=>CIFNFREE01
+group_CIFNNUCL01 = group_Budgets=>CIFNNUCL01
+
+
+vname = "AVEF"
+RI= group_RI->$vname$(0,:,0,:)
+RI@long_name="ice water content"
+RI@units="g/kg"
+
+NI= group_CICE->$vname$(0,:,0,:)
+NI@long_name="ice concentration"
+NI@units="/kg"
+
+RS= group_RS->$vname$(0,:,0,:)
+RS@long_name="snow water content"
+RS@units="g/kg"
+
+RG= group_RG->$vname$(0,:,0,:)
+RG@long_name="graupel water content"
+RG@units="g/kg"
+
+N_IFN_FREE= group_CIFNFREE01->$vname$(0,:,0,:)
+N_IFN_FREE@long_name="concentration of free IFN"
+N_IFN_FREE@units="/kg"
+
+N_IFN_NUCL= group_CIFNNUCL01->$vname$(0,:,0,:)
+N_IFN_NUCL@long_name="concentration of nucleated IFN"
+N_IFN_NUCL@units="/kg"
+
+
+zhat = fichier1->ZHAT(1:50)
+xhat = fichier1->XHAT(1:180)
+zs = fichier1->ZS(1:180) ; ZS
+
+xconf=conform(RI,xhat(0:179),1)
+
+;=================================================;
+; Altitude des niveaux modèles
+;=================================================;
+; Unstagger zhat (from grid 4 to 1)
+ nzh=new(50,double)
+ do k=0,48
+ nzh(k)=(zhat(k)+zhat(k+1))/2.
+ end do
+ nzh(49)=2*zhat(49)-zhat(48)
+
+; Create z (altitude des niveaux modèle)
+ z=new((/50,180/),double)
+ zcoef=new(180,double)
+ zcoef=1.-zs/nzh(49)
+
+ do i=0,179
+ z(:,i) = nzh*zcoef(i)+zs(i)
+ end do
+;=================================================;
+; Set some other basic resources
+;=================================================;
+ resmap = True
+ resmap@gsnFrame = False
+ resmap@gsnDraw = False
+ resmap@gsnMaximize = True
+ resmap@gsnPaperOrientation = "portrait"
+ resmap@gsnSpreadColors = True ; use full range of colormap
+ resmap@tiYAxisString =" "
+ resmap@cnFillOn = True ; turn on color fill
+ resmap@cnLinesOn = False ; turn off contour lines
+
+ resmap@tiXAxisPosition="Left"
+ resmap@tiXAxisFontHeightF=0.015
+; resmap@gsnSpreadColorStart = 0 ; force la première couleur en blanc(= couleur 0 de la palette)
+ resmap@sfYArray = z ; 2D
+ resmap@sfXArray = xconf ; 2D
+ resmap@trGridType = "TriangularMesh"
+;=================================================;
+; TRACE
+;=================================================;
+ gsn_define_colormap(wks,"rainbow") ; Choose colormap
+ res=resmap
+ res@cnLevelSelectionMode = "AutomaticLevels"
+ plot = gsn_csm_contour(wks,RI(:,:),res)
+ draw(plot)
+ frame(wks)
+
+ plot2 = gsn_csm_contour(wks,NI(:,:),res)
+ draw(plot2)
+ frame(wks)
+
+ plot3 = gsn_csm_contour(wks,RS(:,:),res)
+ draw(plot3)
+ frame(wks)
+
+ plot4 = gsn_csm_contour(wks,RG(:,:),res)
+ draw(plot4)
+ frame(wks)
+ plot5 = gsn_csm_contour(wks,N_IFN_FREE(:,:),res)
+ draw(plot5)
+ frame(wks)
+ plot6 = gsn_csm_contour(wks,N_IFN_NUCL(:,:),res)
+ draw(plot6)
+ frame(wks)
+
+end
diff --git a/MY_RUN/KTEST/014_LIMA/003_ncl/run_ncl b/MY_RUN/KTEST/014_LIMA/003_ncl/run_ncl
new file mode 100755
index 0000000000000000000000000000000000000000..2444be340eda2cb7a5a846cb44cd71694626e1c0
--- /dev/null
+++ b/MY_RUN/KTEST/014_LIMA/003_ncl/run_ncl
@@ -0,0 +1,20 @@
+#!/bin/bash
+#MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+#MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+#MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+#MNH_LIC for details. version 1.
+set -x
+
+ln -sf ../002_mesonh/XPREF.1.SEG01.000.nc .
+
+rm -f visu_LIMA.*.png
+
+ncl plot_LIMA.ncl
+
+${POSTRUN} display visu_LIMA.000001.png
+${POSTRUN} display visu_LIMA.000002.png
+${POSTRUN} display visu_LIMA.000003.png
+${POSTRUN} display visu_LIMA.000004.png
+${POSTRUN} display visu_LIMA.000005.png
+${POSTRUN} display visu_LIMA.000006.png
+exit 0
diff --git a/MY_RUN/KTEST/014_LIMA/Makefile b/MY_RUN/KTEST/014_LIMA/Makefile
index c03971716585263f5a1107ad18335ea48c7fc060..d2df0e9b7892909be537addfc14b71a2626254b3 100644
--- a/MY_RUN/KTEST/014_LIMA/Makefile
+++ b/MY_RUN/KTEST/014_LIMA/Makefile
@@ -1,6 +1,7 @@
all:
cd 001_prep_ideal_case && run_prep_ideal_case_xyz
cd 002_mesonh && run_mesonh_xyz
+ cd 003_ncl && run_ncl
ifneq "$(MNH_PYTHON)" "NO"
cd 004_python && run_python
endif
@@ -8,6 +9,7 @@ all:
clean:
cd 001_prep_ideal_case && clean_prep_ideal_case_xyz
cd 002_mesonh && clean_mesonh_xyz
+ cd 003_ncl && clean_ncl
ifneq "$(MNH_PYTHON)" "NO"
cd 004_python && clean_python
endif