Gaelle 2/12/2015 : MAJ 002_3Drelief
authorGaelle Tanguy <gaelle.tanguy@meteo.fr>
Wed, 2 Dec 2015 15:26:01 +0000 (15:26 +0000)
committerPhilippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Thu, 19 May 2016 14:44:49 +0000 (16:44 +0200)
MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/PRE_IDEA1.nam
MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz
MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam
MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH
MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO
MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz
MY_RUN/KTEST/002_3Drelief/002_mesonh/run_mesonh_xyz
MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl [new file with mode: 0755]
MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl [new file with mode: 0644]
MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl [new file with mode: 0755]
MY_RUN/KTEST/002_3Drelief/Makefile

index d11657a..a0b2914 100644 (file)
@@ -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 /
index 289c484..2cea9d3 100644 (file)
@@ -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" ,
index e7fe57c..12538e6 100644 (file)
@@ -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.,
index 4be0a88..c8c1cdf 100644 (file)
@@ -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.,
index f7d857c..faf1e21 100755 (executable)
@@ -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
 
index d77c92c..76f6378 100755 (executable)
@@ -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 (executable)
index 0000000..7b1634b
--- /dev/null
@@ -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 (file)
index 0000000..a86b8dd
--- /dev/null
@@ -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 (executable)
index 0000000..b901502
--- /dev/null
@@ -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
index 19d4b1e..e50e3b3 100755 (executable)
@@ -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