Changeset 89 for palm/trunk/SOURCE


Ignore:
Timestamp:
May 25, 2007 12:08:31 PM (17 years ago)
Author:
raasch
Message:

further changes concerning user-defined profiles

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

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

    r87 r89  
    14051405          CASE ( 'Sw' )
    14061406             dopr_index(i) = 39
     1407             dopr_unit(i)  = 'none'
    14071408             hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )
    14081409
  • palm/trunk/SOURCE/modules.f90

    r87 r89  
    55! Actual revisions:
    66! -----------------
    7 ! +data_output_pr_user, max_pr_user, num_user_profiles, size of data_output_pr,
    8 ! dopr_index, dopr_initial_index and dopr_unit enlarged,
     7! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index,
     8! dopr_initial_index and dopr_unit enlarged,
    99! var_hom and var_sum renamed pr_palm
    1010!
     
    275275                mg_switch_to_pe0_level = 0, ngsrb = 2, nsor = 20, &
    276276                nsor_ini = 100, n_sor, normalizing_region = 0, &
    277                 num_user_profiles = 0, nz_do1d, nz_do3d = -9999, &
    278                 outflow_damping_width = -1, prt_time_count = 0, runnr = 0, &
    279                 skip_do_avs = 0, timestep_count = 0
     277                nz_do1d, nz_do3d = -9999, outflow_damping_width = -1, &
     278                prt_time_count = 0, runnr = 0, skip_do_avs = 0, &
     279                timestep_count = 0
    280280
    281281    INTEGER ::  dist_nxl(0:1), dist_nxr(0:1), dist_nyn(0:1), dist_nys(0:1), &
  • palm/trunk/SOURCE/user_interface.f90

    r87 r89  
    7070
    7171
    72     NAMELIST /userpar/  data_output_pr_user, data_output_user, max_pr_user, &
    73                         region
     72    NAMELIST /userpar/  data_output_pr_user, data_output_user, region
    7473
    7574!
     
    8988    READ ( 11, userpar )
    9089    user_defined_namelist_found = .TRUE.
    91 
    92 !
    93 !-- Determine the maximum number of user-profiles allowed to be calculated.
    94 !-- This number may be (much) larger than num_user_profiles which only gives
    95 !-- the number of user-defined profiles to be output.
    96     IF ( data_output_pr_user(1) /= ' '  .AND.  max_pr_user == 0 )  THEN
    97        max_pr_user = 100
    98     ENDIF
    9990
    10091!
     
    109100       DO  WHILE ( data_output_pr_user(j) /= ' '  .AND.  j <= 100 )
    110101          data_output_pr(i) = data_output_pr_user(j)
    111           num_user_profiles = num_user_profiles + 1
     102          max_pr_user       = max_pr_user + 1
    112103          i = i + 1
    113104          j = j + 1
    114105       ENDDO
    115     ENDIF
    116 
    117     IF ( num_user_profiles > max_pr_user )  THEN
    118        IF ( myid == 0 )  THEN
    119           PRINT*, '+++ user_parin: number of user-defined profiles given by '
    120           PRINT*, '                data_output_pr_user (', num_user_profiles, &
    121                                    ') is larger than the allowed maximum'
    122           PRINT*, '                number of profiles: max_pr_user(', &
    123                                    max_pr_user, ')'
    124        ENDIF
    125        CALL local_stop
    126106    ENDIF
    127107
     
    529509!             DO  k = nzb_s_outer(j,i)+1, nzt
    530510!!
    531 !!--             Sample on how to calculate the profile for vertical velocity
    532 !                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +       &
    533 !                                                 w(k,j,i)               &
     511!!--             Sample on how to calculate the profile of the resolved-scale
     512!!--             horizontal momentum flux u*v*
     513!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +           &
     514!                      ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * &
     515!                      ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) * &
    534516!                                                 * rmask(j,i,sr)
    535517!!
     
    770752!--    the number "pr_palm+1", the second one "pr_palm+2", etc. The respective
    771753!--    user-profile-numbers have also to be used in routine user_statistics!
    772 !       CASE ( 'w_user' )                  ! quantity string as given in
    773 !                                          ! data_output_pr_user
     754!       CASE ( 'u*v*' )                      ! quantity string as given in
     755!                                            ! data_output_pr_user
    774756!          index = pr_palm + 1
    775 !          dopr_index(var_count)  = index  ! the quantities' user-profile-number
    776 !          dopr_unit(var_count)   = 'm/s'  ! quantity unit
    777 !          hom(:,2,index,:)       = SPREAD( zw, 2, statistic_regions+1 )
    778 !                                          ! grid on which the quantity is
    779 !                                          ! defined (use zu or zw)
     757!          dopr_index(var_count)  = index    ! quantities' user-profile-number
     758!          dopr_unit(var_count)   = 'm2/s2'  ! quantity unit
     759!          hom(:,2,index,:)       = SPREAD( zu, 2, statistic_regions+1 )
     760!                                            ! grid on which the quantity is
     761!                                            ! defined (use zu or zw)
    780762
    781763       CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.