I`dont know if the following challenge is possible. I have temperature distributed on 33 levitus-levels and a data set with the values of the inferior limits of isopycnal layers, 22 at most (see below)
currently SET data sets:
1> /dat1/jean/BAS/s.nc
name title I J K L
... ...
S 1:369 1:401 1:22 1:72
2> ./tmp33.nc (default)
name title I J K L
TMP 1:369 1:401 1:33 1:72
For example
yes? list/i=200/j=200/l=1 tmp
VARIABLE : TMP
FILENAME : tmp33.nc
SUBSET : 33 points (DEPTH (m))
LONGITUDE: 43.2W(-43.2)
LATITUDE : 26S
T : 1
43.19W
200
0 / 1: 27.69
10 / 2: 27.77
20 / 3: 26.62
30 / 4: 25.81
50 / 5: 25.30
75 / 6: 24.38
100 / 7: 22.86
125 / 8: 21.91
150 / 9: 21.09
200 / 10: 19.36
250 / 11: 16.83
300 / 12: 15.24
400 / 13: 12.78
500 / 14: 10.12
600 / 15: 7.74
700 / 16: 5.53
800 / 17: 4.25
900 / 18: 3.48
1000 / 19: 3.11
1100 / 20: 2.94
1200 / 21: 2.81
1300 / 22: 2.73
1400 / 23: 2.66
1500 / 24: 2.62
1750 / 25: 2.61
2000 / 26: 2.65
2500 / 27: ....
3000 / 28: ....
3500 / 29: ....
4000 / 30: ....
4500 / 31: ....
5000 / 32: ....
5500 / 33: ....
and...
yes? list/i=200/j=200/l=1 s[d=1]
VARIABLE : S (m)
FILENAME : s.nc
FILEPATH : /dat1/jean/BAS/
SUBSET : 22 points (Z)
LONGITUDE: 43.2W(-43.2)
LATITUDE : 26S
T : 1
43.19W
200
0 / 1: 3.0
1 / 2: 6.4
2 / 3: 10.2
3 / 4: 14.4
4 / 5: 19.2
5 / 6: 24.7
6 / 7: 46.4
7 / 8: 92.1
8 / 9: 190.5
9 / 10: 254.3
10 / 11: 374.9
11 / 12: 457.9
12 / 13: 555.9
13 / 14: 665.3
14 / 15: 827.6
15 / 16: 1007.4
16 / 17: 1215.6
17 / 18: 1517.7
18 / 19: 2283.9
19 / 20: 2283.9
20 / 21: ....
21 / 22: ....
Now, I`d like produce a new Z axis for the temperature with the values of S among the 33 levitus-levels . If it is possible I will use @FLN:n—linear interpolation filler to fill holes .
Best regards
Jean