Changeset 4550 for palm/trunk/SOURCE

Timestamp:

May 29, 2020 3:22:13 PM (4 years ago)

Author:

raasch

Message:

bugfix for reading local restart data

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! bugfix for reading local restart data
+! 
+! 4544 2020-05-21 14:43:05Z raasch
 ! conc_av changed from pointer to allocatable array, array spec_conc_av removed
 ! 
@@ -3090 +3093 @@
 
 
-    DO  lsp = 1, nspec
-
-       IF ( restart_string(1:length) == TRIM( chem_species(lsp)%name) )  THEN
-
-          IF ( k == 1 )  READ ( 13 )  tmp_3d
-          chem_species(lsp)%conc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                      &
+    IF ( ALLOCATED( chem_species ) )  THEN
+
+       DO  lsp = 1, nspec
+
+          IF ( restart_string(1:length) == TRIM( chem_species(lsp)%name) )  THEN
+
+             IF ( k == 1 )  READ ( 13 )  tmp_3d
+             chem_species(lsp)%conc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                   &
                                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
-          found = .TRUE.
-
-       ELSEIF (restart_string(1:length) == TRIM( chem_species(lsp)%name ) // '_av' )  THEN
-
-          IF ( .NOT. ALLOCATED( chem_species(lsp)%conc_av ) )  THEN
-             ALLOCATE( chem_species(lsp)%conc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ) )
+             found = .TRUE.
+
+          ELSEIF (restart_string(1:length) == TRIM( chem_species(lsp)%name ) // '_av' )  THEN
+
+             IF ( .NOT. ALLOCATED( chem_species(lsp)%conc_av ) )  THEN
+                ALLOCATE( chem_species(lsp)%conc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg ) )
+             ENDIF
+             IF ( k == 1 )  READ ( 13 )  tmp_3d
+             chem_species(lsp)%conc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                &
+                                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
+             found = .TRUE.
+
           ENDIF
-          IF ( k == 1 )  READ ( 13 )  tmp_3d
-          chem_species(lsp)%conc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                   &
-                                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
-          found = .TRUE.
-
-       ENDIF
-
-    ENDDO
-
+
+       ENDDO
+
+    ENDIF
 
  END SUBROUTINE chem_rrd_local_ftn