File tree Expand file tree Collapse file tree 2 files changed +3
-7
lines changed
Expand file tree Collapse file tree 2 files changed +3
-7
lines changed Original file line number Diff line number Diff line change @@ -8,7 +8,6 @@ In this page we show how you can download data for a given `NASAPrecipitation` d
88using NASAPrecipitation
99using CairoMakie
1010using DelimitedFiles
11- using Logging
1211
1312download("https://raw.githubusercontent.com/natgeo-wong/GeoPlottingData/main/coastline_resl.txt","coast.cst")
1413coast = readdlm("coast.cst",comments=true)
@@ -47,10 +46,8 @@ Let us download the `IMERGMonthly` Dataset for 2020 over the Caribbean (as defin
4746``` @example download
4847npd = IMERGMonthly(start=Date(2020),stop=Date(2020),path=pwd())
4948geo = GeoRegion("AR6_CAR")
50- hroot = "https://gpm1.gesdisc.eosdis.nasa.gov/opendap/AUXILIARY"
51- npdnc = "GPM_IMERG_LandSeaMask.2/GPM_IMERG_LandSeaMask.2.nc4"
52- npdds = NCDataset(joinpath(hroot,npdnc),"r")
53- @info npdds
49+ lsd = getLandSea(npd,geo)
50+ # download(npd,geo)
5451```
5552
5653## Reading data for a downloaded ` NASAPrecipitationDataset `
Original file line number Diff line number Diff line change @@ -87,14 +87,13 @@ function downloadLandSea(
8787 lon,lat = gpmlonlat ()
8888 nlon = length (lon)
8989 nlat = length (lat)
90- var = zeros (Float32,nlat,nlon)
9190 mask = ones (Int16,nlon,nlat)
9291
9392 @info " $(modulelog ()) - Retrieving the original IMERG Land-Sea Mask data from NASA's EOSDIS OPeNDAP servers"
9493 hroot = " https://gpm1.gesdisc.eosdis.nasa.gov/opendap/AUXILIARY"
9594 npdnc = " GPM_IMERG_LandSeaMask.2/GPM_IMERG_LandSeaMask.2.nc4"
9695 npdds = NCDataset (joinpath (hroot,npdnc)," r" )
97- NCDatasets . load! (npdds[" landseamask" ]. var,var, :,:)
96+ var = nomissing (npdds[" landseamask" ][ :,:], NaN )
9897 close (npdds)
9998
10099 for ilon = 1 : nlon, ilat = 1 : nlat
You can’t perform that action at this time.
0 commit comments