Hi all,
I was trying to calculate the transport using v[x=@din,z=@din], and i
compared these values with v*xbox*zbox followed by x=@sum,z=@sum]. The
results I am getting is not at all comparable. I understood from the
forum that depth axis with irregular dz can cause wrong result. Then
which way I should follow?. Here I attaching the details. I would be
grateful for any suggestions
use Test.nc
let trans_all= vclim*xbox*zbox
let tran =trans_all[y=6,x=82:98@sum,z=0:300@sum]
let transint= vclim[y=6,x=82:98@din,z=0:300@din]
list tran,transint
LIST ZBOXLO[GZ=vclim],ZBOXHI[GZ=vclim],Z[GZ=vclim],ZBOX[GZ=vclim]
Cached data cleared from memory
DATA SET: ./Test.nc
T: -0.5 to 11.5
LONGITUDE: 82E to 98E
LATITUDE: 5.8N
DEPTH (m): 0 to 300
Column 1: TRAN is TRANS_ALL[Y=6,X=82:98@SUM,Z=0:300@SUM]
Column 2: TRANSINT is VCLIM[Y=6,X=82:98@DIN,Z=0:300@DIN]
TRAN TRANSINT
0 / 1: -48.22 -5516444.
1 / 2: -36.44 -4617991.
2 / 3: -1.73 -1081710.
3 / 4: 36.21 3283094.
4 / 5: 37.66 4109916.
5 / 6: -15.87 -1500664.
6 / 7: -76.87 -8493372.
7 / 8: -79.02 -9094106.
8 / 9: -64.63 -7541917.
9 / 10: -17.76 -2386129.
10 / 11: 12.34 1152820.
11 / 12: -6.74 -742074.
DATA SET: ./Test.nc
DEPTH (m): -0.02 to 5500
Column 1: ZBOXLO is ZBOXLO (axis DEPTH)
Column 2: ZBOXHI is ZBOXHI (axis DEPTH)
Column 3: Z is Z (axis DEPTH)
Column 4: ZBOX is ZBOX (axis DEPTH)
ZBOXLO ZBOXHI Z ZBOX
5 / 1: -0. 10. 5. 10.1
15.1 / 2: 10. 20. 15. 10.1
25.3 / 3: 20. 31. 25. 10.3
35.8 / 4: 31. 41. 36. 10.7
46.6 / 5: 41. 52. 47. 11.1
58 / 6: 52. 64. 58. 11.7
70 / 7: 64. 76. 70. 12.5
82.9 / 8: 76. 90. 83. 13.5
96.9 / 9: 90. 105. 97. 14.7
112.3 / 10: 105. 121. 112. 16.3
129.5 / 11: 121. 139. 129. 18.3
149 / 12: 139. 160. 149. 21.0
171.4 / 13: 160. 185. 171. 24.4
197.8 / 14: 185. 214. 198. 29.0
229.5 / 15: 214. 249. 229. 35.3
268.5 / 16: 249. 293. 268. 44.1
317.6 / 17: 293. 350. 318. 56.5
381.4 / 18: 350. 424. 381. 74.1
465.9 / 19: 424. 523. 466. 99.0
579.3 / 20: 523. 654. 579. 131.7
729.4 / 21: 654. 824. 729. 169.5
918.4 / 22: 824. 1029. 918. 204.9
1139.2 / 23: 1029. 1259. 1139. 230.1
1378.6 / 24: 1259. 1502. 1379. 243.3
1625.7 / 25: 1502. 1750. 1626. 248.3
1875.1 / 26: 1750. 2000. 1875. 249.7
2125 / 27: 2000. 2250. 2125. 249.9
2375 / 28: 2250. 2500. 2375. 250.0
2625 / 29: 2500. 2750. 2625. 250.0
2875 / 30: 2750. 3000. 2875. 250.0
3125 / 31: 3000. 3250. 3125. 250.0
3375 / 32: 3250. 3500. 3375. 250.0
3625 / 33: 3500. 3750. 3625. 250.0
3875 / 34: 3750. 4000. 3875. 250.0
4125 / 35: 4000. 4250. 4125. 250.0
4375 / 36: 4250. 4500. 4375. 250.0
4625 / 37: 4500. 4750. 4625. 250.0
4875 / 38: 4750. 5000. 4875. 250.0
5125 / 39: 5000. 5250. 5125. 250.0
5375 / 40: 5250. 5500. 5375. 250.0