GFDL - Geophysical Fluid Dynamics Laboratory

clear all
%input1=[‘ocean_topog_Denmark_width_5_depth_850_sink_near_29w_65n_Gin_new.nc’];
input1=[‘ocean_topog_March_25_2010_Denmark.nc’];
ncload(input1);

vgrid=[‘ocean_vgrid.nc’];
ncload(vgrid);
filename=’ocean_topog_April_11_2010_nolandseachange.nc’ ;
nc=netcdf(filename,’clobber’);
nc.description = ‘1/4 degree topo, changes made on recomendations from Keith Dixon and Whit Anderson. ITF=true TARTAR=true USEAST=true CARIBE=true ANTARCTICA=true ALEUTIANS=true REDSEA=true ‘;
nc.author = ‘Whit Anderson’;
nc.date = date;
% Define dimensions.
nc(‘ntiles’)=1;
nc(‘nx’) = 1440;
nc(‘ny’) = 1070;
epsilon = 1.E-4;

depth_orig=DEPTH;

ITF=true
TARTAR=true
USEAST=true
CARIBE=true
ANTARCTICA=false
ALEUTIANS=true
REDSEA=true
landchanges=true
DRAKE=true
April_landsea=false

if ITF 
%ITF INNER BANDA SEA

DEPTH(460,195)=297.0;
DEPTH(462,198)=297.0;
DEPTH(468,200)=4073.0;
DEPTH(469,200)=4073.0;
DEPTH(468,201)=4073.0;
DEPTH(469,201)=4073.0;
DEPTH(468,202)=4073.0;
DEPTH(469,202)=4073.0;

%ITF BANDA SEA
DEPTH(466,175)=152.0;
DEPTH(467,175)=852.0;
DEPTH(466,177)=280.0;
DEPTH(468,184)=3275.0;
DEPTH(469,184)=3275.0;
DEPTH(468,185)=3275.0;
DEPTH(469,185)=3275.0;
DEPTH(470,187)=1039.0;
DEPTH(470,188)=3275.0;
DEPTH(470,189)=3275.0;
DEPTH(471,189)=3275.0;
DEPTH(475,190)=852.0;
DEPTH(475,191)=852.0;
DEPTH(475,192)=3275.0;
DEPTH(466,191)=616.0;

%ITF LIFAMATOLA
DEPTH(477,185)=3030.0;
DEPTH(478,185)=3030.0;
DEPTH(479,185)=3030.0;
DEPTH(477,186)=3030.0;
DEPTH(481,187)=1916.0;
DEPTH(481,188)=1916.0;
DEPTH(482,189)=570.0;
DEPTH(491,186)=2254.0;
DEPTH(493,186)=2254.0;
DEPTH(491,187)=2254.0;
DEPTH(492,187)=2254.0;
DEPTH(493,187)=2254.0;
DEPTH(494,187)=2254.0;
DEPTH(495,187)=2254.0;
DEPTH(496,187)=2254.0;
DEPTH(492,188)=2254.0;
DEPTH(493,188)=2254.0;
DEPTH(494,188)=2254.0;
DEPTH(495,188)=2254.0;
DEPTH(496,188)=2254.0;
DEPTH(497,188)=2254.0;
DEPTH(498,188)=2254.0;
DEPTH(499,188)=2254.0;
DEPTH(495,189)=2254.0;
DEPTH(496,189)=2254.0;
DEPTH(497,189)=2254.0;

DEPTH(479,184)=    127.0;

DEPTH(482,194)=592.0;
DEPTH(483,191)=592.0;
DEPTH(483,192)=592.0;
DEPTH(483,193)=592.0;
DEPTH(483,194)=592.0;
DEPTH(484,191)=592.0;
DEPTH(484,192)=592.0;
DEPTH(484,193)=592.0;
DEPTH(484,194)=592.0;
DEPTH(485,191)=592.0;
DEPTH(485,192)=592.0;
         
       
if landchanges
DEPTH(480,184)=    0.0;
DEPTH(481,178)=    0.0 ;
DEPTH(481,179)=    0.0;
DEPTH(481,180)=    0.0;
DEPTH(481,181)=    0.0;
DEPTH(481,182)=    0.0;
DEPTH(481,183)=    0.0;
DEPTH(481,184)=    0.0;
DEPTH(481,185)=    0.0;
DEPTH(482,177)=    0.0;
end

%MAKASSAR STRAITS
DEPTH(490,156)=1368.0;
DEPTH(491,156)=977.0;
DEPTH(492,156)=672.0;
DEPTH(487,157)=2413.0;
DEPTH(488,157)=2413.0;
DEPTH(489,157)=2413.0;
DEPTH(490,157)=2413.0;
DEPTH(491,157)=2413.0;
DEPTH(487,158)=1368.0;
DEPTH(488,158)=2413.0;
DEPTH(489,158)=2413.0;
DEPTH(490,158)=2413.0;
DEPTH(489,159)=1368.0;
DEPTH(490,159)=2413.0;

DEPTH(475,154)=2050.0;
DEPTH(476,154)=2050.0;
DEPTH(477,154)=2050.0;
DEPTH(478,154)=2050.0;
DEPTH(479,154)=2050.0;
DEPTH(475,155)=2050.0;
DEPTH(476,155)=2050.0;
DEPTH(477,155)=2050.0;
DEPTH(478,155)=2050.0;
DEPTH(479,155)=2050.0;
DEPTH(480,155)=2050.0;
DEPTH(481,155)=2050.0;

DEPTH(460,147)=503.0;
DEPTH(461,147)=503.0;
DEPTH(462,147)=503.0;
DEPTH(463,147)=503.0;
DEPTH(464,147)=503.0;
DEPTH(460,148)=503.0;
DEPTH(461,148)=503.0;
DEPTH(462,148)=503.0;
DEPTH(463,148)=503.0;
DEPTH(464,148)=503.0;
DEPTH(465,148)=503.0;
DEPTH(466,148)=503.0;
DEPTH(467,148)=503.0;
DEPTH(463,149)=503.0;
DEPTH(464,149)=503.0;
DEPTH(465,149)=503.0;
DEPTH(466,149)=503.0;
DEPTH(467,149)=503.0;
DEPTH(468,149)=503.0;
DEPTH(469,149)=503.0;
DEPTH(456,150)=503.0;
DEPTH(466,150)=503.0;
DEPTH(467,150)=503.0;
DEPTH(466,151)=503.0;

DEPTH(459,152)=61.0;
DEPTH(459,153)=82.0;
DEPTH(460,153)=319.0;
DEPTH(461,153)=319.0;
DEPTH(460,154)=319.0;
DEPTH(461,154)=319.0;
DEPTH(466,154)=483.0;
DEPTH(467,154)=672.0;
DEPTH(468,154)=672.0;
DEPTH(461,155)=319.0;
DEPTH(462,155)=82.0;
DEPTH(464,155)=483.0;
DEPTH(465,155)=672.0;
DEPTH(467,155)=672.0;
DEPTH(468,155)=672.0;
DEPTH(463,156)=483.0;
DEPTH(464,156)=672.0;
DEPTH(465,156)=672.0;
DEPTH(466,156)=672.0;
DEPTH(467,156)=672.0;
DEPTH(468,156)=672.0;
DEPTH(469,156)=672.0;
DEPTH(470,156)=894.0;
DEPTH(464,157)=672.0;
DEPTH(465,157)=672.0;
DEPTH(466,157)=672.0;
DEPTH(467,157)=483.0;
DEPTH(465,158)=672.0;
DEPTH(465,159)=672.0;

%ITF:Sangihe Ridge
DEPTH(498,182)=1368.0;
DEPTH(499,182)=1368.0;
DEPTH(510,182)=1368.0;
DEPTH(499,183)=1368.0;
DEPTH(500,183)=1043.0;
DEPTH(508,183)=1043.0;
DEPTH(509,183)=1368.0;
DEPTH(510,183)=1368.0;

