Changeset 94 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Jun 1, 2007 3:25:22 PM (17 years ago)

Author:

raasch

Message:

preliminary uncomplete changes for ocean version

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (6 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Initial salinity profile is calculated.
+! z_max_do1d is checked only in case of ocean = .f.,
+! +initial temperature profile for the ocean version,
 !
 ! Former revisions:
@@ -153 +155 @@
        ENDIF
     ENDIF
+
+!
+!-- Check ocean setting
+    IF ( ocean )  THEN
+       action = ' '
+       IF ( timestep_scheme(1:8) == 'leapfrog' )  THEN
+          WRITE( action, '(A,A)' )  'timestep_scheme = ', timestep_scheme
+       ENDIF
+       IF ( action /= ' ' )  THEN
+          IF ( myid == 0 )  THEN
+             PRINT*, '+++ check_parameters:'
+             PRINT*, '    ocean = .T. does not allow ', TRIM( action )
+          ENDIF
+          CALL local_stop
+       ENDIF
+    ENDIF
+
 !
 !-- Check whether there are any illegal values
@@ -418 +437 @@
        v_init  = vg_surface
        pt_init = pt_surface
-       IF ( humidity )        q_init = q_surface
-       IF ( passive_scalar )  q_init = s_surface
+       IF ( humidity )        q_init  = q_surface
+       IF ( ocean )           sa_init = sa_surface
+       IF ( passive_scalar )  q_init  = s_surface
 
 !
@@ -504 +524 @@
 
 !
-!--    If required, compute initial temperature profile using the given
-!--    temperature gradient
+!--    Compute initial temperature profile using the given temperature gradients
        i = 1
        gradient = 0.0
-       pt_vertical_gradient_level_ind(1) = 0
-       DO  k = 1, nzt+1
-          IF ( pt_vertical_gradient_level(i) < zu(k)  .AND. &
-               pt_vertical_gradient_level(i) >= 0.0 )  THEN
-             gradient = pt_vertical_gradient(i) / 100.0
-             pt_vertical_gradient_level_ind(i) = k - 1
-             i = i + 1
-             IF ( i > 10 )  THEN
-                IF ( myid == 0 )  THEN
-                   PRINT*, '+++ check_parameters: upper bound 10 of array', &
-                           ' "pt_vertical_gradient_level_ind" exceeded'
-                ENDIF
-                CALL local_stop
-             ENDIF
-          ENDIF
-          IF ( gradient /= 0.0 )  THEN
-             IF ( k /= 1 )  THEN
-                pt_init(k) = pt_init(k-1) + dzu(k) * gradient
+
+       IF ( .NOT. ocean )  THEN
+
+          pt_vertical_gradient_level_ind(1) = 0
+          DO  k = 1, nzt+1
+             IF ( pt_vertical_gradient_level(i) < zu(k)  .AND. &
+                  pt_vertical_gradient_level(i) >= 0.0 )  THEN
+                gradient = pt_vertical_gradient(i) / 100.0
+                pt_vertical_gradient_level_ind(i) = k - 1
+                i = i + 1
+                IF ( i > 10 )  THEN
+                   IF ( myid == 0 )  THEN
+                      PRINT*, '+++ check_parameters: upper bound 10 of array', &
+                              ' "pt_vertical_gradient_level_ind" exceeded'
+                   ENDIF
+                   CALL local_stop
+                ENDIF
+             ENDIF
+             IF ( gradient /= 0.0 )  THEN
+                IF ( k /= 1 )  THEN
+                   pt_init(k) = pt_init(k-1) + dzu(k) * gradient
+                ELSE
+                   pt_init(k) = pt_surface   + 0.5 * dzu(k) * gradient
+                ENDIF
              ELSE
-                pt_init(k) = pt_init(k-1) + 0.5 * dzu(k) * gradient
-             ENDIF
-          ELSE
-             pt_init(k) = pt_init(k-1)
-          ENDIF
-       ENDDO
+                pt_init(k) = pt_init(k-1)
+             ENDIF
+          ENDDO
+
+       ELSE
+
+          pt_vertical_gradient_level_ind(1) = nzt+1
+          DO  k = nzt, 0, -1
+             IF ( pt_vertical_gradient_level(i) > zu(k)  .AND. &
+                  pt_vertical_gradient_level(i) <= 0.0 )  THEN
+                gradient = pt_vertical_gradient(i) / 100.0
+                pt_vertical_gradient_level_ind(i) = k + 1
+                i = i + 1
+                IF ( i > 10 )  THEN
+                   IF ( myid == 0 )  THEN
+                      PRINT*, '+++ check_parameters: upper bound 10 of array', &
+                              ' "pt_vertical_gradient_level_ind" exceeded'
+                   ENDIF
+                   CALL local_stop
+                ENDIF
+             ENDIF
+             IF ( gradient /= 0.0 )  THEN
+                IF ( k /= nzt )  THEN
+                   pt_init(k) = pt_init(k+1) - dzu(k+1) * gradient
+                ELSE
+                   pt_init(k)   = pt_surface - 0.5 * dzu(k+1) * gradient
+                   pt_init(k+1) = pt_surface + 0.5 * dzu(k+1) * gradient
+                ENDIF
+             ELSE
+                pt_init(k) = pt_init(k+1)
+             ENDIF
+          ENDDO
+
+       ENDIF
 
 !
 !--    In case of no given temperature gradients, choose gradient of neutral
 !--    stratification
-       IF ( pt_vertical_gradient_level(1) == -1.0 )  THEN
+       IF ( pt_vertical_gradient_level(1) == -9999999.9 )  THEN
           pt_vertical_gradient_level(1) = 0.0
        ENDIF
 
 !
-!--    Store temperature gradient at the top boundary for possile Neumann
+!--    Store temperature gradient at the top boundary for possible Neumann
 !--    boundary condition
        bc_pt_t_val = ( pt_init(nzt+1) - pt_init(nzt) ) / dzu(nzt+1)
@@ -609 +662 @@
        ENDIF
 
+!
+!--    If required, compute initial salinity profile using the given salinity
+!--    gradients
+       IF ( ocean )  THEN
+
+          i = 1
+          gradient = 0.0
+
+          sa_vertical_gradient_level_ind(1) = nzt+1
+          DO  k = nzt, 0, -1
+             IF ( sa_vertical_gradient_level(i) > zu(k)  .AND. &
+                  sa_vertical_gradient_level(i) <= 0.0 )  THEN
+                gradient = sa_vertical_gradient(i) / 100.0
+                sa_vertical_gradient_level_ind(i) = k + 1
+                i = i + 1
+                IF ( i > 10 )  THEN
+                   IF ( myid == 0 )  THEN
+                      PRINT*, '+++ check_parameters: upper bound 10 of array', &
+                              ' "sa_vertical_gradient_level_ind" exceeded'
+                   ENDIF
+                   CALL local_stop
+                ENDIF
+             ENDIF
+             IF ( gradient /= 0.0 )  THEN
+                IF ( k /= nzt )  THEN
+                   sa_init(k) = sa_init(k+1) - dzu(k+1) * gradient
+                ELSE
+                   sa_init(k)   = sa_surface - 0.5 * dzu(k+1) * gradient
+                   sa_init(k+1) = sa_surface + 0.5 * dzu(k+1) * gradient
+                ENDIF
+             ELSE
+                sa_init(k) = sa_init(k+1)
+             ENDIF
+          ENDDO
+
+       ENDIF
+
     ENDIF
 
@@ -1769 +1859 @@
 !-- y-value range of the coordinate system (PROFIL).
 !-- x-value range determined in plot_1d.
-    cross_uymin = 0.0
-    IF ( z_max_do1d == -1.0 )  THEN
-       cross_uymax = zu(nzt+1)
-    ELSEIF ( z_max_do1d < zu(nzb+1)  .OR.  z_max_do1d > zu(nzt+1) )  THEN
-       IF ( myid == 0 )  PRINT*, '+++ check_parameters:  z_max_do1d=',  &
-                                 z_max_do1d,' must be >= ', zu(nzb+1),  &
-                                 ' or <= ', zu(nzt+1)
-       CALL local_stop
-    ELSE
-       cross_uymax = z_max_do1d
+    IF ( .NOT. ocean )  THEN
+       cross_uymin = 0.0
+       IF ( z_max_do1d == -1.0 )  THEN
+          cross_uymax = zu(nzt+1)
+       ELSEIF ( z_max_do1d < zu(nzb+1)  .OR.  z_max_do1d > zu(nzt+1) )  THEN
+          IF ( myid == 0 )  PRINT*, '+++ check_parameters:  z_max_do1d=',  &
+                                    z_max_do1d,' must be >= ', zu(nzb+1),  &
+                                    ' or <= ', zu(nzt+1)
+          CALL local_stop
+       ELSE
+          cross_uymax = z_max_do1d
+       ENDIF
     ENDIF