Hi,
it may be simpler than this, without waiting for future versions (do the best you can with what you have, no? ;-) )
let's say you have a pressure varying with the coordinate Z of the system CSYS,21 (whatever this coord sys can be).
If you don't want to cope with SFFUN, simply discretize your P-variation law in, say, 10 P values at 10 Z coords.
Then, add a Commands snippet in the Solution and add these commands:
/PREP7
*DIM,P_TABLE,table,10,,,Z,,,21
P_TABLE(1)=p1,p2,...,p10 !pressure values
P_TABLE(1,0)=z1,z2,...,z10 !coord values
cmsel,s,COMPONENT_OF_PRESSURIZED_SURFS
SF,all,pres,%P_TABLE%
allsel
/SOLU
et voilà! The SFFUN method is, by far, more sophisticated, but also more difficult to catch up...
Hope this can help.
Regards