if landchanges
%ITF:Halmahera sea
DEPTH(488,190)=0.0;
DEPTH(489,191)=0.0;
DEPTH(490,191)=0.0;
end

DEPTH(483,195)=592.0;
DEPTH(484,195)=592.0;
DEPTH(485,195)=592.0;
DEPTH(483,196)=592.0;
DEPTH(490,196)=592.0;
DEPTH(491,197)=592.0;
DEPTH(483,197)=249.0;
DEPTH(490,197)=592.0;
DEPTH(491,199)=225.0;
DEPTH(491,200)=249.0;

%ITF:INNER BANDA
DEPTH(460,195)=297.0;
DEPTH(462,198)=297.0;
DEPTH(468,200)=4073.0;
DEPTH(469,200)=4073.0;
DEPTH(469,201)=4073.0;
DEPTH(468,202)=4073.0;
DEPTH(469,202)=4073.0;

%ITF:Leti strait
DEPTH(454,189)=1107.0;
DEPTH(455,189)=1107.0;
DEPTH(454,190)=1107.0;
DEPTH(455,190)=1107.0;
DEPTH(456,196)=1271.0;
DEPTH(455,197)=1271.0;
DEPTH(456,197)=1271.0;
DEPTH(454,198)=1271.0;
DEPTH(455,198)=1271.0;
DEPTH(456,198)=1271.0;
DEPTH(457,198)=1271.0;

%ITF:Lombok strait
DEPTH(453,143)=371.0;
DEPTH(454,143)=298.0;
DEPTH(455,143)=298.0;
DEPTH(453,144)=298.0;
DEPTH(454,144)=298.0;

%ITF:Sumba
DEPTH(445,165)=1132.0;
DEPTH(445,166)=1132.0;
DEPTH(451,157)=896.0;
DEPTH(451,158)=896.0;
DEPTH(451,159)=896.0;
DEPTH(451,160)=896.0;
DEPTH(451,161)=896.0;
DEPTH(451,162)=896.0;
DEPTH(452,157)=896.0;
DEPTH(452,158)=896.0;
DEPTH(452,159)=896.0;
DEPTH(452,160)=896.0;
DEPTH(452,161)=896.0;
DEPTH(452,162)=896.0;
DEPTH(453,157)=896.0;
DEPTH(453,158)=40.0;
DEPTH(453,159)=40.0;
DEPTH(453,160)=40.0;

if landchanges
DEPTH(453,159)=0.0;

%ITF:OMBAI
DEPTH(455,173)=0.0;
DEPTH(455,174)=0.0;
DEPTH(455,175)=0.0;
DEPTH(455,176)=0.0;
DEPTH(455,177)=0.0;

end
DEPTH(452,177)=3095.0;
DEPTH(452,178)=3095.0;
DEPTH(453,173)=3095.0;
DEPTH(453,174)=3095.0;
DEPTH(453,175)=3095.0;
DEPTH(453,176)=3095.0;
DEPTH(453,177)=3095.0;
DEPTH(453,178)=3095.0;
DEPTH(453,179)=3095.0;
DEPTH(453,180)=3095.0;
DEPTH(454,177)=3095.0;
DEPTH(454,178)=3095.0;
DEPTH(454,179)=3095.0;
DEPTH(454,180)=3095.0;
DEPTH(455,179)=3095.0;
DEPTH(455,180)=3095.0;
DEPTH(456,179)=3095.0;
DEPTH(456,180)=3095.0;

DEPTH(445,169)=1111.0;
DEPTH(446,169)=1111.0;
DEPTH(447,169)=1111.0;
DEPTH(444,170)=1111.0;
DEPTH(445,170)=1111.0;
DEPTH(446,170)=1111.0;

DEPTH(441,170)=1890.0;
DEPTH(441,171)=1890.0;
DEPTH(441,172)=1890.0;
DEPTH(442,170)=1890.0;
DEPTH(442,171)=1890.0;
DEPTH(442,172)=1890.0;
DEPTH(442,173)=1890.0;
DEPTH(442,174)=1890.0;
DEPTH(443,171)=1890.0;
DEPTH(443,172)=1890.0;
DEPTH(443,173)=1890.0;
DEPTH(443,174)=1890.0;
DEPTH(443,175)=1890.0;
DEPTH(443,176)=1890.0;
DEPTH(444,173)=1890.0;
DEPTH(444,174)=1890.0;
DEPTH(444,175)=1890.0;
DEPTH(444,176)=1890.0;
DEPTH(444,177)=1630.0;
DEPTH(444,178)=1630.0;
DEPTH(445,179)=1740.0;
DEPTH(445,180)=1740.0;
DEPTH(446,181)=1740.0;
DEPTH(447,182)=1740.0;
DEPTH(448,185)=1628.0;
DEPTH(449,181)=1740.0;
DEPTH(450,184)=1628.0;

end

if CARIBE 

%CARIBE:Cayman Ridge
DEPTH(559,772)=4020.0;
DEPTH(560,772)=4020.0;
DEPTH(561,772)=4020.0;
DEPTH(562,772)=4020.0;
DEPTH(559,773)=4020.0;
DEPTH(560,773)=4020.0;
DEPTH(561,773)=4020.0;
DEPTH(562,773)=4020.0;
DEPTH(563,773)=4020.0;
DEPTH(564,773)=3558.0;
DEPTH(558,774)=4020.0;
DEPTH(559,774)=4020.0;
DEPTH(560,774)=4020.0;
DEPTH(568,774)=3558.0;
DEPTH(569,774)=2491.0;
DEPTH(570,774)=1445.0;
DEPTH(560,775)=4020.0;
DEPTH(571,775)=1650.0;
DEPTH(572,775)=1650.0;
DEPTH(573,776)=1650.0;
DEPTH(575,777)=1892.0;
DEPTH(570,805)=1650.0;
DEPTH(569,806)=1650.0;
DEPTH(572,806)=1650.0;
DEPTH(569,807)=1650.0;

%CARIBE:Yucatan
DEPTH(577,778)=2050.0;
DEPTH(578,778)=2050.0;
DEPTH(579,778)=2050.0;
DEPTH(580,778)=2050.0;
DEPTH(581,778)=2050.0;
DEPTH(577,779)=2050.0;
DEPTH(578,779)=2050.0;
DEPTH(579,779)=2050.0;
DEPTH(583,779)=2050.0;

DEPTH(580,777)=1628.0;
DEPTH(579,780)=1628.0;

%CARIBE:NICARAGUAN RISE
DEPTH(557,801)=1628.0;
DEPTH(558,801)=1628.0;
DEPTH(559,801)=1628.0;
DEPTH(555,802)=1628.0;
DEPTH(556,802)=1628.0;
DEPTH(557,802)=1628.0;
DEPTH(558,802)=1628.0;
DEPTH(559,802)=1628.0;
DEPTH(554,803)=766.0;
DEPTH(555,803)=1628.0;
DEPTH(556,803)=1628.0;
DEPTH(557,803)=1628.0;
DEPTH(558,803)=1628.0;
DEPTH(553,804)=1628.0;
DEPTH(554,804)=1628.0;
DEPTH(555,804)=1628.0;
DEPTH(556,804)=1628.0;
DEPTH(553,805)=1628.0;
DEPTH(554,805)=1628.0;
DEPTH(555,805)=1628.0;
DEPTH(556,805)=1628.0;
DEPTH(553,806)=1628.0;
DEPTH(553,807)=1628.0;
DEPTH(559,807)=646.0;
DEPTH(559,808)=646.0;
DEPTH(560,808)=646.0;
DEPTH(561,808)=646.0;
DEPTH(559,809)=646.0;
DEPTH(560,809)=646.0;
DEPTH(559,810)=646.0;
DEPTH(560,810)=646.0;
DEPTH(559,807)=646.0;
DEPTH(559,811)=646.0;
DEPTH(560,811)=646.0;
DEPTH(560,812)=646.0;
DEPTH(561,819)=1445.0;
DEPTH(562,819)=1445.0;
DEPTH(563,819)=1445.0;
DEPTH(561,820)=1445.0;
DEPTH(562,820)=1445.0;
DEPTH(563,820)=1445.0;

