Different VIMS_SELFE runs have different wave direction variables

[rsignell-usgs]

It appears that some of the Testbed 1 Inundation Tropical VIMS_SELFE runs are using wwm_7 for wave direction, and others are using wwm_15 for wave direction (these wwm_15 output files were renamed pwd at some point). I think wwm_7 is mean_wave_direction while wwm_15 is peak_wave_direction. @drf5n, is this correct?

[testbed@testbed2 VIMS_SELFE]$ pwd
/data/comt_1_archive/inundation_tropical/VIMS_SELFE

[testbed@testbed2 VIMS_SELFE]$ find . -name '1_wwm_7*' -print
./Hurricane_Ike_2D_final_run_with_waves/Output/1_wwm_7.nc
./Hurricane_Rita_2D_final_run_with_waves/Output/1_wwm_7.nc

[testbed@testbed2 VIMS_SELFE]$ find . -name '1_pwd*' -print
./Hurricane_Ike_3D_final_run_with_waves/Output/1_pwd.nc
./Hurricane_Rita_3D_final_run_with_waves/Output/1_pwd.nc
./Hurricane_Rita_2D_final_run_waves_only/Output/1_pwd.nc
./Hurricane_Ike_2D_final_run_wave_only/Output/1_pwd.nc

[drf5n]

Yes, I think that is correct. Some of our old post-processing code was renaming the files, while some was not, and I think there was some confusion about which wave direction folks were using.

Here's a chunk from recent SELFE-WWM code, with wwm_7 being mean direction, and I think wwm_15 was renumbered wwm_16 Peak (dominant) wave direction

!-----------------------------------------------------------------------
! Outputs from WWM (USE_WWM must be on in Makefile)
!-----------------------------------------------------------------------
  wwm_1.61  = 0 !sig. height (m)
  wwm_2.61  = 0 !Mean average period (sec) - TM01
  wwm_3.61  = 0 !Zero down crossing period for comparison with buoy (s) - TM02
  wwm_4.61  = 0 !Average period of wave runup/overtopping - TM10
  wwm_5.61  = 0 !Mean wave number (1/m)
  wwm_6.61  = 0 !Mean wave length (m)
  wwm_7.61  = 0 !Mean average energy transport direction (deg)
  wwm_8.61  = 0 !Mean directional spreading (deg)
  wwm_9.61  = 0 !Discrete peak period (sec)
  wwm_10.61 = 0 !Continues peak period based on higher order moments (sec) 
  wwm_11.61 = 0 !Peak phase vel. (m/s)
  wwm_12.61 = 0 !Peak n-factor [-]
  wwm_13.61 = 0 !Peak group vel. (m/s)
  wwm_14.61 = 0 !Peak wave number (1/m)
  wwm_15.61 = 0 !Peak wave length (m)
  wwm_16.61 = 0 !Peak (dominant) wave direction (degr) ... some buoys record this 
  wwm_17.61 = 0 !Peak directional spreading (deg) ... some buoys record this  
  wwm_18.61 = 0 !Discrete peak direction (deg) ... some buoys record this  
  wwm_19.61 = 0 !Orbital vel. (m/s)
  wwm_20.61 = 0 !RMS orbital vel. (m/s)
  wwm_21.61 = 0 !Bottom excursion period (sec)
  wwm_22.61 = 0 !bottom wave period (sec)
  wwm_23.61 = 0 !Ursell number based on peak period 
  wwm_24.61 = 0 !none 
  wwm_25.62 = 0 !Etot energy vector (m^2)
  wwm_26.62 = 0 !none 

[rsignell-usgs]

@drf5n , okay, thanks for clarifying the mud on this issue. ;-)