Ignore:
Timestamp:
Mar 6, 2007 12:28:36 PM (14 years ago)
Author:
raasch
Message:

preliminary version, several changes to be explained later

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r39 r48  
    801801       ENDIF
    802802
     803!
     804!--    Collect the time series quantities
     805       ts_value(1,sr) = hom(nzb+4,1,var_hom,sr)     ! E
     806       ts_value(2,sr) = hom(nzb+5,1,var_hom,sr)     ! E*
     807       ts_value(3,sr) = dt_3d
     808       ts_value(4,sr) = hom(nzb,1,var_hom,sr)       ! u*
     809       ts_value(5,sr) = hom(nzb+3,1,var_hom,sr)     ! th*
     810       ts_value(6,sr) = u_max
     811       ts_value(7,sr) = v_max
     812       ts_value(8,sr) = w_max
     813       ts_value(9,sr) = hom(nzb+10,1,var_sum,sr)    ! new divergence
     814       ts_value(10,sr) = hom(nzb+9,1,var_hom,sr)    ! old Divergence
     815       ts_value(11,sr) = hom(nzb+6,1,var_hom,sr)    ! z_i(1)
     816       ts_value(12,sr) = hom(nzb+7,1,var_hom,sr)    ! z_i(2)
     817       ts_value(13,sr) = hom(nzb+8,1,var_hom,sr)    ! w*
     818       ts_value(14,sr) = hom(nzb,1,16,sr)           ! w'pt'   at k=0
     819       ts_value(15,sr) = hom(nzb+1,1,16,sr)         ! w'pt'   at k=1
     820       ts_value(16,sr) = hom(nzb+1,1,18,sr)         ! wpt     at k=1
     821       ts_value(17,sr) = hom(nzb,1,4,sr)            ! pt(0)
     822       ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
     823       ts_value(19,sr) = hom(nzb+9,1,var_hom-1,sr)  ! splptx
     824       ts_value(20,sr) = hom(nzb+10,1,var_hom-1,sr) ! splpty
     825       ts_value(21,sr) = hom(nzb+11,1,var_hom-1,sr) ! splptz
     826       IF ( ts_value(5,sr) /= 0.0 )  THEN
     827          ts_value(22,sr) = ts_value(4,sr)**2 / &
     828                            ( kappa * g * ts_value(5,sr) / ts_value(18,sr) ) ! L
     829       ELSE
     830          ts_value(22,sr) = 10000.0
     831       ENDIF
     832
     833!
     834!--    Calculate additional statistics provided by the user interface
     835       CALL user_statistics( sr )
     836
    803837    ENDDO    ! loop of the subregions
    804 
    805 !
    806 !-- Calculate additional statistics provided by the user interface
    807     CALL user_statistics
    808838
    809839!
Note: See TracChangeset for help on using the changeset viewer.