%Caribe:BEATA RIDGE
DEPTH(548,826)=3481.0;
DEPTH(549,826)=4381.0;
DEPTH(550,826)=3864.0;

DEPTH(546,827)=3873.0;
DEPTH(549,827)=3044.0;
DEPTH(550,827)=3044.0;
DEPTH(551,827)=3481.0;

DEPTH(545,828)=4005.0;
DEPTH(546,828)=3873.0;
DEPTH(552,828)=3481.0;
DEPTH(543,829)=4005.0;
DEPTH(544,829)=4005.0;
DEPTH(545,829)=4005.0;
DEPTH(546,829)=3873.0;
DEPTH(553,829)=3284.0;
DEPTH(543,830)=4005.0;
DEPTH(544,830)=4005.0;
DEPTH(544,831)=4005.0;
DEPTH(545,831)=4005.0;
DEPTH(546,831)=3878.0;
DEPTH(552,831)=3284.0;
DEPTH(553,831)=3284.0;
DEPTH(554,831)=3284.0;
DEPTH(555,831)=3284.0;
DEPTH(544,832)=4005.0;
DEPTH(545,832)=4005.0;
DEPTH(546,832)=3873.0;
DEPTH(544,833)=4005.0;
DEPTH(545,833)=4005.0;
DEPTH(546,833)=3873.0;
DEPTH(544,834)=4005.0;
DEPTH(545,834)=4005.0;
DEPTH(546,834)=3873.0;
DEPTH(545,835)=4005.0;
DEPTH(557,833)=2674.0;
DEPTH(558,834)=2674.0;

%CARIBE:AVES RIDGE
DEPTH(538,866)=746.0;
DEPTH(543,868)=446.0;
DEPTH(545,867)=1296.0;
DEPTH(546,866)=1296.0;
DEPTH(547,866)=1586.0;
DEPTH(548,866)=1586.0;
DEPTH(548,869)=1586.0;
DEPTH(549,866)=1479.0;
DEPTH(550,866)=1296.0;
DEPTH(553,866)=1826.0;
DEPTH(554,866)=1826.0;
DEPTH(554,867)=1826.0;
DEPTH(554,868)=1826.0;
DEPTH(554,869)=1826.0;
DEPTH(555,865)=1472.0;
DEPTH(555,868)=1472.0;

%CARIBE:GRENADDA
DEPTH(536,873)=740.0;
DEPTH(536,874)=740.0;
DEPTH(537,874)=740.0;

%CARIBE:St Lucia
DEPTH(542,876)=890.0;
DEPTH(543,877)=890.0;
DEPTH(544,877)=890.0;
DEPTH(542,876)=890.0;
DEPTH(542,876)=890.0;
DEPTH(542,876)=890.0;

%CARIBE:Guadoloupe
DEPTH(556,874)=503.0;
DEPTH(557,874)=503.0;
DEPTH(551,875)=950.0;
DEPTH(552,875)=950.0;
DEPTH(551,876)=950.0;
DEPTH(552,876)=950.0;
DEPTH(549,877)=321.0;

%CARIBE:Anegada
DEPTH(560,859)=1766.0;
DEPTH(561,859)=1766.0;
DEPTH(562,859)=1766.0;
DEPTH(560,860)=1766.0;
DEPTH(560,863)=1861.0;
DEPTH(560,864)=1861.0;
DEPTH(561,862)=1959.0;
DEPTH(562,862)=1959.0;
DEPTH(561,863)=1959.0;
DEPTH(562,863)=1959.0;
DEPTH(561,864)=1959.0;
DEPTH(562,864)=1959.0;
DEPTH(563,864)=1959.0;
DEPTH(564,864)=1959.0;
DEPTH(561,865)=1959.0;
DEPTH(562,865)=1959.0;
DEPTH(563,865)=1959.0;
DEPTH(564,865)=1959.0;
DEPTH(561,866)=1766.0;
DEPTH(562,866)=1959.0;
DEPTH(563,866)=1959.0;
DEPTH(564,866)=1959.0;

%CARIBE:Mona
DEPTH(562,849)=459.0;
DEPTH(562,850)=459.0;

%Caribe:Windward
DEPTH(582,821)=1683.0;
DEPTH(583,821)=1683.0;
DEPTH(584,821)=1683.0;
DEPTH(582,822)=1683.0;
DEPTH(561,819)=1492.0;
DEPTH(562,819)=1492.0;
DEPTH(563,819)=1492.0;
DEPTH(560,820)=1492.0;
DEPTH(561,820)=1492.0;
DEPTH(562,820)=1492.0;
DEPTH(563,820)=1492.0;
DEPTH(561,821)=1492.0;
DEPTH(562,821)=1492.0;
DEPTH(570,823)=1492.0;
DEPTH(567,824)=1683.0;
DEPTH(571,824)=1189.0;
DEPTH(572,824)=1683.0;
DEPTH(568,825)=1683.0;
DEPTH(568,826)=1683.0;
DEPTH(578,825)=1492.0;
DEPTH(578,826)=1683.0;
DEPTH(579,826)=1683.0;
DEPTH(579,827)=1683.0;
DEPTH(580,827)=1683.0;
DEPTH(581,827)=1683.0;
DEPTH(576,828)=1189.0;
DEPTH(579,828)=1683.0;
DEPTH(580,828)=1683.0;
DEPTH(571,829)=1189.0;
DEPTH(580,829)=977.0;
DEPTH(573,834)=1683.0;
DEPTH(573,835)=1683.0;
DEPTH(573,836)=1683.0;
DEPTH(573,837)=1683.0;
DEPTH(573,838)=1683.0;
DEPTH(574,838)=977.0;
DEPTH(570,836)=977.0;
DEPTH(570,837)=1492.0;

%CARIBE:Old Bahamas Channel
DEPTH(584,798)=262.0;
DEPTH(586,799)=262.0;
DEPTH(584,801)=262.0;
DEPTH(585,801)=262.0;
DEPTH(585,802)=262.0;
DEPTH(586,802)=356.0;
DEPTH(582,804)=356.0;
DEPTH(587,804)=356.0;
DEPTH(582,805)=356.0;
DEPTH(582,806)=356.0;
DEPTH(583,806)=356.0;
DEPTH(582,807)=356.0;
DEPTH(583,807)=356.0;
DEPTH(583,808)=356.0;
DEPTH(582,808)=356.0;
DEPTH(584,808)=302.0;
DEPTH(580,809)=356.0;
DEPTH(581,809)=356.0;
DEPTH(582,809)=356.0;
DEPTH(583,809)=356.0;
DEPTH(579,810)=356.0;
DEPTH(580,810)=356.0;
DEPTH(581,810)=356.0;
DEPTH(582,810)=356.0;
DEPTH(583,810)=356.0;
DEPTH(579,811)=356.0;
DEPTH(580,811)=356.0;
DEPTH(581,811)=356.0;

%CARIBE:Providence Channel
DEPTH(597,804)=645.0;
DEPTH(598,804)=645.0;
DEPTH(597,805)=645.0;
DEPTH(598,805)=645.0;
DEPTH(596,807)=645.0;
DEPTH(596,808)=645.0;
DEPTH(594,809)=645.0;
DEPTH(595,809)=645.0;
DEPTH(596,809)=645.0;
DEPTH(597,809)=645.0;
DEPTH(598,809)=645.0;
DEPTH(597,810)=645.0;
DEPTH(591,811)=300.0;
DEPTH(592,811)=300.0;

