Skip to content
Snippets Groups Projects
Commit 4e66f31c authored by Gaelle Tanguy's avatar Gaelle Tanguy Committed by WAUTELET Philippe
Browse files

M.Leriche 2014 : correction décalage d’indice dans le schéma réactionnel ReLACS-AQ

parent 35769266
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 25 additions and 24 deletions
src/MNH/BASIC.f90
+ 25
− 24
View file @ 4e66f31c
......@@ -3,6 +3,7 @@
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!
!
!========================================================================
!
! The following species appear on the left hand side of a reaction,
......@@ -4463,11 +4464,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
CALL SUB0
CALL SUB1
......@@ -5576,11 +5577,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
CALL SUB0
CALL SUB1
......@@ -6479,11 +6480,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
PJAC(:,:,:) = 0.0
CALL SUBJ0
......@@ -24447,11 +24448,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
PJAC(:,:,:) = 0.0
CALL SUBJ0
......@@ -30261,11 +30262,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
PTERMS(:,:,:) = 0.0
CALL SUBT0
......@@ -33612,11 +33613,11 @@ TPK%O1D(:)=(TPK%K002(:)*PCONC(:,JP_O3))/(TPK%K020(:)*TPK%N2(:)+TPK%K021(:)*TPK%O
&TPK%K022(:)*TPK%H2O(:))
TPK%O3P(:)=(TPK%K001(:)*PCONC(:,JP_NO2)+TPK%K003(:)*PCONC(:,JP_O3)+TPK%K008(:)*PCONC(:,JP_NO3)+&
&TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
&0.00000*TPK%K076(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K077(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
&0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
&(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K053(:)*PCONC(:,JP_BIO)+&
&TPK%K054(:)*PCONC(:,JP_CARBO))
&TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
&TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
PTERMS(:,:,:) = 0.0
CALL SUBT0
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment

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