Dear NCL users,<br>I'm using NCL (version 6.2.0) on a Linux cluster.<br><br>I'm basically trying to reproduce with NCL what is explained here using Ferret (<a href="http://ferret.pmel.noaa.gov/Ferret/faq/depth-axis-to-density-axis" target="_BLANK">http://ferret.pmel.noaa.gov/Ferret/faq/depth-axis-to-density-axis</a>), through the function ZAXREPLACE. My aim is to have potential vorticity in density coordinate, in order to plot PV values on isopycnal surfaces.<br>Up to date, I tried to use int2p, but results are not the same I got with Ferret. I'd like to know if there is another method I could use. I attach my script, in which I compute PV, I change vertical axis, and then I put the outcome in a netcdf file. If you need also input and output files, I can upload them in ftp.<br><br>Thank you very much.<br>Kind regards,<br>Giorgio<br>