%CARIBE:Florida-Cuba
DEPTH(588,790)=735.0;
DEPTH(588,791)=735.0;
DEPTH(588,792)=735.0;
DEPTH(588,793)=735.0;
DEPTH(588,794)=735.0;
DEPTH(584,792)=503.0;
DEPTH(584,793)=503.0;
DEPTH(589,796)=356.0;
DEPTH(589,797)=735.0;
DEPTH(590,798)=735.0;
DEPTH(578,799)=735.0;
DEPTH(591,799)=735.0;
DEPTH(588,800)=735.0;
DEPTH(592,800)=735.0;
DEPTH(589,801)=735.0;
DEPTH(589,802)=735.0;

%CARIBE:Yucatan
DEPTH(577,778)=2050.0;
DEPTH(578,778)=2050.0;
DEPTH(579,778)=2050.0;
DEPTH(580,778)=2050.0;
DEPTH(581,778)=2050.0;
DEPTH(577,779)=2050.0;
DEPTH(578,779)=2050.0;
DEPTH(579,779)=2050.0;
DEPTH(583,779)=2050.0;

end

if ALEUTIANS 

DEPTH(741,362)=1960.0;
DEPTH(742,361)=1960.0;
DEPTH(742,362)=1960.0;
DEPTH(743,361)=1960.0;
DEPTH(747,361)=3462.0;
DEPTH(747,366)=3462.0;
DEPTH(747,367)=3462.0;
DEPTH(748,365)=3462.0;
DEPTH(749,358)=3462.0;
DEPTH(749,362)=3462.0;
DEPTH(750,358)=3462.0;
DEPTH(750,359)=3462.0;
DEPTH(751,358)=3462.0;
DEPTH(751,359)=3462.0;

end

if USEAST 
if landchanges
%Cape Cod
DEPTH(672,839)=0.0;
DEPTH(672,840)=0.0;
DEPTH(672,841)=0.0;
DEPTH(672,842)=0.0;
DEPTH(673,842)=0.0;

%Long Island
DEPTH(668,830)=0.0;

%Chesapeake
DEPTH(652,818)=0.0;

%Pamlico Sound
DEPTH(640,814)=0.0;
DEPTH(639,814)=0.0;
DEPTH(638,814)=0.0;
DEPTH(640,815)=0.0;
DEPTH(639,815)=0.0;
DEPTH(640,816)=0.0;

end
end

%Tartar
if landchanges
DEPTH(732,247)=0.0;
end

if REDSEA
DEPTH(538,1293)=138.5;
DEPTH(538,1294)=138.5;
DEPTH(539,1293)=138.5;
DEPTH(539,1294)=138.5;
DEPTH(539,1295)=138.5;
DEPTH(540,1292)=138.5;
DEPTH(540,1293)=138.5;
DEPTH(540,1294)=138.5;
DEPTH(541,1291)=138.5;
DEPTH(541,1292)=138.5;
DEPTH(541,1293)=138.5;
DEPTH(541,1294)=138.5;
DEPTH(542,1291)=138.5;
DEPTH(542,1292)=138.5;
DEPTH(542,1293)=138.5;
DEPTH(543,1290)=138.5;
DEPTH(543,1291)=138.5;
DEPTH(543,1292)=138.5;
DEPTH(543,1293)=138.5;
DEPTH(544,1290)=138.5;
DEPTH(544,1291)=138.5;
DEPTH(544,1292)=138.5;
DEPTH(545,1290)=138.5;
DEPTH(545,1291)=138.5;
DEPTH(546,1290)=138.5;
DEPTH(546,1291)=138.5;
DEPTH(547,1289)=138.5;
DEPTH(548,1288)=138.5;
DEPTH(548,1291)=138.5;

DEPTH(549,1288)=182.7;
DEPTH(551,1287)=182.7;
DEPTH(552,1286)=182.7;
DEPTH(552,1289)=182.7;
DEPTH(556,1287)=182.7;
DEPTH(557,1286)=182.7;
DEPTH(558,1285)=182.7;
DEPTH(555,1287)=428.8;

DEPTH(561,1280)=819.0;
DEPTH(562,1279)=819.0;
DEPTH(562,1282)=819.0;
DEPTH(564,1278)=819.0;
DEPTH(565,1277)=819.0;
DEPTH(567,1275)=859.0;
DEPTH(568,1278)=859.0;
DEPTH(570,1276)=859.0;
DEPTH(574,1272)=867.0;
DEPTH(580,1270)=867.0;
DEPTH(581,1273)=867.0;
DEPTH(585,1271)=867.0;
DEPTH(592,1264)=845.0;
DEPTH(593,1264)=845.0;
DEPTH(598,1260)=845.0;
DEPTH(603,1257)=40.0;
DEPTH(604,1257)=40.0;

DEPTH(535,1328)=876.0;
DEPTH(538,1328)=876.0;
DEPTH(538,1329)=876.0;

end

if ANTARCTICA
if landchanges

DEPTH(121,139:145)=0.0;
DEPTH(122,139:145)=0.0;
DEPTH(123,139:145)=0.0;
DEPTH(124,139:144)=0.0;
DEPTH(125,139:143)=0.0;

DEPTH(116,257:262)=0.0;
DEPTH(117,257:262)=0.0;
DEPTH(118,257:262)=0.0;
DEPTH(119,257:262)=0.0;
DEPTH(120,257:262)=0.0;
DEPTH(121,257:262)=0.0;
DEPTH(122,257:262)=0.0;
DEPTH(123,257:262)=0.0;

end

%land bridge
DEPTH(78,848:851)=40.47;

DEPTH(86,1096:1087)=60.8;
DEPTH(82,1117)=40.47;

DEPTH(98,1182:1183)=40.47;
DEPTH(95,1187:1188)=40.47;

if landchanges
DEPTH(88,1211:1215)=0.0;
DEPTH(89,1211:1215)=0.0;
DEPTH(90,1211:1215)=0.0;
DEPTH(91,1211:1215)=0.0;
DEPTH(92,1211:1214)=0.0;

DEPTH(100,1246:1253)=0.0;
DEPTH(99,1246:1253)=0.0;
DEPTH(98,1246:1253)=0.0;
DEPTH(97,1246:1253)=0.0;
DEPTH(96,1246:1253)=0.0;
DEPTH(95,1246:1253)=0.0;
DEPTH(94,1246:1253)=0.0;

DEPTH(122,1319:1323)=0.0;
DEPTH(123,1319:1323)=0.0;
DEPTH(124,1319:1322)=0.0;
DEPTH(125,1320:1322)=0.0;
DEPTH(126,1322)=0.0;
end

end

if DRAKE
DEPTH(166,863)=3446.0;
DEPTH(166,864)=3446.0;
DEPTH(167,858)=3446.0;
DEPTH(167,862)=3446.0;
DEPTH(168,864)=3446.0;
DEPTH(169,860)=3446.0;
DEPTH(190,858)=3446.0;
DEPTH(190,859)=3446.0;
DEPTH(190,860)=3446.0;
DEPTH(191,859)=3446.0;
DEPTH(191,860)=3446.0;
DEPTH(191,861)=3446.0;
DEPTH(192,855)=3446.0;
DEPTH(192,859)=3446.0;
DEPTH(193,858)=3446.0;
DEPTH(194,856)=3446.0;

DEPTH(191,862)=3446.0;

DEPTH(180,894)=3138.0;
DEPTH(180,895)=3138.0;
DEPTH(181,894)=3138.0;

DEPTH(237,928)=3058.0;
DEPTH(237,929)=3058.0;
DEPTH(238,928)=3058.0;
DEPTH(238,929)=3058.0;

DEPTH(240,948)=3319.0;
DEPTH(242,947)=3319.0;

DEPTH(210,945)=2480.0;
DEPTH(210,946)=3321.0;
DEPTH(213,948)=3321.0;
DEPTH(214,948)=3321.0;
DEPTH(216,954)=3310.0;
DEPTH(219,952)=3310.0;
DEPTH(220,948)=2480.0;
DEPTH(220,952)=3310.0;
DEPTH(223,952)=3310.0;
DEPTH(224,952)=3310.0;

