! jnl file that plots the station location of Argo floats used to make the grid
! requires input arguments for time range, depth range and lat/lon range

! to run the program in ferret, you will need to type the following line with the various values filled in:

!  go fer_psal_station_loc time_start time_end depth_min depth_max lat_min lat_max lon_min lon_max iland 

! If you choose not to give values for the various variables, default ones will be substituted for you.
! If you only want to choose some values, you can, but you cannot skip any values.  This means you can 
! fill in the first 2 and no others, or if you want to specify the lat/lon ranges, you must fill in values
! for the first 4 variables as well as the lat/lon ranges.

! depth values currently correspond to the surface region 0-5db
! time values are entered in corresponding to monthly values starting with Jan, 2004 

! iland  0=hig res, 1=low res go fland
! if no values are entered, the default region is 50S to 25N and 130E to 250E at the surface in Jan 2006
cancel/all data
cancel/all var
cancel/all sym
cancel region
set window/clear
go set_pixel_size 1024 768
set mem/size=200

! assigning variables to arguments passed in the go command as well as default values
let lat_min=`$5%-50%`
let lat_max=`$6%25%`
let lon_min=`$7%130%`
let lon_max=`$8%250%`
let depth_min=`$3%5%`
let depth_max=`$4%5%`
let time_start=`$1%132%`
let time_end=`$2%133%`
let iland=`$9%2%`
let dia=0.2
set mode interpolate
let lonmax=if lon_max gt 360 then 1 else 0 
!let rcheck=`lon_max`-`lon_min`

! setting map scales properly
let beq=if (`lon_min`) eq (`lon_max`) then 1 else 0
if `beq` then
let lon_max=`lon_max`+360
endif
let blm=if (`lon_min`) gt (`lon_max`) then 1 else 0
if `blm` then
let lon_max=`lon_max`+360
endif
define symbol xdeg `lon_max-lon_min`
define symbol ydeg `lat_max-lat_min`
define symbol alasp `($ydeg)/($xdeg)`
set win/asp=($alasp):axis

define axis/x=20.5E:379.5E:1/units=degrees lon
set axis/modulo lon
define axis/y= 89.5S: 89.5N:1/units=degrees lat
!define axis/t=15-jan-2004:15-dec-2009:1/units=months time
define grid/x=lon/y=lat bath


file/order=xy/format="(10f8.3)"/var=count/col=360/skip=0/grid=bath "station_density"
set variable/bad=0 count
shade/set/nokey/nolab/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels="(0,1000,1000)" count+5000
ppl ylab
ppl xlab
ppl shade

define axis/x/x=1:1000000:1 longax
define grid/x=longax biggrid


file/var=fl_month/skip=0 "argo_list_map_sum_psal_2004"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2004"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2005"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2005"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2006"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2006"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2007"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2007"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2008"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2008"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2009"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2009"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2010"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2010"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2011"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2011"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2012"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2012"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2013"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2013"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2014"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2014"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2015"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2015"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2016"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2016"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2017"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2017"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2018"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2018"

! defining constants to help make plots

