[GRASS-SVN] r31580 - grass-addons/gipe/r.soiltex2prop
svn_grass at osgeo.org
svn_grass at osgeo.org
Fri May 30 05:33:34 EDT 2008
Author: ychemin
Date: 2008-05-30 05:33:34 -0400 (Fri, 30 May 2008)
New Revision: 31580
Modified:
grass-addons/gipe/r.soiltex2prop/prct2hf.c
grass-addons/gipe/r.soiltex2prop/prct2ksat.c
Log:
bug fixing
Modified: grass-addons/gipe/r.soiltex2prop/prct2hf.c
===================================================================
--- grass-addons/gipe/r.soiltex2prop/prct2hf.c 2008-05-30 08:34:22 UTC (rev 31579)
+++ grass-addons/gipe/r.soiltex2prop/prct2hf.c 2008-05-30 09:33:34 UTC (rev 31580)
@@ -16,8 +16,6 @@
double temp, hf;
double silt_input=0.0; //Rawls et al (1990)
//do not have silt input
- // set up mark index for inside/outside polygon check
- double mark[POLYGON_DIMENSION]={0.0};
//printf("in prct2hf(), cm\n");
//setup the 3Dvectors and initialize them
struct vector cls[POLYGON_DIMENSION] = {0.0};
@@ -27,9 +25,6 @@
cls[i].clay=0.0;
cls[i].silt=0.0;
}
- //transform input from [0,1] to [0,100]
- sand_input *= 100.0;
- clay_input *= 100.0;
//fill up initial polygon points
cls[0].sand=0.0;
cls[0].clay=100.0;
@@ -38,8 +33,8 @@
cls[2].sand=17.0;
cls[2].clay=55.0;
//Get started
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
hf=175.0;
index=1;
//printf("hf=175.0\n");
@@ -51,8 +46,8 @@
cls[1].clay=55.0;
cls[2].sand=30.0;
cls[2].clay=60.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=175.0\n");
hf=175.0;
@@ -65,8 +60,8 @@
cls[1].clay=66.0;
cls[2].sand=30.0;
cls[2].clay=60.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=175.0\n");
hf=175.0;
@@ -79,8 +74,8 @@
cls[1].clay=0.0;
cls[2].sand=65.0;
cls[2].clay=15.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=5.0\n");
hf=5.0;
@@ -93,8 +88,8 @@
cls[1].clay=30.0;
cls[2].sand=65.0;
cls[2].clay=15.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=5.0\n");
hf=5.0;
@@ -107,8 +102,8 @@
cls[1].clay=30.0;
cls[2].sand=67.0;
cls[2].clay=33.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=5.0\n");
hf=5.0;
@@ -121,8 +116,8 @@
cls[1].clay=30.0;
cls[2].sand=67.0;
cls[2].clay=33.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=15.0\n");
hf=15.0;
@@ -135,8 +130,8 @@
cls[1].clay=30.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=15.0\n");
hf=15.0;
@@ -149,8 +144,8 @@
cls[1].clay=30.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=15.0\n");
hf=15.0;
@@ -163,8 +158,8 @@
cls[1].clay=0.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=15.0\n");
hf=15.0;
@@ -177,8 +172,8 @@
cls[1].clay=14.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=25.0\n");
hf=25.0;
@@ -191,8 +186,8 @@
cls[1].clay=14.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=25.0\n");
hf=25.0;
@@ -205,8 +200,8 @@
cls[1].clay=42.0;
cls[2].sand=30.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=25.0\n");
hf=25.0;
@@ -219,8 +214,8 @@
cls[1].clay=42.0;
cls[2].sand=55.0;
cls[2].clay=45.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=25.0\n");
hf=25.0;
@@ -233,8 +228,8 @@
cls[1].clay=50.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -247,8 +242,8 @@
cls[1].clay=24.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -261,8 +256,8 @@
cls[1].clay=24.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -275,8 +270,8 @@
cls[1].clay=24.0;
cls[2].sand=20.0;
cls[2].clay=14.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -289,8 +284,8 @@
cls[1].clay=14.0;
cls[2].sand=20.0;
cls[2].clay=14.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -303,8 +298,8 @@
cls[1].clay=14.0;
cls[2].sand=20.0;
cls[2].clay=14.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -317,8 +312,8 @@
cls[1].clay=14.0;
cls[2].sand=7.0;
cls[2].clay=3.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -331,8 +326,8 @@
cls[1].clay=0.0;
cls[2].sand=7.0;
cls[2].clay=3.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -345,8 +340,8 @@
cls[1].clay=0.0;
cls[2].sand=7.0;
cls[2].clay=3.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -359,8 +354,8 @@
cls[1].clay=0.0;
cls[2].sand=0.0;
cls[2].clay=9.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=35.0\n");
hf=35.0;
@@ -373,8 +368,8 @@
cls[1].clay=28.0;
cls[2].sand=0.0;
cls[2].clay=9.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -387,8 +382,8 @@
cls[1].clay=28.0;
cls[2].sand=7.0;
cls[2].clay=3.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -401,8 +396,8 @@
cls[1].clay=28.0;
cls[2].sand=7.0;
cls[2].clay=3.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -415,8 +410,8 @@
cls[1].clay=28.0;
cls[2].sand=10.0;
cls[2].clay=27.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -429,8 +424,8 @@
cls[1].clay=28.0;
cls[2].sand=10.0;
cls[2].clay=27.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -443,8 +438,8 @@
cls[1].clay=28.0;
cls[2].sand=10.0;
cls[2].clay=27.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -457,8 +452,8 @@
cls[1].clay=28.0;
cls[2].sand=16.0;
cls[2].clay=35.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -471,8 +466,8 @@
cls[1].clay=28.0;
cls[2].sand=16.0;
cls[2].clay=35.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -485,8 +480,8 @@
cls[1].clay=28.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -499,8 +494,8 @@
cls[1].clay=46.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -513,8 +508,8 @@
cls[1].clay=46.0;
cls[2].sand=35.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -527,8 +522,8 @@
cls[1].clay=46.0;
cls[2].sand=45.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=50.0\n");
hf=50.0;
@@ -541,8 +536,8 @@
cls[1].clay=46.0;
cls[2].sand=45.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -555,8 +550,8 @@
cls[1].clay=58.0;
cls[2].sand=45.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -569,8 +564,8 @@
cls[1].clay=46.0;
cls[2].sand=23.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -583,8 +578,8 @@
cls[1].clay=44.0;
cls[2].sand=23.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -597,8 +592,8 @@
cls[1].clay=44.0;
cls[2].sand=23.0;
cls[2].clay=42.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -611,8 +606,8 @@
cls[1].clay=44.0;
cls[2].sand=14.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -625,8 +620,8 @@
cls[1].clay=44.0;
cls[2].sand=14.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -639,8 +634,8 @@
cls[1].clay=44.0;
cls[2].sand=0.0;
cls[2].clay=28.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=80.0\n");
hf=80.0;
@@ -653,8 +648,8 @@
cls[1].clay=44.0;
cls[2].sand=0.0;
cls[2].clay=54.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=125.0\n");
hf=125.0;
@@ -667,8 +662,8 @@
cls[1].clay=44.0;
cls[2].sand=38.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=125.0\n");
hf=125.0;
@@ -681,8 +676,8 @@
cls[1].clay=60.0;
cls[2].sand=38.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=125.0\n");
hf=125.0;
@@ -695,8 +690,8 @@
cls[1].clay=60.0;
cls[2].sand=38.0;
cls[2].clay=55.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=125.0\n");
hf=125.0;
@@ -709,8 +704,8 @@
cls[1].clay=60.0;
cls[2].sand=34.0;
cls[2].clay=66.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("hf=125.0\n");
hf=125.0;
Modified: grass-addons/gipe/r.soiltex2prop/prct2ksat.c
===================================================================
--- grass-addons/gipe/r.soiltex2prop/prct2ksat.c 2008-05-30 08:34:22 UTC (rev 31579)
+++ grass-addons/gipe/r.soiltex2prop/prct2ksat.c 2008-05-30 09:33:34 UTC (rev 31580)
@@ -16,8 +16,6 @@
double temp,ksat;
double silt_input=0.0; //Rawls et al (1990)
//do not have silt input
- // set up mark index for inside/outside polygon check
- double mark[POLYGON_DIMENSION]={0.0};
//printf("in prct2ksat(), cm/h\n");
//setup the 3Dvectors and initialize them
struct vector cls[POLYGON_DIMENSION] = {0.0};
@@ -27,9 +25,6 @@
cls[i].clay=0.0;
cls[i].silt=0.0;
}
- //transform input from [0,1] to [0,100]
- sand_input *= 100.0;
- clay_input *= 100.0;
//fill up initial polygon points
cls[0].sand=0.0;
cls[0].clay=100.0;
@@ -38,8 +33,8 @@
cls[2].sand=7.0;
cls[2].clay=56.0;
//Get started
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
ksat=0.0025;
index=1;
//printf("Ksat=0.0025\n");
@@ -51,8 +46,8 @@
cls[1].clay=56.0;
cls[2].sand=30.0;
cls[2].clay=56.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0025\n");
ksat=0.0025;
@@ -65,8 +60,8 @@
cls[1].clay=56.0;
cls[2].sand=40.0;
cls[2].clay=63.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0025\n");
ksat=0.0025;
@@ -79,8 +74,8 @@
cls[1].clay=0.0;
cls[2].sand=90.0;
cls[2].clay=10.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=25.0\n");
ksat=25.0;
@@ -93,8 +88,8 @@
cls[1].clay=10.0;
cls[2].sand=80.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=15.0\n");
ksat=15.0;
@@ -107,8 +102,8 @@
cls[1].clay=10.0;
cls[2].sand=85.0;
cls[2].clay=15.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=15.0\n");
ksat=15.0;
@@ -121,8 +116,8 @@
cls[1].clay=0.0;
cls[2].sand=70.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=7.5\n");
ksat=7.5;
@@ -135,8 +130,8 @@
cls[1].clay=15.0;
cls[2].sand=75.0;
cls[2].clay=25.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=7.5\n");
ksat=7.5;
@@ -149,8 +144,8 @@
cls[1].clay=25.0;
cls[2].sand=68.0;
cls[2].clay=23.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=2.5\n");
ksat=2.5;
@@ -163,8 +158,8 @@
cls[1].clay=23.0;
cls[2].sand=70.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=2.5\n");
ksat=2.5;
@@ -177,8 +172,8 @@
cls[1].clay=0.0;
cls[2].sand=68.0;
cls[2].clay=23.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=2.5\n");
ksat=2.5;
@@ -191,8 +186,8 @@
cls[1].clay=0.0;
cls[2].sand=22.0;
cls[2].clay=5.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=1.5\n");
ksat=1.5;
@@ -205,8 +200,8 @@
cls[1].clay=0.0;
cls[2].sand=68.0;
cls[2].clay=23.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=1.5\n");
ksat=1.5;
@@ -219,8 +214,8 @@
cls[1].clay=23.0;
cls[2].sand=60.0;
cls[2].clay=25.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=1.5\n");
ksat=1.5;
@@ -233,8 +228,8 @@
cls[1].clay=23.0;
cls[2].sand=70.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=1.5\n");
ksat=1.5;
@@ -247,8 +242,8 @@
cls[1].clay=35.0;
cls[2].sand=70.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=1.5\n");
ksat=1.5;
@@ -261,8 +256,8 @@
cls[1].clay=0.0;
cls[2].sand=22.0;
cls[2].clay=5.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -275,8 +270,8 @@
cls[1].clay=0.0;
cls[2].sand=22.0;
cls[2].clay=5.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -289,8 +284,8 @@
cls[1].clay=20.0;
cls[2].sand=22.0;
cls[2].clay=5.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -303,8 +298,8 @@
cls[1].clay=20.0;
cls[2].sand=22.0;
cls[2].clay=5.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -317,8 +312,8 @@
cls[1].clay=20.0;
cls[2].sand=57.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -331,8 +326,8 @@
cls[1].clay=35.0;
cls[2].sand=57.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -345,8 +340,8 @@
cls[1].clay=35.0;
cls[2].sand=57.0;
cls[2].clay=30.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.8\n");
ksat=0.8;
@@ -359,8 +354,8 @@
cls[1].clay=40.0;
cls[2].sand=63.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -373,8 +368,8 @@
cls[1].clay=30.0;
cls[2].sand=63.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -387,8 +382,8 @@
cls[1].clay=30.0;
cls[2].sand=38.0;
cls[2].clay=23.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -401,8 +396,8 @@
cls[1].clay=30.0;
cls[2].sand=38.0;
cls[2].clay=23.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -415,8 +410,8 @@
cls[1].clay=20.0;
cls[2].sand=20.0;
cls[2].clay=12.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -429,8 +424,8 @@
cls[1].clay=20.0;
cls[2].sand=20.0;
cls[2].clay=12.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -443,8 +438,8 @@
cls[1].clay=0.0;
cls[2].sand=20.0;
cls[2].clay=12.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -457,8 +452,8 @@
cls[1].clay=0.0;
cls[2].sand=0.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -471,8 +466,8 @@
cls[1].clay=3.0;
cls[2].sand=0.0;
cls[2].clay=0.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.5\n");
ksat=0.5;
@@ -485,8 +480,8 @@
cls[1].clay=3.0;
cls[2].sand=9.0;
cls[2].clay=18.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -499,8 +494,8 @@
cls[1].clay=3.0;
cls[2].sand=9.0;
cls[2].clay=18.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -513,8 +508,8 @@
cls[1].clay=20.0;
cls[2].sand=9.0;
cls[2].clay=18.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -527,8 +522,8 @@
cls[1].clay=20.0;
cls[2].sand=9.0;
cls[2].clay=18.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -541,8 +536,8 @@
cls[1].clay=20.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -555,8 +550,8 @@
cls[1].clay=20.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -569,8 +564,8 @@
cls[1].clay=35.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -583,8 +578,8 @@
cls[1].clay=35.0;
cls[2].sand=55.0;
cls[2].clay=35.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -597,8 +592,8 @@
cls[1].clay=35.0;
cls[2].sand=55.0;
cls[2].clay=35.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -611,8 +606,8 @@
cls[1].clay=35.0;
cls[2].sand=57.0;
cls[2].clay=44.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.3\n");
ksat=0.3;
@@ -625,8 +620,8 @@
cls[1].clay=35.0;
cls[2].sand=57.0;
cls[2].clay=44.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -639,8 +634,8 @@
cls[1].clay=45.0;
cls[2].sand=57.0;
cls[2].clay=44.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -653,8 +648,8 @@
cls[1].clay=35.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -667,8 +662,8 @@
cls[1].clay=38.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -681,8 +676,8 @@
cls[1].clay=38.0;
cls[2].sand=33.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -695,8 +690,8 @@
cls[1].clay=38.0;
cls[2].sand=13.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -709,8 +704,8 @@
cls[1].clay=18.0;
cls[2].sand=13.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -723,8 +718,8 @@
cls[1].clay=18.0;
cls[2].sand=13.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -737,8 +732,8 @@
cls[1].clay=25.0;
cls[2].sand=13.0;
cls[2].clay=29.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.15\n");
ksat=0.15;
@@ -751,8 +746,8 @@
cls[1].clay=25.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -765,8 +760,8 @@
cls[1].clay=25.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -779,8 +774,8 @@
cls[1].clay=38.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -793,8 +788,8 @@
cls[1].clay=38.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -807,8 +802,8 @@
cls[1].clay=38.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -821,8 +816,8 @@
cls[1].clay=38.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -835,8 +830,8 @@
cls[1].clay=45.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.075\n");
ksat=0.075;
@@ -849,8 +844,8 @@
cls[1].clay=41.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -863,8 +858,8 @@
cls[1].clay=53.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -877,8 +872,8 @@
cls[1].clay=53.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -891,8 +886,8 @@
cls[1].clay=58.0;
cls[2].sand=50.0;
cls[2].clay=50.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -905,8 +900,8 @@
cls[1].clay=41.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -919,8 +914,8 @@
cls[1].clay=33.0;
cls[2].sand=15.0;
cls[2].clay=38.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -933,8 +928,8 @@
cls[1].clay=33.0;
cls[2].sand=0.0;
cls[2].clay=51.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.025\n");
ksat=0.025;
@@ -947,8 +942,8 @@
cls[1].clay=60.0;
cls[2].sand=0.0;
cls[2].clay=51.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -961,8 +956,8 @@
cls[1].clay=50.0;
cls[2].sand=0.0;
cls[2].clay=51.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -975,8 +970,8 @@
cls[1].clay=50.0;
cls[2].sand=18.0;
cls[2].clay=53.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -989,8 +984,8 @@
cls[1].clay=56.0;
cls[2].sand=18.0;
cls[2].clay=53.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -1003,8 +998,8 @@
cls[1].clay=56.0;
cls[2].sand=18.0;
cls[2].clay=53.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -1017,8 +1012,8 @@
cls[1].clay=56.0;
cls[2].sand=43.0;
cls[2].clay=58.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
@@ -1031,8 +1026,8 @@
cls[1].clay=56.0;
cls[2].sand=43.0;
cls[2].clay=58.0;
- mark[0]=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
- if(mark[0]==1){
+ index=point_in_triangle(sand_input,clay_input,silt_input,cls[0].sand,cls[0].clay,cls[0].silt,cls[1].sand,cls[1].clay,cls[1].silt,cls[2].sand,cls[2].clay,cls[2].silt);
+ if(index==1){
index=1;
//printf("Ksat=0.0075\n");
ksat=0.0075;
More information about the grass-commit
mailing list