Changeset 221 for palm/trunk/SOURCE

Timestamp:

Jan 12, 2009 3:32:23 PM (16 years ago)

Author:

raasch

Message:

new environment variable mpilib in configuration file, used by mrun, mbuild, and subjob; small bugfix concerning t<0 in init_cloud_physics

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (2 diffs)
• init_cloud_physics.f90 (modified) (2 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -9 +9 @@
 DVRP arguments changed to single precision, mode pathlines added.
 
-User can add additional routines in files user_additional_routines.f90
+User can add additional routines in files user_additional_routines
 
 User can check user parameters and deduce further quantities in user_check_parameters
@@ -76 +76 @@
 of grid points do not match (read_3d_binary)
 
-advec_particles, combine_plot_fields, data_output_2d, header, production_e, read_3d_binary
+Bugfix: abort in case that absolute temperature is below zero (init_cloud_physics)
+
+advec_particles, combine_plot_fields, data_output_2d, header, init_cloud_physics, production_e, read_3d_binary

palm/trunk/SOURCE/init_cloud_physics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: abort in case that absolute temperature is below zero
 !
 ! Former revisions:
@@ -62 +62 @@
     t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
     DO  k = nzb, nzt+1
+!
+!--    Check temperature in case of too large domain height
+       IF ( ( t_surface - g/cp * zu(k) ) < 0.0 )  THEN
+          WRITE( message_string, * )  'absolute temperature < 0.0 at zu(', k, &
+                                      ') = ', zu(k)
+          CALL handle_palm_message( 'init_cloud_physics', 'PA0142', 1, 2, 0, 6, &
+                                    0 )
+       ENDIF
        hyp(k)    = surface_pressure * 100.0 * &
                    ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)