Ocean script mal function

Hi,

I am using ocean scripting for calculating propagation delay and energy per transition of my circuit. There I observed a weird behaviour. When I run the script for one certain value, it generates some value. But when I run a parametric simulation, for the same points the values changed a lot. My script is as follows. The interesting point is, the mismatch is happening only in case of energy calculation not anything else.

simulator( 'spectre )
design( "/homes/schatterjee/subvt_28/SUBVT28CORE_test/simulation/LSLH_LVT/spectre/schematic/netlist/netlist")
resultsDir( "/homes/schatterjee/subvt_28/SUBVT28CORE_test/simulation/LSLH_LVT/spectre/schematic" )
modelFile(
'("/homes/schatterjee/subvt_28/SUBVT28CORE_test/corners.scs" "")
)
definitionFile(
"models.scs"
)
analysis('tran ?stop "4u" )
analysis('dc ?dev "/V1" ?param "dc" ?start "0"
?stop "1" )
desVar( "t" 500n )
desVar( "wp1" 540n )
desVar( "lp1" 600n )
desVar( "wp2" 150n )
desVar( "wp3" 320n )
desVar( "lp3" 60n )
desVar( "wn1" 160n )
desVar( "wn2" 500n )
desVar( "vdd_h" 1 )
desVar( "vdd_l" 300m )
envOption(
'analysisOrder list("tran" "dc")
)
save( 'i "/V0/MINUS" "/V2/MINUS" "/V1/MINUS" )
temp( 27 )
run()
paramAnalysis("t" ?values '(500n 625n 1000n 1250n 2500n 5000n)
; paramAnalysis("vdd_l" ?start 250m ?stop 1000m ?step 50m
; paramAnalysis("wp1" ?start 200n ?stop 800n ?step 6n
; paramAnalysis("lp1" ?start 400n ?stop 800n ?step 4n
; paramAnalysis("wp2" ?start 80n ?stop 480n ?step 4n
; paramAnalysis("wp3" ?start 180n ?stop 580n ?step 4n
; paramAnalysis("lp3" ?start 50n ?stop 65n ?step 1n
; paramAnalysis("wn1" ?start 160n ?stop 660n ?step 5n
; paramAnalysis("ln1" ?start 48n ?stop 60n ?step 1n
; paramAnalysis("wn2" ?start 300n ?stop 800n ?step 5n
; )
; )
; )
; )
; )
; )
; )
; )
; )
)
paramRun()

vddl = ymax(v("in" ?result "tran"))
vddl_50 = vddl/2;
vddh_l = ymax(v("z" ?result "tran"))*.10
vddh_h = ymax(v("z" ?result "tran"))*.90

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;; ;;
;;The following section is for the calculation of delays and time related informations ;;
;; ;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

tpd_r_z = delay(?wf1 v("/in" ?result "tran") ?value1 vddl_50 ?edge1 "rising" ?nth1 1 ?td1 0.0 ?wf2 v("/z" ?result "tran") ?value2 0.5 ?edge2 "rising" ?nth2 1 ?td2r0 nil ?stop nil ?multiple nil)
tpd_f_z = delay(?wf1 v("/in" ?result "tran") ?value1 vddl_50 ?edge1 "falling" ?nth1 1 ?td1 0.0 ?wf2 v("/z" ?result "tran") ?value2 0.5 ?edge2 "falling" ?nth2 1 ?td2r0 nil ?stop nil ?multiple nil)

tpd = ((tpd_r_z + tpd_f_z)*1e9) / 2; The average is taken of the rise and fall propagation delay

plot(tpd ?expr '("tpd"))

;; Energy per transition calculation

tr_l = cross(v("/z" ?result "tran") vddh_l 1 "rising" nil nil )
tr_h = cross(v("/z" ?result "tran") vddh_h 1 "rising" nil nil )

enr_trans_r= vddl * integ((IT("/V2/MINUS") + IT("/V0/MINUS")), value(tr_l, vddl), value(tr_h, vddl)) + integ(IT("/V1/MINUS"), value(tr_l, vddl), value(tr_h, vddl))

tf_l = cross(v("/z" ?result "tran") vddh_l 1 "falling" nil nil )
tf_h = cross(v("/z" ?result "tran") vddh_h 1 "falling" nil nil )

enr_trans_f= vddl * integ((IT("/V2/MINUS") + IT("/V0/MINUS")), value(tf_h, vddl), value(tf_l, vddl)) + integ(IT("/V1/MINUS"), value(tf_h, vddl), value(tf_l, vddl))

enr_trans = 1e15 * (enr_trans_r + enr_trans_f)

plot(enr_trans ?expr '("enr_trans"))

Any help would be appreciated.

Thanks

Saikat