Changeset 2756


Ignore:
Timestamp:
Jan 16, 2018 6:11:14 PM (7 years ago)
Author:
suehring
Message:

Fill values for 3D data output of chemical species introduced

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r2718 r2756  
    2727! -----------------
    2828! $Id$
     29! Fill values in 3D output introduced.
     30!
     31! 2718 2018-01-02 08:49:38Z maronga
    2932! Initial revision
    3033!
     
    12001203
    12011204
    1202    SUBROUTINE chem_data_output_3d( av, variable, found, local_pf )
     1205   SUBROUTINE chem_data_output_3d( av, variable, found, local_pf, fill_value )
    12031206
    12041207
     
    12131216      LOGICAL      ::  found !<
    12141217      INTEGER(iwp) ::  av    !<
    1215       REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !<
     1218
     1219      REAL(wp) ::  fill_value !<
     1220      REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf
     1221
    12161222
    12171223      !-- local variables
     
    12341240                   DO  j = nys, nyn
    12351241                      DO  k = nzb, nzt+1
    1236                           local_pf(i,j,k) = chem_species(lsp)%conc(k,j,i)
     1242                          local_pf(i,j,k) = MERGE(                             &
     1243                                              chem_species(lsp)%conc(k,j,i),   &
     1244                                              REAL( fill_value, KIND = wp ),   &
     1245                                              BTEST( wall_flags_0(k,j,i), 0 ) )
    12371246                      ENDDO
    12381247                   ENDDO
     
    12431252                   DO  j = nys, nyn
    12441253                      DO  k = nzb, nzt+1
    1245                           local_pf(i,j,k) = chem_species(lsp)%conc_av(k,j,i)
     1254                          local_pf(i,j,k) = MERGE(                             &
     1255                                              chem_species(lsp)%conc_av(k,j,i),&
     1256                                              REAL( fill_value, KIND = wp ),   &
     1257                                              BTEST( wall_flags_0(k,j,i), 0 ) )
    12461258                      ENDDO
    12471259                   ENDDO
  • palm/trunk/SOURCE/data_output_3d.f90

    r2746 r2756  
    2525! -----------------
    2626! $Id$
     27! Fill values for 3D output of chemical species introduced.
     28!
     29! 2746 2018-01-15 12:06:04Z suehring
    2730! Move flag plant canopy to modules
    2831!
     
    680683             IF ( .NOT. found  .AND.  air_chemistry )  THEN
    681684                CALL chem_data_output_3d( av, do3d(av,if), found,              &
    682                                           local_pf )
     685                                          local_pf, fill_value )
    683686                resorted = .TRUE.
    684687             ENDIF
Note: See TracChangeset for help on using the changeset viewer.