0

i am working on netcdf file and interpolating data for geopotential height 850hpa , 700hpa, 500hpa and 3oohpa

    height =getvar(ncfile, 'z', timeidx=i)
    pressure=getvar(ncfile,'pressure',timeidx=i)
    temp=interplevel(height,pressure,850)

while inspecting data,height and pressure are completely numeric value, and while interpolating data for 500 and 300 , it is working well

what might be the reason behind this and solution for this

Improve this question
3
  • The 850hPa level corresponds to about 1,450masl: your NaNs are probably low mountains. 500hPa is about 5,500masl, above most terrain worldwide Commented Sep 11, 2024 at 7:30
  • yes, actually this region is himalyan region where value is NaN, i got the reason, but what should i do ,leave it blank, this will create a white region at nan value , which i think is not good appraoch Commented Sep 11, 2024 at 8:33
  • It all depends on your intended use. If there is no 850hPa atmospheric level then leaving it blank seems like a sensible approach Commented Sep 11, 2024 at 9:46
Related questions
1003
Peak detection in a 2D array
2
Interpolating 3D WRF Data against a 2D array in python
4
Merging multiple Netcdf files according to the vertical level
Load 2 more related questions Show fewer related questions

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.