The clavrx_merra2 from joleenf

clavrx_merra2's People

Contributors

clavrx_merra2's Issues

fill value incorrect/sets wrong attribute fill

clavrx_MERRA2/merra2/merra24clavrx_brett.py

Lines 500 to 510 in 86ac237

 elif self.out_name == 'rh': 

 temp_sds = in_sds.variables['T'] # temperature in [K] (Time, Height, Y, X) 

 temp_k = temp_sds[time_ind] 

 temp_k = _trim_toa(temp_k) 

 fill = temp_sds._FillValue 

 temp_k[temp_k == fill] = np.nan 

 fill = in_sds.variables[self.in_name]._FillValue 

 data[data == fill] = np.nan 

 rh = qv_to_rh(data, temp_k) 

 rh[np.isnan(rh)] = fill 

 out_sds.set(_refill(_reshape(rh, self.ndims_out, fill), fill))

In this block, the fill value later written as an attribute to the hdf is set to 1e+15. However, the calculated field (data) is a masked_array with fill_value of 1e+20. This creates a mismatch when the mask array is filled before writing to the hdf.

Also, trim_toa is truncating the levels without checking range.

Recommend Projects

React

A declarative, efficient, and flexible JavaScript library for building user interfaces.

Vue.js

🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

Typescript

TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

TensorFlow

An Open Source Machine Learning Framework for Everyone

Django

The Web framework for perfectionists with deadlines.

Laravel

A PHP framework for web artisans

D3

Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

javascript

JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

web

Some thing interesting about web. New door for the world.

server

A server is a program made to process requests and deliver data to clients.

Machine learning

Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

Visualization

Some thing interesting about visualization, use data art

Game

Some thing interesting about game, make everyone happy.

Recommend Org

Facebook

We are working to build community through open source technology. NB: members must have two-factor auth.

Microsoft

Open source projects and samples from Microsoft.

Google

Google ❤️ Open Source for everyone.

Alibaba

Alibaba Open Source for everyone

D3

Data-Driven Documents codes.

Tencent

China tencent open source team.

	elif self.out_name == 'rh':
	temp_sds = in_sds.variables['T'] # temperature in [K] (Time, Height, Y, X)
	temp_k = temp_sds[time_ind]
	temp_k = _trim_toa(temp_k)
	fill = temp_sds._FillValue
	temp_k[temp_k == fill] = np.nan
	fill = in_sds.variables[self.in_name]._FillValue
	data[data == fill] = np.nan
	rh = qv_to_rh(data, temp_k)
	rh[np.isnan(rh)] = fill
	out_sds.set(_refill(_reshape(rh, self.ndims_out, fill), fill))