NADCON5-ng  0.0.2
NADCON5 Next Generation Documentation
plotcoast.f
Go to the documentation of this file.
1 c> \ingroup core
2 c> \if MANPAGE
3 c> \page plotcoast
4 c> \endif
5 c>
6 c> Subroutine to write GMT-based commands to create a shoreline
7 c> Write GMT-based commands to create a shoreline
8 c> based on region.
9 c>
10 c> Use GMT-default coastline for:
11 c> - `conus`
12 c> - `alaska`
13 c>
14 c> Use Dru's custome coastline for:
15 c> - hawaii
16 c> - prvi
17 c> - as
18 c> - guamcnmi
19 c> - stlawrence
20 c> - stmatthew
21 c> - stgeorge
22 c> - stpaul
23 c>
24 c> \param[in] region The Region to create coastline for
25 c> \param[in] ifnum the file descriptor of the output file to write `GMT` commands to
26 c>
27 c>
28 c> ## Changelog
29 c>
30 c> ### 2016 01 07:
31 c> Forced the Alaska region to plot the
32 c> islands of St. George, St. Matthew and St. Paul
33 c> (St. Lawrence is already plotted), as well
34 c> as 35 missing Aleutian Islands
35 c>
36 c> ### 2015 09 23:
37 c> Added four new regions:
38 c> - St. Lawrence Island, Alaska
39 c> - St. Matthew Island, Alaska
40 c> - St. George Island, Alaska
41 c> - St. Paul Island, Alaska
42 c>
43  subroutine plotcoast(region,ifnum)
44 
45 c - 2016 01 07
46 c Forced the Alaska region to plot the
47 c islands of St. George, St. Matthew and St. Paul
48 c (St. Lawrence is already plotted), as well
49 c as 35 missing Aleutian Islands
50 
51 c - 2015 09 23: Added four new regions:
52 c St. Lawrence Island, Alaska
53 c St. Matthew Island, Alaska
54 c St. George Island, Alaska
55 c St. Paul Island, Alaska
56 
57  character*10 region
58 
59 c - Write GMT-based commands to create a shoreline
60 c - based on region.
61 c - Use GMT-default coastline for:
62 c conus
63 c alaska
64 c - Use Dru's custome coastline for:
65 c hawaii
66 c prvi
67 c as
68 c guamcnmi
69 c stlawrence
70 c stmatthew
71 c stgeorge
72 c stpaul
73 
74 
75  if (trim(region).eq.'conus')then
76  write(ifnum,100)
77  elseif(trim(region).eq.'alaska')then
78  write(ifnum,109)
79  elseif(trim(region).eq.'hawaii')then
80  write(ifnum,101)
81  elseif(trim(region).eq.'prvi')then
82  write(ifnum,102)
83  elseif(trim(region).eq.'guamcnmi')then
84  write(ifnum,103)
85  elseif(trim(region).eq.'as')then
86  write(ifnum,104)
87  elseif(trim(region).eq.'stlawrence')then
88  write(ifnum,105)
89  elseif(trim(region).eq.'stmatthew')then
90  write(ifnum,106)
91  elseif(trim(region).eq.'stgeorge')then
92  write(ifnum,107)
93  elseif(trim(region).eq.'stpaul')then
94  write(ifnum,108)
95 
96  else
97  write(6,200)region
98  stop 10001
99  endif
100  return
101 
102  200 format('FATAL in plotcoast.f: Unknown region: ',a10)
103 
104 c - CONUS
105  100 format(
106  *'pscoast -Df -R -JM -W0.25p -N1 -N2 -A1200 -O >> plot.ps')
107 
108 c - HAWAII
109  101 format(
110  *'psxy Boundaries/HI_Hawaii.gmt -R -JM -B -O -K',
111  *' >> plot.ps',/,
112  *'psxy Boundaries/HI_Kahoolawe.gmt -R -JM -B -O -K ',
113  *' >> plot.ps',/,
114  *'psxy Boundaries/HI_Kauai.gmt -R -JM -B -O -K ',
115  *' >> plot.ps',/,
116  *'psxy Boundaries/HI_Lanai.gmt -R -JM -B -O -K ',
117  *' >> plot.ps',/,
118  *'psxy Boundaries/HI_Maui.gmt -R -JM -B -O -K ',
119  *' >> plot.ps',/,
120  *'psxy Boundaries/HI_Molokai.gmt -R -JM -B -O -K ',
121  *' >> plot.ps',/,
122  *'psxy Boundaries/HI_Nihau.gmt -R -JM -B -O -K ',
123  *' >> plot.ps',/,
124  *'psxy Boundaries/HI_Oahu.gmt -R -JM -B -O ',
125  *' >> plot.ps')
126 
127 c - PRVI
128  102 format(
129  *'psxy Boundaries/VQ_StCroix.gmt -R -JM -B -O -K',
130  *' >> plot.ps',/,
131  *'psxy Boundaries/VQ_StJohn.gmt -R -JM -B -O -K ',
132  *' >> plot.ps',/,
133  *'psxy Boundaries/VQ_StThomas.gmt -R -JM -B -O -K ',
134  *' >> plot.ps',/,
135  *'psxy Boundaries/PR_Culebra.gmt -R -JM -B -O -K ',
136  *' >> plot.ps',/,
137  *'psxy Boundaries/PR_Desecheo.gmt -R -JM -B -O -K ',
138  *' >> plot.ps',/,
139  *'psxy Boundaries/PR_Mona.gmt -R -JM -B -O -K ',
140  *' >> plot.ps',/,
141  *'psxy Boundaries/PR_Viequez.gmt -R -JM -B -O -K ',
142  *' >> plot.ps',/,
143  *'psxy Boundaries/PR_PuertoRico.gmt -R -JM -B -O ',
144  *' >> plot.ps')
145 
146 c - GUAM/CNMI
147  103 format(
148  *'psxy Boundaries/GQ_Guam.gmt -R -JM -B -O -K',
149  *' >> plot.ps',/,
150  *'psxy Boundaries/CQ_Rota.gmt -R -JM -B -O -K ',
151  *' >> plot.ps',/,
152  *'psxy Boundaries/CQ_Saipan.gmt -R -JM -B -O -K ',
153  *' >> plot.ps',/,
154  *'psxy Boundaries/CQ_Tinian.gmt -R -JM -B -O -K ',
155  *' >> plot.ps',/,
156  *'psxy Boundaries/CQ_Anatahan.gmt -R -JM -B -O -K ',
157  *' >> plot.ps',/,
158  *'psxy Boundaries/CQ_Sarigan.gmt -R -JM -B -O -K ',
159  *' >> plot.ps',/,
160  *'psxy Boundaries/CQ_Guguan.gmt -R -JM -B -O -K ',
161  *' >> plot.ps',/,
162  *'psxy Boundaries/CQ_Alamagan.gmt -R -JM -B -O -K ',
163  *' >> plot.ps',/,
164  *'psxy Boundaries/CQ_Pagan.gmt -R -JM -B -O -K ',
165  *' >> plot.ps',/,
166  *'psxy Boundaries/CQ_Agrihan.gmt -R -JM -B -O -K ',
167  *' >> plot.ps',/,
168  *'psxy Boundaries/CQ_Asuncion.gmt -R -JM -B -O -K ',
169  *' >> plot.ps',/,
170  *'psxy Boundaries/CQ_MaugWest.gmt -R -JM -B -O -K ',
171  *' >> plot.ps',/,
172  *'psxy Boundaries/CQ_MaugEast.gmt -R -JM -B -O -K ',
173  *' >> plot.ps',/,
174  *'psxy Boundaries/CQ_MaugNorth.gmt -R -JM -B -O -K ',
175  *' >> plot.ps',/,
176  *'psxy Boundaries/CQ_Pajaros.gmt -R -JM -B -O',
177  *' >> plot.ps')
178 
179 c - AS
180  104 format(
181  *'psxy Boundaries/AS_Aunuu.gmt -R -JM -B -O -K',
182  *' >> plot.ps',/,
183  *'psxy Boundaries/AS_Nuutele.gmt -R -JM -B -O -K ',
184  *' >> plot.ps',/,
185  *'psxy Boundaries/AS_Ofu.gmt -R -JM -B -O -K ',
186  *' >> plot.ps',/,
187  *'psxy Boundaries/AS_Olosega.gmt -R -JM -B -O -K ',
188  *' >> plot.ps',/,
189  *'psxy Boundaries/AS_Rose.gmt -R -JM -B -O -K ',
190  *' >> plot.ps',/,
191  *'psxy Boundaries/AS_Swains.gmt -R -JM -B -O -K ',
192  *' >> plot.ps',/,
193  *'psxy Boundaries/AS_Tau.gmt -R -JM -B -O -K ',
194  *' >> plot.ps',/,
195  *'psxy Boundaries/AS_Tutila.gmt -R -JM -B -O ',
196  *' >> plot.ps')
197 
198 c - St. Lawrence Island, Alaska
199  105 format(
200  *'psxy Boundaries/AK_StLawrence.gmt -R -JM -B -O -K',
201  *' >> plot.ps',/,
202  *'psxy Boundaries/AK_NorthPunuk.gmt -R -JM -B -O -K ',
203  *' >> plot.ps',/,
204  *'psxy Boundaries/AK_MiddlePunuk.gmt -R -JM -B -O -K ',
205  *' >> plot.ps',/,
206  *'psxy Boundaries/AK_SouthPunuk.gmt -R -JM -B -O ',
207  *' >> plot.ps')
208 
209 c - St. Matthew Island, Alaska
210  106 format(
211  *'psxy Boundaries/AK_StMatthew.gmt -R -JM -B -O -K',
212  *' >> plot.ps',/,
213  *'psxy Boundaries/AK_Hall.gmt -R -JM -B -O -K ',
214  *' >> plot.ps',/,
215  *'psxy Boundaries/AK_Pinnacle.gmt -R -JM -B -O ',
216  *' >> plot.ps')
217 
218 c - St. George Island, Alaska
219  107 format(
220  *'psxy Boundaries/AK_StGeorge.gmt -R -JM -B -O',
221  *' >> plot.ps')
222 
223 c - St. Paul Island, Alaska
224  108 format(
225  *'psxy Boundaries/AK_StPaul.gmt -R -JM -B -O -K',
226  *' >> plot.ps',/,
227  *'psxy Boundaries/AK_Otter.gmt -R -JM -B -O -K ',
228  *' >> plot.ps',/,
229  *'psxy Boundaries/AK_Walrus.gmt -R -JM -B -O ',
230  *' >> plot.ps')
231 
232 c - ALASKA
233  109 format(
234  *'pscoast -Df -R -JM -W0.25p -N1 -N2 -A1200 -O -K ',
235  *' >> plot.ps',/,
236  *'psxy Boundaries/AK_NorthPunuk.gmt -R -JM -B -O -K ',
237  *' >> plot.ps',/,
238  *'psxy Boundaries/AK_MiddlePunuk.gmt -R -JM -B -O -K ',
239  *' >> plot.ps',/,
240  *'psxy Boundaries/AK_SouthPunuk.gmt -R -JM -B -O -K ',
241  *' >> plot.ps',/,
242  *'psxy Boundaries/AK_StMatthew.gmt -R -JM -B -O -K',
243  *' >> plot.ps',/,
244  *'psxy Boundaries/AK_Hall.gmt -R -JM -B -O -K ',
245  *' >> plot.ps',/,
246  *'psxy Boundaries/AK_Pinnacle.gmt -R -JM -B -O -K ',
247  *' >> plot.ps',/,
248  *'psxy Boundaries/AK_StGeorge.gmt -R -JM -B -O -K',
249  *' >> plot.ps',/,
250  *'psxy Boundaries/AK_StPaul.gmt -R -JM -B -O -K',
251  *' >> plot.ps',/,
252  *'psxy Boundaries/AK_Otter.gmt -R -JM -B -O -K ',
253  *' >> plot.ps',/,
254  *'psxy Boundaries/AK_Walrus.gmt -R -JM -B -O -K ',
255  *' >> plot.ps',/,
256  *'psxy Boundaries/AK_Adak.gmt -R -JM -B -O -K ',
257  *' >> plot.ps',/,
258  *'psxy Boundaries/AK_Agattu.gmt -R -JM -B -O -K ',
259  *' >> plot.ps',/,
260  *'psxy Boundaries/AK_Amchitka.gmt -R -JM -B -O -K ',
261  *' >> plot.ps',/,
262  *'psxy Boundaries/AK_Amlia.gmt -R -JM -B -O -K ',
263  *' >> plot.ps',/,
264  *'psxy Boundaries/AK_Amukta.gmt -R -JM -B -O -K ',
265  *' >> plot.ps',/,
266  *'psxy Boundaries/AK_Andreanof.gmt -R -JM -B -O -K ',
267  *' >> plot.ps',/,
268  *'psxy Boundaries/AK_Atka.gmt -R -JM -B -O -K ',
269  *' >> plot.ps',/,
270  *'psxy Boundaries/AK_Attu.gmt -R -JM -B -O -K ',
271  *' >> plot.ps',/,
272  *'psxy Boundaries/AK_Buldir.gmt -R -JM -B -O -K ',
273  *' >> plot.ps',/,
274  *'psxy Boundaries/AK_Carlisle.gmt -R -JM -B -O -K ',
275  *' >> plot.ps',/,
276  *'psxy Boundaries/AK_Chagulak.gmt -R -JM -B -O -K ',
277  *' >> plot.ps',/,
278  *'psxy Boundaries/AK_Chugul.gmt -R -JM -B -O -K ',
279  *' >> plot.ps',/,
280  *'psxy Boundaries/AK_Fenimore.gmt -R -JM -B -O -K ',
281  *' >> plot.ps',/,
282  *'psxy Boundaries/AK_FourMountains.gmt -R -JM -B -O -K ',
283  *' >> plot.ps',/,
284  *'psxy Boundaries/AK_GreatSitkin.gmt -R -JM -B -O -K ',
285  *' >> plot.ps',/,
286  *'psxy Boundaries/AK_Herbert.gmt -R -JM -B -O -K ',
287  *' >> plot.ps',/,
288  *'psxy Boundaries/AK_Igitkin.gmt -R -JM -B -O -K ',
289  *' >> plot.ps',/,
290  *'psxy Boundaries/AK_Ikiginak.gmt -R -JM -B -O -K ',
291  *' >> plot.ps',/,
292  *'psxy Boundaries/AK_Kagalaska.gmt -R -JM -B -O -K ',
293  *' >> plot.ps',/,
294  *'psxy Boundaries/AK_Kagami.gmt -R -JM -B -O -K ',
295  *' >> plot.ps',/,
296  *'psxy Boundaries/AK_Kanaga.gmt -R -JM -B -O -K ',
297  *' >> plot.ps',/,
298  *'psxy Boundaries/AK_Kasatochi.gmt -R -JM -B -O -K ',
299  *' >> plot.ps',/,
300  *'psxy Boundaries/AK_Kiska.gmt -R -JM -B -O -K ',
301  *' >> plot.ps',/,
302  *'psxy Boundaries/AK_Koniuji.gmt -R -JM -B -O -K ',
303  *' >> plot.ps',/,
304  *'psxy Boundaries/AK_Oglodak.gmt -R -JM -B -O -K ',
305  *' >> plot.ps',/,
306  *'psxy Boundaries/AK_Seguam.gmt -R -JM -B -O -K ',
307  *' >> plot.ps',/,
308  *'psxy Boundaries/AK_Semichi.gmt -R -JM -B -O -K ',
309  *' >> plot.ps',/,
310  *'psxy Boundaries/AK_Semisopochnoi.gmt -R -JM -B -O -K ',
311  *' >> plot.ps',/,
312  *'psxy Boundaries/AK_Shemya.gmt -R -JM -B -O -K ',
313  *' >> plot.ps',/,
314  *'psxy Boundaries/AK_Tagalak.gmt -R -JM -B -O -K ',
315  *' >> plot.ps',/,
316  *'psxy Boundaries/AK_Tanaga.gmt -R -JM -B -O -K ',
317  *' >> plot.ps',/,
318  *'psxy Boundaries/AK_Ulak.gmt -R -JM -B -O -K ',
319  *' >> plot.ps',/,
320  *'psxy Boundaries/AK_Uliaga.gmt -R -JM -B -O -K ',
321  *' >> plot.ps',/,
322  *'psxy Boundaries/AK_Umak.gmt -R -JM -B -O -K ',
323  *' >> plot.ps',/,
324  *'psxy Boundaries/AK_Yunaska.gmt -R -JM -B -O ',
325  *' >> plot.ps')
326 
327  end
subroutine plotcoast(region, ifnum)
Subroutine to write GMT-based commands to create a shoreline Write GMT-based commands to create a sho...
Definition: plotcoast.f:44