标签：angel dataset depth values attrs 瞎搞 vel delft3d

def add_dep(dataset):
    if 'LAYER_INTERFACE' in dataset:
        print("LAYER_INTERFACE property alread in DataSet, renaming it to depth...")
    else:
        print("Calculating depth of interfaces...")
        depth = dataset.S1
        depth.values = dataset.DP0.values[:] + dataset.S1.values[:]
        # depth_center =depth_center.assign_attrs({"unit": "m", "long_name": "Depth at Sigma-layer centers"})
        depth = depth.assign_attrs({"unit": "m", "long_name": "Depth at zeta"})
        dataset['dep'] = depth
    return dataset



def makevel(dataset, transpose=True):
    vel = dataset.U1
    angel = dataset.V1
    u = dataset.U1.values
    v =  dataset.V1.values
    vel.values = np.sqrt(u**2+v**2)
    angel.values = np.degrees(np.arctan(u/v))
    vel = vel.assign_attrs({"unit": "m/s", "long_name": "Depth at zeta"})
    angel = angel.assign_attrs({"unit": "°", "long_name": "Depth at zeta"})
    dataset['vel'] = vel
    dataset['angel'] = angel
    return dataset

标签：angel,dataset,depth,values,attrs,瞎搞,vel,delft3d
来源： https://www.cnblogs.com/chenyun-delft3d/p/14736567.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。