Changeset 4550


Ignore:
Timestamp:
May 29, 2020 3:22:13 PM (4 years ago)
Author:
raasch
Message:

bugfix for reading local restart data

File:
1 edited

Legend:

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

    r4544 r4550  
    2727! -----------------
    2828! $Id$
     29! bugfix for reading local restart data
     30!
     31! 4544 2020-05-21 14:43:05Z raasch
    2932! conc_av changed from pointer to allocatable array, array spec_conc_av removed
    3033!
     
    30903093
    30913094
    3092     DO  lsp = 1, nspec
    3093 
    3094        IF ( restart_string(1:length) == TRIM( chem_species(lsp)%name) )  THEN
    3095 
    3096           IF ( k == 1 )  READ ( 13 )  tmp_3d
    3097           chem_species(lsp)%conc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                      &
     3095    IF ( ALLOCATED( chem_species ) )  THEN
     3096
     3097       DO  lsp = 1, nspec
     3098
     3099          IF ( restart_string(1:length) == TRIM( chem_species(lsp)%name) )  THEN
     3100
     3101             IF ( k == 1 )  READ ( 13 )  tmp_3d
     3102             chem_species(lsp)%conc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                   &
    30983103                                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    3099           found = .TRUE.
    3100 
    3101        ELSEIF (restart_string(1:length) == TRIM( chem_species(lsp)%name ) // '_av' )  THEN
    3102 
    3103           IF ( .NOT. ALLOCATED( chem_species(lsp)%conc_av ) )  THEN
    3104              ALLOCATE( chem_species(lsp)%conc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ) )
     3104             found = .TRUE.
     3105
     3106          ELSEIF (restart_string(1:length) == TRIM( chem_species(lsp)%name ) // '_av' )  THEN
     3107
     3108             IF ( .NOT. ALLOCATED( chem_species(lsp)%conc_av ) )  THEN
     3109                ALLOCATE( chem_species(lsp)%conc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ) )
     3110             ENDIF
     3111             IF ( k == 1 )  READ ( 13 )  tmp_3d
     3112             chem_species(lsp)%conc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                &
     3113                                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     3114             found = .TRUE.
     3115
    31053116          ENDIF
    3106           IF ( k == 1 )  READ ( 13 )  tmp_3d
    3107           chem_species(lsp)%conc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                   &
    3108                                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    3109           found = .TRUE.
    3110 
    3111        ENDIF
    3112 
    3113     ENDDO
    3114 
     3117
     3118       ENDDO
     3119
     3120    ENDIF
    31153121
    31163122 END SUBROUTINE chem_rrd_local_ftn
Note: See TracChangeset for help on using the changeset viewer.