%Orkne      y
DEPTH(169,971)=3446.0;
DEPTH(169,972)=3446.0;
DEPTH(169,973)=3446.0;
DEPTH(170,971)=3446.0;
DEPTH(170,972)=3446.0;
DEPTH(170,973)=3446.0;
DEPTH(171,971)=3446.0;
DEPTH(171,972)=3446.0;
DEPTH(173,968)=3446.0;
DEPTH(174,967)=3446.0;
DEPTH(174,968)=3446.0;
DEPTH(175,966)=3446.0;
DEPTH(175,967)=3446.0;
DEPTH(175,968)=3446.0;
DEPTH(175,969)=3446.0;
DEPTH(176,966)=3446.0;
DEPTH(176,967)=3446.0;
DEPTH(176,968)=3446.0;

DEPTH(178,962)=3446.0;
DEPTH(179,962)=3446.0;
DEPTH(180,952)=3446.0;
DEPTH(180,961)=3446.0;
DEPTH(180,962)=3446.0;
DEPTH(180,963)=3446.0;
DEPTH(181,952)=3446.0;
DEPTH(181,955)=3446.0;
DEPTH(181,956)=3446.0;
DEPTH(181,961)=3446.0;
DEPTH(181,962)=3446.0;
DEPTH(181,963)=3446.0;
DEPTH(182,952)=3446.0;
DEPTH(182,953)=3446.0;
DEPTH(183,952)=3446.0;
DEPTH(183,953)=3446.0;
DEPTH(184,952)=3446.0;
DEPTH(184,953)=3446.0;
DEPTH(185,951)=3446.0;
DEPTH(185,952)=3446.0;

%Georgasagia
DEPTH(201,967)=3238.0;
DEPTH(203,965)=3238.0;
DEPTH(204,966)=3238.0;
DEPTH(204,967)=3238.0;
DEPTH(205,966)=3238.0;
DEPTH(205,967)=3238.0;
DEPTH(206,966)=3238.0;
DEPTH(206,967)=3238.0;
DEPTH(207,966)=3238.0;
DEPTH(207,967)=3238.0;
DEPTH(208,966)=3238.0;
DEPTH(208,967)=3238.0;
DEPTH(209,966)=3238.0;
DEPTH(210,966)=3238.0;
DEPTH(222,957)=3238.0;
DEPTH(222,958)=3238.0;
DEPTH(222,959)=3238.0;
DEPTH(223,957)=3238.0;
DEPTH(223,958)=3238.0;
DEPTH(223,959)=3238.0;
DEPTH(223,960)=3238.0;
DEPTH(223,961)=3238.0;
DEPTH(223,962)=3238.0;
DEPTH(224,959)=3238.0;
DEPTH(224,960)=3238.0;
DEPTH(224,961)=3238.0;
DEPTH(224,962)=3238.0;

DEPTH(205,983)=3253.0;
DEPTH(205,984)=3253.0;
DEPTH(205,985)=3253.0;
DEPTH(206,984)=3253.0;
DEPTH(206,985)=3253.0;

DEPTH(213,983)=3238.0;
DEPTH(213,986)=3238.0;
DEPTH(214,986)=3238.0;
DEPTH(214,990)=3238.0;
DEPTH(214,991)=3238.0;
DEPTH(214,992)=3238.0;
DEPTH(219,992)=3238.0;
DEPTH(219,993)=3238.0;
DEPTH(220,992)=3238.0;
DEPTH(220,993)=3238.0;
DEPTH(221,992)=3238.0;
DEPTH(221,993)=3238.0;
DEPTH(222,992)=3238.0;
DEPTH(222,993)=3238.0;

DEPTH(229,982)=3238.0;
DEPTH(229,983)=3238.0;
DEPTH(230,981)=3238.0;
DEPTH(230,982)=3238.0;
DEPTH(230,983)=3238.0;
DEPTH(230,984)=3238.0;
DEPTH(231,981)=3238.0;
DEPTH(231,982)=3238.0;
DEPTH(231,983)=3238.0;

DEPTH(174,962)=3030.0;
DEPTH(174,963)=3030.0;
DEPTH(174,964)=3030.0;
DEPTH(175,962)=3030.0;
DEPTH(182,968)=2888.0;
DEPTH(183,966)=2888.0;

%Sandwich
DEPTH(169,1031)=4716.0;
DEPTH(169,1032)=4716.0;
DEPTH(170,1030)=4547.0;
DEPTH(170,1031)=4547.0;
DEPTH(170,1031)=4547.0;
DEPTH(171,1030)=4547.0;
DEPTH(171,1031)=4547.0;
DEPTH(171,1032)=4547.0;
DEPTH(172,1030)=4547.0;
DEPTH(172,1031)=4547.0;
DEPTH(172,1032)=4547.0;
DEPTH(173,1030)=4547.0;
DEPTH(173,1031)=4547.0;
DEPTH(173,1032)=4547.0;
DEPTH(174,1030)=4547.0;
DEPTH(174,1031)=4547.0;
DEPTH(174,1032)=4547.0;

DEPTH(179,1016)=4537.0;
DEPTH(179,1017)=4537.0;
DEPTH(179,1023)=4537.0;
DEPTH(179,1024)=4537.0;
DEPTH(180,1023)=4537.0;
DEPTH(180,1024)=4537.0;
DEPTH(180,1025)=4912.0;
DEPTH(179,1017)=4537.0;
DEPTH(181,1018)=4537.0;
DEPTH(182,1020)=4537.0;
DEPTH(183,1017)=4537.0;

DEPTH(246,1008)=4549.0;
DEPTH(246,1009)=4549.0;
DEPTH(246,1010)=4549.0;
DEPTH(247,1008)=4549.0;
DEPTH(247,1009)=4549.0;
DEPTH(247,1010)=4549.0;
DEPTH(248,1008)=4549.0;
DEPTH(248,1009)=4549.0;

DEPTH(338,1005)=4383.0;
DEPTH(338,1006)=4383.0;
DEPTH(342,1005)=4383.0;
DEPTH(342,1006)=4383.0;
DEPTH(343,1005)=4383.0;
DEPTH(343,1006)=4383.0;
DEPTH(344,1005)=4383.0;
DEPTH(344,1006)=4383.0;
DEPTH(347,1015)=4218.0;
DEPTH(348,1015)=4218.0;
DEPTH(349,1014)=4218.0;
DEPTH(349,1015)=4218.0;
DEPTH(350,1014)=4218.0;
DEPTH(350,1015)=4218.0;

%Vema
DEPTH(346,961)=4600.0;
DEPTH(349,962)=4600.0;
DEPTH(349,963)=4600.0;
DEPTH(350,962)=4600.0;
DEPTH(350,963)=4600.0;
DEPTH(351,962)=4600.0;
DEPTH(351,963)=4600.0;
DEPTH(352,962)=4600.0;
DEPTH(352,963)=4600.0;
DEPTH(353,962)=4600.0;
DEPTH(354,962)=4600.0;
DEPTH(354,963)=4600.0;
DEPTH(355,962)=4600.0;
DEPTH(355,963)=4600.0;
DEPTH(356,962)=4600.0;
DEPTH(356,963)=4600.0;
DEPTH(357,962)=4600.0;
DEPTH(357,963)=4600.0;
DEPTH(357,964)=4600.0;
DEPTH(358,962)=4600.0;
DEPTH(358,963)=4600.0;
DEPTH(358,964)=4600.0;
DEPTH(359,963)=4600.0;
DEPTH(359,964)=4600.0;
DEPTH(360,963)=4600.0;
DEPTH(360,964)=4600.0;
DEPTH(361,963)=4600.0;
DEPTH(361,964)=4600.0;
DEPTH(362,963)=4600.0;
DEPTH(362,964)=4600.0;
DEPTH(363,963)=4600.0;
DEPTH(363,964)=4600.0;
DEPTH(364,963)=4600.0;
DEPTH(364,964)=4600.0;

