Can't output convection source in global energy source array

Dear developers,
I try to use des_usr_var to output convection source in global energy source array.According to the description, Qcv is the heat transfer between the particle and the fluid, but the output value during my calculation is 0.I don’t know why,I use the dmp solver.
solar_03_2022-08-12T045253.996079.zip (297.1 KB)

1660305228089

Adding

WRITE(*,*) "Qcv=", Qcv

to your code shows that Qcv is in fact 0.

! Calculate the heat transfer coefficient.
            CALL CALC_GAMMA_DES(NP, GAMMA)
            if (gamma .ne. 0.0) write(*,*) "********GAMMA", gamma
! Calculate heat transfer for a real particle in the parcel, then scale W times
            if(CGDEM) then
               GAMMAxSA(NP) = GAMMA* l4Pi*DES_CGP_RPR(NP)*DES_CGP_RPR(NP) *DES_CGP_STW(NP)
            else
               GAMMAxSA(NP) = GAMMA* l4Pi*DES_RADIUS(NP)*DES_RADIUS(NP)
            endif

! Calculate the rate of heat transfer to the particle
            Qcv = GAMMAxSA(NP)*(lTg - DES_T_s(NP))

the value of GAMMA cming from CALC_GAMMA_DES is 0, which causes GAMMAxSA and Qcv to both be 0.

So the question is, why does CALC_GAMMA_DES always return 0?

image
image
When I use the udf kg in usr_properties.f,the gamma is zero.Is a bug?