Changeset 221 for palm/trunk/SOURCE/init_cloud_physics.f90

Timestamp:

Jan 12, 2009 3:32:23 PM (15 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

File:

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

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)