[ferret_users] integration from surface pressure level at each grid point

[Pratik Kad <pratikkad17@xxxxxxxxx>]

Dear Ferret users,

I would like to take vertical integration from surface pressure level to top level. I have estimated height with respect to pressure level; I need to take surface pressure where we have a first positive value.

For example, at one gird point (here positive values of "H" is Height above the surface)

 

yes? list/x=83/y=28  PRES,H
Z (hPa): 0 to 1013
LONGITUDE: 83
LATITUDE: 28
 Column  1: PRES is Z[GZ=H]
 Column  2: H is Height (m)
              PRES     H
10     /  1:    10.  29854.
30     /  2:    30.  22620.
50     /  3:    50.  19412.
70     /  4:    70.  17367.
100    /  5:   100.  15271.
150    /  6:   150.  12855.
175    /  7:   175.  11897.
200    /  8:   200.  11043.
225    /  9:   225.  10277.
250    / 10:   250.   9592.
300    / 11:   300.   8351.
350    / 12:   350.   7276.
400    / 13:   400.   6325.
450    / 14:   450.   5468.
500    / 15:   500.   4674.
550    / 16:   550.   3937.
600    / 17:   600.   3261.
650    / 18:   650.   2623.
700    / 19:   700.   2023.
750    / 20:   750.   1456.
775    / 21:   775.   1185.
800    / 22:   800.    922.
825    / 23:   825.    666.
850    / 24:   850.    416.
875    / 25:   875.    171.
900    / 26:   900.    -70.
925    / 27:   925.   -311.
950    / 28:   950.   -551.
975    / 29:   975.   -791.
1000   / 30:  1000.  -1031.

yes? let integrate = any_para[z=875:10@din] 

How can I define the surface pressure level at each grid point? Please suggest. 

Best Regards

Pratik