Changeset 3302 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Oct 3, 2018 2:39:40 AM (6 years ago)

Author:

raasch

Message:

Craik-Leibovich force and wave breaking effects added to ocean mode, header output for ocean mode

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (12 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Stokes drift + wave breaking term added
+! 
+! 3298 2018-10-02 12:21:11Z kanani
 ! Code added for decycling chemistry (basit)
 !
@@ -392 +395 @@
 
     USE ocean_mod,                                                             &
-        ONLY:  eqn_state_seawater, ocean_prognostic_equations
+        ONLY:  ocean_prognostic_equations, stokes_drift_terms, stokes_force,   &
+               wave_breaking, wave_breaking_term
 
     USE plant_canopy_model_mod,                                                &
@@ -595 +599 @@
 
 !
+!--          Effect of Stokes drift (in ocean mode only)
+             IF ( stokes_force )  CALL stokes_drift_terms( i, j, 1 )
+
+!
 !--          Forces by wind turbines
              IF ( wind_turbine )  CALL wtm_tendencies( i, j, 1 )
@@ -612 +620 @@
                                                  )
              ENDDO
+
+!
+!--          Add turbulence generated by wave breaking (in ocean mode only)
+             IF ( wave_breaking  .AND.                                         &
+               intermediate_timestep_count == intermediate_timestep_count_max )&
+             THEN
+                CALL wave_breaking_term( i, j, 1 )
+             ENDIF
 
 !
@@ -666 +682 @@
 
 !
+!--          Effect of Stokes drift (in ocean mode only)
+             IF ( stokes_force )  CALL stokes_drift_terms( i, j, 2 )
+
+!
 !--          Forces by wind turbines
              IF ( wind_turbine )  CALL wtm_tendencies( i, j, 2 )
@@ -682 +702 @@
                                                  )
              ENDDO
+
+!
+!--          Add turbulence generated by wave breaking (in ocean mode only)
+             IF ( wave_breaking  .AND.                                         &
+               intermediate_timestep_count == intermediate_timestep_count_max )&
+             THEN
+                CALL wave_breaking_term( i, j, 2 )
+             ENDIF
 
 !
@@ -731 +759 @@
 !--       Drag by plant canopy
           IF ( plant_canopy )  CALL pcm_tendency( i, j, 3 )
+
+!
+!--       Effect of Stokes drift (in ocean mode only)
+          IF ( stokes_force )  CALL stokes_drift_terms( i, j, 3 )
 
 !
@@ -1383 +1415 @@
 
 !
+!-- Effect of Stokes drift (in ocean mode only)
+    IF ( stokes_force )  CALL stokes_drift_terms( 1 )
+
+!
 !-- Forces by wind turbines
     IF ( wind_turbine )  CALL wtm_tendencies( 1 )
@@ -1405 +1441 @@
 
 !
+!-- Add turbulence generated by wave breaking (in ocean mode only)
+    IF ( wave_breaking  .AND.                                                  &
+         intermediate_timestep_count == intermediate_timestep_count_max )      &
+    THEN
+       CALL wave_breaking_term( 1 )
+    ENDIF
+
+!
 !-- Calculate tendencies for the next Runge-Kutta step
     IF ( timestep_scheme(1:5) == 'runge' )  THEN
@@ -1466 +1510 @@
 !-- Nudging
     IF ( nudging )  CALL nudge( simulated_time, 'v' )
+
+!
+!-- Effect of Stokes drift (in ocean mode only)
+    IF ( stokes_force )  CALL stokes_drift_terms( 2 )
 
 !
@@ -1490 +1538 @@
 
 !
+!-- Add turbulence generated by wave breaking (in ocean mode only)
+    IF ( wave_breaking  .AND.                                                  &
+         intermediate_timestep_count == intermediate_timestep_count_max )      &
+    THEN
+       CALL wave_breaking_term( 2 )
+    ENDIF
+
+!
 !-- Calculate tendencies for the next Runge-Kutta step
     IF ( timestep_scheme(1:5) == 'runge' )  THEN
@@ -1547 +1603 @@
 !-- Drag by plant canopy
     IF ( plant_canopy )  CALL pcm_tendency( 3 )
+
+!
+!-- Effect of Stokes drift (in ocean mode only)
+    IF ( stokes_force )  CALL stokes_drift_terms( 3 )
 
 !