Changeset 334 in lmdz_wrf


Ignore:
Timestamp:
Feb 27, 2015, 6:05:11 PM (10 years ago)
Author:
lfita
Message:

Fixing problem on the statisitics variable on radius close to borders in 'compute_tevolboxtraj'

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/tools/nc_var_tools.py

    r331 r334  
    1097110971                    maskedvals2 = maskedvals*maskedvals
    1097210972                    for iz in range(dimz):
    10973                         rvarvals[it,iz,0] = varvalst[iz,box2,box2]
    10974                         rvarvals[it,iz,1] = np.min(varvalst[iz,:,:])
    10975                         rvarvals[it,iz,2] = np.max(varvalst[iz,:,:])
    10976                         rvarvals[it,iz,3] = np.mean(varvalst[iz,:,:])
    10977                         rvarvals[it,iz,4] = maskedvals2[iz,:,:].mean()
    10978                         rvarvals[it,iz,5] = np.sqrt(rvarvals[it,iz,4] -              \
    10979                           rvarvals[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])
    1098010980
    1098110981                else:
Note: See TracChangeset for help on using the changeset viewer.