Gaelle 2/12/2015 : MAJ 002_3Drelief

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 /

MY_RUN/KTEST/002_3Drelief/001_prep_ideal_case/clean_prep_ideal_case_xyz

 #!/bin/sh
 set -x
 rm -f RELIEF3D* OUTPUT_LISTING* pipe* *.tex
+rm -f file_for_xtransfer

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" ,

MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_CEN4TH

+&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.,

MY_RUN/KTEST/002_3Drelief/002_mesonh/EXSEG1.nam_WENO

+&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.,

MY_RUN/KTEST/002_3Drelief/002_mesonh/clean_mesonh_xyz

 set -x
-rm -f PRESSURE    
+rm -f PRESSURE file_for_xtransfer
 rm -f RELIEF3D.*  
 rm -f REL3D.1* OUT* pipe_name
 

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}

MY_RUN/KTEST/002_3Drelief/005_ncl/clean_ncl

+
+rm -f  *.nc4 *.ps 

MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl

+;================================================;
+  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

MY_RUN/KTEST/002_3Drelief/005_ncl/run_ncl

+#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

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
+	cd 004_diaprog         && clean_diaprog
\ No newline at end of file
+	cd 005_ncl             && clean_ncl