ncview can not show density output

Hi @jvegas,

I have been able to compute the density diagnostic without issues, nevertheless I would like to ask if it was possible to add a mask to the final output.

I added a mask with cdo or python, but it would be nice to have the final diagnostic already masked.

You can see an example of original output here:

/esarchive/exp/ecearth/a2er_sixth/diags/OMIP/EC-Earth-Consortium/EC-Earth3/omip2/r1i1p1f1/Omon/sigma0/gn/v20200301/

you can check that ncview gives this error:

Note: no Ncview app-defaults file found, using internal defaults
calculating min and maxes for sigma0...
*** genlevs: step =-1, failed
*** genlevs: Input data: mindat=-inf maxdat=1029.308105 nlevels=9
*** trialstep=inf, mant=inf, exp=-2147483648.000000

after masking it, there is no error.

Here you can find a masked file I created as test:

/esarchive/scratch/Earth/vsicardi/OMIP/a2er_sixth/test_mask.nc

thank you