DEPTH(365,963)=4600.0;
DEPTH(365,964)=4600.0;
DEPTH(365,965)=4600.0;
DEPTH(365,966)=4600.0;
DEPTH(366,963)=4600.0;
DEPTH(366,964)=4600.0;
DEPTH(366,965)=4600.0;
DEPTH(366,966)=4600.0;
DEPTH(367,965)=4600.0;
DEPTH(367,966)=4600.0;
DEPTH(367,967)=4600.0;
DEPTH(367,968)=4600.0;
DEPTH(368,965)=4600.0;
DEPTH(368,966)=4600.0;
DEPTH(368,967)=4600.0;
DEPTH(368,968)=4600.0;
DEPTH(369,967)=4600.0;
DEPTH(369,968)=4600.0;
DEPTH(370,967)=4600.0;
DEPTH(370,968)=4600.0;
DEPTH(371,967)=4600.0;
DEPTH(371,968)=4600.0;
DEPTH(372,968)=4600.0;
DEPTH(372,967)=4600.0;
DEPTH(372,968)=4600.0;
DEPTH(373,967)=4600.0;
DEPTH(373,968)=4600.0;

DEPTH(374,967)=4600.0;
DEPTH(374,968)=4600.0;
DEPTH(374,969)=4600.0;
DEPTH(374,970)=4600.0;
DEPTH(374,974)=4600.0;
DEPTH(375,967)=4600.0;
DEPTH(375,968)=4600.0;
DEPTH(375,969)=4600.0;
DEPTH(375,976)=4600.0;
DEPTH(376,973)=4600.0;
DEPTH(376,976)=4600.0;
DEPTH(376,977)=4600.0;
DEPTH(376,978)=4600.0;
DEPTH(376,979)=4600.0;
DEPTH(376,980)=4600.0;
DEPTH(377,975)=4600.0;
DEPTH(377,976)=4600.0;
DEPTH(377,977)=4600.0;
DEPTH(377,978)=4600.0;
DEPTH(377,979)=4600.0;
DEPTH(377,980)=4600.0;
DEPTH(377,981)=4600.0;
DEPTH(377,982)=4600.0;
DEPTH(377,983)=4600.0;
DEPTH(378,979)=4600.0;
DEPTH(378,980)=4600.0;
DEPTH(378,981)=4600.0;
DEPTH(378,982)=4600.0;
DEPTH(378,983)=4600.0;
DEPTH(379,982)=4600.0;
DEPTH(379,983)=4600.0;
DEPTH(380,982)=4600.0;
DEPTH(380,983)=4600.0;
DEPTH(380,984)=4600.0;
DEPTH(381,982)=4600.0;
DEPTH(381,983)=4600.0;

%Romanche
DEPTH(487,1040)=4338.0;
DEPTH(487,1041)=4338.0;
DEPTH(488,1040)=4338.0;
DEPTH(488,1041)=4338.0;
DEPTH(489,1040)=4338.0;
DEPTH(489,1041)=4338.0;

%Brazil Basin
DEPTH(484,986)=4403.0;
DEPTH(484,987)=4403.0;
DEPTH(487,984)=4403.0;
DEPTH(488,969)=4411.0;
DEPTH(489,967)=4411.0;
DEPTH(489,968)=4411.0;
DEPTH(490,967)=4411.0;
DEPTH(490,968)=4411.0;
DEPTH(491,962)=4411.0;
DEPTH(491,963)=4411.0;
DEPTH(491,964)=4411.0;
DEPTH(491,965)=4411.0;
DEPTH(491,968)=4411.0;
DEPTH(492,962)=4411.0;
DEPTH(492,963)=4411.0;
DEPTH(492,964)=4411.0;
DEPTH(492,965)=4411.0;
DEPTH(492,968)=4411.0;
DEPTH(493,963)=4411.0;
DEPTH(494,963)=4411.0;
DEPTH(503,965)=4432.0;
DEPTH(503,966)=4432.0;
DEPTH(504,965)=4432.0;
DEPTH(504,966)=4432.0;

%Vema
DEPTH(527,946)=4450.0;
DEPTH(527,947)=4450.0;
DEPTH(527,948)=4450.0;
DEPTH(527,949)=4450.0;
DEPTH(527,950)=4450.0;
DEPTH(528,951)=4450.0;
DEPTH(528,945)=4450.0;
DEPTH(528,946)=4450.0;
DEPTH(528,947)=4450.0;
DEPTH(528,948)=4450.0;
DEPTH(528,949)=4450.0;
DEPTH(528,950)=4450.0;
DEPTH(528,951)=4450.0;
DEPTH(529,950)=4450.0;
DEPTH(529,951)=4450.0;
DEPTH(530,950)=4450.0;
DEPTH(530,951)=4450.0;
DEPTH(530,952)=4450.0;
DEPTH(531,951)=4450.0;
DEPTH(531,952)=4450.0;
DEPTH(531,953)=4450.0;
DEPTH(532,951)=4450.0;
DEPTH(532,952)=4450.0;
DEPTH(532,953)=4450.0;
DEPTH(532,954)=4450.0;

DEPTH(532,956)=4450.0;
DEPTH(532,957)=4450.0;
DEPTH(532,958)=4450.0;
DEPTH(532,959)=4450.0;
DEPTH(532,960)=4450.0;
DEPTH(532,961)=4450.0;
DEPTH(532,962)=4450.0;
DEPTH(533,956)=4450.0;
DEPTH(533,957)=4450.0;
DEPTH(533,958)=4450.0;
DEPTH(533,959)=4450.0;
DEPTH(533,960)=4450.0;
DEPTH(533,961)=4450.0;
DEPTH(533,962)=4450.0;

%Gibbs
DEPTH(729,967)=3550.0;
DEPTH(730,971)=3550.0;
DEPTH(730,977)=3550.0;
DEPTH(731,965)=3550.0;
DEPTH(731,966)=3550.0;
DEPTH(731,971)=3550.0;
DEPTH(731,977)=3550.0;
DEPTH(731,978)=3550.0;
DEPTH(731,979)=3550.0;
DEPTH(732,970)=3550.0;
DEPTH(732,983)=3550.0;
DEPTH(732,1021)=3550.0;
DEPTH(732,1025)=3550.0;
DEPTH(733,1019)=3550.0;
DEPTH(733,1020)=3550.0;
DEPTH(733,1021)=3550.0;
DEPTH(733,1022)=3550.0;
DEPTH(733,1023)=3550.0;

DEPTH(734,972)=3550.0;
DEPTH(734,973)=3550.0;
DEPTH(734,977)=3550.0;
DEPTH(734,978)=3550.0;
DEPTH(734,979)=3550.0;
DEPTH(734,980)=3550.0;
DEPTH(734,981)=3550.0;
DEPTH(734,984)=3550.0;
DEPTH(734,985)=3550.0;
DEPTH(734,986)=3550.0;
DEPTH(734,987)=3550.0;
DEPTH(734,988)=3550.0;
DEPTH(734,989)=3550.0;
DEPTH(734,990)=3550.0;
DEPTH(734,991)=3550.0;
DEPTH(734,1010)=3550.0;
DEPTH(734,1021)=3550.0;
DEPTH(735,979)=3550.0;
DEPTH(735,980)=3550.0;
DEPTH(735,981)=3550.0;
DEPTH(735,984)=3550.0;
DEPTH(735,985)=3550.0;
DEPTH(735,986)=3550.0;
DEPTH(735,987)=3550.0;
DEPTH(735,988)=3550.0;
DEPTH(735,989)=3550.0;
DEPTH(735,990)=3550.0;
DEPTH(735,991)=3550.0;
DEPTH(735,1021)=3550.0;