let s2004=if `time_start` le 12 then 1 else 0
let e2004=if `time_end` le 12 then 1 else 0
let s2005=if `time_start` le 24 and `time_start` gt 12 then 1 else 0
let e2005=if `time_end` le 24 and `time_end` gt 12 then 1 else 0
let s2006=if `time_start` le 36 and `time_start` gt 24 then 1 else 0
let e2006=if `time_end` le 36 and `time_end` gt 24 then 1 else 0
let s2007=if `time_start` le 48 and `time_start` gt 36 then 1 else 0
let e2007=if `time_end` le 48 and `time_end` gt 36 then 1 else 0
let s2008=if `time_start` le 60 and `time_start` gt 48 then 1 else 0
let e2008=if `time_end` le 60 and `time_end` gt 48 then 1 else 0
let s2009=if `time_start` le 72 and `time_start` gt 60 then 1 else 0
let e2009=if `time_end` le 72 and `time_end` gt 60 then 1 else 0
let s2010=if `time_start` le 84 and `time_start` gt 72 then 1 else 0
let e2010=if `time_end` le 84 and `time_end` gt 72 then 1 else 0
let s2011=if `time_start` le 96 and `time_start` gt 84 then 1 else 0
let e2011f=if `time_end` le 90 and `time_start` gt 84 then 1 else 0
let e2011s=if `time_end` le 96 and `time_start` gt 90 then 1 else 0
let s2012=if `time_start` le 108 and `time_start` gt 96 then 1 else 0
let e2012f=if `time_end` le 102 and `time_end` gt 96 then 1 else 0
let e2012s=if `time_end` le 108 and `time_end` gt 102 then 1 else 0
let s2013=if `time_start` le 120 and `time_start` gt 108 then 1 else 0
let e2013f=if `time_end` le 114 and `time_end` gt 108 then 1 else 0
let e2013s=if `time_end` le 120 and `time_end` gt 114 then 1 else 0
let s2014=if `time_start` le 132 and `time_start` gt 120 then 1 else 0
let e2014f=if `time_end` le 126 and `time_end` gt 120 then 1 else 0
let e2014s=if `time_end` le 132 and `time_end` gt 126 then 1 else 0
let s2015=if `time_start` le 144 and `time_start` gt 132 then 1 else 0
let e2015f=if `time_end` le 138 and `time_end` gt 132 then 1 else 0
let e2015s=if `time_end` le 144 and `time_end` gt 138 then 1 else 0
let s2016=if `time_start` le 156 and `time_start` gt 144 then 1 else 0
let e2016f=if `time_end` le 150 and `time_end` gt 144 then 1 else 0
let e2016s=if `time_end` le 156 and `time_end` gt 150 then 1 else 0
let s2017=if `time_start` le 168 and `time_start` gt 156 then 1 else 0
let e2017f=if `time_end` le 162 and `time_end` gt 156 then 1 else 0
let e2017s=if `time_end` le 168 and `time_end` gt 162 then 1 else 0
let s2018=if `time_start` le 180 and `time_start` gt 168 then 1 else 0
let e2018f=if `time_end` le 174 and `time_end` gt 168 then 1 else 0
let e2018s=if `time_end` le 180 and `time_end` gt 174 then 1 else 0


if `s2004` then 
if `e2004` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2004,i=`time_start`]
let time_end_s=`time_end`+1
let sumend=fl_month[d=argo_list_map_sum_psal_2004,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2004,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2004,i=`sumstart`:`sumend`]
say 2004
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2004 eq 0` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2004,i=`time_start`]
let sumend=fl_month[d=argo_list_map_sum_psal_2004,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2004,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2004,i=`sumstart`:`sumend`]
say 2004plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !e2004

endif ! s2004



if `s2005` then
! starts in 2005 and ends in 2005
if `e2005` then
let time_start_s=`time_start`-12
let sumstart=fl_month[d=argo_list_map_sum_psal_2005,i=`time_start_s`]
let time_end_s=`time_end`-11
let sumend=fl_month[d=argo_list_map_sum_psal_2005,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
say 2005
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2004 eq 0` then
! starts in 2005 and ends later

let time_start_s=`time_start`-12
let sumstart=fl_month[d=argo_list_map_sum_psal_2005,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2005,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
say 2005plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !e2005

! starts before 2005 and ends in 2005
elif `e2005` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2005,i=1]
let time_end_s=`time_end`-11
let sumend=fl_month[d=argo_list_map_sum_psal_2005,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
say 2005end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !s2005
 
! starts before 2005 and ends after 2005

let t5= if (`time_start` lt 12 and `time_end` gt 24) then 1 else 0
if `t5` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2005,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2005,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2005,i=`sumstart`:`sumend`]
say all2005
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t5



if `s2006` then
! starts and ends in 2006
if `e2006` then
let time_start_s=`time_start`-24
let sumstart=fl_month[d=argo_list_map_sum_psal_2006,i=`time_start_s`]
let time_end_s=`time_end`-23
let sumend=fl_month[d=argo_list_map_sum_psal_2006,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
say 2006
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2006 eq 0` then

let time_start_s=`time_start`-24
let sumstart=fl_month[d=argo_list_map_sum_psal_2006,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2006,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
say 2006plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2006

! starts before 2006 and ends in 2006
elif `e2006` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2006,i=1]
let time_end_s=`time_end`-23
let sumend=fl_month[d=argo_list_map_sum_psal_2006,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
say 2006end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif !s2006

! starts before 2006 and ends after 2006

let t6= if (`time_start` lt 25 and `time_end` gt 36) then 1 else 0
if `t6` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2006,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2006,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2006,i=`sumstart`:`sumend`]
say all2006
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t6


if `s2007` then
! start and end in 2007
if `e2007` then
let time_start_s=`time_start`-36
let sumstart=fl_month[d=argo_list_map_sum_psal_2007,i=`time_start_s`]
let time_end_s=`time_end`-35
let sumend=fl_month[d=argo_list_map_sum_psal_2007,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
say 2007
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2007 and ends later
elif `e2007 eq 0` then

let time_start_s=`time_start`-36
let sumstart=fl_month[d=argo_list_map_sum_psal_2007,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2007,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
say 2007plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2007

! starts before 2007 and ends in 2007
elif `e2007` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2007,i=1]
let time_end_s=`time_end`-35
let sumend=fl_month[d=argo_list_map_sum_psal_2007,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
say 2007end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif !s2007

! starts before 2007 and ends after 2007

let t7= if (`time_start` lt 37 and `time_end` gt 48) then 1 else 0
if `t7` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2007,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2007,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2007,i=`sumstart`:`sumend`]
say all2007
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t7



if `s2008` then
! start and end in 2008
if `e2008` then
let time_start_s=`time_start`-48
let sumstart=fl_month[d=argo_list_map_sum_psal_2008,i=`time_start_s`]
let time_end_s=`time_end`-47
let sumend=fl_month[d=argo_list_map_sum_psal_2008,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
say 2008
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2008 and ends later
elif `e2008 eq 0` then

let time_start_s=`time_start`-48
let sumstart=fl_month[d=argo_list_map_sum_psal_2008,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2008,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
say 2008plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2008

! starts before 2008 and ends in 2008
elif `e2008` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2008,i=1]
let time_end_s=`time_end`-47
let sumend=fl_month[d=argo_list_map_sum_psal_2008,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
say 2008end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2008

! starts before 2008 and ends after 2008

let t8= if (`time_start` lt 49 and `time_end` gt 60) then 1 else 0
if `t8` then
let sumstart=fl_month[d=argo_list_map_sum_psal_2008,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2008,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2008,i=`sumstart`:`sumend`]
say all2008
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t8



if `s2009` then
! start and end in 2009
if `e2009` then
let time_start_s=`time_start`-60
let sumstart=fl_month[d=argo_list_map_sum_psal_2009,i=`time_start_s`]
let time_end_s=`time_end`-59
let sumend=fl_month[d=argo_list_map_sum_psal_2009,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
say 2009
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2009 and ends later
elif `e2009 eq 0` then

let time_start_s=`time_start`-60
let sumstart=fl_month[d=argo_list_map_sum_psal_2009,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2009,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
say 2009plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2009

! starts before 2009 and ends in 2009
elif `e2009` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2009,i=1]
let time_end_s=`time_end`-59
let sumend=fl_month[d=argo_list_map_sum_psal_2009,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
say 2009end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2009

! starts before 2009 and ends after 2009

let t9=if (`time_start lt 61` and `time_end gt 72`) then 1 else 0
if `t9` then 
let sumstart=fl_month[d=argo_list_map_sum_psal_2009,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2009,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2009,i=`sumstart`:`sumend`]
say all2009
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t9

if `s2010` then
! start and end in 2010
if `e2010` then
let time_start_s=`time_start`-72
let sumstart=fl_month[d=argo_list_map_sum_psal_2010,i=`time_start_s`]
let time_end_s=`time_end`-71
let sumend=fl_month[d=argo_list_map_sum_psal_2010,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
say 2010
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2010 and ends later
elif `e2010 eq 0` then

let time_start_s=`time_start`-72
let sumstart=fl_month[d=argo_list_map_sum_psal_2010,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2010,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
say 2010plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2010

! starts before 2010 and ends in 2010
elif `e2010` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2010,i=1]
let time_end_s=`time_end`-71
let sumend=fl_month[d=argo_list_map_sum_psal_2010,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
say 2010end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2010

! starts before 2010 and ends after 2010

let t10=if (`time_start lt 73` and `time_end gt 84`) then 1 else 0
if `t10` then 
let sumstart=fl_month[d=argo_list_map_sum_psal_2010,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2010,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2010,i=`sumstart`:`sumend`]
say all2010
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t10


if `s2011` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up

! start and end in 2011
	if `e2011f`  then
let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2011s` then
! find out if it starts in first or second half
	if `time_start le 90` then
! do first half first
let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	elif `time_start gt 90` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-84 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 90 loop


! starts in 2011 and ends later
	elif `e2011f eq 0` then
		if `time_start le 90` then 
let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=7] ! only doing first half of year
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plusfirst
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 90` then

let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=13] ! only doing first half of year
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 90 loop


	elif `e2011s eq 0` then

		if `time_start le 90` then 
let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=7] ! only doing first half of year
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plusfirst
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 90` then

let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=13] ! only doing first half of year
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 90 loop
	endif ! starts in 2011 


! starts before 2011 and ends in 2011
elif `e2011f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=1]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2011s` then

! do first have first
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say 2011ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2011

! starts before 2011 and ends after 2011

let t11=if (`time_start lt 85` and `time_end gt 96`) then 1 else 0
if `t11` then 
let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say all2011
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2011,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2011,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2011,i=`sumstart`:`sumend`]
say all2011
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t11

if `s2012` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2012
	if `e2012f` then
let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2012s` then
! find out if it starts in first or second half
	if `time_start le 102` then
! do first half first
let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `time_start gt 102` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-96 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 102 loop





! starts in 2012 and ends later - same for either e2012f or e2012s
	elif `e2012f eq 0` then


		if `time_start le 102` then 

let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 102` then 

let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 102 loop


	elif `e2012s eq 0` then

		if `time_start le 102` then 

let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 102` then 

let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 102 loop


	endif ! starts in 2012

! starts before 2012 and ends in 2012
elif `e2012f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=1]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2012s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say 2012ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2012

! starts before 2012 and ends after 2012

let t12=if (`time_start lt 97` and `time_end gt 108`) then 1 else 0
if `t12` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say all2012first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2012,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2012,i=`sumstart`:`sumend`]
say all2012second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t12




if `s2013` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2013
	if `e2013f` then
let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let time_end_s=`time_end`-107
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2013s` then
! find out if it starts in first or second half 
	if `time_start le 114` then 
! do first half first
let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let time_end_s=`time_end`-107
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia



	elif `time_start gt 114` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-108 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let time_end_s=`time_end`-107
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 114 loop


! starts in 2013 and ends later - same for either e2013f or e2013s
	elif `e2013f eq 0` then


		if `time_start le 114` then 

let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 114` then 

let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 114 loop


	elif `e2013s eq 0` then

		if `time_start le 114` then 

let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 114` then 

let time_start_s=`time_start`-108
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 114 loop


	endif ! starts in 2013

! starts before 2013 and ends in 2013
elif `e2013f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=1]
let time_end_s=`time_end`-107
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2013s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let time_end_s=`time_end`-107
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say 2013ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2013

! starts before 2013 and ends after 2013

let t13=if (`time_start lt 107` and `time_end gt 120`) then 1 else 0
if `t13` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say all2013first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2013,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2013,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2013,i=`sumstart`:`sumend`]
say all2013second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t13




if `s2014` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2014
	if `e2014f` then
let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let time_end_s=`time_end`-119
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2014s` then
! find out if it starts in first or second half 
	if `time_start le 126` then
! do first half first
let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 ! start at second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let time_end_s=`time_end`-119
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	elif `time_start gt 126` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-120 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let time_end_s=`time_end`-119
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 126 loop



! starts in 2014 and ends later - same for either e2014f or e2014s
	elif `e2014f eq 0` then


		if `time_start le 126` then 

let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 126` then 

let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 126 loop


	elif `e2014s eq 0` then

		if `time_start le 126` then 

let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 126` then 

let time_start_s=`time_start`-120
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 126 loop


	endif ! starts in 2014

! starts before 2014 and ends in 2014
elif `e2014f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=1]
let time_end_s=`time_end`-119
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2014s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let time_end_s=`time_end`-119
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say 2014ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2014

! starts before 2014 and ends after 2014

let t14=if (`time_start lt 119` and `time_end gt 132`) then 1 else 0
if `t14` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say all2014first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2014,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2014,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2014,i=`sumstart`:`sumend`]
say all2014second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t14


if `s2015` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2015
	if `e2015f` then
let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let time_end_s=`time_end`-131
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2015s` then
! find out if it starts in first or second half 
	if `time_start le 138` then
! do first half first
let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 ! start at second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let time_end_s=`time_end`-131
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `time_start gt 138` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-132 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let time_end_s=`time_end`-131
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 201ds
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 138 loop



! starts in 2015 and ends later - same for either e2015f or e2015s
	elif `e2015f eq 0` then


		if `time_start le 138` then 

let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 138` then 

let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 138 loop


	elif `e2015s eq 0` then

		if `time_start le 138` then 

let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 138` then 

let time_start_s=`time_start`-132
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 138 loop


	endif ! starts in 2015

! starts before 2015 and ends in 2015
elif `e2015f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=1]
let time_end_s=`time_end`-131
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2015s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let time_end_s=`time_end`-131
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say 2015ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2015

! starts before 2015 and ends after 2015

let t15=if (`time_start lt 131` and `time_end gt 144`) then 1 else 0
if `t15` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say all2015first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2015,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2015,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2015,i=`sumstart`:`sumend`]
say all2015second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t15


if `s2016` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2016
	if `e2016f` then
let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let time_end_s=`time_end`-143
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2016s` then
! find out if it starts in first or second half 
	if `time_start le 150` then
! do first half first
let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 ! start at second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let time_end_s=`time_end`-143
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `time_start gt 150` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-144 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let time_end_s=`time_end`-143
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016ds
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 150 loop



! starts in 2016 and ends later - same for either e2016f or e2016s
	elif `e2016f eq 0` then


		if `time_start le 150` then 

let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 150` then 

let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 150 loop


	elif `e2016s eq 0` then

		if `time_start le 150` then 

let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 150` then 

let time_start_s=`time_start`-144
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 150 loop


	endif ! starts in 2016

! starts before 2016 and ends in 2016
elif `e2016f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=1]
let time_end_s=`time_end`-143
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2016s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let time_end_s=`time_end`-143
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say 2016ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2016

! starts before 2016 and ends after 2016

let t16=if (`time_start lt 143` and `time_end gt 156`) then 1 else 0
if `t16` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]

let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say all2016first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2016,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2016,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2016,i=`sumstart`:`sumend`]
say all2016second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t16


if `s2017` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2017
	if `e2017f` then
let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let time_end_s=`time_end`-155
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2017s` then
! find out if it starts in first or second half 
	if `time_start le 162` then
! do first half first
let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 
! start at second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let time_end_s=`time_end`-155
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `time_start gt 162` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-156 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let time_end_s=`time_end`-155
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017ds
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 162 loop



! starts in 2017 and ends later - same for either e2017f or e2017s
	elif `e2017f eq 0` then


		if `time_start le 162` then 

let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 162` then 

let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 162 loop


	elif `e2017s eq 0` then

		if `time_start le 162` then 

let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 162` then 

let time_start_s=`time_start`-156
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 162 loop


	endif ! starts in 2017

! starts before 2017 and ends in 2017
elif `e2017f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=1]
let time_end_s=`time_end`-155
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2017s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let time_end_s=`time_end`-155
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say 2017ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2017

! starts before 2017 and ends after 2017

let t17=if (`time_start lt 155` and `time_end gt 168`) then 1 else 0
if `t17` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say all2017first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2017,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2017,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2017,i=`sumstart`:`sumend`]
say all2017second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t17





if `s2018` then
! now there are too many measurements to get a whole year in one plot, so must 
! split it up


! start and end in 2018
	if `e2018f` then
let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let time_end_s=`time_end`-167
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `e2018s` then
! find out if it starts in first or second half 
	if `time_start le 174` then
! do first half first
let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let time_end_s=7 ! do entire first half of year
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018f
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! do second half
let time_start_s=7 
! start at second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let time_end_s=`time_end`-167
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018s
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

	elif `time_start gt 174` then
! starts and ends in second half.  only have to do second half
let time_start_s=`time_start`-168 ! subtract 6 to get to second half of year
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let time_end_s=`time_end`-167
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018ds
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


	endif ! time_start 174 loop



! starts in 2018 and ends later - same for either e2018f or e2018s
	elif `e2018f eq 0` then


		if `time_start le 174` then 

let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 174` then 

let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 162 loop


	elif `e2018s eq 0` then

		if `time_start le 174` then 

let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plusfirst 
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

		elif `time_start gt 174` then 

let time_start_s=`time_start`-168
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018plussecond
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
		endif ! time start le 174 loop


	endif ! starts in 2018

! starts before 2018 and ends in 2018
elif `e2018f` then

let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=1]
let time_end_s=`time_end`-167
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

elif `e2018s` then 

! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018endf
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! then do second half
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let time_end_s=`time_end`-167
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say 2018ends
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


endif !s2018

! starts before 2018 and ends after 2018

let t18=if (`time_start lt 167` and `time_end gt 174`) then 1 else 0
if `t18` then 
! do first half first
let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=1]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say all2018first
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

let sumstart=fl_month[d=argo_list_map_sum_psal_2018,i=7]
let sumend=fl_month[d=argo_list_map_sum_psal_2018,i=13]
let lonx=xlon[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_psal_2018,i=`sumstart`:`sumend`]
say all2018second
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t18





! now use etopo6 to add bathymetry
! check zonal range to see if using high or low land resoluion
let rcheck=`lon_max`-`lon_min`
if `rcheck gt 120` then
if `iland eq 0` then
use etopo12_grid.nc
shade/over/nokey/nolab/pal=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/lev=(0,10000,10000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
shade/over/nokey/nolab/palette=lightgreyscale/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-1000,0,1000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
shade/over/nokey/nolab/pal=grey/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/lev=(-200,0,200) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif
if `iland eq 1` then
use etopo12_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif
endif
if `rcheck le 120` then
if `iland eq 0` then
use etopo6_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
fill/over/nokey/nolab/palette=lightgreyscale/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-1000,0,1000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
fill/over/nokey/nolab/palette=grey/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-200,0,200) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
else
use etopo6_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif

endif




say Plot is finished
say To adjust parameters, use the command line to pass arguments
say For this file, the form to pass arguments is
say go fer_station_loc time_start time_end depth_min depth_max lat_min lat_max lon_min lon_max