Changeset 1927 in lmdz_wrf for trunk/tools


Ignore:
Timestamp:
Jul 12, 2018, 11:50:47 PM (7 years ago)
Author:
lfita
Message:

Adding

  • `uavaFROMwswd': ua,va from ws,wd
Location:
trunk/tools
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/tools/diagnostics.inf

    r1909 r1927  
    4141ua, WRFua, U@V@SINALPHA@COSALPHA
    4242va, WRFva, U@V@SINALPHA@COSALPHA
     43uavaFwswd, uavaFROMwswd, ws@wd
    4344uavaz, WRFzwind, U@V@WRFz@U10@V10@SINALPHA@COSALPHA@z=100.
    4445uavaz, WRFzwind_log, U@V@WRFz@U10@V10@SINALPHA@COSALPHA@z=100.
  • trunk/tools/diagnostics.py

    r1909 r1927  
    8585availdiags = ['ACRAINTOT', 'accum', 'clt', 'cllmh', 'deaccum', 'fog_K84', 'fog_RUC', \
    8686  'LMDZrh', 'mslp', 'OMEGAw', 'RAINTOT',                                             \
    87   'rvors', 'td', 'turbulence', 'WRFcape_afwa', 'WRFclivi', 'WRFclwvi',               \
    88   'WRF_denszint', 'WRFgeop',                                                         \
     87  'rvors', 'td', 'turbulence', 'uavaFROMwswd', 'WRFcape_afwa', 'WRFclivi',           \
     88  'WRFclwvi', 'WRF_denszint', 'WRFgeop',                                             \
    8989  'WRFmrso', 'WRFpotevap_orPM', 'WRFp', 'WRFpsl_ecmwf',                              \
    9090  'WRFpsl_ptarget', 'WRFrvors', 'WRFslw', 'ws', 'wds', 'wss', 'WRFheight',           \
     
    876876          " Taylor decomposition turbulence term")
    877877
     878# ua va from ws wd (deg)
     879    elif diagn == 'uavaFROMwswd':
     880           
     881        var0 = ncobj.variables[depvars[0]][:]
     882        var1 = ncobj.variables[depvars[1]][:]
     883
     884        ua = var0*np.cos(var1*np.pi/180.)
     885        va = var0*np.sin(var1*np.pi/180.)
     886
     887        dnamesvar = ncobj.variables[depvars[0]].dimensions
     888        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
     889
     890        ncvar.insert_variable(ncobj, 'ua', ua, dnamesvar, dvnamesvar, newnc)
     891        ncvar.insert_variable(ncobj, 'va', va, dnamesvar, dvnamesvar, newnc)
     892
     893
    878894# WRFbils fom WRF as HFX + LH
    879895    elif diagn == 'WRFbils':
Note: See TracChangeset for help on using the changeset viewer.