Changeset 4029 for palm/trunk/SOURCE/pmc_interface_mod.f90

Timestamp:

Jun 14, 2019 2:04:35 PM (5 years ago)

Author:

raasch

Message:

bugfix: decycling of chemistry species after nesting data transfer, exchange of ghost points and boundary conditions separated for chemical species and SALSA module, nest_chemistry option removed, netcdf variable NF90_NOFILL is used as argument instead of 1 in calls to NF90_DEF_VAR_FILL

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (13 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! nest_chemistry switch removed
+! 
+! 4026 2019-06-12 16:50:15Z suehring
 ! Masked topography at boundary grid points in mass conservation, in order to 
 ! avoid that mean velocities within topography are imposed
@@ -465 +468 @@
 
     USE chemistry_model_mod,                                                   &
-        ONLY:  nest_chemistry, spec_conc_2
+        ONLY:  spec_conc_2
 
     USE cpulog,                                                                &
@@ -1429 +1432 @@
        ENDIF
        
-       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+       IF ( air_chemistry )  THEN
           DO n = 1, nspec
              CALL pmc_set_dataarray_name( 'parent', 'chem_' // TRIM( chem_species(n)%name ),        &
@@ -2546 +2549 @@
 !
 !-- Chemistry, depends on number of species
-    IF ( air_chemistry  .AND.  nest_chemistry )  pmc_max_array = pmc_max_array + nspec
+    IF ( air_chemistry )  pmc_max_array = pmc_max_array + nspec
 !
 !-- SALSA, depens on the number aerosol size bins and chemical components + 
@@ -2623 +2626 @@
     IF ( TRIM(name) == "s"    )  p_3d_sec => s_2
     IF ( TRIM(name) == "diss" )  p_3d_sec => diss_2
-    IF ( INDEX( TRIM(name), "chem_" ) /= 0      )  p_3d_sec => spec_conc_2(:,:,:,n)
-    IF ( INDEX( TRIM(name), "an_" ) /= 0  )  p_3d_sec => nconc_2(:,:,:,n)
-    IF ( INDEX( TRIM(name), "am_" ) /= 0 )  p_3d_sec => mconc_2(:,:,:,n)
-    IF ( INDEX( TRIM(name), "sg_" ) /= 0  .AND.  .NOT. salsa_gases_from_chem ) &
-       p_3d_sec => gconc_2(:,:,:,n)
+    IF ( INDEX( TRIM(name), "chem_" ) /= 0 )  p_3d_sec => spec_conc_2(:,:,:,n)
+    IF ( INDEX( TRIM(name), "an_" )   /= 0 )  p_3d_sec => nconc_2(:,:,:,n)
+    IF ( INDEX( TRIM(name), "am_" )   /= 0 )  p_3d_sec => mconc_2(:,:,:,n)
+    IF ( INDEX( TRIM(name), "sg_" )   /= 0  .AND.  .NOT. salsa_gases_from_chem ) &
+                                 p_3d_sec => gconc_2(:,:,:,n)
 
     IF ( ASSOCIATED( p_3d ) )  THEN
@@ -2953 +2956 @@
        ENDIF
 
-       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+       IF ( air_chemistry )  THEN
           DO  n = 1, nspec
              CALL pmci_interp_1sto_all ( chem_species(n)%conc, chem_spec_c(:,:,:,n),                &
@@ -3630 +3633 @@
              ENDIF
 
-             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+             IF ( air_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_1sto_lr( chem_species(n)%conc, chem_spec_c(:,:,:,n),            &
@@ -3701 +3704 @@
              ENDIF
 
-             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+             IF ( air_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_1sto_lr( chem_species(n)%conc, chem_spec_c(:,:,:,n),            &
@@ -3772 +3775 @@
              ENDIF
 
-             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+             IF ( air_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_1sto_sn( chem_species(n)%conc, chem_spec_c(:,:,:,n),            &
@@ -3843 +3846 @@
              ENDIF
 
-             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+             IF ( air_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_1sto_sn( chem_species(n)%conc, chem_spec_c(:,:,:,n),            &
@@ -3909 +3912 @@
        ENDIF
 
-       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+       IF ( air_chemistry )  THEN
           DO  n = 1, nspec
              CALL pmci_interp_1sto_t( chem_species(n)%conc, chem_spec_c(:,:,:,n),                   &
@@ -4000 +4003 @@
       ENDIF
 
-      IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+      IF ( air_chemistry )  THEN
          DO  n = 1, nspec
             CALL pmci_anterp_tophat( chem_species(n)%conc, chem_spec_c(:,:,:,n),                    &
@@ -5052 +5055 @@
 !
 !-- Set Neumann boundary conditions for chemical species
-    IF ( air_chemistry  .AND.  nest_chemistry )  THEN
+    IF ( air_chemistry )  THEN
        IF ( ibc_cs_b == 1 )  THEN
           DO  n = 1, nspec