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