Changeset 48 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Mar 6, 2007 12:28:36 PM (18 years ago)

Author:

raasch

Message:

preliminary version, several changes to be explained later

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (1 diff)

palm/trunk/SOURCE/flow_statistics.f90

@@ -801 +801 @@
        ENDIF
 
+!
+!--    Collect the time series quantities
+       ts_value(1,sr) = hom(nzb+4,1,var_hom,sr)     ! E
+       ts_value(2,sr) = hom(nzb+5,1,var_hom,sr)     ! E*
+       ts_value(3,sr) = dt_3d
+       ts_value(4,sr) = hom(nzb,1,var_hom,sr)       ! u*
+       ts_value(5,sr) = hom(nzb+3,1,var_hom,sr)     ! th*
+       ts_value(6,sr) = u_max
+       ts_value(7,sr) = v_max
+       ts_value(8,sr) = w_max
+       ts_value(9,sr) = hom(nzb+10,1,var_sum,sr)    ! new divergence
+       ts_value(10,sr) = hom(nzb+9,1,var_hom,sr)    ! old Divergence
+       ts_value(11,sr) = hom(nzb+6,1,var_hom,sr)    ! z_i(1)
+       ts_value(12,sr) = hom(nzb+7,1,var_hom,sr)    ! z_i(2)
+       ts_value(13,sr) = hom(nzb+8,1,var_hom,sr)    ! w*
+       ts_value(14,sr) = hom(nzb,1,16,sr)           ! w'pt'   at k=0
+       ts_value(15,sr) = hom(nzb+1,1,16,sr)         ! w'pt'   at k=1
+       ts_value(16,sr) = hom(nzb+1,1,18,sr)         ! wpt     at k=1
+       ts_value(17,sr) = hom(nzb,1,4,sr)            ! pt(0)
+       ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
+       ts_value(19,sr) = hom(nzb+9,1,var_hom-1,sr)  ! splptx
+       ts_value(20,sr) = hom(nzb+10,1,var_hom-1,sr) ! splpty
+       ts_value(21,sr) = hom(nzb+11,1,var_hom-1,sr) ! splptz
+       IF ( ts_value(5,sr) /= 0.0 )  THEN
+          ts_value(22,sr) = ts_value(4,sr)**2 / &
+                            ( kappa * g * ts_value(5,sr) / ts_value(18,sr) ) ! L
+       ELSE
+          ts_value(22,sr) = 10000.0
+       ENDIF
+
+!
+!--    Calculate additional statistics provided by the user interface
+       CALL user_statistics( sr )
+
     ENDDO    ! loop of the subregions
-
-!
-!-- Calculate additional statistics provided by the user interface
-    CALL user_statistics
 
 !