Ignore:
Timestamp:
Jan 9, 2008 8:17:38 AM (14 years ago)
Author:
raasch
Message:

second preliminary update for turbulent inflow

File:
1 edited

Legend:

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

    r143 r145  
    55! -----------------
    66! After reading the global variables, unit 13 is closed for all PEs except PE0
    7 ! -cross_ts_*,
    8 ! hom_sum moved, volume_flow_area, volume_flow_initial moved from
     7! +numprocs_previous_run, hor_index_bounds_previous_run,
     8! -cross_ts_*, npex, npey,
     9! hom_sum, volume_flow_area, volume_flow_initial moved from
    910! read_3d_binary to here
    1011!
     
    7273!-- Make version number check first
    7374    READ ( 13 )  version_on_file
    74     binary_version = '3.1'
     75    binary_version = '3.2'
    7576    IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
    7677       IF ( myid == 0 )  THEN
     
    8485       CALL local_stop
    8586    ENDIF
     87
     88!
     89!-- Read number of PEs and horizontal index bounds of all PEs used in previous
     90!-- run
     91    READ ( 13 )  variable_chr
     92    IF ( TRIM( variable_chr ) /= 'numprocs' )  THEN
     93       PRINT*, '+++ read_var_list: numprocs not found in data from prior run', &
     94                    ' on PE ', myid
     95       CALL local_stop
     96    ENDIF
     97    READ ( 13 )  numprocs_previous_run
     98
     99    ALLOCATE( hor_index_bounds_previous_run(4,0:numprocs_previous_run-1) )
     100
     101    READ ( 13 )  variable_chr
     102    IF ( TRIM( variable_chr ) /= 'hor_index_bounds' )  THEN
     103       PRINT*, '+++ read_var_list: hor_index_bounds not found in data from ', &
     104                    'prior run on PE ', myid
     105       CALL local_stop
     106    ENDIF
     107    READ ( 13 )  hor_index_bounds_previous_run
    86108
    87109!
     
    124146    READ ( 13 )  statistic_regions
    125147    ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1), v_init(0:nz+1), &
    126               pt_init(0:nz+1), q_init(0:nz+1), sa_init(0:nz+1), &
    127               hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions) )
     148              pt_init(0:nz+1), q_init(0:nz+1), sa_init(0:nz+1),       &
     149              hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions),  &
     150              hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
    128151
    129152!
     
    283306          CASE ( 'netcdf_precision' )
    284307             READ ( 13 )  netcdf_precision
    285           CASE ( 'npex' )
    286              READ ( 13 )  npex
    287           CASE ( 'npey' )
    288              READ ( 13 )  npey
    289308          CASE ( 'nsor_ini' )
    290309             READ ( 13 )  nsor_ini
Note: See TracChangeset for help on using the changeset viewer.