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

Maud et Gaelle 26/04/2018 : TUV 531

parent 653fc5b5
No related branches found
Tags PACK-MNH-V5-4--0-12_CHIMIE
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
MY_RUN/KTEST/009_ICARTT/001_pgd1/get_chimie_files
+ 1
− 1
View file @ 5780a42a
...@@ -20,7 +20,7 @@ cd ${CHIMIE_FILES} ...@@ -20,7 +20,7 @@ cd ${CHIMIE_FILES}
CHIMIE_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_CHIMIE" CHIMIE_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_CHIMIE"
WGET="wget" WGET="wget"
export CHIMIE_URL export CHIMIE_URL
for dir in EMISSIONS tuv50 for dir in EMISSIONS tuv531
do do
[ -d $dir ] || ( ${WGET} -c -nd $CHIMIE_URL/$dir.tar.gz ; tar xvfz $dir.tar.gz ; rm -f $dir.tar.gz ) [ -d $dir ] || ( ${WGET} -c -nd $CHIMIE_URL/$dir.tar.gz ; tar xvfz $dir.tar.gz ; rm -f $dir.tar.gz )
done done
......
MY_RUN/KTEST/009_ICARTT/003_mesonh/run_mesonh_xyz
+ 2
− 3
View file @ 5780a42a
...@@ -9,9 +9,8 @@ rm -f ICART* OUT* EXSEG?.nam ...@@ -9,9 +9,8 @@ rm -f ICART* OUT* EXSEG?.nam
export CHIMIE_FILES=${CHIMIE_FILES:-"$HOME/CHIMIE_FILES"} export CHIMIE_FILES=${CHIMIE_FILES:-"$HOME/CHIMIE_FILES"}
ln -sf ${CHIMIE_FILES}/tuv50/DATAE1 . ln -sf ${CHIMIE_FILES}/tuv531/DATAE1 .
ln -sf ${CHIMIE_FILES}/tuv50/DATAJ1 . ln -sf ${CHIMIE_FILES}/tuv531/DATAJ1 .
ln -sf ${CHIMIE_FILES}/tuv50/DATAS1 .
cp EXSEG1.nam.src EXSEG1.nam cp EXSEG1.nam.src EXSEG1.nam
ln -sf ../002_arp2lfi/CPLCH20040810.18.{des,lfi,nc} . ln -sf ../002_arp2lfi/CPLCH20040810.18.{des,lfi,nc} .
ln -sf ../001_pgd1/ICARTT1008_PGD_15km.{des,lfi,nc} . ln -sf ../001_pgd1/ICARTT1008_PGD_15km.{des,lfi,nc} .
......
MY_RUN/KTEST/011_KW78CHEM/002_mesonh/run_mesonh_xyz
+ 2
− 3
View file @ 5780a42a
...@@ -5,9 +5,8 @@ ...@@ -5,9 +5,8 @@
set -x set -x
set -e set -e
export CHIMIE_FILES=${CHIMIE_FILES:-"$HOME/CHIMIE_FILES"} export CHIMIE_FILES=${CHIMIE_FILES:-"$HOME/CHIMIE_FILES"}
ln -sf ${CHIMIE_FILES}/tuv50/DATAE1 . ln -sf ${CHIMIE_FILES}/tuv531/DATAE1 .
ln -sf ${CHIMIE_FILES}/tuv50/DATAJ1 . ln -sf ${CHIMIE_FILES}/tuv531/DATAJ1 .
ln -sf ${CHIMIE_FILES}/tuv50/DATAS1 .
ln -fs ../001_prep_ideal_case/KWRAIN.{des,lfi,nc} . ln -fs ../001_prep_ideal_case/KWRAIN.{des,lfi,nc} .
rm -f KWRAI.1.SEGCH OUT* rm -f KWRAI.1.SEGCH OUT*
time ${MPIRUN} MESONH${XYZ} time ${MPIRUN} MESONH${XYZ}
src/MNH/BASIC.f90
+ 46
− 43
View file @ 5780a42a
!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.
! !
!======================================================================== !========================================================================
! !
...@@ -250,7 +246,7 @@ ...@@ -250,7 +246,7 @@
! KTR38=@KT(0.03,60.,TPK%T,TPK%RADR,KVECNPT)/(@HEFFA(4.1e3,-6200.,1.74e-5,0.,0.,0.,TPK%PHR,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_ORA2-->ORA2 ! KTR38=@KT(0.03,60.,TPK%T,TPK%RADR,KVECNPT)/(@HEFFA(4.1e3,-6200.,1.74e-5,0.,0.,0.,TPK%PHR,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_ORA2-->ORA2
! KTR39=@KT(0.05,47.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(2.7e0,-2030.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_MO2-->MO2 ! KTR39=@KT(0.05,47.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(2.7e0,-2030.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_MO2-->MO2
! KTR40=@KT(0.007,48.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(3.e2,-5280.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_OP1-->OP1 ! KTR40=@KT(0.007,48.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(3.e2,-5280.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_OP1-->OP1
! KC1=7.e-6::WC_H2O2-->WC_OH+WC_OH ! KC1=!ZRATES(:,018)::WC_H2O2-->WC_OH+WC_OH
! KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD::WC_OH+WC_OH-->WC_H2O2 ! KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD::WC_OH+WC_OH-->WC_H2O2
! KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e-5)/(1.6e-5+10.**(-TPK%PHC)))/TPK%MOL2MOLECCLOUD::WC_OH+WC_HO2--> ! KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e-5)/(1.6e-5+10.**(-TPK%PHC)))/TPK%MOL2MOLECCLOUD::WC_OH+WC_HO2-->
! KC4=(3.2E+7*exp(-1700.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD::WC_H2O2+WC_OH-->WC_HO2 ! KC4=(3.2E+7*exp(-1700.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD::WC_H2O2+WC_OH-->WC_HO2
...@@ -262,7 +258,7 @@ ...@@ -262,7 +258,7 @@
! KC10=2.6E-2*10.**(-TPK%PHC)/(1.26e-6+10.**(-TPK%PHC))::WC_HNO4-->WC_HO2+WC_NO2 ! KC10=2.6E-2*10.**(-TPK%PHC)/(1.26e-6+10.**(-TPK%PHC))::WC_HNO4-->WC_HO2+WC_NO2
! KC11=1.1E0*1.26e-6/(1.26e-6+10.**(-TPK%PHC))::WC_HNO4-->WC_HONO ! KC11=1.1E0*1.26e-6/(1.26e-6+10.**(-TPK%PHC))::WC_HNO4-->WC_HONO
! KC12=(3.3E+5*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*(10.**(-TPK%PHC))**2./((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD::WC_HNO4+WC_SO2-->WC_SULF+WC_HNO3 ! KC12=(3.3E+5*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*(10.**(-TPK%PHC))**2./((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD::WC_HNO4+WC_SO2-->WC_SULF+WC_HNO3
! KC13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHC))::WC_HNO3-->WC_NO2+WC_OH ! KC13=!ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHC))::WC_HNO3-->WC_NO2+WC_OH
! KC14=1.0E+10::WC_N2O5-->WC_HNO3+WC_HNO3 ! KC14=1.0E+10::WC_N2O5-->WC_HNO3+WC_HNO3
! KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SULF-->WC_HNO3+WC_ASO4 ! KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SULF-->WC_HNO3+WC_ASO4
! KC16=(1.3E+9*exp(-2200.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SO2-->WC_HNO3+WC_ASO3 ! KC16=(1.3E+9*exp(-2200.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SO2-->WC_HNO3+WC_ASO3
...@@ -280,7 +276,7 @@ ...@@ -280,7 +276,7 @@
! KC28=4.6E+2*exp(-1100.*(1./TPK%T-1./298.15))::WC_ASO4-->WC_SULF+WC_OH ! KC28=4.6E+2*exp(-1100.*(1./TPK%T-1./298.15))::WC_ASO4-->WC_SULF+WC_OH
! KC29=((3.7E+5*exp(-5500.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+1.5E+9*exp(-5300.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)))/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_SO2+WC_O3-->WC_SULF ! KC29=((3.7E+5*exp(-5500.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+1.5E+9*exp(-5300.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)))/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_SO2+WC_O3-->WC_SULF
! KC30=(9.1E+7*exp(-3600.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)*10.**(-TPK%PHC)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_SO2+WC_H2O2-->WC_SULF ! KC30=(9.1E+7*exp(-3600.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)*10.**(-TPK%PHC)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_SO2+WC_H2O2-->WC_SULF
! KR1=7.e-6::WR_H2O2-->WR_OH+WR_OH ! KR1=!ZRATES(:,018)::WR_H2O2-->WR_OH+WR_OH
! KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN::WR_OH+WR_OH-->WR_H2O2 ! KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN::WR_OH+WR_OH-->WR_H2O2
! KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e-5)/(1.6e-5+10.**(-TPK%PHR)))/TPK%MOL2MOLECRAIN::WR_OH+WR_HO2--> ! KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e-5)/(1.6e-5+10.**(-TPK%PHR)))/TPK%MOL2MOLECRAIN::WR_OH+WR_HO2-->
! KR4=(3.2E+7*exp(-1700.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN::WR_H2O2+WR_OH-->WR_HO2 ! KR4=(3.2E+7*exp(-1700.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN::WR_H2O2+WR_OH-->WR_HO2
...@@ -292,7 +288,7 @@ ...@@ -292,7 +288,7 @@
! KR10=2.6E-2*10.**(-TPK%PHR)/(1.26e-6+10.**(-TPK%PHR))::WR_HNO4-->WR_HO2+WR_NO2 ! KR10=2.6E-2*10.**(-TPK%PHR)/(1.26e-6+10.**(-TPK%PHR))::WR_HNO4-->WR_HO2+WR_NO2
! KR11=1.1E0*1.26e-6/(1.26e-6+10.**(-TPK%PHR))::WR_HNO4-->WR_HONO ! KR11=1.1E0*1.26e-6/(1.26e-6+10.**(-TPK%PHR))::WR_HNO4-->WR_HONO
! KR12=(3.3E+5*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*(10.**(-TPK%PHR))**2./((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN::WR_HNO4+WR_SO2-->WR_SULF+WR_HNO3 ! KR12=(3.3E+5*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*(10.**(-TPK%PHR))**2./((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN::WR_HNO4+WR_SO2-->WR_SULF+WR_HNO3
! KR13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHR))::WR_HNO3-->WR_NO2+WR_OH ! KR13=!ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHR))::WR_HNO3-->WR_NO2+WR_OH
! KR14=1.0E+10::WR_N2O5-->WR_HNO3+WR_HNO3 ! KR14=1.0E+10::WR_N2O5-->WR_HNO3+WR_HNO3
! KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SULF-->WR_HNO3+WR_ASO4 ! KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SULF-->WR_HNO3+WR_ASO4
! KR16=(1.3E+9*exp(-2200.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SO2-->WR_HNO3+WR_ASO3 ! KR16=(1.3E+9*exp(-2200.*(1./TPK%T-1./298.15))*1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)/(1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SO2-->WR_HNO3+WR_ASO3
...@@ -3572,7 +3568,7 @@ IF (GFIRSTCALL) THEN ...@@ -3572,7 +3568,7 @@ IF (GFIRSTCALL) THEN
&030.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_MO2-->MO2' &030.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_MO2-->MO2'
CFULLREACS(212) = 'KTR40=@KT(0.007,48.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(3.e2,-5& CFULLREACS(212) = 'KTR40=@KT(0.007,48.,TPK%T,TPK%RADR,KVECNPT)/(@HENRY(3.e2,-5&
&280.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_OP1-->OP1' &280.,TPK%T,KVECNPT)*TPK%RCH*TPK%T)::WR_OP1-->OP1'
CFULLREACS(213) = 'KC1=7.e-6::WC_H2O2-->WC_OH+WC_OH' CFULLREACS(213) = 'KC1=!ZRATES(:,018)::WC_H2O2-->WC_OH+WC_OH'
CFULLREACS(214) = 'KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECC& CFULLREACS(214) = 'KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECC&
&LOUD::WC_OH+WC_OH-->WC_H2O2' &LOUD::WC_OH+WC_OH-->WC_H2O2'
CFULLREACS(215) = 'KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1& CFULLREACS(215) = 'KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1&
...@@ -3602,8 +3598,8 @@ IF (GFIRSTCALL) THEN ...@@ -3602,8 +3598,8 @@ IF (GFIRSTCALL) THEN
&PK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**& &PK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**&
&(-TPK%PHC))**2.)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD::WC_HNO4+WC_SO& &(-TPK%PHC))**2.)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD::WC_HNO4+WC_SO&
&2-->WC_SULF+WC_HNO3' &2-->WC_SULF+WC_HNO3'
CFULLREACS(225) = 'KC13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHC))::WC_HNO3-->WC_NO& CFULLREACS(225) = 'KC13=!ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHC))::WC_HNO&
&2+WC_OH' &3-->WC_NO2+WC_OH'
CFULLREACS(226) = 'KC14=1.0E+10::WC_N2O5-->WC_HNO3+WC_HNO3' CFULLREACS(226) = 'KC14=1.0E+10::WC_N2O5-->WC_HNO3+WC_HNO3'
CFULLREACS(227) = 'KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK& CFULLREACS(227) = 'KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK&
&%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SULF-->WC_HNO3+WC_& &%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD::WC_NO3+WC_SULF-->WC_HNO3+WC_&
...@@ -3659,7 +3655,7 @@ IF (GFIRSTCALL) THEN ...@@ -3659,7 +3655,7 @@ IF (GFIRSTCALL) THEN
&./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1& &./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15))+1.3e-2*exp(1965.*(1&
&./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD& &./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.))/TPK%MOL2MOLECCLOUD&
&::WC_SO2+WC_H2O2-->WC_SULF' &::WC_SO2+WC_H2O2-->WC_SULF'
CFULLREACS(243) = 'KR1=7.e-6::WR_H2O2-->WR_OH+WR_OH' CFULLREACS(243) = 'KR1=!ZRATES(:,018)::WR_H2O2-->WR_OH+WR_OH'
CFULLREACS(244) = 'KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECR& CFULLREACS(244) = 'KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECR&
&AIN::WR_OH+WR_OH-->WR_H2O2' &AIN::WR_OH+WR_OH-->WR_H2O2'
CFULLREACS(245) = 'KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1& CFULLREACS(245) = 'KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1&
...@@ -3689,8 +3685,8 @@ IF (GFIRSTCALL) THEN ...@@ -3689,8 +3685,8 @@ IF (GFIRSTCALL) THEN
&PK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**& &PK%T-1./298.15))+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**&
&(-TPK%PHR))**2.)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN::WR_HNO4+WR_SO2& &(-TPK%PHR))**2.)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN::WR_HNO4+WR_SO2&
&-->WR_SULF+WR_HNO3' &-->WR_SULF+WR_HNO3'
CFULLREACS(255) = 'KR13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHR))::WR_HNO3-->WR_NO& CFULLREACS(255) = 'KR13=!ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHR))::WR_HNO&
&2+WR_OH' &3-->WR_NO2+WR_OH'
CFULLREACS(256) = 'KR14=1.0E+10::WR_N2O5-->WR_HNO3+WR_HNO3' CFULLREACS(256) = 'KR14=1.0E+10::WR_N2O5-->WR_HNO3+WR_HNO3'
CFULLREACS(257) = 'KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK& CFULLREACS(257) = 'KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK&
&%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SULF-->WR_HNO3+WR_A& &%PHR)+(10.**(-TPK%PHR))**2.))/TPK%MOL2MOLECRAIN::WR_NO3+WR_SULF-->WR_HNO3+WR_A&
...@@ -28919,7 +28915,6 @@ SUBROUTINE SUBSRW8 ...@@ -28919,7 +28915,6 @@ SUBROUTINE SUBSRW8
! !
!Indices 213 a 222 !Indices 213 a 222
! !
TPK%KC1=7.e-6
TPK%KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD TPK%KC2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECCLOUD
TPK%KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e& TPK%KC3=((2.8E+10*10.**(-TPK%PHC)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e&
&-5)/(1.6e-5+10.**(-TPK%PHC)))/TPK%MOL2MOLECCLOUD &-5)/(1.6e-5+10.**(-TPK%PHC)))/TPK%MOL2MOLECCLOUD
...@@ -28951,7 +28946,6 @@ SUBROUTINE SUBSRW9 ...@@ -28951,7 +28946,6 @@ SUBROUTINE SUBSRW9
&((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)& &((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)&
&)+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.& &)+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHC)+(10.**(-TPK%PHC))**2.&
&)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD &)*(1.26e-6+10.**(-TPK%PHC))))/TPK%MOL2MOLECCLOUD
TPK%KC13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHC))
TPK%KC14=1.0E+10 TPK%KC14=1.0E+10
TPK%KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHC)+(10.**(-T& TPK%KC15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHC)+(10.**(-T&
&PK%PHC))**2.))/TPK%MOL2MOLECCLOUD &PK%PHC))**2.))/TPK%MOL2MOLECCLOUD
...@@ -29012,7 +29006,6 @@ SUBROUTINE SUBSRW11 ...@@ -29012,7 +29006,6 @@ SUBROUTINE SUBSRW11
! !
!Indices 243 a 252 !Indices 243 a 252
! !
TPK%KR1=7.e-6
TPK%KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN TPK%KR2=(3.6E+9*exp(-930.*(1./TPK%T-1./298.15)))/TPK%MOL2MOLECRAIN
TPK%KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e& TPK%KR3=((2.8E+10*10.**(-TPK%PHR)+3.5E+10*exp(-720.*(1./TPK%T-1./298.15))*1.6e&
&-5)/(1.6e-5+10.**(-TPK%PHR)))/TPK%MOL2MOLECRAIN &-5)/(1.6e-5+10.**(-TPK%PHR)))/TPK%MOL2MOLECRAIN
...@@ -29044,7 +29037,6 @@ SUBROUTINE SUBSRW12 ...@@ -29044,7 +29037,6 @@ SUBROUTINE SUBSRW12
&((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)& &((1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*6.4e-8*exp(1430.*(1./TPK%T-1./298.15)&
&)+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.& &)+1.3e-2*exp(1965.*(1./TPK%T-1./298.15))*10.**(-TPK%PHR)+(10.**(-TPK%PHR))**2.&
&)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN &)*(1.26e-6+10.**(-TPK%PHR))))/TPK%MOL2MOLECRAIN
TPK%KR13=4.e-7*2.2e+1/(2.2e+1+10.**(-TPK%PHR))
TPK%KR14=1.0E+10 TPK%KR14=1.0E+10
TPK%KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHR)+(10.**(-T& TPK%KR15=(1.0E+5*1.0e+3*1.0e-2/(1.0e+3*1.0e-2+1.0e+3*10.**(-TPK%PHR)+(10.**(-T&
&PK%PHR))**2.))/TPK%MOL2MOLECRAIN &PK%PHR))**2.))/TPK%MOL2MOLECRAIN
...@@ -29193,8 +29185,8 @@ REAL,DIMENSION(:,:,:,:), INTENT(IN) :: PJVALUES ! Tuv coefficient ...@@ -29193,8 +29185,8 @@ REAL,DIMENSION(:,:,:,:), INTENT(IN) :: PJVALUES ! Tuv coefficient
! /BEGIN_SET_PHOTO_RATES/ ! /BEGIN_SET_PHOTO_RATES/
! parameter for use by subroutine JVALUES, ! parameter for use by subroutine JVALUES,
! contains the actual photolysis rates ! contains the actual photolysis rates
REAL, DIMENSION(KVECNPT,21) :: ZRATESIO ! TUV photolysis rates at one level REAL, DIMENSION(KVECNPT,41) :: ZRATESIO ! TUV photolysis rates at one level
REAL, DIMENSION(KVECNPT,17) :: ZRATES ! photolysis rates of RACM (vector) REAL, DIMENSION(KVECNPT,19) :: ZRATES ! photolysis rates of RACM (vector)
INTEGER :: JITPK ! loop counter for J-Value transfer INTEGER :: JITPK ! loop counter for J-Value transfer
INTEGER :: IDTI,IDTJ INTEGER :: IDTI,IDTJ
INTEGER :: JITPKPLUS INTEGER :: JITPKPLUS
...@@ -29229,29 +29221,34 @@ DO JITPK = 0, KVECNPT-1 ...@@ -29229,29 +29221,34 @@ DO JITPK = 0, KVECNPT-1
! !
! associate TUV J-Values to ReLACS J-Values ! associate TUV J-Values to ReLACS J-Values
! !
ZRATES(JITPK+1, 1) = ZRATESIO(JITPK+1,2) ! change according to original coefficients and modified RACM
ZRATES(JITPK+1, 2) = ZRATESIO(JITPK+1,3) !
ZRATES(JITPK+1, 3) = ZRATESIO(JITPK+1,4) ZRATES(JITPK+1, 1) = ZRATESIO(JITPK+1,5)
ZRATES(JITPK+1, 4) = ZRATESIO(JITPK+1,9) ZRATES(JITPK+1, 2) = ZRATESIO(JITPK+1,2)
ZRATES(JITPK+1, 5) = ZRATESIO(JITPK+1,10) ZRATES(JITPK+1, 3) = ZRATESIO(JITPK+1,3)
ZRATES(JITPK+1, 6) = ZRATESIO(JITPK+1,11) ZRATES(JITPK+1, 4) = ZRATESIO(JITPK+1,8)
ZRATES(JITPK+1, 7) = ZRATESIO(JITPK+1,5) ZRATES(JITPK+1, 5) = ZRATESIO(JITPK+1,9)
ZRATES(JITPK+1, 8) = ZRATESIO(JITPK+1,6) ZRATES(JITPK+1, 6) = ZRATESIO(JITPK+1,10)
ZRATES(JITPK+1, 9) = ZRATESIO(JITPK+1,12) ZRATES(JITPK+1, 7) = ZRATESIO(JITPK+1,6)
ZRATES(JITPK+1, 10) = ZRATESIO(JITPK+1,14) ZRATES(JITPK+1, 8) = ZRATESIO(JITPK+1,7)
ZRATES(JITPK+1, 11) = ZRATESIO(JITPK+1,13) ZRATES(JITPK+1, 9) = ZRATESIO(JITPK+1,4)
ZRATES(JITPK+1, 12) = ZRATESIO(JITPK+1,17) ZRATES(JITPK+1, 10) = ZRATESIO(JITPK+1,12)
ZRATES(JITPK+1, 13) = ZRATESIO(JITPK+1,15) ZRATES(JITPK+1, 11) = ZRATESIO(JITPK+1,11)
ZRATES(JITPK+1, 14) = 0.962055 *ZRATESIO(JITPK+1,15)+& ZRATES(JITPK+1, 12) = ZRATESIO(JITPK+1,13)
& 1.06247E-02 *ZRATESIO(JITPK+1,12) ZRATES(JITPK+1, 13) = ZRATESIO(JITPK+1,17)
ZRATES(JITPK+1, 15) = ZRATESIO(JITPK+1,20) ZRATES(JITPK+1, 14) = 0.962055 *ZRATESIO(JITPK+1,17)+&
ZRATES(JITPK+1, 16) = 3.16657 *ZRATESIO(JITPK+1,15)& & 3.79454E-02 *ZRATESIO(JITPK+1,38)
&+ 0.372446 *ZRATESIO(JITPK+1,15)& ZRATES(JITPK+1, 15) = ZRATESIO(JITPK+1,33)
&+ 8.42257 *ZRATESIO(JITPK+1,15)& ZRATES(JITPK+1, 16) = 0.20842 *ZRATESIO(JITPK+1,35)&
&+ 207.5913 *ZRATESIO(JITPK+1,20)& &+ 6.43207E-02 *ZRATESIO(JITPK+1,36)&
&+ 0.0 *ZRATESIO(JITPK+1,20)& &+ 3.10372E-02 *ZRATESIO(JITPK+1,34)&
&+ 8.44837E-04 *ZRATESIO(JITPK+1,20) &+ 0.376 *ZRATESIO(JITPK+1,37)&
ZRATES(JITPK+1, 17) = ZRATESIO(JITPK+1,16) &+ 0.31937 *ZRATESIO(JITPK+1,26)
ZRATES(JITPK+1, 17) = 0.20*ZRATESIO(JITPK+1,20)&
&+ 0.80*ZRATESIO(JITPK+1,21)
! aqueous phase photolysis
ZRATES(JITPK+1, 18) = ZRATESIO(JITPK+1,40)
ZRATES(JITPK+1, 19) = ZRATESIO(JITPK+1,41)
! !
END DO END DO
! !
...@@ -29273,6 +29270,12 @@ END DO ...@@ -29273,6 +29270,12 @@ END DO
TPK%K015=ZRATES(:,015) TPK%K015=ZRATES(:,015)
TPK%K016=ZRATES(:,016) TPK%K016=ZRATES(:,016)
TPK%K017=ZRATES(:,017) TPK%K017=ZRATES(:,017)
IF (TPK%LUSECHAQ) THEN
TPK%KC1=ZRATES(:,018)
TPK%KC13=ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHC))
TPK%KR1=ZRATES(:,018)
TPK%KR13=ZRATES(:,019)*2.2e+1/(2.2e+1+10.**(-TPK%PHR))
END IF
TPK%NOUT = KOUT TPK%NOUT = KOUT
TPK%NVERB = KVERB TPK%NVERB = KVERB
RETURN RETURN
......
src/MNH/ch_f77.fx90
+ 8590
− 3817
View file @ 5780a42a
This diff is collapsed.
src/MNH/ch_init_jvalues.f90
+ 3
− 3
View file @ 5780a42a
...@@ -72,7 +72,7 @@ END MODULE MODI_CH_INIT_JVALUES ...@@ -72,7 +72,7 @@ END MODULE MODI_CH_INIT_JVALUES
!! !!
!! EXTERNAL !! EXTERNAL
!! -------- !! --------
!! TUV39.f (Fortran 77 code from S. Madronich) !! TUV version 5.3.1 (Fortran 77 code from S. Madronich)
!! !!
!! IMPLICIT ARGUMENTS !! IMPLICIT ARGUMENTS
!! ------------------ !! ------------------
...@@ -148,8 +148,8 @@ IF (.NOT.ALLOCATED(XJDATA)) ALLOCATE(XJDATA(NSZA_INCR,NZZ_JVAL,JPJVMAX,NBALB)) ...@@ -148,8 +148,8 @@ IF (.NOT.ALLOCATED(XJDATA)) ALLOCATE(XJDATA(NSZA_INCR,NZZ_JVAL,JPJVMAX,NBALB))
ZLWC(JKLEV)= 0.0 ZLWC(JKLEV)= 0.0
ENDDO ENDDO
! !
!* CALL TUV 5.0 !* CALL TUV 5.3.1
! ------------ ! --------------
! !
DO JALB=1,NBALB DO JALB=1,NBALB
ZALBLOOP=0.02+0.20*FLOAT(JALB-1)/FLOAT(NBALB-1) ZALBLOOP=0.02+0.20*FLOAT(JALB-1)/FLOAT(NBALB-1)
......
src/MNH/modd_ch_init_jvalues.f90
+ 1
− 1
View file @ 5780a42a
...@@ -33,7 +33,7 @@ INTEGER :: NSZA_INCR = 99 + 1 ...@@ -33,7 +33,7 @@ INTEGER :: NSZA_INCR = 99 + 1
REAL, ALLOCATABLE, DIMENSION(:) :: XSZA_JVAL REAL, ALLOCATABLE, DIMENSION(:) :: XSZA_JVAL
INTEGER, PARAMETER :: NZZ_JVAL = 30 + 1 INTEGER, PARAMETER :: NZZ_JVAL = 30 + 1
REAL, ALLOCATABLE, DIMENSION(:) :: XZZ_JVAL REAL, ALLOCATABLE, DIMENSION(:) :: XZZ_JVAL
INTEGER, PARAMETER :: JPJVMAX = 21 INTEGER, PARAMETER :: JPJVMAX = 41
INTEGER :: NBALB = 10 INTEGER :: NBALB = 10
! !
END MODULE MODD_CH_INIT_JVALUES END MODULE MODD_CH_INIT_JVALUES
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