DEPTH(736,974)=3550.0;
DEPTH(736,989)=3550.0;
DEPTH(736,990)=3550.0;
DEPTH(736,991)=3550.0;
DEPTH(736,997)=3550.0;
DEPTH(736,998)=3550.0;
DEPTH(736,999)=3550.0;
DEPTH(736,1000)=3550.0;
DEPTH(736,1004)=3550.0;
DEPTH(736,1013)=3550.0;
DEPTH(736,1025)=3550.0;
DEPTH(737,961)=3550.0;
DEPTH(737,962)=3550.0;
DEPTH(737,963)=3550.0;
DEPTH(737,972)=3550.0;
DEPTH(737,973)=3550.0;
DEPTH(737,989)=3550.0;
DEPTH(737,990)=3550.0;
DEPTH(737,991)=3550.0;
DEPTH(737,992)=3550.0;
DEPTH(737,993)=3550.0;
DEPTH(737,994)=3550.0;
DEPTH(737,995)=3550.0;
DEPTH(737,996)=3550.0;
DEPTH(737,997)=3550.0;
DEPTH(737,998)=3550.0;
DEPTH(737,999)=3550.0;
DEPTH(737,1000)=3550.0;
DEPTH(737,1004)=3550.0;
DEPTH(737,1019)=3550.0;
DEPTH(737,1020)=3550.0;

DEPTH(738,989)=3550.0;
DEPTH(738,990)=3550.0;
DEPTH(738,991)=3550.0;
DEPTH(738,992)=3550.0;
DEPTH(738,993)=3550.0;
DEPTH(738,994)=3550.0;
DEPTH(738,995)=3550.0;
DEPTH(738,996)=3550.0;
DEPTH(738,997)=3550.0;
DEPTH(738,998)=3550.0;
DEPTH(740,1018)=3550.0;
DEPTH(740,1019)=3550.0;
DEPTH(741,1019)=3550.0;
DEPTH(741,1020)=3550.0;
DEPTH(741,1021)=3550.0;
DEPTH(741,1022)=3550.0;
DEPTH(741,1025)=3550.0;
DEPTH(742,1020)=3550.0;
DEPTH(742,1025)=3550.0;

end

if landchanges
%Galapagos
DEPTH(487,756)=0.0;
DEPTH(488,756)=0.0;
end

%dzt fixes
DEPTH(883,1091)=DEPTH(883,1091)-10.0;
DEPTH(882,1092)=DEPTH(882,1092)-10.0;
DEPTH(882,1097)=DEPTH(882,1097)-10.0;
DEPTH(882,1100)=DEPTH(882,1100)-10.0;
DEPTH(877,1097)=DEPTH(877,1097)-10.0;
DEPTH(876,1090)=DEPTH(876,1090)-10.0;
DEPTH(872,1097)=DEPTH(872,1097)-10.0;
DEPTH(871,1099)=DEPTH(871,1099)-10.0;
DEPTH(869,1091)=DEPTH(869,1091)-10.0;
DEPTH(869,1097)=DEPTH(869,1097)-10.0;
DEPTH(866,1096)=DEPTH(866,1096)-10.0;
DEPTH(864,1091)=DEPTH(864,1091)-10.0;
DEPTH(864,1096)=DEPTH(864,1096)-10.0;
DEPTH(862,1097)=DEPTH(862,1097)-10.0;
DEPTH(856,1097)=DEPTH(856,1097)-10.0;
DEPTH(855,1090)=DEPTH(855,1090)-10.0;
DEPTH(855,1097)=DEPTH(855,1097)-10.0;
DEPTH(854,1090)=DEPTH(854,1090)-10.0;
DEPTH(854,1097)=DEPTH(854,1097)-10.0;
DEPTH(853,1094)=DEPTH(853,1094)-10.0;
DEPTH(852,1094)=DEPTH(852,1094)-10.0;
DEPTH(852,1100)=DEPTH(852,1100)-10.0;
DEPTH(851,1094)=DEPTH(851,1094)-10.0;
DEPTH(851,1100)=DEPTH(851,1100)-10.0;
DEPTH(847,1100)=DEPTH(847,1100)-10.0;
DEPTH(844,1093)=DEPTH(844,1093)-10.0;
DEPTH(843,1097)=DEPTH(843,1097)-10.0;
DEPTH(841,1095)=DEPTH(841,1095)-10.0;
DEPTH(839,1094)=DEPTH(839,1094)-10.0;
DEPTH(836,1090)=DEPTH(836,1090)-10.0;

% March additions from Keith
DEPTH(448,165)= 1111.0;
DEPTH(459,165)=  305.0;
DEPTH(460,165)=   93.0;
DEPTH(455,166)=  593.0;
DEPTH(455,167)=  305.0;
DEPTH(459,167)=  475.0;
DEPTH(455,168)=  153.0;
DEPTH(459,168)=  593.0;
DEPTH(455,169)=   40.0;
DEPTH(459,169)=  907.0;
DEPTH(454,170)=   40.0;
DEPTH(455,170)=   40.0;
DEPTH(454,171)=   40.0;
DEPTH(455,171)=   40.0;
DEPTH(453,180)= 2480.0;

DEPTH(464,155)= 672.0;
DEPTH(466,155)= 672.0;
DEPTH(456,149)= 999.0;
DEPTH(457,149)= 999.0;
DEPTH(456,150)= 999.0;
DEPTH(457,150)= 999.0;
DEPTH(456,151)= 999.0;
DEPTH(457,151)= 999.0;
DEPTH(456,152)= 999.0;
DEPTH(457,152)= 999.0;
DEPTH(458,152)= 999.0;
DEPTH(456,153)= 999.0;

DEPTH(463,157)= 672.0;

DEPTH(491,196)= 40.0;
DEPTH(488,198)= 40.0;
DEPTH(489,198)=565.0;
DEPTH(490,200)=190.0;
DEPTH(490,201)= 71.0;

DEPTH(585,787)=1959.0;
DEPTH(586,787)=1959.0;
DEPTH(587,787)=1339.0;
DEPTH(588,787)=1107.0;
DEPTH(590,787)= 190.0;

DEPTH(585,788 )=1865.0;
DEPTH(586,788 )=1865.0;
DEPTH(587,788 )=1339.0;
DEPTH(588,788 )=1107.0;
DEPTH(589,788 )= 390.0;
DEPTH(590,788 )=  71.0;
DEPTH(585,789 )=1766.0;
DEPTH(586,789 )=1766.0;
DEPTH(587,789 )=1271.0;
DEPTH(588,789 )=1107.0;
DEPTH(589,789 )= 344.0;
DEPTH(590,789 )=  71.0;
DEPTH(584,790 )=1271.0;
DEPTH(585,790 )=1766.0;
DEPTH(586,790 )=1766.0;
DEPTH(587,790 )=1271.0;
DEPTH(588,790 )=1107.0;
DEPTH(589,790 )= 344.0;
DEPTH(590,790 )=  71.0;
DEPTH(584,791 )=1271.0;
DEPTH(585,791 )=1766.0;
DEPTH(586,791 )=1766.0;
DEPTH(587,791 )=1271.0;
DEPTH(588,791 )=1107.0;
DEPTH(589,791 )= 344.0;
DEPTH(590,791 )=  71.0;
DEPTH(584,792 )=1271.0;
DEPTH(585,792 )=1646.0;
DEPTH(586,792 )=1646.0;
DEPTH(588,792 )=1107.0;
DEPTH(589,792 )= 344.0;
DEPTH(590,792 )=  71.0;
DEPTH(584,793 )=1107.0;
DEPTH(588,793 )=1107.0;
DEPTH(589,793 )= 344.0;
DEPTH(590,793 )=  71.0;
DEPTH(584,794 )=1107.0;
DEPTH(588,794 )=1107.0;
DEPTH(589,794 )= 391.0;
DEPTH(590,794 )=  71.0;
DEPTH(584,795 )= 965.0;
DEPTH(588,795 )=1107.0;
DEPTH(589,795 )= 455.0;
DEPTH(590,795 )= 190.0;
DEPTH(584,796 )= 620.0;
DEPTH(588,796 )=1107.0;
DEPTH(589,796 )= 721.0;
DEPTH(590,796 )= 268.0;
DEPTH(584,797 )= 398.0;
DEPTH(588,797 )=1107.0;
DEPTH(589,797 )= 976.0;
DEPTH(590,797 )= 721.0;
DEPTH(591,797 )= 268.0;
DEPTH(584,798 )= 536.0;
DEPTH(585,798 )= 976.0;
DEPTH(586,798 )=1107.0;
DEPTH(587,798 )=1107.0;
DEPTH(588,798 )=1107.0;
DEPTH(589,798 )=1107.0;
DEPTH(590,798 )= 976.0;
DEPTH(591,798 )= 721.0;

