Skip to content
Snippets Groups Projects
Commit ce1a28f6 authored by Gaelle DELAUTIER's avatar Gaelle DELAUTIER
Browse files

Gaelle 17/5/2018 : new ktest 005_ARM

parent ef5c6e53
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 968 additions and 1 deletion
MY_RUN/KTEST/005_ARM/001_prep_ideal/PRE_IDEA1.nam 0 → 100755
+ 104
− 0
View file @ ce1a28f6
&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F/
&NAM_REAL_PGD /
&NAM_DIMn_PRE NIMAX=1, NJMAX=1 /
&NAM_CONF_PRE LCARTESIAN=.TRUE., NVERB=5,
CIDEAL='RSOU', CZS='FLAT', LFORCING=.TRUE., LPACK=.FALSE.,
LBOUSS=.FALSE., CEQNSYS='DUR', LPERTURB=.FALSE.,
JPHEXT=1,NHALO=1 /
&NAM_PERT_PRE /
&NAM_CONFn LUSERV=.TRUE. /
&NAM_GRID_PRE XLAT0=35.762 XLON0=-96. /
&NAM_GRIDH_PRE XDELTAX=40000., XDELTAY=40000. /
&NAM_VER_GRID LTHINSHELL=.TRUE., NKMAX=100,
ZDZGRD=40., ZDZTOP=40.,
ZZMAX_STRGRD=1000. , ZSTRGRD=0., ZSTRTOP=0. /
&NAM_LUNITn CINIFILE='eurocs_noground', CINIFILEPGD='eurocs_noground_pgd' /
&NAM_POST_PRE /
&NAM_DYNn_PRE /
&NAM_LBCn_PRE CLBCX=2*"CYCL", CLBCY=2*"CYCL" /
&NAM_VPROF_PRE /
&NAM_GRn_PRE CSURF='EXTE'/
&NAM_CH_MNHCn_PRE /
&NAM_BLANK /
&NAM_PGD_SCHEMES CSEA='FLUX ' /
&NAM_COVER XUNIF_COVER(1)=1. /
RSOU
1997 6 21 41400.
'ZUVTHDMR'
0.
97000.
299.
15.20E-03
2
0. 10. 0.
5500. 10. 0.
8
50. 301.50 15.17E-03
350. 302.50 14.98E-03
650. 303.53 14.80E-03
700. 303.70 14.70E-03
1300. 307.13 13.50E-03
2500. 314.00 3.00E-03
5500. 343.20 3.00E-03
ZFRC
6
1997 6 21 41400.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
1000. 10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
1997 6 21 52200.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
1000. 10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
1997 6 21 63000.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
1000. 10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
1997 6 21 73800.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
1000. 10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
1997 6 21 84600.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
1000. 10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
1997 6 22 9000.
0.
97000.
299.
1.51913E-02
4
0. 10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
1000. 10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
3000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
5000. 10. 0. 0. 0. 0. 0. 0. 0. 0.
MY_RUN/KTEST/005_ARM/001_prep_ideal/clean_prep_ideal_case_xyz 0 → 100755
+ 1
− 0
View file @ ce1a28f6
rm -f eurocs* OUTPUT_LISTING* pipe* *.tex
MY_RUN/KTEST/005_ARM/001_prep_ideal/run_prep_ideal_case_xyz 0 → 100755
+ 6
− 0
View file @ ce1a28f6
#!/bin/sh
set -x
set -e
rm -f eurocs* OUTPUT_LISTING* pipe* *.tex
time ${MONORUN} PREP_IDEAL_CASE${XYZ}
#ddd --directory=~/DEV/MNH.V4.6.2/src/dir_obj_bug2 PREP_PGD${XYZ}
MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam 0 → 100755
+ 198
− 0
View file @ ce1a28f6
&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F /
&NAM_LUNITn CINIFILE = "eurocs_noground", CINIFILEPGD = "eurocs_noground_pgd" /
&NAM_CONFn LUSERV = .TRUE./
&NAM_DYNn XTSTEP = 60. /
&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" ,CTEMP_SCHEME='LEFR'/
&NAM_PARAMn CTURB = "TKEL", CRAD = "NONE",
CCLOUD = "ICE3", CDCONV="NONE", CSCONV="EDKF" /
&NAM_PARAM_RADn /
&NAM_PARAM_KAFRn XDTCONV= 60., LDIAGCONV=.TRUE./
&NAM_PARAM_GROUNDN /
&NAM_SEAFLUXn /
&NAM_ISBAn /
&NAM_DIAG_SURFn /
&NAM_LBCn /
&NAM_TURBn XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
LTURB_DIAG = .TRUE., LRMC01=.TRUE.,
LTURB_FLX = .TRUE., LSUBG_COND =.TRUE., LSIGMAS=.TRUE.,
LSIG_CONV=.FALSE., CSUBG_AUCV='CLFR' /
&NAM_PARAM_MFSHALLn CMF_UPDRAFT='EDKF',CMF_CLOUD='DIRE', LMIXUV=T, LMF_FLX=T /
&NAM_CH_MNHCn /
&NAM_CONF CCONF = "START", LFLAT = .TRUE., NMODEL = 1,
NVERB = 4, CEQNSYS='DUR', lforcing=.true.,
CEXP = "ARM__", CSEG = "CEN4T" /
&NAM_DYN XSEGLEN = 54000. , XASSELIN = 0.1, LCORIO = .true.,
LNUMDIFU = .FALSE.,
XALKTOP = 0.01, XALZBOT = 3500. /
&NAM_NESTING NDAD(1) = 1 , NDAD(2) = 1 , NDTRATIO(1) = 1, NDTRATIO(2) = 2 ,
XWAY(1) = 1., XWAY(2) = 2. /
&NAM_BACKUP XBAK_TIME(1,1)= 36000. /
&NAM_BUDGET CBUTYPE = 'NONE' /
&NAM_BU_RU LBU_RU = .FALSE. /
&NAM_BU_RV LBU_RV = .FALSE. /
&NAM_BU_RW LBU_RW = .FALSE. /
&NAM_BU_RTH LBU_RTH = .FALSE. /
&NAM_BU_RTKE LBU_RTKE = .FALSE. /
&NAM_BU_RRV LBU_RRV = .FALSE. /
&NAM_BU_RRC LBU_RRC = .FALSE. /
&NAM_BU_RRR LBU_RRR = .FALSE. /
&NAM_BU_RRI LBU_RRI = .FALSE. /
&NAM_BU_RRS LBU_RRS = .FALSE. /
&NAM_BU_RRG LBU_RRG = .FALSE. /
&NAM_BU_RRH LBU_RRH = .FALSE. /
&NAM_BU_RSV LBU_RSV = .FALSE. /
&NAM_LES LLES_MEAN=.TRUE., LLES_RESOLVED=.TRUE., LLES_SUBGRID=.TRUE.,
LLES_UPDRAFT=.FALSE., LLES_DOWNDRAFT=.FALSE.,
CBL_HEIGHT_DEF='WTV',
CLES_NORM_TYPE='CONV',
XLES_TEMP_SAMPLING=300.,
LLES_NEB_MASK = .TRUE.,
XLES_TEMP_MEAN_STEP=3600.,
XLES_TEMP_MEAN_START=3600., XLES_TEMP_MEAN_END=32400. /
&NAM_BLANK /
&NAM_FRC LGEOST_UV_FRC=.TRUE., LTEND_THRV_FRC=.TRUE. /
&NAM_CH_SOLVER /
&NAM_PARAM_ICE /
&NAM_IDEAL_FLUX
NFORCT = 2,
NFORCF = 31,
XTIMET(1) = 0.,
XTIMET(2) = 54000.,
XTIMEF(1) = 0.,
XTIMEF(2) = 1800.,
XTIMEF(3) = 3600.,
XTIMEF(4) = 5400.,
XTIMEF(5) = 7200.,
XTIMEF(6) = 9000.,
XTIMEF(7) = 10800.,
XTIMEF(8) = 12600.,
XTIMEF(9) = 14400.,
XTIMEF(10) = 16200.,
XTIMEF(11) = 18000.,
XTIMEF(12) = 19800.,
XTIMEF(13) = 21600.,
XTIMEF(14) = 23400.,
XTIMEF(15) = 25200.,
XTIMEF(16) = 27000.,
XTIMEF(17) = 28800.,
XTIMEF(18) = 30600.,
XTIMEF(19) = 32400.,
XTIMEF(20) = 34200.,
XTIMEF(21) = 36000.,
XTIMEF(22) = 37800.,
XTIMEF(23) = 39600.,
XTIMEF(24) = 41400.,
XTIMEF(25) = 43200.,
XTIMEF(26) = 45000.,
XTIMEF(27) = 46800.,
XTIMEF(28) = 48600.,
XTIMEF(29) = 50400.,
XTIMEF(30) = 52200.,
XTIMEF(31) = 54000.,
XSFTH(1) = -30.,
XSFTH(2) = -15.0,
XSFTH(3) = 0.,
XSFTH(4) = 15.0,
XSFTH(5) = 30.0,
XSFTH(6) = 45.0,
XSFTH(7) = 60.0,
XSFTH(8) = 75.0,
XSFTH(9) = 90.0,
XSFTH(10) = 100.0,
XSFTH(11) = 110.0,
XSFTH(12) = 120.0,
XSFTH(13) = 130.0,
XSFTH(14) = 140.0,
XSFTH(15) = 140.0,
XSFTH(16) = 140.0,
XSFTH(17) = 132.0,
XSFTH(18) = 124.0,
XSFTH(19) = 116.0,
XSFTH(20) = 108.0,
XSFTH(21) = 100.0,
XSFTH(22) = 78.0,
XSFTH(23) = 56.0,
XSFTH(24) = 34.0,
XSFTH(25) = 12.0,
XSFTH(26) = -10.0,
XSFTH(27) = -10.0,
XSFTH(28) = -10.0,
XSFTH(29) = -10.0,
XSFTH(30) = -10.0,
XSFTH(31) = -10.0,
CSFTQ='kg/m2/s'
XSFTQ(1)=1.99936020473448506E-006
XSFTQ(2)=1.42454414587332047E-005
XSFTQ(3)=2.64915227127319252E-005
XSFTQ(4)=3.87376039667306440E-005
XSFTQ(5)=5.09836852207293695E-005
XSFTQ(6)=6.32297664747280883E-005
XSFTQ(7)=7.54758477287268071E-005
XSFTQ(8)=8.77219289827255259E-005
XSFTQ(9)=9.99680102367242446E-005
XSFTQ(10)=1.15962891874600132E-004
XSFTQ(11)=1.31957773512476006E-004
XSFTQ(12)=1.47952655150351879E-004
XSFTQ(13)=1.63947536788227780E-004
XSFTQ(14)=1.79942418426103654E-004
XSFTQ(15)=1.89939219449776085E-004
XSFTQ(16)=1.99936020473448489E-004
XSFTQ(17)=1.93538067818298140E-004
XSFTQ(18)=1.87140115163147790E-004
XSFTQ(19)=1.80742162507997441E-004
XSFTQ(20)=1.74344209852847091E-004
XSFTQ(21)=1.67946257197696742E-004
XSFTQ(22)=1.48752399232245693E-004
XSFTQ(23)=1.29558541266794618E-004
XSFTQ(24)=1.10364683301343569E-004
XSFTQ(25)=9.11708253358925209E-005
XSFTQ(26)=7.19769673704414589E-005
XSFTQ(27)=5.39827255278310908E-005
XSFTQ(28)=3.59884836852207294E-005
XSFTQ(29)=1.79942418426103647E-005
XSFTQ(30)=0.0000000000000000
XSFTQ(31)=-1.79942418426103647E-005
XSFCO2(1) = 0.,
XSFCO2(2) = 0.,
XSFCO2(3) = 0.,
XSFCO2(4) = 0.,
XSFCO2(5) = 0.,
XSFCO2(6) = 0.,
XSFCO2(7) = 0.,
XSFCO2(8) = 0.,
XSFCO2(9) = 0.,
XSFCO2(10) = 0.,
XSFCO2(11) = 0.,
XSFCO2(12) = 0.,
XSFCO2(13) = 0.,
XSFCO2(14) = 0.,
XSFCO2(15) = 0.,
XSFCO2(16) = 0.,
XSFCO2(17) = 0.,
XSFCO2(18) = 0.,
XSFCO2(19) = 0.,
XSFCO2(20) = 0.,
XSFCO2(21) = 0.,
XSFCO2(22) = 0.,
XSFCO2(23) = 0.,
XSFCO2(24) = 0.,
XSFCO2(25) = 0.,
XSFCO2(26) = 0.,
XSFCO2(27) = 0.,
XSFCO2(28) = 0.,
XSFCO2(29) = 0.,
XSFCO2(30) = 0.,
XSFCO2(31) = 0.,
CUSTARTYPE = 'Z0 ',
XZ0=0.035,
XALB = 0.,
XEMIS = 1.,
XTSRAD(1) = 273.16,
XTSRAD(2) = 273.16
/
MY_RUN/KTEST/005_ARM/002_mesonh/clean_mesonh_xyz 0 → 100755
+ 2
− 0
View file @ ce1a28f6
rm -f ARM* pipe_name
rm -f eurocs* OUT*
MY_RUN/KTEST/005_ARM/002_mesonh/run_mesonh_xyz 0 → 100755
+ 6
− 0
View file @ ce1a28f6
set -x
set -e
ln -fs ../001_prep_ideal/eurocs_noground.* .
ln -fs ../001_prep_ideal/eurocs_noground_pgd.* .
rm -f ARM* OUT*
time ${MONORUN} MESONH${XYZ}
MY_RUN/KTEST/005_ARM/003_diaprog/clean_diaprog 0 → 100755
+ 3
− 0
View file @ ce1a28f6
rm -f ARM**
rm -f gmeta OUT*
rm -f rm dir.*:*
MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm 0 → 100755
+ 77
− 0
View file @ ce1a28f6
visu
_file1_'ARM__.1.CEN4T.000'
print groups
lcolarea=t
lminmax=t
liso=f
xhmax=3000
LCOLZERO=t
NCOLZERO=1
LMINMAX=t
LVPTUSER=T
XVPTVL=0.20
XVPTR=0.80
NIMNMX=1
LTITDEF=F
LVPTUSER=T
XVPTVL=0.20
XVPTR=0.80
CTITB1='white'
CTITB2='white'
CTITB3='white'
CTITB4='white'
CTITT1='white'
CTITT2='white'
CTITVAR2='white'
LDATFILE=.FALSE.
XISOMIN_SBG_TKE=0
XISOMAX_SBG_TKE=2.5
XDIAINT_SBG_TKE=0.25
XISOMIN_MEAN_TH=298
XISOMAX_MEAN_TH=316
XDIAINT_MEAN_TH=1
XISOMIN_MEAN_U=3
XISOMAX_MEAN_U=11.
XDIAINT_MEAN_U=0.5
XISOMIN_MEAN_V=0
XISOMAX_MEAN_V=3.2
XDIAINT_MEAN_V=0.2
XISOMIN_MEAN_RC=0
XISOMAX_MEAN_RC=0.2E-03
XDIAINT_MEAN_RC=0.1E-04
XISOMIN_MEAN_RR=0
XISOMAX_MEAN_RR=0.2E-04
XDIAINT_MEAN_RR=0.1E-05
XISOMIN_SBG_WTHL=-0.3
XISOMAX_SBG_WTHL=0.12
XDIAINT_SBG_WTHL=0.02
XISOMIN_SBG_WRT=-0.1E-03
XISOMAX_SBG_WRT=0.4E-03
XDIAINT_SBG_WRT=0.05E-03
SBG_TKE_pvt_
MEAN_TH_pvt_
MEAN_U_pvt_
MEAN_V_pvt_
MEAN_RC_pvt_
MEAN_RR_pvt_
SBG_WTHL_pvt_
SBG_WRT_pvt_
THLUP_MF_pvt_
RTUP_MF_pvt_
RVUP_MF_pvt_
RCUP_MF_pvt_
RIUP_MF_pvt_
WUP_MF_pvt_
MAFLX_MF_pvt_
DETR_MF_pvt_
ENTR_MF_pvt_
FRCUP_MF_pvt_
THVUP_MF_pvt_
WTHL_MF_pvt_
WRT_MF_pvt_
WTHV_MF_pvt_
WU_MF_pvt_
WV_MF_pvt_
quit
MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog 0 → 100755
+ 7
− 0
View file @ ce1a28f6
set -x
set -e
ln -sf ../002_mesonh/ARM__.1.CEN4T.000.* .
rm -f dir.0*
${POSTRUN} diaprog < dir_arm
MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl 0 → 100755
+ 2
− 0
View file @ ce1a28f6
rm -f ARM__.1.CEN4T.000* *.png
MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl 0 → 100644
+ 513
− 0
View file @ ce1a28f6
;================================================;
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
vname="SBG_TKE___PROC1"
SBGTKE = fichier1->$vname$(0,:,:,0,0)
SBGTKE@long_name="SBG_TKE"
SBGTKE!1="nzh"
SBGTKE!0="TIME"
SBG_TKE=transpose(SBGTKE)
vname="MEAN_TH___PROC1"
MEANTH = fichier1->$vname$(0,:,:,0,0) ; MEAN_TH
MEANTH@long_name="MEAN_TH"
MEANTH@units="K"
MEANTH!1="nzh"
MEANTH!0="TIME"
MEAN_TH=transpose(MEANTH)
vname="MEAN_U___PROC1"
MEANU = fichier1->$vname$(0,:,:,0,0)
MEANU!1="nzh"
MEANU!0="TIME"
MEAN_U=transpose(MEANU)
MEAN_U@long_name="MEAN_U"
vname="MEAN_V___PROC1"
MEANV = fichier1->$vname$(0,:,:,0,0)
MEANV!1="nzh"
MEANV!0="TIME"
MEAN_V=transpose(MEANV)
MEAN_V@long_name="MEAN_V"
vname="MEAN_RC___PROC1"
MEANRC = fichier1->$vname$(0,:,:,0,0)
MEANRC!1="nzh"
MEANRC!0="TIME"
MEAN_RC=transpose(MEANRC)
MEAN_RC@long_name="MEAN_RC"
vname="MEAN_RR___PROC1"
MEANRR = fichier1->$vname$(0,:,:,0,0)
MEANRR!1="nzh"
MEANRR!0="TIME"
MEAN_RR=transpose(MEANRR)
MEAN_RR@long_name="MEAN_RR"
vname="SBG_WTHL___PROC1"
SBGWTHL = fichier1->$vname$(0,:,:,0,0)
SBGWTHL!1="nzh"
SBGWTHL!0="TIME"
SBG_WTHL=transpose(SBGWTHL)
SBG_WTHL@long_name="SBG_WTHL"
vname="SBG_WRT___PROC1"
SBGWRT = fichier1->$vname$(0,:,:,0,0)
SBGWRT!1="nzh"
SBGWRT!0="TIME"
SBG_WRT=transpose(SBGWRT)
SBG_WRT@long_name="SBG_WRT"
vname="THLUP_MF___PROC1"
THLUPMF= fichier1->$vname$(0,:,:,0,0)
THLUPMF!1="nzh"
THLUPMF!0="TIME"
THLUP_MF=transpose(THLUPMF)
THLUP_MF@long_name="THLUP_MF"
vname="RTUP_MF___PROC1"
RTUPMF = fichier1->$vname$(0,:,:,0,0)
RTUPMF!1="nzh"
RTUPMF!0="TIME"
RTUP_MF=transpose(RTUPMF)
RTUP_MF@long_name="RTUP_MF"
vname="RVUP_MF___PROC1"
RVUPMF = fichier1->$vname$(0,:,:,0,0)
RVUPMF!1="nzh"
RVUPMF!0="TIME"
RVUP_MF=transpose(RVUPMF)
RVUP_MF@long_name="RVUP_MF"
vname="RCUP_MF___PROC1"
RCUPMF = fichier1->$vname$(0,:,:,0,0)
RCUPMF!1="nzh"
RCUPMF!0="TIME"
RCUP_MF=transpose(RCUPMF)
RCUP_MF@long_name="RCUP_MF"
vname="RIUP_MF___PROC1"
RIUPMF = fichier1->$vname$(0,:,:,0,0)
RIUPMF!1="nzh"
RIUPMF!0="TIME"
RIUP_MF=transpose(RIUPMF)
RIUP_MF@long_name="RIUP_MF"
vname="WUP_MF___PROC1"
WUPMF = fichier1->$vname$(0,:,:,0,0)
WUPMF!1="nzh"
WUPMF!0="TIME"
WUP_MF=transpose(WUPMF)
WUP_MF@long_name="WUP_MF"
vname="MAFLX_MF___PROC1"
MAFLXMF = fichier1->$vname$(0,:,:,0,0)
MAFLXMF!1="nzh"
MAFLXMF!0="TIME"
MAFLX_MF=transpose(MAFLXMF)
MAFLX_MF@long_name="MAFLX_MF"
vname="DETR_MF___PROC1"
DETRMF = fichier1->$vname$(0,:,:,0,0)
DETRMF!1="nzh"
DETRMF!0="TIME"
DETR_MF=transpose(DETRMF)
DETR_MF@long_name="DETR_MF"
vname="ENTR_MF___PROC1"
ENTRMF = fichier1->$vname$(0,:,:,0,0)
ENTRMF!1="nzh"
ENTRMF!0="TIME"
ENTR_MF=transpose(ENTRMF)
ENTR_MF@long_name="ENTR_MF"
vname="FRCUP_MF___PROC1"
FRCUPMF = fichier1->$vname$(0,:,:,0,0)
FRCUPMF!1="nzh"
FRCUPMF!0="TIME"
FRCUP_MF=transpose(FRCUPMF)
FRCUP_MF@long_name="FRCUP_MF"
vname="THVUP_MF___PROC1"
THVUPMF = fichier1->$vname$(0,:,:,0,0)
THVUPMF!1="nzh"
THVUPMF!0="TIME"
THVUP_MF=transpose(THVUPMF)
THVUP_MF@long_name="THVUP_MF"
vname="WTHL_MF___PROC1"
WTHLMF = fichier1->$vname$(0,:,:,0,0)
WTHLMF!1="nzh"
WTHLMF!0="TIME"
WTHL_MF=transpose(WTHLMF)
WTHL_MF@long_name="WTHL_MF"
vname="WRT_MF___PROC1"
WRTMF = fichier1->$vname$(0,:,:,0,0)
WRTMF!1="nzh"
WRTMF!0="TIME"
WRT_MF=transpose(WRTMF)
WRT_MF@long_name="WRT_MF"
vname="WTHV_MF___PROC1"
WTHVMF = fichier1->$vname$(0,:,:,0,0)
WTHVMF!1="nzh"
WTHVMF!0="TIME"
WTHV_MF=transpose(WTHVMF)
WTHV_MF@long_name="WTHV_MF"
vname="WU_MF___PROC1"
WUMF = fichier1->$vname$(0,:,:,0,0)
WUMF!1="nzh"
WUMF!0="TIME"
WU_MF=transpose(WUMF)
WU_MF@long_name="WU_MF"
vname="WV_MF___PROC1"
WVMF = fichier1->$vname$(0,:,:,0,0)
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
MY_RUN/KTEST/005_ARM/004_ncl/run_ncl 0 → 100755
+ 34
− 0
View file @ ce1a28f6
#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.* .
ncl plot_arm.ncl
display visu_ARM.000001.png
display visu_ARM.000002.png
display visu_ARM.000003.png
display visu_ARM.000004.png
display visu_ARM.000005.png
display visu_ARM.000006.png
display visu_ARM.000007.png
display visu_ARM.000008.png
display visu_ARM.000009.png
display visu_ARM.000010.png
display visu_ARM.000011.png
display visu_ARM.000012.png
display visu_ARM.000013.png
display visu_ARM.000014.png
display visu_ARM.000015.png
display visu_ARM.000016.png
display visu_ARM.000017.png
display visu_ARM.000018.png
display visu_ARM.000019.png
display visu_ARM.000020.png
display visu_ARM.000021.png
display visu_ARM.000022.png
display visu_ARM.000023.png
display visu_ARM.000024.png
exit 0
MY_RUN/KTEST/005_ARM/Makefile 0 → 100755
+ 11
− 0
View file @ ce1a28f6
all:
cd 001_prep_ideal && run_prep_ideal_case_xyz
cd 002_mesonh && run_mesonh_xyz
cd 003_diaprog && run_diaprog
cd 004_ncl && run_ncl
clean:
cd 001_prep_ideal && clean_prep_ideal_case_xyz
cd 002_mesonh && clean_mesonh_xyz
cd 003_diaprog && clean_diaprog
cd 004_ncl && clean_ncl
MY_RUN/KTEST/005_ARM/run_005_ARM 0 → 100755
+ 3
− 0
View file @ ce1a28f6
set -e
make clean
make
\ No newline at end of file
src/Makefile
+ 1
− 1
View file @ ce1a28f6
...@@ -420,7 +420,7 @@ endif ...@@ -420,7 +420,7 @@ endif
PWD_KTEST=$(SRC_MESONH)/MY_RUN/KTEST PWD_KTEST=$(SRC_MESONH)/MY_RUN/KTEST
EXAMPLES=001_2Drelief 002_3Drelief 003_KW78 004_Reunion 007_16janvier 014_LIMA EXAMPLES=001_2Drelief 002_3Drelief 003_KW78 004_Reunion 005_ARM 007_16janvier 014_LIMA
examples : $(EXAMPLES) examples : $(EXAMPLES)
EXAMPLES_CHIMIE=009_ICARTT 011_KW78CHEM EXAMPLES_CHIMIE=009_ICARTT 011_KW78CHEM
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme