Downdraft Mass Flux

I am trying to output the downdraft mass flux at the surface globally. In the tphysbc routine, after the call to the zm_convr subrourinte (Zhang-McFarlane (1996) convection parameterization), I have added the following line:

call outfld('DNMSSFLX',md2(:,pver) ,pcols ,lchnk )


Additionally, in the history.F90 module, I have added the following line to the h_default subroutine:

call add_default ('DNMSSFLX', 1, ' ')

and, in the history.F90..bldfld routine:

call addfld ('DNMSSFLX','kg/m2/s ',1,'A','Downdraft Mass Flux ',phys_decomp)

The output looks strange and I wanted to make sure that my output method is correct.

Thank you.
 
Back
Top