DEPTH(584,799)= 536.0;
DEPTH(585,799)= 536.0;
DEPTH(586,799)= 777.0;
DEPTH(587,799)=1025.0;
DEPTH(588,799)=1107.0;
DEPTH(589,799)=1107.0;
DEPTH(590,799)= 976.0;
DEPTH(591,799)= 976.0;
DEPTH(592,799)= 464.0;
DEPTH(584,800)= 536.0;
DEPTH(585,800)= 536.0;
DEPTH(588,800)= 826.0;
DEPTH(592,800)= 790.0;
DEPTH(593,800)= 299.0;
DEPTH(594,800)= 191.0;
DEPTH(595,800)= 127.0;
DEPTH(598,800)= 115.0;
DEPTH(599,800)=  92.0;
DEPTH(600,800)=  61.0;
DEPTH(584,801)= 536.0;
DEPTH(585,801)= 536.0;
DEPTH(588,801)= 826.0;
DEPTH(589,801)= 826.0;
DEPTH(604,801)= 826.0;
DEPTH(605,801)= 721.0;
DEPTH(606,801)= 620.0;
DEPTH(607,801)= 536.0;
DEPTH(584,802)= 536.0;
DEPTH(585,802)= 536.0;
DEPTH(586,802)= 536.0;
DEPTH(587,802)= 536.0;
DEPTH(588,802)= 721.0;
DEPTH(589,802)= 826.0;
DEPTH(606,802)= 826.0;
DEPTH(607,802)= 721.0;
DEPTH(584,803)= 536.0;
DEPTH(585,803)= 536.0;
DEPTH(586,803)= 536.0;
DEPTH(587,803)= 536.0;
DEPTH(588,803)= 620.0;
DEPTH(589,803)= 721.0;
DEPTH(590,803)= 826.0;
DEPTH(591,803)= 826.0;
DEPTH(588,804)= 190.0;
DEPTH(591,804)= 190.0;
DEPTH(601,804)= 269.0;
DEPTH(602,804)= 475.0;
DEPTH(603,804)= 826.0;
DEPTH(604,805)= 826.0;
DEPTH(599,806)=  92.0;

DEPTH(581,807)= 356.0;
DEPTH(581,808)= 356.0;
DEPTH(585,808)=  61.0;
DEPTH(586,808)=  40.0;
DEPTH(584,809)= 103.0;
DEPTH(585,809)=  40.0;
DEPTH(586,809)=  40.0;
DEPTH(583,810)= 103.0;
DEPTH(584,810)=  40.0;
DEPTH(583,811)=  40.0;
DEPTH(597,811)= 305.0;
DEPTH(582,812)=  40.0;
DEPTH(583,812)=  40.0;
DEPTH(582,813)=  40.0;
DEPTH(583,813)=  40.0;

if April_landsea
%First the land/sea mask changes for the Philippines:

DEPTH(545,168)= 0.0;
DEPTH(545,169)= 0.0;
DEPTH(541,176)= 0.0;
DEPTH(544,177)= 0.0;
DEPTH(529,180)= 0.0;
DEPTH(527,183)= 0.0;

%…and for the Lesser Sunda Islands
DEPTH( 454,169)= 0.0;
DEPTH( 455,170)= 0.0;
DEPTH( 455,171)= 0.0;

end
%A set of mods to better establish an observed gap in the ridge between Jamaica and Nicaragua.

DEPTH(551,796)= 349.0;
DEPTH(552,796)= 540.0;
DEPTH(553,796)= 540.0;
DEPTH(554,796)= 540.0;
DEPTH(555,796)= 540.0;
DEPTH(556,796)= 540.0;
DEPTH(557,796)= 540.0;
DEPTH(553,797)= 540.0;
DEPTH(554,797)= 540.0;
DEPTH(555,797)= 540.0;
DEPTH(556,797)= 540.0;
DEPTH(557,797)= 540.0;
DEPTH(558,797)= 540.0;
DEPTH(554,798)= 540.0;
DEPTH(557,805)= 721.0;

%Our earlier mod got the sill depth correct at the western end of the channel, but we could have done a better job setting up the contours of the eastern entrance. ÊHere’s the mods to do so.
DEPTH(596,808 )=956.0;
DEPTH(595,809 )=956.0;
DEPTH(596,809 )=956.0;
DEPTH(597,809 )=956.0;
DEPTH(598,809 )=956.0;
DEPTH(595,810 )=1492.0;
DEPTH(596,810 )=1492.0;
DEPTH(597,810 )=956.0;
DEPTH(596,811 )=1766.0;
DEPTH(597,811 )=481.0;
DEPTH(596,812 )=1766.0;
DEPTH(594,814 )=645.0;
DEPTH(593,815 )=481.0;

%Florida Shave

DEPTH(590,786)= 346.0;
DEPTH(591,786)= 114.0;
DEPTH(589,787)= 618.0;
DEPTH(589,788)= 541.0;
DEPTH(590,788)= 155.0;
DEPTH(589,789)= 541.0;
DEPTH(590,789)= 139.0;
DEPTH(589,790)= 541.0;
DEPTH(590,790)= 171.0;
DEPTH(589,791)= 541.0;
DEPTH(590,791)= 190.0;
DEPTH(589,792)= 541.0;
DEPTH(590,792)= 190.0;
DEPTH(589,793)= 541.0;
DEPTH(590,793)= 190.0;
DEPTH(589,794)= 541.0;
DEPTH(590,794)= 190.0;
DEPTH(589,795)= 670.0;
DEPTH(590,795)= 268.0;
DEPTH(590,796)= 555.0;
DEPTH(591,796)= 169.0;
DEPTH(591,797)= 474.0;
DEPTH(592,797)= 169.0;

DEPTH(592,798)=190.0;
DEPTH(592,799)=646.0;
DEPTH(593,800)=616.0;
DEPTH(594,800)=616.0;
DEPTH(595,800)=548.0;
DEPTH(596,800)=548.0;
DEPTH(597,800)=464.0;
DEPTH(598,800)=464.0;
DEPTH(599,800)=464.0;
DEPTH(600,800)=413.0;
DEPTH(601,800)=397.0;
DEPTH(602,800)=397.0;
DEPTH(603,800)=349.0;
DEPTH(604,800)=268.0;
DEPTH(605,800)=155.0;
DEPTH(606,800)=115.0;
DEPTH(607,800)=93.0;
DEPTH(586,804)=429.0;
DEPTH(587,804)=429.0;
DEPTH(588,804)=429.0;
DEPTH(589,804)=429.0;
DEPTH(590,804)=429.0;
DEPTH(591,804)=429.0;
DEPTH(592,804)=397.0;
DEPTH(593,804)=397.0;
DEPTH(594,804)=349.0;
DEPTH(595,804)=349.0;
DEPTH(596,804)=349.0;
DEPTH(600,804)=516.0;
DEPTH(601,804)=516.0;
DEPTH(602,804)=562.0;

 % Define variables.

nc{‘depth’}= ncdouble(‘ny’,’nx’);

nc{‘depth’}.units= ‘meters’;
nc{‘depth’}.standard_name=’topographic depth at T-cell centers’;

% Store the data
nc{‘depth’}(:) = DEPTH;
nc = close(nc);