d.legend GLOBCOVER22_CA_lab use=30,40,50,60,70,90,100,110,120,130,140,150,160,170,180,190,200,210,200 -m

r.category GLOBCOVER22_CA_lab rules=-<<EOF
  14:Rainfed cropland
  20:Mosaic cropland / grassland/shrubland/forest
? 30:Mosaic grassland/shrubland/forest / cropland
? 40:Closed to open (>15%) broadleaved evergreen forest (>5m)
? 50:Closed (>40%) broadleaved deciduous forest (>5m)
? 60:Open (15-40%) broadleaved deciduous forest/woodland (>5m)
? 70:Closed (>40%) needleleaved evergreen forest (>5m)
? 90:Open (15-40%) needleleaved deciduous or evergreen forest (>5m)
? 100:Closed to open (>15%) mixed forest (>5m)
? 110:Mosaic forest or shrubland / grassland
? 120:Mosaic grassland / forest or shrubland
? 130:Closed to open (>15%) shrubland (<5m)
? 140:Closed to open (>15%) herbaceous vegetation
? 150:Sparse (<15%) vegetation
? 160:Closed to open (>15%) broadleaved forest regularly flooded
? 170:Closed (>40%) broadleaved forest or shrubland permanently flooded
? 180:Closed to open (>15%) grassland or woody vegetation on regularly flooded
? 190:Artificial surfaces and associated areas (Urban areas >50%)
? 200:Bare areas
? 210:Water bodies
? 220:Permanent snow and ice
? EOF

r.recode GLOBCOVER22_CA_lab out=Landcov_forest rules=lcforest.txt
14:30:1:1
40:110:2:2
120:150:1:1
160:170:2:2
180:220:1:1

r.category Landcov_forest rules=-<<EOF
1 other
2 forest

d.vect bstreamsacc50t co=aqua
nsfall
d.rast srtm90ctr_shade
d.vect bstreamsacc50t 
d.vect srtmstream90med3t co=blue width=2
d.vect nsfpts siz=18 icon=basic/marker co=red fcol=red
d.vect pacorabasbigBW type=boundary co=yellow width=2
d.vect chagbasinBW type=boundary co=green width=2
d.vect pacorabasbigBW type=boundary co=yellow

nsf_basins only
v.overlay ain=bstreamsacc50t atype=line bin=pacorabasbigBW out=pacora_bstreamsacc50 op=and
v.overlay ain=srtmstream90med3t atype=line bin=pacorabasbigBW out=pacora_stream90med op=and
v.overlay ain=srtmstream90med3t atype=line bin=basin_chagres out=chag_stream90med op=and
v.overlay ainput=bstreamsacc50t atype=line binput=basin_chagres out=chag_bstreamsacc50 operator=and --o

d.vect chagbasinBW fcol=220:220:220 col=green 
d.vect pacorabasbigBW fcol=230:230:230 col=yellow
d.vect chag_bstreamsacc50 co=aqua
d.vect pacora_bstreamsacc50 co=aqua 
d.vect pacora_stream90med co=blue width=2
d.vect chag_stream90med co=blue width=2

d.vect nsfpts siz=18 icon=basic/marker co=red fcol=red 

d.rgb b=LSp012r054_20020528z17_30 g=LSp012r054_20020528z17_40 r=LSp012r054_20020528z17_50

awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt>Transpanam_map2009basrec1.txt

r.watershed elev="srtm90panamalin" threshold=5000 accum="accum200" drain="drain200" bas="basin200
r.watershed elev="panama_srtm3_90m" threshold=5000 accum="accumv3_200" drain="drainiv3_200" bas="basiniv3_200
r.watershed elev="SRTM_v2_90m" threshold=20000 accum="accumv2_90" drain="drainiv2_90" bas="basiniv2_90"
r.mapcalc "streams_srtm3_200m=if(abs(accumv3_200)>50,1,null())"
r.mapcalc "streams_srtm3_200med=if(abs(accumv3_200)>500,1,null())"
r.mapcalc "streams_srtm3_200maj=if(abs(accumv3_200)>5000,1,null())"

