Bring coupling between ocean and wave components using vortex formulation

.gitmodules

@@ -4,8 +4,8 @@
   branch = develop
 [submodule "WW3"]
   path = WW3
-  url = https://github.com/oceanmodeling/WW3
-  branch = dev/ufs-coastal
+  url = https://github.com/DeniseWorthen/WW3
+  branch = feature/addflds4schism
 [submodule "stochastic_physics"]
   path = stochastic_physics
   url = https://github.com/NOAA-PSL/stochastic_physics
@@ -17,7 +17,8 @@
 [submodule "CMEPS"]
   path = CMEPS-interface/CMEPS
   url = https://github.com/oceanmodeling/CMEPS
-  branch = feature/coastalapp
+  #branch = feature/coastalapp
+  branch = feature/schism_3d
 [submodule "HYCOM"]
   path = HYCOM-interface/HYCOM
   url = https://github.com/NOAA-EMC/HYCOM-src
@@ -53,11 +54,13 @@
 [submodule "SCHISM"]
   path = SCHISM-interface/SCHISM
   url = https://github.com/oceanmodeling/schism
-  branch = master
+  #branch = master
+  branch = feature/schism_3d
 [submodule "SCHISM-ESMF"]
   path = SCHISM-interface/SCHISM-ESMF
   url = https://github.com/oceanmodeling/schism-esmf
-  branch = master
+  #branch = master
+  branch = feature/schism_3d
 [submodule "FVCOM"]
   path = FVCOM-interface/FVCOM
   url = https://github.com/oceanmodeling/FVCOM

CMakeLists.txt

@@ -81,7 +81,8 @@ set(OPENMP          ON  CACHE BOOL "Enable OpenMP threading")
 set(PARALLEL_NETCDF OFF CACHE BOOL "Enable parallel NetCDF")
 set(JEDI_DRIVER     OFF CACHE BOOL "Enable JEDI as top level driver")
 set(CMEPS_AOFLUX    OFF CACHE BOOL "Enable atmosphere-ocean flux calculation in mediator")
-set(PDLIB           OFF CACHE BOOL "Enable Domain Decomposition in WW3 via PDLIB")
+set(PDLIB           OFF CACHE BOOL "Enable Domain Decomposition in WW3 via PDLIB with BT1")
+set(PDLIB_BT4       OFF CACHE BOOL "Enable Domain Decomposition in WW3 via PDLIB with BT4")
 set(CDEPS_INLINE    OFF CACHE BOOL "Enable CDEPS inline capability")
 set(HYDRO           OFF CACHE BOOL "Enable hydrostatic set")
 
@@ -227,8 +228,12 @@ if(WW3)
     set(UFS_CAP "MULTI_ESMF" CACHE STRING "Build with MULTI_ESMF cap")
     list(APPEND _ufs_defs_private FRONT_WW3=WMESMFMD)
   else()
-    if(PDLIB)
-      set(SWITCH "meshcap_pdlib" CACHE STRING "NUOPC mesh cap")
+    if(PDLIB OR PDLIB_BT4)
+      if(PDLIB)
+        set(SWITCH "meshcap_pdlib" CACHE STRING "NUOPC mesh cap")
+      else()
+	set(SWITCH "meshcap_pdlib_bt4" CACHE STRING "NUOPC mesh cap")
+      endif()
     else()
       set(SWITCH "meshcap" CACHE STRING "NUOPC mesh cap")
     endif()

SCHISM-interface/CMakeLists.txt

@@ -83,10 +83,6 @@ endif()
 set(CMAKE_MODULE_PATH ${CMAKE_MODULE_PATH} "${CMAKE_CURRENT_SOURCE_DIR}/SCHISM/cmake/modules")
 add_subdirectory(SCHISM/src)
 
