Changeset 4752 for palm


Ignore:
Timestamp:
Oct 21, 2020 1:54:51 PM (4 years ago)
Author:
suehring
Message:

virtual measurements: Remove unnecessary output and merge output for quantities that are represented by the same variable in PALM, e.g. surface temperature and brightness temperature

File:
1 edited

Legend:

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

    r4707 r4752  
    2525! -----------------
    2626! $Id$
     27! Remove unnecessary output and merge output for quantities that are represented by the same
     28! variable in PALM, e.g. surface temperature and brightness temperature
     29!
     30! 4707 2020-09-28 14:25:28Z suehring
    2731! Bugfix for previous commit + minor formatting adjustments
    2832!
     
    151155!> 'real-world' observation coordinates. More precisely, coordinates and measured quanties will be
    152156!> read from a NetCDF file which contains all required information. In the model, the same
    153 !> quantities (as long as all the required components are switched-on) will be sampled at the
     157!> quantities (as long as all the required PALM-modules are switched-on) will be sampled at the
    154158!> respective positions and output into an extra file, which allows for straight-forward comparison
    155159!> of model results with observations.
     
    368372                                                        'smp   ' /)
    369373
    370     CHARACTER(LEN=10), DIMENSION(0:1,1:8) ::  chem_vars = RESHAPE( (/ 'mcpm1     ', 'PM1       ',  &
     374    CHARACTER(LEN=10), DIMENSION(0:1,1:9) ::  chem_vars = RESHAPE( (/ 'mcpm1     ', 'PM1       ',  &
    371375                                                                      'mcpm2p5   ', 'PM2.5     ',  &
    372376                                                                      'mcpm10    ', 'PM10      ',  &
     
    375379                                                                      'mcno2     ', 'NO2       ',  &
    376380                                                                      'mcno      ', 'NO        ',  &
    377                                                                       'tro3      ', 'O3        '   &
    378                                                                     /), (/ 2, 8 /) )
     381                                                                      'tro3      ', 'O3        ',  &
     382                                                                      'ncaa      ', 'PM10      '   & !simply assume ncaa to be PM10
     383                                                                    /), (/ 2, 9 /) )
    379384
    380385    INTEGER(iwp) ::  maximum_name_length = 32  !< maximum name length of station names
     
    23792384                 ENDIF
    23802385
    2381               CASE ( 't_va' )
    2382 
    23832386              CASE ( 'hus' ) ! mixing ratio
    23842387                 IF ( humidity )  THEN
     
    25592562!--           PALM are in ppm and no distinction is made between mole-fraction and concentration
    25602563!--           quantities (all are output in ppm so far).
    2561               CASE ( 'mcpm1', 'mcpm2p5', 'mcpm10', 'mfno', 'mfno2', 'mcno', 'mcno2', 'tro3' )
     2564              CASE ( 'mcpm1', 'mcpm2p5', 'mcpm10', 'mfno', 'mfno2', 'mcno', 'mcno2', 'tro3', 'ncaa' )
    25622565                 IF ( air_chemistry )  THEN
    25632566!
     
    26812684                    ENDDO
    26822685                 ENDDO
    2683 
    2684               CASE ( 'lwcs' )  ! liquid water of soil layer
    2685 !                  DO  m = 1, vmea(l)%ns
    2686 ! !
    2687 ! !--                 Surface data is only available on inner subdomains, not on ghost points. Hence,
    2688 ! !--                 limit the indices.
    2689 !                     j = MERGE( vmea(l)%j(m), nys, vmea(l)%j(m) > nys )
    2690 !                     j = MERGE( j           , nyn, j            < nyn )
    2691 !                     i = MERGE( vmea(l)%i(m), nxl, vmea(l)%i(m) > nxl )
    2692 !                     i = MERGE( i           , nxr, i            < nxr )
    2693 !
    2694 !                     DO  mm = surf_lsm_h(0)%start_index(j,i), surf_lsm_h(0)%end_index(j,i)
    2695 !                        vmea(l)%measured_vars(m,n) = ?
    2696 !                     ENDDO
    2697 !                  ENDDO
    26982686
    26992687              CASE ( 'rnds' ) ! surface net radiation
     
    29112899                 ENDDO
    29122900
    2913               CASE ( 'm_soil' ) ! soil moisture
     2901              CASE ( 'm_soil', 'lwcs' ) ! soil moisture
    29142902                 DO  m = 1, vmea(l)%ns_soil
    29152903                    j = MERGE( vmea(l)%j_soil(m), nys, vmea(l)%j_soil(m) > nys )
     
    29282916                 ENDDO
    29292917
    2930               CASE ( 'ts' ) ! surface temperature
     2918              CASE ( 'ts', 'tb' ) ! surface temperature and brighness temperature
    29312919                 DO  m = 1, vmea(l)%ns
    29322920!
     
    29632951                 vmea(l)%measured_vars(:,n) = surface_pressure
    29642952
    2965               CASE ( 'pswrtg' ) ! platform speed above ground
    2966                  vmea(l)%measured_vars(:,n) = 0.0_wp
    2967 
    2968               CASE ( 'pswrta' ) ! platform speed in air
    2969                  vmea(l)%measured_vars(:,n) = 0.0_wp
    2970 
    29712953              CASE ( 't_lw' ) ! water temperature
    29722954                 DO  m = 1, vmea(l)%ns
     
    29862968                 ENDDO
    29872969!
    2988 !--           More will follow ...
    2989               CASE ( 'ncaa' )
    2990 !
    29912970!--           No match found - just set a fill value
    29922971              CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.