I am facing a problem about the conversion of z-axis from density to z-axis of meters depth. For this purpose I am using ZAXREPLACE function. I am going to explain what I did and where I am facing the problem. I replace the `m` depth axis of my data by using zaxreplace function to the depth axis of density.
sh gr heat[d=1]
GRID GTS1
name axis # pts start end subset
LONGITUDE LONGITUDE 360mr 20.5E 19.5E(379.5) full
LATITUDE LATITUDE 145 r 64.5S 79.5N full
DEPTH_AXIS DEPTH (m) 58 i- 2.4811 1950.9 full
TIME TIME 180 i 16-JAN-2004 05:14 16-DEC-2018 10:03 full
sh gr rho[d=2]
GRID GFP1
name axis # pts start end subset
LON LONGITUDE 360 r 20.5E 19.5E(379.5) full
LAT LATITUDE 145 r 64.5S 79.5N full
LEV DEPTH (dbar) 58 i- 2.4811 1950.9 full
TIME TIME 180 i 16-JAN-2004 05:14 16-DEC-2018 10:03 full
define axis/z=`RHO[x=@min,y=@min,l=@min,z=@min,d=2]`:`RHO[x=@max,y=@max,l=@max,z=@max,d=2]`:.1 axden
let heatsig=zaxreplace(HEAT[d=1],NEUTRALRHO[d=2,gx=heat[d=1],gy=heat[d=1],gt=heat[d=1]],z[gz=axden])
then I computed the time slope of heatsig which is shown like this.
sh gr slope
GRID GRC1
name axis # pts start end subset
LONGITUDE LONGITUDE 360mr 20.5E 19.5E(379.5) full
LATITUDE LATITUDE 145 r 64.5S 79.5N full
AXDEN Z 118 r 16.81 28.51 full
normal T
now I want to map it back to the `m` axis of depth which is define by `
ohc.nc` so how I should do it ? or Do I map it first to `m` axis and after that I compute the slope ?