Commit 29ba5082 authored by Ronan Trépos's avatar Ronan Trépos
Browse files

ContrainteEau: constrain STAOCKC* not to be zero

parent fcb9e453
......@@ -260,7 +260,7 @@ public :
double STOCKC1_tmp = 0.0;
STOCKC1_tmp = ( ( ps.Hpf_C1 / 100.0 ) * ps.da_C1 * psi.zC1 ) + C1();
if ( STOCKC1_tmp < 0.0 ){
if ( STOCKC1_tmp <= 0.0 ){
STOCKC1_tmp = 0.1;
}
STOCKC1 = STOCKC1_tmp;
......@@ -269,7 +269,7 @@ public :
double STOCKC2_tmp = 0.0;
STOCKC2_tmp = ( ( ps.Hpf_C2 / 100.0 ) * ps.da_C2 * zC2() ) + C2();
if ( STOCKC2_tmp < 0.0 ){
if ( STOCKC2_tmp <= 0.0 ){
STOCKC2_tmp = 0.1;
}
STOCKC2 = STOCKC2_tmp;
......@@ -277,7 +277,7 @@ public :
{
double STOCKC3_tmp = 0.0;
STOCKC3_tmp = ( ( ps.Hpf_C2 / 100.0 ) * ps.da_C2 * zC3() ) + C3();
if ( STOCKC3_tmp < 0.0 ){
if ( STOCKC3_tmp <= 0.0 ){
STOCKC3_tmp = 0.1;
}
STOCKC3 = STOCKC3_tmp;
......@@ -495,7 +495,7 @@ public :
double STOCKC1_tmp = 0.0;
STOCKC1_tmp = ( ( ps.Hpf_C1 / 100.0 ) * ps.da_C1 * psi.zC1 ) + C1();
if ( STOCKC1_tmp < 0.0 ){
if ( STOCKC1_tmp <= 0.0 ){
STOCKC1_tmp = 0.1;
}
......@@ -506,7 +506,7 @@ public :
double STOCKC2_tmp = 0.0;
STOCKC2_tmp = ( ( ps.Hpf_C2 / 100.0 ) * ps.da_C2 * zC2() ) + C2();
if ( STOCKC2_tmp < 0.0 ){
if ( STOCKC2_tmp <= 0.0 ){
STOCKC2_tmp = 0.1;
}
STOCKC2 = STOCKC2_tmp;
......@@ -516,7 +516,7 @@ public :
double STOCKC3_tmp = 0.0;
STOCKC3_tmp = ( ( ps.Hpf_C2 / 100.0 ) * ps.da_C2 * zC3() ) + C3();
if ( STOCKC3_tmp < 0.0 ){
if ( STOCKC3_tmp <= 0.0 ){
STOCKC3_tmp = 0.1;
}
STOCKC3 = STOCKC3_tmp;
......
Supports Markdown
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