Ignore:
Timestamp:
May 30, 2017 5:47:52 PM (4 years ago)
Author:
suehring
Message:

Adjustments according new topography and surface-modelling concept implemented

File:
1 edited

Legend:

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

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Example code added for accessing output quantities stored on surface-data
     23! types
    2324!
    2425! Former revisions:
     
    6162    USE kinds
    6263
     64    USE surface_mod
     65
    6366    USE user
    6467
     
    6871    CHARACTER (LEN=*) ::  variable !<
    6972
    70     INTEGER(iwp) ::  av !<
    71     INTEGER(iwp) ::  i  !<
    72     INTEGER(iwp) ::  j  !<
    73     INTEGER(iwp) ::  k  !<
     73    INTEGER(iwp) ::  av     !< flag to control data output of instantaneous or time-averaged data
     74    INTEGER(iwp) ::  i      !< grid index along x-direction
     75    INTEGER(iwp) ::  j      !< grid index along y-direction
     76    INTEGER(iwp) ::  k      !< grid index along z-direction
     77    INTEGER(iwp) ::  m      !< running index surface elements
    7478    INTEGER(iwp) ::  nzb_do !< lower limit of the domain (usually nzb)
    7579    INTEGER(iwp) ::  nzt_do !< upper limit of the domain (usually nzt+1)
     
    110114!
    111115!          grid = 'zu'
     116!
     117!--    In case two-dimensional surface variables are outputted, the user
     118!--    has to access related surface-type. Uncomment and extend following lines
     119!--    appropriately (example output of vertical surface momentum flux of u-
     120!--    component). Please note, surface elements can be distributed over
     121!--    several data type, depending on their respective surface properties.
     122!       CASE ( 'usws_xy' )
     123!          IF ( av == 0 )  THEN
     124!
     125!--           Horizontal default-type surfaces
     126!             DO  m = 1, surf_def_h(0)%ns
     127!                i = surf_def_h(0)%i(m)
     128!                j = surf_def_h(0)%j(m)
     129!                local_pf(i,j,1) = surf_def_h(0)%usws(m)
     130!             ENDDO
     131!
     132!--           Horizontal natural-type surfaces
     133!             DO  m = 1, surf_lsm_h%ns
     134!                i = surf_lsm_h%i(m)
     135!                j = surf_lsm_h%j(m)
     136!                local_pf(i,j,1) = surf_lsm_h%usws(m)
     137!             ENDDO
     138!
     139!--           Horizontal urban-type surfaces
     140!             DO  m = 1, surf_usm_h%ns
     141!                i = surf_usm_h%i(m)
     142!                j = surf_usm_h%j(m)
     143!                local_pf(i,j,1) = surf_usm_h%usws(m)
     144!             ENDDO
     145!          ENDIF
     146!
     147!          grid = 'zu'
     148!--       
    112149
    113150
Note: See TracChangeset for help on using the changeset viewer.