Changeset 4017 for palm/trunk/SOURCE/header.f90
- Timestamp:
- Jun 6, 2019 12:16:46 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/header.f90
r3761 r4017 1 ! > @file header.f901 ! !> @file header.f90 2 2 !------------------------------------------------------------------------------! 3 3 ! This file is part of the PALM model system. … … 437 437 ONLY: day_of_year_init, time_utc_init 438 438 439 USE dvrp_variables, &440 ONLY: use_seperate_pe_for_dvrp_output441 442 439 USE grid_variables, & 443 440 ONLY: dx, dy … … 462 459 sa_vertical_gradient, sa_vertical_gradient_level, & 463 460 sa_vertical_gradient_level_ind 464 465 USE particle_attributes, &466 ONLY: bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, collision_kernel, &467 curvature_solution_effects, &468 density_ratio, dissipation_classes, dt_min_part, dt_prel, &469 dt_write_particle_data, end_time_prel, &470 number_of_particle_groups, particle_advection, &471 particle_advection_start, &472 particles_per_point, pdx, pdy, pdz, psb, psl, psn, psr, pss, &473 pst, radius, radius_classes, random_start_position, &474 seed_follows_topography, &475 total_number_of_particles, use_sgs_for_particles, &476 vertical_particle_advection, write_particle_statistics477 461 478 462 USE pegrid … … 509 493 510 494 CHARACTER (LEN=40) :: output_format !< netcdf format 511 495 512 496 CHARACTER (LEN=70) :: char1 !< dummy varialbe used for various strings 513 497 CHARACTER (LEN=70) :: char2 !< string containing informating about the advected distance in case of Galilei transformation … … 651 635 WRITE ( io, 107 ) 'y' 652 636 ENDIF 653 IF ( use_seperate_pe_for_dvrp_output ) WRITE ( io, 105 )654 637 IF ( numprocs /= maximum_parallel_io_streams ) THEN 655 638 WRITE ( io, 108 ) maximum_parallel_io_streams … … 1836 1819 ENDIF 1837 1820 1838 #if defined( __dvrp_graphics )1839 !1840 !-- Dvrp-output1841 IF ( dt_dvrp /= 9999999.9_wp ) THEN1842 WRITE ( io, 360 ) dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &1843 TRIM( dvrp_username ), TRIM( dvrp_directory )1844 i = 11845 l = 01846 m = 01847 DO WHILE ( mode_dvrp(i) /= ' ' )1848 IF ( mode_dvrp(i)(1:10) == 'isosurface' ) THEN1849 READ ( mode_dvrp(i), '(10X,I2)' ) j1850 l = l + 11851 IF ( do3d(0,j) /= ' ' ) THEN1852 WRITE ( io, 361 ) TRIM( do3d(0,j) ), threshold(l), &1853 isosurface_color(:,l)1854 ENDIF1855 ELSEIF ( mode_dvrp(i)(1:6) == 'slicer' ) THEN1856 READ ( mode_dvrp(i), '(6X,I2)' ) j1857 m = m + 11858 IF ( do2d(0,j) /= ' ' ) THEN1859 WRITE ( io, 362 ) TRIM( do2d(0,j) ), &1860 slicer_range_limits_dvrp(:,m)1861 ENDIF1862 ENDIF1863 i = i + 11864 ENDDO1865 1866 WRITE ( io, 365 ) groundplate_color, superelevation_x, &1867 superelevation_y, superelevation, clip_dvrp_l, &1868 clip_dvrp_r, clip_dvrp_s, clip_dvrp_n1869 1870 IF ( TRIM( topography ) /= 'flat' ) THEN1871 WRITE ( io, 366 ) topography_color1872 IF ( cluster_size > 1 ) THEN1873 WRITE ( io, 367 ) cluster_size1874 ENDIF1875 ENDIF1876 1877 ENDIF1878 #endif1879 1880 1881 1821 WRITE ( io, 99 ) 1882 1822 … … 1912 1852 WRITE ( io, 431 ) 1913 1853 ENDIF 1914 IF ( humidity .AND. cloud_droplets ) THEN1915 WRITE ( io, 433 )1916 IF ( curvature_solution_effects ) WRITE ( io, 434 )1917 IF ( collision_kernel /= 'none' ) THEN1918 WRITE ( io, 435 ) TRIM( collision_kernel )1919 IF ( collision_kernel(6:9) == 'fast' ) THEN1920 WRITE ( io, 436 ) radius_classes, dissipation_classes1921 ENDIF1922 ELSE1923 WRITE ( io, 437 )1924 ENDIF1925 ENDIF1926 1927 1854 ! 1928 1855 !-- LES / turbulence parameters … … 1967 1894 WRITE ( io, 477 ) q_surface_initial_change 1968 1895 ENDIF 1969 1970 IF ( particle_advection ) THEN1971 !1972 !-- Particle attributes1973 WRITE ( io, 480 ) particle_advection_start, dt_prel, bc_par_lr, &1974 bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &1975 end_time_prel1976 IF ( use_sgs_for_particles ) WRITE ( io, 488 ) dt_min_part1977 IF ( random_start_position ) WRITE ( io, 481 )1978 IF ( seed_follows_topography ) WRITE ( io, 496 )1979 IF ( particles_per_point > 1 ) WRITE ( io, 489 ) particles_per_point1980 WRITE ( io, 495 ) total_number_of_particles1981 IF ( dt_write_particle_data /= 9999999.9_wp ) THEN1982 WRITE ( io, 485 ) dt_write_particle_data1983 IF ( netcdf_data_format > 1 ) THEN1984 output_format = 'netcdf (64 bit offset) and binary'1985 ELSE1986 output_format = 'netcdf and binary'1987 ENDIF1988 IF ( netcdf_deflate == 0 ) THEN1989 WRITE ( io, 344 ) output_format1990 ELSE1991 WRITE ( io, 354 ) TRIM( output_format ), netcdf_deflate1992 ENDIF1993 ENDIF1994 IF ( dt_dopts /= 9999999.9_wp ) WRITE ( io, 494 ) dt_dopts1995 IF ( write_particle_statistics ) WRITE ( io, 486 )1996 1997 WRITE ( io, 487 ) number_of_particle_groups1998 1999 DO i = 1, number_of_particle_groups2000 IF ( i == 1 .AND. density_ratio(i) == 9999999.9_wp ) THEN2001 WRITE ( io, 490 ) i, 0.0_wp2002 WRITE ( io, 492 )2003 ELSE2004 WRITE ( io, 490 ) i, radius(i)2005 IF ( density_ratio(i) /= 0.0_wp ) THEN2006 WRITE ( io, 491 ) density_ratio(i)2007 ELSE2008 WRITE ( io, 492 )2009 ENDIF2010 ENDIF2011 WRITE ( io, 493 ) psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &2012 pdx(i), pdy(i), pdz(i)2013 IF ( .NOT. vertical_particle_advection(i) ) WRITE ( io, 482 )2014 ENDDO2015 2016 ENDIF2017 2018 1896 2019 1897 ! … … 2055 1933 104 FORMAT (' Number of PEs:',10X,I6,4X,'Tasks:',I4,' threads per task:',I4/ & 2056 1934 35X,'Processor grid (x,y): (',I4,',',I4,')',1X,A) 2057 105 FORMAT (35X,'One additional PE is used to handle'/37X,'the dvrp output!')2058 1935 107 FORMAT (35X,'A 1d-decomposition along ',A,' is used') 2059 1936 108 FORMAT (35X,'Max. # of parallel I/O streams is ',I5) … … 2278 2155 353 FORMAT (/' Number of output time levels allowed: unlimited' /) 2279 2156 354 FORMAT (' Output format: ',A, ' compressed with level: ',I1/) 2280 #if defined( __dvrp_graphics )2281 360 FORMAT (' Plot-Sequence with dvrp-software:'/ &2282 ' Output every ',F7.1,' s'/ &2283 ' Output mode: ',A/ &2284 ' Host / User: ',A,' / ',A/ &2285 ' Directory: ',A// &2286 ' The sequence contains:')2287 361 FORMAT (/' Isosurface of "',A,'" Threshold value: ', E12.3/ &2288 ' Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')2289 362 FORMAT (/' Slicer plane ',A/ &2290 ' Slicer limits: [',F6.2,',',F6.2,']')2291 365 FORMAT (/' Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &2292 ' Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &2293 ')'/ &2294 ' Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &2295 ' from y = ',F9.1,' m to y = ',F9.1,' m')2296 366 FORMAT (/' Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')2297 367 FORMAT (' Polygon reduction for topography: cluster_size = ', I1)2298 #endif2299 2157 400 FORMAT (//' Physical quantities:'/ & 2300 2158 ' -------------------'/) … … 2350 2208 ' ----------------------------------'/) 2351 2209 431 FORMAT (' Humidity is considered, bu no condensation') 2352 433 FORMAT (' Cloud droplets treated explicitly using the Lagrangian part', &2353 'icle model')2354 434 FORMAT (' Curvature and solution effecs are considered for growth of', &2355 ' droplets < 1.0E-6 m')2356 435 FORMAT (' Droplet collision is handled by ',A,'-kernel')2357 436 FORMAT (' Fast kernel with fixed radius- and dissipation classes ', &2358 'are used'/ &2359 ' number of radius classes: ',I3,' interval ', &2360 '[1.0E-6,2.0E-4] m'/ &2361 ' number of dissipation classes: ',I2,' interval ', &2362 '[0,1000] cm**2/s**3')2363 437 FORMAT (' Droplet collision is switched off')2364 2210 450 FORMAT (//' LES / Turbulence quantities:'/ & 2365 2211 ' ---------------------------'/) … … 2401 2247 ' value is negative) by ',E8.1,' kg/m**3 at the beginning of', & 2402 2248 ' the 3D-simulation'/) 2403 480 FORMAT (' Particles:'/ &2404 ' ---------'// &2405 ' Particle advection is active (switched on at t = ', F7.1, &2406 ' s)'/ &2407 ' Start of new particle generations every ',F6.1,' s'/ &2408 ' Boundary conditions: left/right: ', A, ' north/south: ', A/&2409 ' bottom: ', A, ' top: ', A/&2410 ' Maximum particle age: ',F9.1,' s'/ &2411 ' Advection stopped at t = ',F9.1,' s'/)2412 481 FORMAT (' Particles have random start positions'/)2413 482 FORMAT (' Particles are advected only horizontally'/)2414 485 FORMAT (' Particle data are written on file every ', F9.1, ' s')2415 486 FORMAT (' Particle statistics are written on file'/)2416 487 FORMAT (' Number of particle groups: ',I2/)2417 488 FORMAT (' SGS velocity components are used for particle advection'/ &2418 ' minimum timestep for advection:', F8.5/)2419 489 FORMAT (' Number of particles simultaneously released at each ', &2420 'point: ', I5/)2421 490 FORMAT (' Particle group ',I2,':'/ &2422 ' Particle radius: ',E10.3, 'm')2423 491 FORMAT (' Particle inertia is activated'/ &2424 ' density_ratio (rho_fluid/rho_particle) =',F6.3/)2425 492 FORMAT (' Particles are advected only passively (no inertia)'/)2426 493 FORMAT (' Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&2427 ' y:',F8.1,' - ',F8.1,' m'/&2428 ' z:',F8.1,' - ',F8.1,' m'/&2429 ' Particle distances: dx = ',F8.1,' m dy = ',F8.1, &2430 ' m dz = ',F8.1,' m'/)2431 494 FORMAT (' Output of particle time series in NetCDF format every ', &2432 F8.2,' s'/)2433 495 FORMAT (' Number of particles in total domain: ',I10/)2434 496 FORMAT (' Initial vertical particle positions are interpreted ', &2435 'as relative to the given topography')2436 2249 500 FORMAT (//' 1D-Model parameters:'/ & 2437 2250 ' -------------------'// &
Note: See TracChangeset
for help on using the changeset viewer.