-if(BUILD_UTILS)
-  add_subdirectory(SCHISM/src/Utility)
-endif()
-
 list(APPEND schism_nuopc_src_files
   SCHISM-ESMF/src/schism/schism_nuopc_util.F90
   SCHISM-ESMF/src/schism/schism_nuopc_cap.F90

tests/default_vars.sh

@@ -844,6 +844,7 @@ export RST_BEG=${RUN_BEG}
 export RST_2_BEG=${RUN_BEG}
 export RST_END=${RUN_END}
 export RST_2_END=${RUN_END}
+export WW3_WLEV='F'
 export WW3_CUR='F'
 export WW3_ICE='F'
 export WW3_IC1='F'
@@ -1019,6 +1020,7 @@ export_ww3() {
   WW3_DTPNT="$(printf "%02d" $(( WW3_OUTDTHR*3600 )) )"
   export WW3_DTFLD
   export WW3_DTPNT
+  export WW3_WLEV='F'
   export WW3_CUR='C'
   export WW3_ICE='C'
   export WW3_IC1='F'
@@ -1430,6 +1432,7 @@ export_hafs_regional ()
   # default hafs with no ice
   export WW3_DOMAIN=natl_6m
   export WW3_MODDEF=mod_def.${WW3_DOMAIN}
+  export WW3_WLEV='F'
   export WW3_ICE='F'
   export WW3_OUTPARS="WND HS T01 T02 DIR FP DP PHS PTP PDIR UST CHA USP"
 

tests/parm/fd_ufs.yaml

@@ -919,6 +919,61 @@
      - standard_name: Sw_z0
        alias: z0rlen
      #
+     - standard_name: Sw_hs
+       alias: sea_surface_wave_significant_height
+       canonical_units: m
+       description: ww3 export - significant wave height
+     #
+     - standard_name: Sw_bhd
+       alias: sea_water_waves_effect_on_currents_bernoulli_head_adjustment
+       canonical_units: m2 s-2
+       description: ww3 export - Bernoulli head (J term)
+     #
+     - standard_name: Sw_tauox
+       alias: sea_surface_x_stress_due_to_waves
+       canonical_units: m2 s-2
+       description: ww3 export - wave to ocean momentum flux x
+     #
+     - standard_name: Sw_tauoy
+       alias: sea_surface_y_stress_due_to_waves
+       canonical_units: m2 s-2
+       description: ww3 export - wave to ocean momentum flux y
+     #
+     - standard_name: Sw_taubblx
+       alias: sea_bottom_upward_x_stress_due_to_waves
+       canonical_units: m2 s-2
+       description: ww3 export - momentum flux due to bottom friction x
+     #
+     - standard_name: Sw_taubbly
+       alias: sea_bottom_upward_y_stress_due_to_waves
+       canonical_units: m2 s-2
+       description: ww3 export - momentum flux due to bottom friction y
+     #
+     - standard_name: Sw_ubrx
+       alias: sea_bed_orbital_x_velocity_due_to_waves
+       canonical_units: m s-1
+       description: ww3 export - near bottom rms wave velocities x
+     #
+     - standard_name: Sw_ubry
+       alias: sea_bed_orbital_y_velocity_due_to_waves
+       canonical_units: m s-1
+       description: ww3 export - near bottom rms wave velocities y
+     #
+     - standard_name: Sw_thm
+       alias: sea_surface_wave_mean_direction
+       canonical_units: deg
+       description: ww3 export - mean wave direction
+     #
+     - standard_name: Sw_t0m1
+       alias: sea_surface_wave_mean_period
+       canonical_units: s
+       description: ww3 export - mean wave period
+     #
+     - standard_name: Sw_wnmean
+       alias: sea_surface_wave_mean_number
+       canonical_units: m-1
+       description: ww3 export - mean wave number
+     #
      - standard_name: Sw_pstokes_x
        alias: eastward_partitioned_stokes_drift_current
        canonical_units: m s-1
@@ -967,27 +1022,6 @@
      - standard_name: Sw_vstokes
        alias: ussy
      #
-     - standard_name: Sw_wbcuru
-       alias: eastward_wave_bottom_current
-       canonical_units: m s-1
-       description: ww3 export
-     - standard_name: Sw_wbcuru
-       alias: wbcuru
-     #
-     - standard_name: Sw_wbcurv
-       alias: northward_wave_bottom_current
-       canonical_units: m s-1
-       description: ww3 export
-     - standard_name: Sw_wbcurv
-       alias: wbcurv
-     #
-     - standard_name: Sw_wbcurp
-       alias: wave_bottom_current_period
-       canonical_units: s
-       description: ww3 export
-     - standard_name: Sw_wbcurp
-       alias: wbcurp
-     #
      - standard_name: Sw_wavsuu
        alias: eastward_wave_radiation_stress
        canonical_units: N m-1

tests/parm/ww3_shel.nml.IN

@@ -1,4 +1,5 @@
 &input_nml
+  input%forcing%water_levels = '@[WW3_WLEV]'
   input%forcing%winds      = 'C'
   input%forcing%currents   = '@[WW3_CUR]'
   input%forcing%ice_conc   = '@[WW3_ICE]'

tests/rt_coastal.conf

@@ -30,10 +30,10 @@ COMPILE | atm2fvc | intel | -DAPP=CSTLF -DCOORDINATE_TYPE=SPHERICAL -DWET_DRY=ON
 RUN | coastal_scituateharbor_atm2fvc | | baseline |
 
 ### SCHISM tests ###
-COMPILE | atm2sch | intel | -DAPP=CSTLS -DUSE_ATMOS=ON -DNO_PARMETIS=OFF -DOLDIO=ON -DBUILD_UTILS=ON | | fv3 |
+COMPILE | atm2sch | intel | -DAPP=CSTLS -DUSE_ATMOS=ON -DNO_PARMETIS=OFF -DOLDIO=ON | | fv3 |
 RUN | coastal_ike_shinnecock_atm2sch | | baseline |
 
-COMPILE | atm2sch2ww3 | intel | -DAPP=CSTLSW -DUSE_ATMOS=ON -DUSE_WW3=ON -DNO_PARMETIS=OFF -DOLDIO=ON -DPDLIB=ON -DBUILD_UTILS=ON | | fv3 |
+COMPILE | atm2sch2ww3 | intel | -DAPP=CSTLSW -DUSE_ATMOS=ON -DUSE_WW3=ON -DNO_PARMETIS=OFF -DOLDIO=ON -DPDLIB=ON | | fv3 |
 RUN | coastal_ike_shinnecock_atm2sch2ww3 | | baseline |
 
 ### ROMS tests ###
@@ -47,7 +47,7 @@ RUN | coastal_ike_shinnecock_atm2ww3 | | baseline |
 
 ### GNU Tests ###
 ### SCHISM tests ###
-COMPILE | atm2sch2ww3 | gnu | -DAPP=CSTLSW -DUSE_ATMOS=ON -DUSE_WW3=ON -DNO_PARMETIS=OFF -DOLDIO=ON -DPDLIB=ON -DBUILD_UTILS=ON | | fv3 |
+COMPILE | atm2sch2ww3 | gnu | -DAPP=CSTLSW -DUSE_ATMOS=ON -DUSE_WW3=ON -DNO_PARMETIS=OFF -DOLDIO=ON -DPDLIB=ON | | fv3 |
 RUN | coastal_ike_shinnecock_atm2sch2ww3 | | baseline |
 
 ### ROMS tests ###