r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90m basin=testbasinsrtm3_90m
r.mapcalc "teststreams_srtm3_90=if(abs(testaccumsrtm3_90m)>100,1,null())"
r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90ms basin=testbasinsrtm3_90ms
r.mapcalc "teststreams_srtm3_90s=if(abs(testaccumsrtm3_90ms)>100,1,null())"

r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_10m basin=testbasinifs_10m
r.mapcalc "teststreams_ifs_10m=if(abs(testaccumifs_10m)>500,1,null())"
r.thin teststreams_ifs_10m out=teststreams_ifs_10mt
r.to.vect teststreams_ifs_10mt out=teststreams_ifs_10mt

r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_30mnn basin=testbasinifs_30mnn
r.mapcalc "teststreams_ifs_30mnn=if(abs(testaccumifs_30mnn)>100,1,null())"

r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_30mnn basin=testbasinsrtm3_30mnn
r.mapcalc "teststreams_srtm3_30mnn=if(abs(testaccumsrtm3_30mnn)>100,1,null())"

r.terraflow -s ifsall_10m_patch filled=ifsall_10m_fill dir=ifsall_10m_dir swater=ifsall_10m_swatersh accum=ifsall_10m_accum tci=ifsall_10m_tci

r.terraflow -s srtm90.8378 filled=srtm90.fill dir=srtm90.dir swater=srtm90.swatersh accum=srtm90.accum tci=srtm90.tci
r.terraflow [-sq] elev=string filled=string direction=string
   swatershed=string accumulation=string tci=string [d8cut=value]
   [memory=value] [STREAM_DIR=string] [stats=string]
r.mapcalc "srtmstream90maj=if(abs(srtm90.accum)>5000 && srtm90.8378 > 0,1,null())"
r.mapcalc "srtmstream90med=if(abs(srtm90.accum)>1000 && srtm90.8378 > 0,1,null())"
r.mapcalc "srtmstream90med2=if(abs(srtm90.accum)>2000 && srtm90.8378 > 0,1,null())"
r.mapcalc "srtmstream90=if(abs(srtm90.accum)>100 && srtm90.8378 > 0,1,null())"
r.mapcalc "srtmstream90med3=if(abs(accum90ctr) > 1000 && srtm90ctr > 0, 1, null())" 

r.mapcalc="if(bbasin.50m30K == 34 || bbasin.50m30K == 36 || bbasin.50m30K == 38 || bbasin.50m30K == 32, 1, null())"
r.mapcalc "pacorabasin_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118, 1, null())"
r.mapcalc "pacorabasino_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 || basinsctr90mask==116 || basinsctr90mask==110, 1, null())"
r.mapcalc "chagresbasin=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98, 1, null())"
r.mapcalc "basinsmask=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98 ||basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 ||basinsctr90mask==116 || basinsctr90mask==110 || basinsctr90mask==92 || basinsctr90mask==34, 1, null())"
r.mapcalc "basinsmask_inv=if(isnull(basinsmask), 1,null())"
r.mapcalc "pequinibasin90m=if(basinsctr90mask==92, 1, null())"
r.mapcalc "cuangobasin90m=if(basinsctr90mask==34, 1, null())"
r.mapcalc "basinscanal=if(basin200mask==62 || basin200mask==60 || basin200mask==50 || basin200mask==44 ||basin200mask==46 || basin200mask==48 || basin200mask==52 || basin200mask==56 ||basin200mask==54 || basin200mask==366, 1, null())"

r.mapcalc "boundaryctr=if(panama_srtm3_90m,1,null())"
r.to.vect boundaryctr out=boundaryctr feature=area --o

d.rast -o pacorabasin_90m
d.rast -o cuangobasin90m
d.rast -o pequinibasin90m
d.rast -o chagresbasin
d.rast -o pacorabasino_90m
d.vect srtmstream90med3t
d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 fcolor=black
d.rast -o basinsmask_inv
d.vect boundaryctr type=boundary width=2
d.barscale -mt

