Changeset 343 for palm/trunk/SOURCE
- Timestamp:
- Jun 24, 2009 12:59:09 PM (16 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r336 r343 37 37 write_var_list) 38 38 39 check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, modules, mrun, package_parin, parin, prognostic_equations, read_var_list, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_dvrp, user_header, user_init_grid, write_var_list 39 Option for user defined 2D data output in xy cross sections at z=nzb+1 40 (data_output_2d, user_data_output_2d) 41 42 check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, modules, mrun, package_parin, parin, prognostic_equations, read_var_list, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_2d, user_data_output_dvrp, user_header, user_init_grid, write_var_list 40 43 41 44 New: set_particle_attributes, set_slicer_attributes_dvrp … … 91 94 Default value of grid_matching changed to strict. 92 95 93 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list 96 Adjustments for runs on lcxt4 (necessary due to an software update on CRAY) and 97 for coupled runs on ibmy (mrun, subjob) 98 99 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, mrun, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, subjob, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list 94 100 95 101 -
palm/trunk/SOURCE/check_parameters.f90
r332 r343 391 391 ELSEIF ( TRIM( coupling_mode ) == 'uncoupled' .AND. & 392 392 TRIM( coupling_char ) == '_O' ) THEN 393 393 394 ! 394 395 !-- Check whether an (uncoupled) atmospheric run has been declared as an 395 396 !-- ocean run (this setting is done via mrun-option -y) 397 396 398 message_string = 'ocean = .F. does not allow coupling_char = "' // & 397 399 TRIM( coupling_char ) // '" set by mrun-option "-y"' -
palm/trunk/SOURCE/data_output_2d.f90
r336 r343 8 8 ! Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0) 9 9 ! Output of messages replaced by message handling routine. 10 ! Output of user defined 2D (XY) arrays at z=nzb+1 is now possible 10 11 ! 11 12 ! … … 512 513 !-- User defined quantity 513 514 CALL user_data_output_2d( av, do2d(av,if), found, grid, & 514 local_pf )515 local_pf, two_d ) 515 516 resorted = .TRUE. 516 517 … … 519 520 ELSEIF ( grid == 'zw' ) THEN 520 521 IF ( mode == 'xy' ) level_z = zw 522 ELSEIF ( grid == 'zu1' ) THEN 523 IF ( mode == 'xy' ) level_z(nzb+1) = zu(nzb+1) 521 524 ENDIF 522 525 -
palm/trunk/SOURCE/flow_statistics.f90
r305 r343 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! Added missing timeseries calculation of w"q"(0), moved timeseries q* to the 7 ! end. 6 8 ! Temperature gradient criterion for estimating the boundary layer height 7 9 ! replaced by the gradient criterion of Sullivan et al. (1998). … … 1025 1027 ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr) ! u'w' at k=0 1026 1028 ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr) ! v'w' at k=0 1027 ts_value(21,sr) = hom(nzb +12,1,pr_palm,sr) ! q*1029 ts_value(21,sr) = hom(nzb,1,48,sr) ! w"q" at k=0 1028 1030 1029 1031 IF ( ts_value(5,sr) /= 0.0 ) THEN … … 1034 1036 ENDIF 1035 1037 1038 ts_value(23,sr) = hom(nzb+12,1,pr_palm,sr) ! q* 1036 1039 ! 1037 1040 !-- Calculate additional statistics provided by the user interface -
palm/trunk/SOURCE/modules.f90
r336 r343 3 3 4 4 !------------------------------------------------------------------------------! 5 ! Actualrevisions:5 ! Current revisions: 6 6 ! ----------------- 7 ! translation error of actual -> current revisions fixed 8 ! +q* in dots_label, dots_unit. increased dots_num respectively 7 9 ! typographical error in dots_unit fixed 8 10 ! +clip_dvrp_*, cluster_size, color_interval, dvrpsize_interval, dvrp_overlap, … … 824 826 replace_num = 6 825 827 826 INTEGER :: dots_num = 2 2828 INTEGER :: dots_num = 23 827 829 828 830 CHARACTER, DIMENSION( replace_num ) :: & … … 856 858 'vmax ', 'wmax ', 'div_new', 'div_old', 'z_i_wpt', 'z_i_pt ', & 857 859 'w* ', 'w"pt"0 ', 'w"pt" ', 'wpt ', 'pt(0) ', 'pt(zp) ', & 858 'w"u"0 ', 'w"v"0 ', 'w"q"0 ', 'mo_L ', 859 ( 'unknown', i9 = 1, 7 8) /)860 'w"u"0 ', 'w"v"0 ', 'w"q"0 ', 'mo_L ', 'q* ', & 861 ( 'unknown', i9 = 1, 77) /) 860 862 861 863 CHARACTER (LEN=7), DIMENSION(dots_max) :: dots_unit = & … … 863 865 'm/s ', 'm/s ', 's-1 ', 's-1 ', 'm ', 'm ', & 864 866 'm/s ', 'K m/s ', 'K m/s ', 'K m/s ', 'K ', 'K ', & 865 'm2/s2 ', 'm2/s2 ', 'kg m/s ', 'm ', 866 ( 'unknown', i9 = 1, 7 8) /)867 'm2/s2 ', 'm2/s2 ', 'kg m/s ', 'm ', 'kg/kg ', & 868 ( 'unknown', i9 = 1, 77 ) /) 867 869 868 870 CHARACTER (LEN=7), DIMENSION(300) :: dopr_unit = 'unknown' -
palm/trunk/SOURCE/time_integration.f90
r291 r343 116 116 !-- Determine size of next time step 117 117 IF ( simulated_time /= 0.0 ) CALL timestep 118 119 118 ! 120 119 !-- Execute the user-defined actions … … 334 333 IF ( coupling_mode /= 'uncoupled' .AND. run_coupled ) THEN 335 334 time_coupling = time_coupling + dt_3d 335 336 336 ! 337 337 !-- In case of model termination initiated by the local model -
palm/trunk/SOURCE/timestep.f90
r320 r343 55 55 REAL, DIMENSION(2) :: uv_gtrans, uv_gtrans_l 56 56 REAL, DIMENSION(nzb+1:nzt) :: dxyz2_min 57 58 57 59 58 CALL cpu_log( log_point(12), 'calculate_timestep', 'start' ) … … 385 384 386 385 ENDIF 387 388 386 CALL cpu_log( log_point(12), 'calculate_timestep', 'stop' ) 389 387 390 391 388 END SUBROUTINE timestep -
palm/trunk/SOURCE/user_data_output_2d.f90
r226 r343 1 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf )1 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf, two_d ) 2 2 3 3 !------------------------------------------------------------------------------! … … 29 29 INTEGER :: av, i, j, k 30 30 31 LOGICAL :: found 31 LOGICAL :: found, two_d 32 32 33 33 REAL, DIMENSION(nxl-1:nxr+1,nys-1:nyn+1,nzb:nzt+1) :: local_pf … … 64 64 ! grid = 'zu' 65 65 66 66 67 CASE DEFAULT 67 68 found = .FALSE.
Note: See TracChangeset
for help on using the changeset viewer.