Changeset 96 for palm/trunk/SOURCE/buoyancy.f90

Timestamp:

Jun 4, 2007 8:07:41 AM (17 years ago)

Author:

raasch

Message:

more preliminary uncomplete changes for ocean version

File:

• palm/trunk/SOURCE/buoyancy.f90 (modified) (7 diffs)

palm/trunk/SOURCE/buoyancy.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! calc_mean_pt_profile renamed calc_mean_profile
 !
 ! Former revisions:
@@ -29 +29 @@
 
     PRIVATE
-    PUBLIC buoyancy, calc_mean_pt_profile
+    PUBLIC buoyancy, calc_mean_profile
 
     INTERFACE buoyancy
@@ -36 +36 @@
     END INTERFACE buoyancy
 
-    INTERFACE calc_mean_pt_profile
-       MODULE PROCEDURE calc_mean_pt_profile
-    END INTERFACE calc_mean_pt_profile
+    INTERFACE calc_mean_profile
+       MODULE PROCEDURE calc_mean_profile
+    END INTERFACE calc_mean_profile
 
  CONTAINS
@@ -207 +207 @@
 
 
-    SUBROUTINE calc_mean_pt_profile( theta, pr )
+    SUBROUTINE calc_mean_profile( var, pr )
 
 !------------------------------------------------------------------------------!
@@ -224 +224 @@
 
        INTEGER ::  i, j, k, omp_get_thread_num, pr, tn
-       REAL, DIMENSION(:,:,:), POINTER  ::  theta
-
-!
-!--    Computation of the horizontally averaged temperature profile, unless
+       REAL, DIMENSION(:,:,:), POINTER  ::  var
+
+!
+!--    Computation of the horizontally averaged profile of variable var, unless
 !--    already done by the relevant call from flow_statistics. The calculation
 !--    is done only for the first respective intermediate timestep in order to
 !--    spare communication time and to produce identical model results with jobs
 !--    which are calling flow_statistics at different time intervals.
-!--    Although this calculation is not required for model runs with a slope,
-!--    it is nevertheless also computed.
        IF ( .NOT. flow_statistics_called  .AND. &
             intermediate_timestep_count == 1 )  THEN
 
 !
-!--       Horizontal average of the potential temperature
+!--       Horizontal average of variable var
           tn           =   0  ! Default thread number in case of one thread
           !$OMP PARALLEL PRIVATE( i, j, k, tn )
@@ -247 +245 @@
              DO  j =  nys, nyn
                 DO  k = nzb_s_outer(j,i), nzt+1
-                   sums_l(k,pr,tn) = sums_l(k,pr,tn) + theta(k,j,i)
+                   sums_l(k,pr,tn) = sums_l(k,pr,tn) + var(k,j,i)
                 ENDDO
              ENDDO
@@ -272 +270 @@
        ENDIF
 
-    END SUBROUTINE calc_mean_pt_profile
+    END SUBROUTINE calc_mean_profile
 
  END MODULE buoyancy_mod