Changeset 4196 for palm/trunk/SOURCE/ocean_mod.f90

Timestamp:

Aug 29, 2019 11:02:06 AM (5 years ago)

Author:

gronemeier

Message:

Consider rotation of model domain for claculation of Coriolis force:

• check_parameters: Overwrite rotation_angle from namelist by value from static driver;
• coriolis: Consider rotation of model domain;
• header: Write information about rotation angle;
• modules: Added rotation_angle;
• netcdf_interface_mod: replaced rotation angle from input-netCDF file by namelist parameter 'rotation_angle';
• ocean_mod: Consider rotation of model domain for calculating the Stokes drift;
• parin: added rotation_angle to initialization_parameters;

File:

• palm/trunk/SOURCE/ocean_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/ocean_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Consider rotation of model domain for calculating the Stokes drift
+! 
+! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
 ! 
@@ -2126 +2129 @@
                v_stokes_zw, w, tend
 
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  pi
+
     USE control_parameters,                                                    &
-        ONLY:  f, fs, message_string
+        ONLY:  f, fs, message_string, rotation_angle
 
     USE grid_variables,                                                        &
@@ -2137 +2143 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  component  !< component of momentum equation
-    INTEGER(iwp) ::  i          !< loop index along x
-    INTEGER(iwp) ::  j          !< loop index along y
-    INTEGER(iwp) ::  k          !< loop index along z
-
+    INTEGER(iwp) ::  component      !< component of momentum equation
+    INTEGER(iwp) ::  i              !< loop index along x
+    INTEGER(iwp) ::  j              !< loop index along y
+    INTEGER(iwp) ::  k              !< loop index along z
+
+    REAL(wp)     ::  cos_rot_angle  !< cosine of model rotation angle
+    REAL(wp)     ::  sin_rot_angle  !< sine of model rotation angle
 
 !
@@ -2182 +2190 @@
 !--    w-component
        CASE ( 3 )
+
+!
+!--       Precalculate cosine and sine of rotation angle
+          cos_rot_angle = COS( rotation_angle * pi / 180.0_wp )
+          sin_rot_angle = SIN( rotation_angle * pi / 180.0_wp )
+
           DO  i = nxl, nxr
              DO  j = nys, nyn
@@ -2197 +2211 @@
                                                    + v(k+1,j+1,i) - v(k,j+1,i) &
                                                    ) * ddzu(k)  )              &
-                                           + fs * u_stokes_zw(k)
+                                           + fs * (                            &
+                                               sin_rot_angle * v_stokes_zw(k)  &
+                                             + cos_rot_angle * u_stokes_zw(k)  &
+                                                  )
                 ENDDO
              ENDDO
@@ -2226 +2243 @@
                v_stokes_zw, w, tend
 
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  pi
+
     USE control_parameters,                                                    &
-        ONLY:  f, fs, message_string
+        ONLY:  f, fs, message_string, rotation_angle
 
     USE grid_variables,                                                        &
@@ -2237 +2257 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  component  !< component of momentum equation
-    INTEGER(iwp) ::  i          !< loop index along x
-    INTEGER(iwp) ::  j          !< loop index along y
-    INTEGER(iwp) ::  k          !< loop incex along z
-
+    INTEGER(iwp) ::  component      !< component of momentum equation
+    INTEGER(iwp) ::  i              !< loop index along x
+    INTEGER(iwp) ::  j              !< loop index along y
+    INTEGER(iwp) ::  k              !< loop incex along z
+
+    REAL(wp)     ::  cos_rot_angle  !< cosine of model rotation angle
+    REAL(wp)     ::  sin_rot_angle  !< sine of model rotation angle
 
 !
@@ -2273 +2295 @@
 !--    w-component
        CASE ( 3 )
+
+!
+!--       Precalculate cosine and sine of rotation angle
+          cos_rot_angle = COS( rotation_angle * pi / 180.0_wp )
+          sin_rot_angle = SIN( rotation_angle * pi / 180.0_wp )
+
           DO  k = nzb+1, nzt
              tend(k,j,i) = tend(k,j,i) + u_stokes_zw(k) * (              &
@@ -2286 +2314 @@
                                                    + v(k+1,j+1,i) - v(k,j+1,i) &
                                                    ) * ddzu(k)  )              &
-                                       + fs * u_stokes_zw(k)
+                                       + fs * ( sin_rot_angle * v_stokes_zw(k) &
+                                              + cos_rot_angle * u_stokes_zw(k) &
+                                              )
           ENDDO