Changeset 334 in lmdz_wrf
- Timestamp:
- Feb 27, 2015, 6:05:11 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/tools/nc_var_tools.py
r331 r334 10971 10971 maskedvals2 = maskedvals*maskedvals 10972 10972 for iz in range(dimz): 10973 r varvals[it,iz,0] = varvalst[iz,box2,box2]10974 r varvals[it,iz,1] = np.min(varvalst[iz,:,:])10975 r varvals[it,iz,2] = np.max(varvalst[iz,:,:])10976 r varvals[it,iz,3] = np.mean(varvalst[iz,:,:])10977 r varvals[it,iz,4] = maskedvals2[iz,:,:].mean()10978 r varvals[it,iz,5] = np.sqrt(rvarvals[it,iz,4] -\10979 r varvals[it,iz,3]*rvarvals[it,iz,3])10973 rstatvarvals[it,iz,0] = varvalst[iz,box2,box2] 10974 rstatvarvals[it,iz,1] = np.min(varvalst[iz,:,:]) 10975 rstatvarvals[it,iz,2] = np.max(varvalst[iz,:,:]) 10976 rstatvarvals[it,iz,3] = np.mean(varvalst[iz,:,:]) 10977 rstatvarvals[it,iz,4] = maskedvals2[iz,:,:].mean() 10978 rstatvarvals[it,iz,5] = np.sqrt(rstatvarvals[it,iz,4] - \ 10979 rstatvarvals[it,iz,3]*rsattvarvals[it,iz,3]) 10980 10980 10981 10981 else:
Note: See TracChangeset
for help on using the changeset viewer.