[Thread Prev][Thread Next][Index]

[ferret_users] zaxreplace - convert z-axis of density to the z-axis of `m`



Dear Ferretrs,
G'day

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])

sh gr heatsig
GRID GBL1
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
TIME               TIME             180 i   16-JAN-2004 05:14   16-DEC-2018 10:03   full

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 ?

Eagerly waiting for the reply.

Cheers, Saurabh

--


REGARDS

Saurabh Rathore
Research Scholar (PhD.)
Centre For Oceans, Rivers, Atmosphere & Land Science Technology
Indian Institute Of Technology, Kharagpur
contact :- 91- 8345984434

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement