Changeset 751 in lmdz_wrf


Ignore:
Timestamp:
May 4, 2016, 11:46:34 AM (9 years ago)
Author:
lfita
Message:

Addig 'gen' on `SpatialWeightedMean?'

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/tools/nc_var_tools.py

    r750 r751  
    1194711947    idim = 0
    1194811948    for dim in list(ovar.dimensions):
    11949         if not searchInlist(seldims,dim):
     11949        if not gen.searchInlist(seldims,dim):
    1195011950            newshape.append(ovar.shape[idim])
    1195111951            noloopdims.append(dim)
     
    1196711967    idim = 0
    1196811968    for dim in ovar.dimensions:
    11969         if searchInlist(dloop,dim):
     11969        if gen.searchInlist(dloop,dim):
    1197011970            slicevals.append(dloopindex[idim])
    1197111971        else:
     
    1226312263              + varname2 + ' as space weights'
    1226412264            if len(loopshape) == 1:
    12265                 newvals = np.ones((loopshape[0]), dtype=np.float)*fillValueF
     12265                newvals = np.ones((loopshape[0]), dtype=np.float)*gen.fillValueF
    1226612266                for id1 in range(loopshape[0]):
    1226712267                    slicevalues = SliceVar(iovar,dimsloop,[id1])
     
    1227212272            elif len(loopshape) == 2:
    1227312273                newvals = np.ones((loopshape[0],loopshape[1]),dtype=np.float)*       \
    12274                   fillValueF
     12274                  gen.fillValueF
    1227512275                for id1 in range(loopshape[0]):
    1227612276                    for id2 in range(loopshape[1]):
     
    1228212282            elif len(loopshape) == 3:
    1228312283                newvals = np.ones((loopshape[0],loopshape[1],loopshape[2]),          \
    12284                   dtype=np.float)*fillValueF
     12284                  dtype=np.float)*gen.fillValueF
    1228512285                for id1 in range(loopshape[0]):
    1228612286                    for id2 in range(loopshape[1]):
     
    1231412314        ivarwgtv = np.abs(np.cos(latv*np.pi/180.))
    1231512315        TOTsumwgt = np.sum(ivarwgtv)
     12316
     12317        print 'Lluis len(looopshape):', len(loopshape)
    1231612318        if len(loopshape) == 1:
    1231712319            newvals = np.ones((loopshape[0]), dtype=np.float)*fillValueF
     
    1236612368##
    1236712369    newvar = onewnc.createVariable(varn + 'spaceweightsum', 'f4', tuple(dimsloop),   \
    12368       fill_value=fillValueF)
     12370      fill_value=gen.fillValueF)
    1236912371    basicvardef(newvar, varn + 'spaceweightsum', 'space summed ' + varn + ' ' +      \
    1237012372      longvarname, iovar.getncattr('units'))
     
    1237412376
    1237512377    newvar = onewnc.createVariable(varn + 'spaceweightmean', 'f4', tuple(dimsloop),  \
    12376       fill_value=fillValueF)
     12378      fill_value=gen.fillValueF)
    1237712379    basicvardef(newvar, varn + 'spaceweightmean', 'space weighted ' + varn + ' ' +   \
    1237812380      longvarname, iovar.getncattr('units'))
     
    1238812390    newdim = onewnc.createDimension(xdimname, outweightvals.shape[1])
    1238912391    newvar = onewnc.createVariable('spatialweight', 'f4',tuple([ydimname, xdimname]),\
    12390       fill_value=fillValueF)
     12392      fill_value=gen.fillValueF)
    1239112393    basicvardef(newvar, 'spatialweight', 'space weight ' + longvarname, '-')
    1239212394    newvar[:] = outweightvals
     
    1242112423    onewnc.sync()
    1242212424    onewnc.close()
    12423 
    1242412425
    1242512426    print fname + ": Successfull written of file: '" + ofile + "' !!"
Note: See TracChangeset for help on using the changeset viewer.