d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 
d.vect PanamaMaster07_Feb05_GIS icon=basic/marker size=16 width=2
d.vect sstreamctr90t
d.vect srtmstream90medt
d.vect sstreamctr90majt co=black width=2

d.rast pan_srtm3_300msh
d.rast lake1_90m -o
d.rast lake2 -o
d.vect canal width=3
d.vect basinscanal type=boundary width=2
d.vect basin_cuango type=boundary
d.vect basin_pacorabig type=boundary
d.vect stream90m_canal
d.vect stream90m_cuango
d.vect stream90m_pacora

v.overlay ain=srtmstream90majt bin=basinscanal out=stream90m_canal atype=line operator=and
v.overlay ain=srtmstream90majt bin=basin_cuango out=stream90m_cuango atype=line operator=and
v.overlay ain=srtmstream90majt bin=basin_pacorabig out=stream90m_pacora atype=line operator=and

v.in.ascii feb05pts.txt 
r.mapcalc "lsfac=exp(fccum90ctr*90.,1.4)*exp(sin(slopectr_srtm3_90m),1.3)"
r.mapcalc "lsfac2=exp(accum90ctr,1.4)*exp(sin(slope.30m),1.3)"
r.average base=basinsctr90mask cover=slope30m_percc out=basins_slopeperc
r.average base=basinsctr90mask cover=lsfac out=basins_lsfac
r.average base=basinsctr90mask cover=lsfac2 out=basins_lsfac2

d.his i=srtm90ctr_shade h=basins_lsfac
r.mapcalc "basins5ctrBW=if(isnull(pacorabasbigBW),0,1)+if(isnull(chagbasinBW),0,2) \
+ if(isnull(pequinibasBW),0,3)+if(isnull(cuangobasBW),0,4)+if(isnull(diosbasinBW),0,5)"

0 white
1 170 170 170
2 110 110 110
3 210 210 210
4 150 150 150
5 130 130 130

551  20:53   d.rast -o pacorabasbigBW
   552  20:53   d.rast -o pequinibasBW
   553  20:54   d.rast -o chagbasinBW
   554  20:54   d.rast -o cuangobasBW
   555  20:55   history
   556  20:55   r.colors diosbasinBW co=rules

MAPS
d.his i=pan_srtm3_300msh h=panama_srtm3_90m bright=20
d.vect Transpanama_map2009_pts attrcol=int_1 disp=attr lcol=blue lsiz=12
d.vect Transpanama_map2009_pts co=red siz=15 icon=basic/marker fcol=red width=2 
d.vect srtmstream90majt co=blue
d.vect sstreamsouth90maj co=blue
d.grid
d.vect Transpanama_pubJuly09xy attrcol=str_1 disp=attr lcol=blue lsiz=12

    r.colors basinsv3rec2_200m co=random
    d.rast -o basinsv3rec2_200m
    d.vect streams_srtm3_200medt
    d.vect Transpanama_map2009_pts co=black siz=20 icon=basic/marker fcol=red

v.in.ascii Transpanama_pubJuly09xy.txt out=Transpanama_pubJuly09xy x=3 y=4 fs=, skip=1
v.out.ascii Transpanama_map2009_pts fs=space | r.what basinsv3_200m
v.out.ascii Transpanama_pub2009xy fs=space | r.what basinsv3_200m
v.out.ascii Transpanama_pubJuly09xy fs=space | r.what accumv3_200ma
  vi Transpanam_map2009_coordbasins.txt
  cp Transpanam_map2009_coordbasins.txt Transpanam_map2009basrec.txt
  vi Transpanam_map2009basrec.txt
  cp Transpanam_map2009basrec.txt Transpanam_map2009basrec.txtbak
  awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt > Transpanam_map2009basrec1.txt
  awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt > Transpanam_map2009basrec1.txt
