Changeset 221 for palm/trunk/SOURCE


Ignore:
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

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r217 r221  
    99DVRP arguments changed to single precision, mode pathlines added.
    1010
    11 User can add additional routines in files user_additional_routines.f90
     11User can add additional routines in files user_additional_routines
    1212
    1313User can check user parameters and deduce further quantities in user_check_parameters
     
    7676of grid points do not match (read_3d_binary)
    7777
    78 advec_particles, combine_plot_fields, data_output_2d, header, production_e, read_3d_binary
     78Bugfix: abort in case that absolute temperature is below zero (init_cloud_physics)
     79
     80advec_particles, combine_plot_fields, data_output_2d, header, init_cloud_physics, production_e, read_3d_binary
  • 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.