From 97a69b7a705884dd908daf295694f58dab464b0a Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Wed, 2 Dec 2015 15:26:01 +0000
Subject: [PATCH] Gaelle 2/12/2015 : MAJ 002_3Drelief
---
.../001_prep_ideal_case/PRE_IDEA1.nam | 2 +-
.../clean_prep_ideal_case_xyz | 1 +
.../KTEST/002_3Drelief/002_mesonh/EXSEG1.nam | 2 +-
.../002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH | 4 +-
.../002_3Drelief/002_mesonh/EXSEG1.nam_WENO | 6 +-
.../002_3Drelief/002_mesonh/clean_mesonh_xyz | 2 +-
.../002_3Drelief/002_mesonh/run_mesonh_xyz | 2 +-
MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl | 2 +
.../002_3Drelief/005_ncl/plot_3Drelief.ncl | 71 +++++++++++++++++++
MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl | 9 +++
MY_RUN/KTEST/002_3Drelief/Makefile | 4 +-
11 files changed, 97 insertions(+), 8 deletions(-)
create mode 100755 MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl
create mode 100644 MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl
create mode 100755 MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl
diff --git a/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/PRE_IDEA1.nam b/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/PRE_IDEA1.nam
index d11657a5e..a0b291440 100644
--- a/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/PRE_IDEA1.nam
+++ b/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/PRE_IDEA1.nam
@@ -17,7 +17,7 @@
&NAM_GRIDH_PRE XDELTAX=4.E3 XDELTAY=4.E3
XHMAX=1000. XAX=10.E3 XAY=10.E3
NIZS=16 NJZS=16 /
-&NAM_LUNITn CINIFILE='RELIEF3D.520' /
+&NAM_LUNITn CINIFILE='RELIEF3D' /
&NAM_DYNn_PRE
! CPRESOPT ='ZRESI' ,
NITR=4 /
diff --git a/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz b/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz
index 05c1c35ba..3dd534c9e 100755
--- a/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz
+++ b/MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz
@@ -1,3 +1,4 @@
#!/bin/sh
set -x
rm -f RELIEF3D* OUTPUT_LISTING* pipe* *.tex
+rm -f file_for_xtransfer
diff --git a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam
index 289c484be..2cea9d3e0 100644
--- a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam
+++ b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam
@@ -2,7 +2,7 @@
&NAM_CONFZ
! NZ_VERB = 5 , NZ_PROC = 0 , NB_PROCIO_R = 8 , NB_PROCIO_W = 8
/
-&NAM_LUNITn CINIFILE = "RELIEF3D.520" /
+&NAM_LUNITn CINIFILE = "RELIEF3D" /
&NAM_CONFn LUSERV = F /
&NAM_DYNn XTSTEP = 180., LITRADJ=T,
! CPRESOPT = "ZRESI" ,
diff --git a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH
index e7fe57c9c..12538e69b 100644
--- a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH
+++ b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH
@@ -1,7 +1,8 @@
+&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F /
&NAM_CONFZ
! NZ_VERB = 5 , NZ_PROC = 0 , NB_PROCIO_R = 8 , NB_PROCIO_W = 1
/
-&NAM_LUNITn CINIFILE = "RELIEF3D.410" /
+&NAM_LUNITn CINIFILE = "RELIEF3D" /
&NAM_CONFn LUSERV = F /
&NAM_DYNn XTSTEP = 45., LITRADJ=T,
! CPRESOPT = "ZRESI" ,
@@ -12,6 +13,7 @@
&NAM_CONF CCONF = "START", NMODEL = 1,
CSPLIT='BSPLITTING' ,
NVERB = 5, CEXP = "REL3D", CSEG = "EXP01" ,
+ !JPHEXT=3 NHALO=3
/
&NAM_DYN XSEGLEN = 1000., XASSELIN = 0.2, LCORIO = F,
XALKTOP = 0.005, XALZBOT = 12570.,
diff --git a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO
index 4be0a8813..c8c1cdf51 100644
--- a/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO
+++ b/MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO
@@ -1,17 +1,19 @@
+&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F /
&NAM_CONFZ
! NZ_VERB = 5 , NZ_PROC = 0 , NB_PROCIO_R = 8 , NB_PROCIO_W = 1
/
-&NAM_LUNITn CINIFILE = "RELIEF3D.410" /
+&NAM_LUNITn CINIFILE = "RELIEF3D" /
&NAM_CONFn LUSERV = F /
&NAM_DYNn XTSTEP = 180., LITRADJ=T,
! CPRESOPT = "ZRESI" ,
LHORELAX_UVWTH = T, LHORELAX_RV = F, LVE_RELAX = T,
NRIMX = 6, NRIMY = 3, XRIMKMAX = .001666 /
-&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",CMET_ADV_SCHEME="PPM_01", CSV_ADV_SCHEME="PPM_01" /
+&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",CMET_ADV_SCHEME="PPM_01", CSV_ADV_SCHEME="PPM_01" NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_CONF CCONF = "START", NMODEL = 1,
CSPLIT='BSPLITTING' ,
NVERB = 5, CEXP = "REL3D", CSEG = "EXP01" ,
+ !JPHEXT=3 NHALO=3
/
&NAM_DYN XSEGLEN = 1000., XASSELIN = 0.2, LCORIO = F,
XALKTOP = 0.005, XALZBOT = 12570.,
diff --git a/MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz b/MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz
index f7d857c17..faf1e21a6 100755
--- a/MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz
+++ b/MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz
@@ -1,5 +1,5 @@
set -x
-rm -f PRESSURE
+rm -f PRESSURE file_for_xtransfer
rm -f RELIEF3D.*
rm -f REL3D.1* OUT* pipe_name
diff --git a/MY_RUN/KTEST/002_3Drelief/002_mesonh/run_mesonh_xyz b/MY_RUN/KTEST/002_3Drelief/002_mesonh/run_mesonh_xyz
index d77c92c5c..76f63781e 100755
--- a/MY_RUN/KTEST/002_3Drelief/002_mesonh/run_mesonh_xyz
+++ b/MY_RUN/KTEST/002_3Drelief/002_mesonh/run_mesonh_xyz
@@ -4,6 +4,6 @@
#MNH_LIC for details. version 1.
set -x
set -e
-ln -fs ../001_prep_ideal_case/RELIEF3D.520.*??? .
+ln -fs ../001_prep_ideal_case/RELIEF3D.*??? .
rm -f REL3D.1* OUT*
time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl b/MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl
new file mode 100755
index 000000000..7b1634bda
--- /dev/null
+++ b/MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl
@@ -0,0 +1,2 @@
+
+rm -f *.nc4 *.ps
diff --git a/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl b/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl
new file mode 100644
index 000000000..a86b8dd7d
--- /dev/null
+++ b/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl
@@ -0,0 +1,71 @@
+;================================================;
+ 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("REL3D.1.EXP01.002.nc4", "r")
+;==================================================;
+; Open the workstation
+;==================================================;
+ type = "ps"
+ wks = gsn_open_wks(type,"visu_3Drelief")
+;=================================================;
+; lecture des différents champs
+;=================================================;
+ zs = fichier1->ZS(0,:) ; ZS
+ zhat = fichier1->ZHAT(:) ; ZHAT
+ xhat = fichier1->XHAT(:) ; XHAT
+
+ ut= fichier1->UT(:,:,:) ; UT
+ ut@long_name="Zonal wind"
+ ut@units="m/s"
+
+
+;=================================================;
+; Récupération des dimensions des tableaux
+;=================================================;
+ mdims = dimsizes(ut)
+ nd=dimsizes(mdims)
+ imax=mdims(nd-1)
+ jmax=mdims(nd-2)
+ kmax=mdims(nd-3)
+
+;=================================================;
+; Unstagger U
+;=================================================;
+ ut!0="Z"
+ ut!1="Y"
+ ut!2="X"
+ u = wrf_user_unstagger(ut,"X")
+
+;=================================================;
+; Set some other basic resources
+;=================================================;
+ resmap = True
+ resmap@gsnFrame = False
+ resmap@gsnDraw = False
+ resmap@gsnMaximize = True
+ resmap@gsnPaperOrientation = "portrait"
+ resmap@gsnSpreadColors= True
+ resmap@tiMainString ="UT level 2 "
+ resmap@cnFillOn= True
+ resmap@cnLinesOn= False
+
+;=================================================;
+; TRACE
+;=================================================;
+ gsn_define_colormap(wks,"rainbow")
+
+ resmap@cnLevelSelectionMode = "ManualLevels"
+ resmap@cnLevelSpacingF = 0.25
+ resmap@cnMinLevelValF = 8
+ resmap@cnMaxLevelValF = 14
+ plot_ut = gsn_csm_contour(wks,u(1,:,1:),resmap)
+ draw(plot_ut)
+ frame(wks)
+
+end
diff --git a/MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl b/MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl
new file mode 100755
index 000000000..b9015024b
--- /dev/null
+++ b/MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl
@@ -0,0 +1,9 @@
+#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
+set -e
+ln -sf ../002_mesonh/REL3D.1.EXP01.002.nc4 .
+
+ncl plot_3Drelief.ncl
diff --git a/MY_RUN/KTEST/002_3Drelief/Makefile b/MY_RUN/KTEST/002_3Drelief/Makefile
index 19d4b1eaa..e50e3b333 100755
--- a/MY_RUN/KTEST/002_3Drelief/Makefile
+++ b/MY_RUN/KTEST/002_3Drelief/Makefile
@@ -3,9 +3,11 @@ all:
cd 002_mesonh && run_mesonh_xyz
cd 003_convdia && run_conv2dia
cd 004_diaprog && run_diaprog
+ cd 005_ncl && run_ncl
clean:
cd 001_prep_ideal_case && clean_prep_ideal_case_xyz
cd 002_mesonh && clean_mesonh_xyz
cd 003_convdia && clean_conv2dia
- cd 004_diaprog && clean_diaprog
\ No newline at end of file
+ cd 004_diaprog && clean_diaprog
+ cd 005_ncl && clean_ncl
--
GitLab