Skip to main content
13 votes
Accepted

Regarding compression, shuffle filter of netcdf4

I made some pictures on my tablet to explain everything a bit better. Sorry that they are not that nice. Shuffle If we do not us shuffle, we store on value (float, double, integer, short integer, ......
daniel.heydebreck's user avatar
12 votes
Accepted

Language to modernize a scientific Fortran codebase

The modern language you are looking for is called modern Fortran 2018. Fortran 2008 and 2018 have everything a numerical computing project would need and so many features that many other languages ...
Scientist's user avatar
  • 236
8 votes

Language to modernize a scientific Fortran codebase

As a former Fortran programmer I did a small online review several weeks ago into current trends for scientific programming. To begin with, despite its age and its sometimes archaic style of ...
Fred's user avatar
  • 24.7k
7 votes

Python - How to write interpolated raster to netcdf file

I changed your model solution a little bit; but it's like Ingvar Lukas wrote in his answer: you redefined xi and yi in the ...
Erik's user avatar
  • 1,999
6 votes

Python - How to write interpolated raster to netcdf file

You have a shape mismatch, as you are overwritingxi and yi using np.meshgrid, then assigning ...
Ingvar Lukas's user avatar
4 votes
Accepted

Digitalizing precipitation data from the image of weather map

You can convert the image to NetCDF data using gdal_translate. The command line looks something like: ...
Gemechu Fanta Garuma's user avatar
3 votes

Efficient way to parse lat-lon grid at certain depth in tripolar model from nc file?

First of all: NEVER, EVER, EVER (I repeat, EVER) loop in Python over an array to subset some values :) xarray was created exactly to avoid having to write something ...
Droid's user avatar
  • 421
3 votes

How to determine the orientation of coastline from NetCDF file

It should be noted that NetCDF data just describes the format that the data is in. One dataset that I know of (IBTrACS) contains data that is not gridded data, and also contains a landmask. But this ...
BarocliniCplusplus's user avatar
3 votes

Huge netcdf files handling

For processing large data it is a good practice to call data into RAM in slices(Either by spliting time axis or spatial domain). In the interest of earth sciences python packages Xarray, iris, netCDF4 ...
Jithu's user avatar
  • 31
2 votes

How to correct PSMSL data for GIA?

The simple answer is that you just subtract the nearest value from the glacial isostatic adjustment (GIA) model. The values provided by Peltier in your netCDF file are in a grid that is 1 degree x 1 ...
arkaia's user avatar
  • 15.5k
2 votes
Accepted

Storing borehole interval data (logs) in netcdf

tl;dr I would suggest considering the borehole measurements as a vertical profile. In terms of the CF Convention it would be a featureType ...
daniel.heydebreck's user avatar
2 votes
Accepted

Temperature at 2 m MERRA2

When I visit that link, the first option at the top ("Data Product") defaults to inst3_3d_asm_Np, which is atmospheric variables in pressure levels. I suspect you ...
Deditos's user avatar
  • 3,950
2 votes
Accepted

WRF/Chem: Editing NetCDF (wrfchemi files) for own emissions

You can try the R package eixport, with wrf_put ...
Sergio's user avatar
  • 412
2 votes
Accepted

How to prepare GrADS's .ctl file to convert grd file into .nc file in batch mode?

After some trial n error I figured it out using the following shell script ...
Amit Ghosh's user avatar
1 vote
Accepted

Correlation calculation of masked array:

I'm going to hazard a guess that you're using the NOAA OISST data, although I don't think that makes a difference for the answer here. Note that your loops, ...
Deditos's user avatar
  • 3,950
1 vote

Language to modernize a scientific Fortran codebase

Right now, I think that Julia makes the cut, because it provides a great tradeoff between "fast scripting" and "high-performance scientific computing&...
Mike76's user avatar
  • 213
1 vote

Huge netcdf files handling

One option is to not load the entire file at a time. You can use ncgeodataset. The routine allows for the extraction of a subset of data without having to load the entire file or even an array into ...
arkaia's user avatar
  • 15.5k
1 vote

Adding time dimension and corresponding new variable to a netcdf file

There are two issues (and one optimization) with/for the code above: Array indexing starts with 0 (not 1) in Python. Therefore, ...

Only top scored, non community-wiki answers of a minimum length are eligible