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

Timestamp:

Aug 27, 2018 4:58:37 PM (6 years ago)

Author:

suehring

Message:

Additional namelist parameter to switch on/off the nesting of chemical species; Enable the input of soil data from dynamic input files independent on atmosphere in order to initialize soil properties in nested child domains from dynamic input; Revise error message number for static/dynamic input; Revise and add checks for static/dynamic input; Bugfix, add netcdf4_parallel directive for collective read operation

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Apply an additional switch controlling the nesting of chemical species
 ! 
 ! Former revisions:
@@ -327 +327 @@
 
     USE chemistry_model_mod,                                                   &
-        ONLY:  chem_species, spec_conc_2
+        ONLY:  chem_species, nest_chemistry, spec_conc_2
 
     USE cpulog,                                                                &
@@ -1291 +1291 @@
        ENDIF
        
-       IF ( air_chemistry )  THEN
+       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
           DO  n = 1, nspec
              CALL pmc_set_dataarray_name( 'coarse',                            &
@@ -3515 +3515 @@
     IF ( TRIM(name) == "s"          )  p_3d => s
     IF ( TRIM(name) == "diss"       )  p_3d => diss   
-    IF ( TRIM(name) == "nr_part"    )   i_2d => nr_part
+    IF ( TRIM(name) == "nr_part"    )  i_2d => nr_part
     IF ( TRIM(name) == "part_adr"   )  i_2d => part_adr
     IF ( INDEX( TRIM(name), "chem_" ) /= 0 )  p_3d => chem_species(n)%conc
@@ -3873 +3873 @@
        ENDIF
 
-       IF ( air_chemistry )  THEN
+       IF ( air_chemistry  .AND.  nest_chemistry)  THEN
           DO  n = 1, nspec
              CALL pmci_interp_tril_all ( chem_species(n)%conc,                 &
@@ -4632 +4632 @@
              ENDIF
 
-             IF ( air_chemistry )  THEN
+             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_tril_lr( chem_species(n)%conc,             &
@@ -4746 +4746 @@
              ENDIF
 
-             IF ( air_chemistry )  THEN
+             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_tril_lr( chem_species(n)%conc,             &
@@ -4857 +4857 @@
              ENDIF
 
-             IF ( air_chemistry )  THEN
+             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_tril_sn( chem_species(n)%conc,             &
@@ -4968 +4968 @@
              ENDIF
 
-             IF ( air_chemistry )  THEN
+             IF ( air_chemistry  .AND.  nest_chemistry )  THEN
                 DO  n = 1, nspec
                    CALL pmci_interp_tril_sn( chem_species(n)%conc,             &
@@ -5041 +5041 @@
        ENDIF
 
-       IF ( air_chemistry )  THEN
+       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
           DO  n = 1, nspec
              CALL pmci_interp_tril_t( chem_species(n)%conc,                    &
@@ -5124 +5124 @@
       ENDIF
 
-      IF ( air_chemistry )  THEN
+      IF ( air_chemistry  .AND.  nest_chemistry )  THEN
          DO  n = 1, nspec
             CALL pmci_anterp_tophat( chem_species(n)%conc,                     &
@@ -5987 +5987 @@
 !
 !-- Set Neumann boundary conditions for chemical species
-    IF ( air_chemistry )  THEN
+    IF ( air_chemistry  .AND.  nest_chemistry )  THEN
        IF ( ibc_cs_b == 1 )  THEN
           DO  n = 1, nspec