You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: src/get.jl
+12Lines changed: 12 additions & 0 deletions
Original file line number
Diff line number
Diff line change
@@ -47,10 +47,12 @@ function getLandSea(
47
47
if!isfile(lsmfnc)
48
48
49
49
@info "$(modulelog()) - The ETOPO $(uppercase(type)) Relief dataset for the \"$(geo.ID)\" GeoRegion is not available, extracting from Global ETOPO $(uppercase(type)) Relief dataset ..."
50
+
flush(stderr)
50
51
51
52
area = (geo.E-geo.W)/360* (geo.N-geo.S)/180
52
53
if area >0.1* (resolution/60)^2
53
54
@info "$(modulelog()) - The total specified area covers more than $(10* (resolution/60)^2)% of global points, set `downloadglb = true` ..."
@info "$(modulelog()) - The Global ETOPO $(uppercase(type)) Relief dataset has already been downloaded, now using it for extraction ..."
64
+
flush(stderr)
62
65
else
63
66
@info "$(modulelog()) - Downloading the Global ETOPO $(uppercase(type)) Relief dataset ..."
67
+
flush(stderr)
64
68
setup(type,etd.path,resolution)
65
69
end
66
70
@@ -75,6 +79,7 @@ function getLandSea(
75
79
nlat = length(ggrd.ilat)
76
80
77
81
@info "$(modulelog()) - Extracting regional ETOPO $(uppercase(type)) Relief data for the \"$(geo.ID)\" GeoRegion from the Global ETOPO Relief dataset ..."
82
+
flush(stderr)
78
83
79
84
roro = extract(goro,ggrd)
80
85
rlsm = deepcopy(roro)
@@ -83,12 +88,14 @@ function getLandSea(
83
88
rlsm[isnan.(roro)] .=NaN
84
89
85
90
@info "$(modulelog()) - Saving the regional ETOPO Relief data for \"$(geo.ID)\" GeoRegion ..."
@info "$(modulelog()) - Extracting regional ETOPO $(uppercase(type)) Land-Sea mask for the \"GLB\" GeoRegion from the Global ETOPO Land-Sea mask dataset ..."
0 commit comments