Ignore:
Timestamp:
Dec 28, 2007 9:44:53 AM (14 years ago)
Author:
raasch
Message:

first preliminary update for turbulent inflow

File:
1 edited

Legend:

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

    r110 r143  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Restart file has to be re-opened on all PEs except PE0
     7! -dopr_time_count,
     8! hom_sum moved, volume_flow_area, volume_flow_initial moved to read_var_list,
     9! reading of old profil parameters (cross_..., dopr_crossindex, profile_***)
     10! removed
    711!
    812! Former revisions:
     
    5458    CHARACTER (LEN=10) ::  binary_version, version_on_file
    5559    CHARACTER (LEN=20) ::  field_chr
    56     CHARACTER (LEN=10),  DIMENSION(:), ALLOCATABLE ::  chdum10
    57     CHARACTER (LEN=40),  DIMENSION(:), ALLOCATABLE ::  chdum40
    58     CHARACTER (LEN=100), DIMENSION(:), ALLOCATABLE ::  chdum100
    5960
    6061    INTEGER ::  idum1, myid_on_file, numprocs_on_file, nxl_on_file, &
    6162                nxr_on_file, nyn_on_file, nys_on_file, nzb_on_file, nzt_on_file
    6263
    63     INTEGER, DIMENSION(:), ALLOCATABLE ::  idum
    64 
    65     REAL, DIMENSION(:), ALLOCATABLE ::  rdum
    66 
    6764!
    6865!-- Read data from previous model run. unit 13 already opened in parin
    6966    CALL cpu_log( log_point_s(14), 'read_3d_binary', 'start' )
     67
     68!
     69!-- Restart file has to be re-opened on all PEs except PE0
     70    CALL check_open( 13 )
    7071
    7172!
     
    173174       CALL local_stop
    174175    ENDIF
    175 
    176 !
    177 !-- Local arrays that may be required for possible temporary information
    178 !-- storage in the following
    179     ALLOCATE( chdum10(crmax), chdum40(crmax), chdum100(crmax), &
    180               idum(100*crmax), rdum(100*crmax) )
    181176
    182177!
     
    304299          CASE ( 'shf_m' )
    305300             READ ( 13 )  shf_m
    306           CASE ( 'tswst' )
    307              READ ( 13 )  tswst
    308           CASE ( 'tswst_m' )
    309              READ ( 13 )  tswst_m
    310301          CASE ( 'spectrum_x' )
    311302             READ ( 13 )  spectrum_x
     
    317308             ALLOCATE( ts_av(nys-1:nyn+1,nxl-1:nxr+1) )
    318309             READ ( 13 )  ts_av
     310          CASE ( 'tswst' )
     311             READ ( 13 )  tswst
     312          CASE ( 'tswst_m' )
     313             READ ( 13 )  tswst_m
    319314          CASE ( 'u' )
    320315             READ ( 13 )  u
     
    347342          CASE ( 'v' )
    348343             READ ( 13 )  v
    349           CASE ( 'volume_flow_area' )
    350              READ ( 13 )  volume_flow_area
    351           CASE ( 'volume_flow_initial' )
    352              READ ( 13 )  volume_flow_initial
    353344          CASE ( 'v_av' )
    354345             ALLOCATE( v_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     
    400391             READ ( 13 )  z0_av
    401392
    402           CASE ( 'cross_linecolors' )
    403              IF ( use_prior_plot1d_parameters )  THEN
    404                 READ ( 13 )  cross_linecolors
    405              ELSE
    406                 READ ( 13 )  idum
    407              ENDIF
    408           CASE ( 'cross_linestyles' )
    409              IF ( use_prior_plot1d_parameters )  THEN
    410                 READ ( 13 )  cross_linestyles
    411              ELSE
    412                 READ ( 13 )  idum
    413              ENDIF
    414           CASE ( 'cross_normalized_x' )
    415              IF ( use_prior_plot1d_parameters )  THEN
    416                 READ ( 13 )  cross_normalized_x
    417              ELSE
    418                 READ ( 13 )  chdum10
    419              ENDIF
    420           CASE ( 'cross_normalized_y' )
    421              IF ( use_prior_plot1d_parameters )  THEN
    422                 READ ( 13 )  cross_normalized_y
    423              ELSE
    424                 READ ( 13 )  chdum10
    425              ENDIF
    426           CASE ( 'cross_normx_factor' )
    427              IF ( use_prior_plot1d_parameters )  THEN
    428                 READ ( 13 )  cross_normx_factor
    429              ELSE
    430                 READ ( 13 )  rdum
    431              ENDIF
    432           CASE ( 'cross_normy_factor' )
    433              IF ( use_prior_plot1d_parameters )  THEN
    434                 READ ( 13 )  cross_normy_factor
    435              ELSE
    436                 READ ( 13 )  rdum
    437              ENDIF
    438           CASE ( 'cross_profiles' )
    439              IF ( use_prior_plot1d_parameters )  THEN
    440                 READ ( 13 )  cross_profiles
    441              ELSE
    442                 READ ( 13 )  chdum100
    443              ENDIF
    444           CASE ( 'cross_profile_n_coun' )
    445              IF ( use_prior_plot1d_parameters )  THEN
    446                 READ ( 13 )  cross_profile_number_count
    447              ELSE
    448                 READ ( 13 )  idum(1:crmax)
    449              ENDIF
    450           CASE ( 'cross_profile_number' )
    451              IF ( use_prior_plot1d_parameters )  THEN
    452                 READ ( 13 )  cross_profile_numbers
    453              ELSE
    454                 READ ( 13 )  idum
    455              ENDIF
    456           CASE ( 'cross_uxmax' )
    457              IF ( use_prior_plot1d_parameters )  THEN
    458                 READ ( 13 )  cross_uxmax
    459              ELSE
    460                 READ ( 13 )  rdum(1:crmax)
    461              ENDIF
    462           CASE ( 'cross_uxmax_computed' )
    463              IF ( use_prior_plot1d_parameters )  THEN
    464                 READ ( 13 )  cross_uxmax_computed
    465              ELSE
    466                 READ ( 13 )  rdum(1:crmax)
    467              ENDIF
    468           CASE ( 'cross_uxmax_normaliz' )
    469              IF ( use_prior_plot1d_parameters )  THEN
    470                 READ ( 13 )  cross_uxmax_normalized
    471              ELSE
    472                 READ ( 13 )  rdum(1:crmax)
    473              ENDIF
    474           CASE ( 'cross_uxmax_norm_com' )
    475              IF ( use_prior_plot1d_parameters )  THEN
    476                 READ ( 13 )  cross_uxmax_normalized_computed
    477              ELSE
    478                 READ ( 13 )  rdum(1:crmax)
    479              ENDIF
    480           CASE ( 'cross_uxmin' )
    481              IF ( use_prior_plot1d_parameters )  THEN
    482                 READ ( 13 )  cross_uxmin
    483              ELSE
    484                 READ ( 13 )  rdum(1:crmax)
    485              ENDIF
    486           CASE ( 'cross_uxmin_computed' )
    487              IF ( use_prior_plot1d_parameters )  THEN
    488                 READ ( 13 )  cross_uxmin_computed
    489              ELSE
    490                 READ ( 13 )  rdum(1:crmax)
    491              ENDIF
    492           CASE ( 'cross_uxmin_normaliz' )
    493              IF ( use_prior_plot1d_parameters )  THEN
    494                 READ ( 13 )  cross_uxmin_normalized
    495              ELSE
    496                 READ ( 13 )  rdum(1:crmax)
    497              ENDIF
    498           CASE ( 'cross_uxmin_norm_com' )
    499              IF ( use_prior_plot1d_parameters )  THEN
    500                 READ ( 13 )  cross_uxmin_normalized_computed
    501              ELSE
    502                 READ ( 13 )  rdum(1:crmax)
    503              ENDIF
    504           CASE ( 'cross_uymax' )
    505              IF ( use_prior_plot1d_parameters )  THEN
    506                 READ ( 13 )  cross_uymax
    507              ELSE
    508                 READ ( 13 )  rdum(1:crmax)
    509              ENDIF
    510           CASE ( 'cross_uymin' )
    511              IF ( use_prior_plot1d_parameters )  THEN
    512                 READ ( 13 )  cross_uymin
    513              ELSE
    514                 READ ( 13 )  rdum(1:crmax)
    515              ENDIF
    516           CASE ( 'cross_xtext' )
    517              IF ( use_prior_plot1d_parameters )  THEN
    518                 READ ( 13 )  cross_xtext
    519              ELSE
    520                 READ ( 13 )  chdum40
    521              ENDIF
    522           CASE ( 'dopr_crossindex' )
    523              IF ( use_prior_plot1d_parameters )  THEN
    524                 READ ( 13 )  dopr_crossindex
    525              ELSE
    526                 READ ( 13 )  idum(1:100)
    527              ENDIF
    528           CASE ( 'dopr_time_count' )
    529              IF ( use_prior_plot1d_parameters )  THEN
    530                 READ ( 13 )  dopr_time_count
    531              ELSE
    532                 READ ( 13 )  idum1
    533              ENDIF
    534           CASE ( 'hom_sum' )
    535              READ ( 13 )  hom_sum
    536           CASE ( 'profile_columns' )
    537              IF ( use_prior_plot1d_parameters )  THEN
    538                 READ ( 13 )  profile_columns
    539              ELSE
    540                 READ ( 13 )  idum1
    541              ENDIF
    542           CASE ( 'profile_number' )
    543              IF ( use_prior_plot1d_parameters )  THEN
    544                 READ ( 13 )  profile_number
    545              ELSE
    546                 READ ( 13 )  idum1
    547              ENDIF
    548           CASE ( 'profile_rows' )
    549              IF ( use_prior_plot1d_parameters )  THEN
    550                 READ ( 13 )  profile_rows
    551              ELSE
    552                 READ ( 13 )  idum1
    553              ENDIF
    554 
    555393          CASE DEFAULT
    556394             PRINT*, '+++ init_3d_model: unknown field named "', &
     
    566404    ENDDO
    567405
    568     DEALLOCATE( chdum10, chdum40, chdum100, idum, rdum )
    569 
    570406!
    571407!-- End of time measuring for reading binary data
Note: See TracChangeset for help on using the changeset viewer.