awk 'BEGIN {FS="|"; OFS=":"} {print $4,$4,-1*$4,-1*$4}' Transpanama_pub2009_coordbasins.txt >Transpanama_pub2009basrec.txt
awk 'BEGIN {FS="|"; OFS=" "} {print $4,$5,$2,$1}' Transpanam_names_xy.txt >Transpanam_xy_names.txt
r.what basinsv3_200m < Transpanam_xy_names2.txt>Transpanam_xy_names_bas.txt
r.what basinsv3recpub_200m4 < Transpanam_xy_names2.txt > Transpanam_xy_names_bas4.txt
awk 'BEGIN {FS="|"; OFS=":"} {print -1.*$4,$3}' Transpanam_xy_names_bas.txt >Transpanam_baslabelsneg.txt
awk 'BEGIN {FS="|"; OFS=":"} {print $4,$3}' Transpanam_xy_names_bas4.txt >Transpanam_baslabel.txt
r.category basinsv3rec_200m rules=Transpanam_baslabels.txt
r.category basinsv3recpub_200m4 rules=Transpanam_baslabel.txt
r.report basinsv3recpub_200m4,GLOBCOVER22_CA_lab unit=p,k -n
  1470  23:43   vi Transpanam_map2009basrec1.txt
  1471  23:43   r.recode help
  1472  23:44   r.recode basinsv3_200m out=basinsv3_200m rules=Transpanam_map2009basrec1.txt
  1473  23:44   r.recode basinsv3_200m out=basinsv3rec_200m rules=Transpanam_map2009basrec1.txt
  1473  23:44   r.recode basinsv3_200m2 out=basinsv3rec_200m rules=Transpanam_map2009basrec2.txt
  1474  23:45   d.rast basinsv3rec_200m
  1475  23:47   d.his h=basinsv3rec_200m i=pan_srtm3_300msh
  1476  23:48   d.vect Transpanama_map2009_pts co=black siz=15 icon=basic/marker fcol=red

r.watershed elev="panama_srtm3_90m" threshold=3000 accum="accumv3w2_90" drain="drainiv3w2_90" bas="basiniv3w2_90"
r.watershed elev="panama_srtm3_90m" threshold=700 accum="accumv3s1_90" drain="drainiv3s1_90" bas="basiniv3s1_90"
r.recode basinsv3_200m out=basinsv3recpub_200m3 rules=Transpanama_pub2009basrec3.txt
r.water.outlet drain=drainiv3w2_90 basin=basinsmwest e=293515 n=913393.
r.water.outlet drain=drainiv3n_90 basin=RioNobreDios e=666131 n=1057687
v.out.ascii Transpanama_pub2009xy fs=space | r.what basinsv3_200m
vi Transpanama_pub2009_coordbasins.txt
awk 'BEGIN {FS="|"; OFS=":"} {print $4,$4,-1*$4,-1*$4}' Transpanama_pub2009_coordbasins.txt > Transpanama_pub2009basrec.txt
r.recode basinsv3_200m out=basinsv3recpub_200m rules=Transpanama_pub2009basrec.txt

r.recode basinsv3_200m out=basinsv3recpub_200m3 rules=Transpanama_pub2009basrec3.txt
r.recode basinsv3_200m out=basinsv3recpub_200m4 rules=Transpanama_pub2009basrec4.txt
r.recode basiniv3s2_90 out=LajasHato_90m rules=LajasHatorecode.txt

Data range is -458 to -34
r.report basinsv3recpub_200m,GLOBCOVER22_CA_lab unit=p,k -n

r.mapcalc "lancovpub2009=if(basinsv3recpub_200m,GLOBCOVER22_CA_lab)"
  r.colors lancovpub2009 rast=GLOBCOVER22_CA_lab
  d.rast -o lancovpub2009
  d.vect Transpanama_pub2009xy siz=9 icon=basic/circle fco=black
  d.vect streams_srtm3_200medt
  d.out.file transpan_basinlandcov

> -458 20 20 20
> -362 80 80 80
> -200 120 120 120
> -100 150 150 150
> -34 200 200 200
> end

 -458 160 160 160
 -440 20 20 20
 -406 60 60 60
 -394 80 80 80
 -362 100 100 100
 -243 120 120 120
 -238 140 140 140
 -196 150 150 150
 -62 170 170 170
 -34 190 190 190
