! NOAA/PMEL TMAP
 ! FERRET v7.6 (optimized)
 ! Linux 4.18.0-372.32.1.el8_6.x86_64 64-bi - 07/06/24
 ! 10-Feb-25 15:06     

let ver="0.1"
use drainage+antarctica-v`ver`-vec.nc
let landmask=if (dcount ge 0) then 1 else 0

palette blue_red_nowhite_centered

let regnames={ "world", \
"baffingreen", "northamerica", "southamerica", \
"europe", "mideurope", "scandinavia", \
"mediteran", "antarctica", "africa", "china", \
"arabia", "siberia", \
"amazon", "stlawrence", "mackenzie", \
"lena", "yenisei", "ob", "mississippi", \
"mesopotamia", "yseawest", "chinasouth", \
"bassiberia", "baspolaramerica", "basgreenland", \
"seasia", "indiaeast", "indiawest" }

let regtitles={ "World", \
"Baffin+Greenland", "North America", "South America", \
"Europe", "Mid Europe", "Scandinavia", \
"Mediteran", "Antarctica", "Africa", "China", \
"Arabia", "Siberia", \
"Amazon", "St Lawrence", "MacKenzie", \
"Lena", "Yenisei", "Ob", "Mississippi", \
"Mesopotamia", "Yellow Sea West", "South China", \
"N Pole/Siberia", "N Pole/America", "Greenland",
"SouthEast Asia", "India East", "India West" }

let regbounds={ "X=0:360,Y=-90:90", \
"X=-93:-15,Y=58:85", "X=-168:-54,Y=22:80", "X=-85:-35,Y=-60:15", \
"X=-10:40,Y=36:72", "X=-10:30,Y=36:60", "X=4:40,Y=55:72", \
"X=-10:40,Y=30:45", "X=0:360,Y=-90:-60", "X=-20:62,Y=-37:38", "X=60:150,Y=20:70", \
"X=30:65,Y=10:35", "X=60:192,Y=60:80", \
"X=80W:40W,Y=10S:12N", "X=-85:-60,Y=40:55", "X=-136:-110,Y=60:72", \
"X=120:135,Y=60:75", "X=80:90,Y=60:75", "X=60:75,Y=60:75", "X=-100:-80,Y=20:40" , \
"X=X=35:55,Y=25:45", "X=110:125,Y=28:45", "X=100:120,Y=8:25", \
"X=20:180,Y=70:90", "X=-180:-60,Y=65:90", "X=-72:10,Y=58:90", \
"X=90:110,Y=8:25", "X=75:93,Y=7:27", "X=65:85,Y=7:25" \
}

let nregs=`regnames[i=@ngd]`

repeat/name=rr/range=1:`nregs` ( \
  let regname=regnames[i=`rr`] ; \
  let regtitle=regtitles[i=`rr`] ; \
  let regbound=regbounds[i=`rr`] ; \
  let dreg=dcount[`regbound`] ; \
  let lmreg=landmask[`regbound`] ; \
  shade/levels=(0,1000,10)(inf)/title="`regtitle`" if dreg gt 10 then dreg ; \
  contour/over/nolabel/nokey/levels=2/size=0.00001/color=black lmreg ; \
  plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy ; \
  frame/file=drainage+antarctica-v`ver`-`regname`.gif ; \
)



def sym world=        X=-180:180/Y=-90:90
def sym baffingreen=  X=-93:-12/Y=59:85
def sym northamerica= X=-168:-54/Y=22:80
def sym southamerica= X=-85:-35/Y=-60:15
def sym mideurope=    X=-10:30/Y=36:60
def sym scandinavia=  X=4:40/Y=55:72
def sym europe=       X=-10:40/Y=36:72
def sym mediteran=    X=-10:40/Y=30:45
def sym antarctica=   X=0:360/Y=-90:-60
def sym africa=       X=-20:62/Y=-37:38
def sym china=        X=60:150/Y=20:70
def sym arabia=       X=30:65/Y=10:35
def sym siberia=      X=60:192/Y=60:80

def sym amazon=       X=80W:40W/Y=10S:12N
def sym stlawrence=   X=-85:-60/Y=40:55
def sym mackenzie=    X=-136:-110/Y=60:72
def sym lena=         X=120:135/Y=60:75
def sym yenisei=      X=80:90/Y=60:75
def sym ob=           X=60:75/Y=60:75
def sym mississippi=  X=-100:-80/Y=20:40
def sym mesopotamia=  X=35:55/Y=25:45
def sym yseawest=     X=110:125/Y=28:45

! ocean basins for checking runoff into ocean
def sym bassiberia=X=20:180/Y=70:90
def sym baspolaramerica=X=-180:-60/Y=65:90
def sym basgreenland=    X=-72:-10/Y=58:90
def sym seasia=       X=90:110/Y=8:25
def sym indiaeast=    X=75:93/Y=7:27
def sym indiawest=    X=65:85/Y=7:25

!shade/levels=(0,1000,10)(inf)/title="DDM30+antarctica" if dcount gt 10 then dcount
!contour/over/nolabel/nokey/levels=2/size=0.00001/color=black landmask
!!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-world.gif


!shade/($amazon)/levels=(0,1000,10)(inf)/title="Amazon" if dcount gt 10 then dcount
!contour/($amazon)/over/nolabel/nokey/levels=2/size=0.00001/color=black landmask
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-amazon.gif

!shade/($stlawrence)/title="St Lawrence" dcount
!contour/($stlawrence)/over/nolabel/nokey/levels=2/size=0.00001/color=black landmask
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-stlawrence.gif

!shade/($mackenzie)/title="MacKenzie" dcount
!contour/($mackenzie)/over/nolabel/nokey/levels=2/size=0.00001/color=black landmask
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-mackenzie.gif

!shade/($northamerica)/title="North America" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-northamerica.gif

!shade/($africa)/title="Africa" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-africa.gif

!shade/($siberia)/title="Siberia" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-siberia.gif

!shade/($ob)/title="Ob" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-ob.gif

!shade/($yenisei)/title="Yenisei" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-yenisei.gif

!shade/($lena)/title="Lena" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-lena.gif

!shade/($china)/title="China" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-china.gif

!shade/($mediteran)/title="Mediteran" dcount
!plot/vs/line/over/color=red/nolabel basinboundsx,basinboundsy
!frame/file=drainage+antarctica-v`ver`-mediteran.gif