Ignore:
Timestamp:
Jan 12, 2009 3:32:23 PM (13 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_cloud_physics.f90

    r98 r221  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: abort in case that absolute temperature is below zero
    77!
    88! Former revisions:
     
    6262    t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
    6363    DO  k = nzb, nzt+1
     64!
     65!--    Check temperature in case of too large domain height
     66       IF ( ( t_surface - g/cp * zu(k) ) < 0.0 )  THEN
     67          WRITE( message_string, * )  'absolute temperature < 0.0 at zu(', k, &
     68                                      ') = ', zu(k)
     69          CALL handle_palm_message( 'init_cloud_physics', 'PA0142', 1, 2, 0, 6, &
     70                                    0 )
     71       ENDIF
    6472       hyp(k)    = surface_pressure * 100.0 * &
    6573                   ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
Note: See TracChangeset for help on using the changeset viewer.