Changeset 4842 for palm/trunk/SOURCE/bulk_cloud_model_mod.f90

Timestamp:

Jan 14, 2021 10:42:28 AM (4 years ago)

Author:

raasch

Message:

reading of namelist file and actions in case of namelist errors revised so that statement labels and goto statements are not required any more, deprecated namelists removed

File:

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

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! reading of namelist file and actions in case of namelist errors revised so that statement labels
+! and goto statements are not required any more
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! Enable 3D data output also with 64-bit precision
 !
@@ -664 +668 @@
 !> Parin for &bulk_cloud_parameters for the bulk cloud module
 !--------------------------------------------------------------------------------------------------!
-    SUBROUTINE bcm_parin
-
-
-       IMPLICIT NONE
-
-       CHARACTER (LEN=80)  ::  line  !< dummy string that contains the current line of the parameter
-                                     !< file
-
-       NAMELIST /bulk_cloud_parameters/                                                            &
-          aerosol_bulk,                                                                            &
-          bulk_cloud_model,                                                                        &
-          c_sedimentation,                                                                         &
-          call_microphysics_at_all_substeps,                                                       &
-          cloud_scheme,                                                                            &
-          cloud_water_sedimentation,                                                               &
-          collision_turbulence,                                                                    &
-          curvature_solution_effects_bulk,                                                         &
-          dry_aerosol_radius,                                                                      &
-          graupel,                                                                                 &
-          ice_crystal_sedimentation,                                                               &
-          in_init,                                                                                 &
-          limiter_sedimentation,                                                                   &
-          microphysics_ice_phase,                                                                  &
-          na_init,                                                                                 &
-          nc_const,                                                                                &
-          sigma_bulk,                                                                              &
-          snow,                                                                                    &
-          start_ice_microphysics,                                                                  &
-          ventilation_effect
-
-       line = ' '
-!
-!--    Try to find bulk cloud module namelist
-       REWIND ( 11 )
-       line = ' '
-       DO   WHILE ( INDEX( line, '&bulk_cloud_parameters' ) == 0 )
-          READ ( 11, '(A)', END=10 )  line
-       ENDDO
-       BACKSPACE ( 11 )
-!
-!--    Read user-defined namelist
-       READ ( 11, bulk_cloud_parameters )
-!
-!--    Set flag that indicates that the bulk cloud module is switched on
-       !bulk_cloud_model = .TRUE.
-
-10     CONTINUE
-
-
-    END SUBROUTINE bcm_parin
+ SUBROUTINE bcm_parin
+
+
+    IMPLICIT NONE
+
+    CHARACTER(LEN=100)  ::  line  !< Dummy string that contains the current line of the parameter
+                                  !< file
+
+    INTEGER(iwp)  ::  io_status   !< Status after reading the namelist file
+
+
+    NAMELIST /bulk_cloud_parameters/                                                               &
+       aerosol_bulk,                                                                               &
+       bulk_cloud_model,                                                                           &
+       c_sedimentation,                                                                            &
+       call_microphysics_at_all_substeps,                                                          &
+       cloud_scheme,                                                                               &
+       cloud_water_sedimentation,                                                                  &
+       collision_turbulence,                                                                       &
+       curvature_solution_effects_bulk,                                                            &
+       dry_aerosol_radius,                                                                         &
+       graupel,                                                                                    &
+       ice_crystal_sedimentation,                                                                  &
+       in_init,                                                                                    &
+       limiter_sedimentation,                                                                      &
+       microphysics_ice_phase,                                                                     &
+       na_init,                                                                                    &
+       nc_const,                                                                                   &
+       sigma_bulk,                                                                                 &
+       snow,                                                                                       &
+       start_ice_microphysics,                                                                     &
+       ventilation_effect
+
+!
+!-- Move to the beginning of the namelist file and try to find and read the namelist.
+    REWIND( 11 )
+    READ( 11, bulk_cloud_parameters, IOSTAT=io_status )
+!
+!-- Action depending on the READ status
+    IF ( io_status == 0 )  THEN
+!
+!--    bulk_cloud_parameters namelist was found and read correctly. Set flag that
+!--    bulk_cloud_model_mod is switched on.
+       bulk_cloud_model = .TRUE.
+
+    ELSEIF ( io_status > 0 )  THEN
+!
+!--    bulk_cloud_parameters namelist was found, but contained errors. Print an error message
+!--    containing the line that caused the problem.
+       BACKSPACE( 11 )
+       READ( 11 , '(A)') line
+       CALL parin_fail_message( 'bulk_cloud_parameters', line )
+
+    ENDIF
+
+ END SUBROUTINE bcm_parin