Cloud map plot and tz-plot produce much too low values

[GoogleCodeExporter]

The cloud map plot and the final plot produced by tzWrf produce much too low 
values for the sum of cloud water etc in the column. The reason for this is 
that the code fails to multiply the values by the layer thickness.

In mapWRF.mapCloud the following code needs to be added:
    zfull = (nc.variables['PH'][time,:,:,:] + nc.variables['PHB'][time,:,:,:])/g
    dz = np.zeros((Nz,Ny,Nx))
    for k in arange(Nz):
        dz[k,:,:] = zfull[k+1,:,:]-zfull[k,:,:]

and then use dz when summing the cloud water in the column:
    qcloud = 1000.0*np.sum(dz*rho*nc.variables['QCLOUD'][time,:,:,:],axis=0) 

Similarly for qrain.

Original issue reported on code.google.com by gwaa...@gmail.com on 14 Feb 2013 at 12:44