diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam index fd46b178fb402351cb83c779b3713c83158dc704..f588e40759202d49a658a076f0fa7961f505a9ef 100644 --- a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam +++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam @@ -3,12 +3,12 @@ / &NAM_LUNITn CINIFILE = "REUNION_IDEA_410", CINIFILEPGD="REUNION_PGD_1km5_410" / -&NAM_DYNn XTSTEP = 5., +&NAM_DYNn XTSTEP = 20., !CPRESOPT = "ZRESI", LITRADJ = T, LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T, NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. / -&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" , CMET_ADV_SCHEME = "FCT2ND" / +&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K" , CMET_ADV_SCHEME = "PPM_01" / &NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE", CCLOUD = "NONE", CDCONV = "NONE" / &NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" / @@ -20,7 +20,7 @@ NVERB = 5, CEXP = "REUNI", CSEG = "00A20" , CSPLIT="BSPLITTING" / &NAM_DYN XSEGLEN = 40., - LCORIO = F, LNUMDIFU = T, + LCORIO = F, LNUMDIFU = F, XALKTOP = 0.01, XALZBOT = 14000. / &NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20., XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. / diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH new file mode 100644 index 0000000000000000000000000000000000000000..cf5d50f09dac499c5022b46f638ff7ccab153e73 --- /dev/null +++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH @@ -0,0 +1,38 @@ +&NAM_CONFZ + ! NZ_VERB=5 , NB_PROCIO_R=8 , NB_PROCIO_W=1 +/ +&NAM_LUNITn CINIFILE = "REUNION_IDEA_410", + CINIFILEPGD="REUNION_PGD_1km5_410" / +&NAM_DYNn XTSTEP = 5., + !CPRESOPT = "ZRESI", + LITRADJ = T, + LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T, + NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. / +&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" , CMET_ADV_SCHEME = "PPM_01" / +&NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE", + CCLOUD = "NONE", CDCONV = "NONE" / +&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" / +&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F / +&NAM_SERIESn NIBOXL=25 NIBOXH=50 NJBOXL=20 NJBOXH=40 + NKCLS=3 NKLOW=12 NKUP=18 NKMID=14 + NBJSLICE=1 NJSLICEL(1)=30 NJSLICEH(1)=35 / +&NAM_CONF CCONF = "START", NMODEL = 1, + NVERB = 5, CEXP = "REUNI", CSEG = "00A20" , + CSPLIT="BSPLITTING" / +&NAM_DYN XSEGLEN = 40., + LCORIO = F, LNUMDIFU = T, + XALKTOP = 0.01, XALZBOT = 14000. / +&NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20., + XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. / +&NAM_BUDGET CBUTYPE='CART', XBULEN=200., + NBUKL=2, NBUKH=51, NBUIL=15, NBUIH=64, NBUJL=15, NBUJH=68, + LBU_KCP=.TRUE., LBU_ICP=.FALSE., LBU_JCP=.FALSE. / +&NAM_BU_RTH LBU_RTH=.FALSE. / +&NAM_BU_RTKE LBU_RTKE=.FALSE., + NDPTKE=1, NTPTKE=1, NDISSTKE=0, NTRTKE=0 / +&NAM_SERIES LSERIES=.TRUE. LMASKLANDSEA=.TRUE. LWMINMAX=.FALSE. / +&NAM_SEAFLUXn CSEA_ALB = 'UNIF' / +&NAM_ISBAn / +&NAM_DIAG_SURFn LSURF_BUDGET=F / +&NAM_DIAG_ISBAn / + diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO new file mode 100644 index 0000000000000000000000000000000000000000..f588e40759202d49a658a076f0fa7961f505a9ef --- /dev/null +++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO @@ -0,0 +1,38 @@ +&NAM_CONFZ + ! NZ_VERB=5 , NB_PROCIO_R=8 , NB_PROCIO_W=1 +/ +&NAM_LUNITn CINIFILE = "REUNION_IDEA_410", + CINIFILEPGD="REUNION_PGD_1km5_410" / +&NAM_DYNn XTSTEP = 20., + !CPRESOPT = "ZRESI", + LITRADJ = T, + LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T, + NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. / +&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K" , CMET_ADV_SCHEME = "PPM_01" / +&NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE", + CCLOUD = "NONE", CDCONV = "NONE" / +&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" / +&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F / +&NAM_SERIESn NIBOXL=25 NIBOXH=50 NJBOXL=20 NJBOXH=40 + NKCLS=3 NKLOW=12 NKUP=18 NKMID=14 + NBJSLICE=1 NJSLICEL(1)=30 NJSLICEH(1)=35 / +&NAM_CONF CCONF = "START", NMODEL = 1, + NVERB = 5, CEXP = "REUNI", CSEG = "00A20" , + CSPLIT="BSPLITTING" / +&NAM_DYN XSEGLEN = 40., + LCORIO = F, LNUMDIFU = F, + XALKTOP = 0.01, XALZBOT = 14000. / +&NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20., + XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. / +&NAM_BUDGET CBUTYPE='CART', XBULEN=200., + NBUKL=2, NBUKH=51, NBUIL=15, NBUIH=64, NBUJL=15, NBUJH=68, + LBU_KCP=.TRUE., LBU_ICP=.FALSE., LBU_JCP=.FALSE. / +&NAM_BU_RTH LBU_RTH=.FALSE. / +&NAM_BU_RTKE LBU_RTKE=.FALSE., + NDPTKE=1, NTPTKE=1, NDISSTKE=0, NTRTKE=0 / +&NAM_SERIES LSERIES=.TRUE. LMASKLANDSEA=.TRUE. LWMINMAX=.FALSE. / +&NAM_SEAFLUXn CSEA_ALB = 'UNIF' / +&NAM_ISBAn / +&NAM_DIAG_SURFn LSURF_BUDGET=F / +&NAM_DIAG_ISBAn / + diff --git a/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion b/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion index 67da61d2a5be81a9c982021ce8681929499c0dad..fa9ef813d9b7f69396c15427438519454d2396b8 100644 --- a/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion +++ b/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion @@ -7,25 +7,25 @@ _file_'REUNI.1.00A20.004dg' ! ouverture de la fenetre d affichage visu ! module du vent et vecteurs du niveau K=2 -MUMVM_K_2_ON_UMVM_K_2 +MUTVT_K_2_ON_UTVT_K_2 ! temperature potentielle a 1500m -THM_Z_1500 +THT_Z_1500 pression et vent en vecteurs a 1500m -PABSM_Z_1500_ON_UMVM_Z_1500 +PABST_Z_1500_ON_UTVT_Z_1500 ! definition d une coupe verticale NIDEBCOU=35 NJDEBCOU=2 NLANGLE=90 NLMAX=80 ! trace de la coupe verticale dans le plan horizontal LTRACECV=T ! temperature potentielle -THM_CV_ +THT_CV_ ! plus de trace de la coupe verticale LTRACECV=F ! difference pour la temperature potentielle avec le champ initial -THM_CV__MINUS_LSTHM_CV_ +THT_CV__MINUS_LSTHM_CV_ ! difference pour la composante du vent selon Y avec le champ initial -VM_CV__MINUS_LSVM_CV_ +VT_CV__MINUS_LSVM_CV_ ! difference pour la vitesse verticale avec le champ initial -WM_CV__ON_ULMWM_CV_ +WT_CV__ON_ULMWM_CV_ ! ! a vous d indiquer vos directives quit diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf index 020304e70b5116c4380776343cab75ba6cbc13f3..f39dba49c772930ca52215ac0a49211adfdbc27e 100644 --- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf +++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf @@ -14,16 +14,16 @@ KCDL 0 0 CONF - UM - VM - WM + UT + VT + WT ALT -THM -PABSM -RVM -RCM -RRM -RIM +THT +PABST +RVT +RCT +RRT +RIT ACPRR ZSBIS END diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl index 0517e0dfdf7793ef78453a7b5e8add67ae19a826..38e12ff80915470a1caa2b445ec7256769bf4042 100644 --- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl +++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl @@ -43,8 +43,8 @@ begin ;=================================================; ; Read water wapor mixing ratio if it exists ; ----------------------------------------- - if(isfilevar(a,"RVM")) - qv = a->RVM(0,:,:,:) + if(isfilevar(a,"RVT")) + qv = a->RVT(0,:,:,:) qv = qv*1000. qv@units = "g/kg" qv@description="Water vapor mixing ratio" @@ -70,8 +70,8 @@ begin ; Read ice mixing ratio if it exists ; ----------------------------------------- - if(isfilevar(a,"RIM")) - qi = a->RIM(0,:,:,:) + if(isfilevar(a,"RIT")) + qi = a->RIT(0,:,:,:) qi = qi*1000. qi@units = "g/kg" qi@description="Ice mixing ratio" diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl index 90f709bdbbd10f12b0fc195c5ff912b4780af95f..2667ac88b12016b8865bc56de197dc0f428166f0 100644 --- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl +++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl @@ -34,17 +34,17 @@ begin ;=================================================; ; Read the variables we need ;=================================================; - tk = a->THM(0,:,:,:) ; potential temperature + tk = a->THT(0,:,:,:) ; potential temperature tk@description="Potential Temperature" - u = a->UM(0,:,:,:) ; u + u = a->UT(0,:,:,:) ; u - v = a->VM(0,:,:,:) ; v + v = a->VT(0,:,:,:) ; v - p = a->PABSM(0,:,:,:) ; pressure + p = a->PABST(0,:,:,:) ; pressure p@description="Pressure" - rv = a->RVM(0,:,:,:)*1000 ; water vapor mixing ratio + rv = a->RVT(0,:,:,:)*1000 ; water vapor mixing ratio rv@description="Water vapor mixing ratio" ZS = a->ZSBIS(0,:,:) ; terrain