#d.vect streams_srtm3_200t
#d.vect streams_srtm3_200majt co=blue width=2

MAP publ
g.region rast=basinsv3recpub_200m4 -p
d.rast pan_srtm3_300msh
d.rast -o basinsv3recpub_200m4col2
d.rast -o basins_smallw3
d.rast -o cuangobasin90m
d.rast -o RioNobreDios
d.rast -o LajasHato_90m
d.vect streams_srtm3_200medt
d.vect Transpanama_pub2009xy siz=7 icon=basic/circle fco=black
d.barscale -m
d.out.file transpan_basinsall_fix

r.mapcalc "TranspanBasinpub=basinsv3recpub_200m4+basins_smallw3+cuangobasin90m+RioNobreDios+LajasHato_90m"
r.mapcalc "Tbasins_small=basins_smallw3+10*cuangobasin90m+15*RioNobreDios+LajasHato_90m+chagresbasin"
#r.mapcalc "Tbasins_test=if(cuangobasin90m,11,0)+if(chagresbasin,13,0)" does not work
r.mapcalc "Tbasins_test=if(isnull(cuangobasin90m),0,11)+if(isnull(chagresbasin),0,15)" 
r.mapcalc "Transbasins_all=if(isnull(cuangobasin90m),0,11)+if(isnull(chagresbasin),0,15)+if(isnull(basinsv3recpub_200m4),0,basinsv3recpub_200m4)+if(isnull(basins_smallw3),0,basins_smallw3)+if(isnull(RioNobreDios),0,18)+if(isnull(LajasHato_90m),0,LajasHato_90m)" 
15*RioNobreDios+LajasHato_90m+chagresbasin"
r.mapcalc "lancovpub2009f=if(basinsv3recpub_200m4col2,GLOBCOVER22_CA_lab)"
r.mapcalc "lancovbas_smallw=if(basins_smallw3,GLOBCOVER22_CA_lab)"
r.mapcalc "lancovcuango=if(cuangobasin90m,GLOBCOVER22_CA_lab)"
r.colors lancovcuango rast=GLOBCOVER22_CA_lab
r.mapcalc "lancovRioNobreDios=if(RioNobreDios,GLOBCOVER22_CA_lab)"
r.colors lancovRioNobreDios rast=GLOBCOVER22_CA_lab
r.mapcalc "lancovLajasHato=if(LajasHato_90m,GLOBCOVER22_CA_lab)"
r.colors lancovLajasHato rast=GLOBCOVER22_CA_lab

g.region rast=basinsv3recpub_200m4 -p
d.rast pan_srtm3_300msh
d.rast lancovpub2009f -o
d.vect basinsv3recpub_200m4 type=boundary co=brown width=2
d.rast -o lancovcuango
d.vect basin_cuango co=brown type=boundary width=2
d.rast -o lancovRioNobreDios
d.vect RioNobreDios co=brown width=2 type=boundary
d.rast -o lancovbas_smallw
d.vect basins_smallw3 co=brown type=boundary 
d.rast -o lancovLajasHato
#d.vect lancovLajasHato co=brown type=boundary 
d.vect streams_srtm3_200medt
d.vect Transpanama_pub2009xy siz=7 icon=basic/circle fco=black
d.legend GLOBCOVER22_CA_lab use=30,40,50,60,70,90,100,110,120,130,140,150,160,170,180,190,200,210,200 -m
d.barscale -m

g.region rast=basinsv3recpub_200m4 -p
d.rast pan_srtm3_300msh
d.rast -o Transbasins_alln_co3
d.vect streams_srtm3_200medt
d.vect Transpanama_pubJuly09xy siz=9 icon=basic/circle fco=black

r.mapcalc "lancovforpub2009f=if(Transbasins_alln_co3,Landcov_forest)"

d.rast pan_srtm3_300msh
d.rast -o lancovforpub2009f
d.vect Transbasins_alln_co3 type=boundary co=brown width=2
d.vect streams_srtm3_200medt
d.vect Transpanama_pubJuly09xy siz=9 icon=basic/circle fco=black