Ignore:
Timestamp:
Sep 13, 2012 2:08:46 PM (9 years ago)
Author:
raasch
Message:

leapfrog timestep scheme and upstream-spline advection scheme completely removed from the code,
reading of dt_fixed from restart file removed

File:
1 edited

Legend:

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

    r979 r1001  
    44! Current revisions:
    55! -----------------
    6 !
     6! all actions concerning leapfrog scheme removed
    77!
    88! Former revisions:
     
    376376                            tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    377377
    378                 CASE ( 'e_m' )
    379                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    380                    e_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    381                              tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    382 
    383378                CASE ( 'iran' ) ! matching random numbers is still unresolved
    384379                                ! issue
     
    390385                             tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    391386
    392                 CASE ( 'kh_m' )
    393                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    394                    kh_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    395                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    396 
    397387                CASE ( 'km' )
    398388                   IF ( k == 1 )  READ ( 13 )  tmp_3d
    399389                   km(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    400390                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    401 
    402                 CASE ( 'km_m' )
    403                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    404                    km_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    405                                 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    406391
    407392                CASE ( 'lpt_av' )
     
    476461                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    477462
    478                 CASE ( 'pt_m' )
    479                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    480                    pt_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    481                                     tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    482 
    483463                CASE ( 'q' )
    484464                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    494474                                     tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    495475
    496                 CASE ( 'q_m' )
    497                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    498                    q_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    499                                       tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    500 
    501476                CASE ( 'ql' )
    502477                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    544519                   IF ( k == 1 )  READ ( 13 )  tmp_2d
    545520                   qsws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    546                                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    547 
    548                 CASE ( 'qsws_m' )
    549                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    550                    qsws_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    551521                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    552522
     
    564534                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    565535
    566                 CASE ( 'qswst_m' )
    567                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    568                    qswst_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    569                                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    570 
    571536                CASE ( 'qv_av' )
    572537                   IF ( .NOT. ALLOCATED( qv_av ) )  THEN
     
    594559                                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    595560
    596                 CASE ( 'rif_m' )
    597                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    598                    rif_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    599                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    600 
    601561                CASE ( 'rif_wall' )
    602562                   IF ( k == 1 )  THEN
     
    644604                      tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    645605
    646                 CASE ( 'shf_m' )
    647                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    648                    shf_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    649                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    650606                CASE ( 'shf_av' )
    651607                   IF ( .NOT. ALLOCATED( shf_av ) )  THEN
     
    655611                   shf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    656612                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     613
    657614                CASE ( 'spectrum_x' )
    658615                   IF ( k == 1 )  THEN
     
    703660                         tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    704661
    705                 CASE ( 'tswst_m' )
    706                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    707                    tswst_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    708                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    709 
    710662                CASE ( 'u' )
    711663                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    721673                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    722674
    723                 CASE ( 'u_m' )
    724                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    725                    u_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    726                                  tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    727 
    728675                CASE ( 'u_m_l' )
    729676                   IF ( k == 1 )  THEN
     
    781728                        tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    782729
    783                 CASE ( 'usws_m' )
    784                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    785                    usws_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    786                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    787 
    788                 CASE ( 'uswst_m' )
    789                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    790                    uswst_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    791                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    792 
    793730                CASE ( 'us_av' )
    794731                   IF ( .NOT. ALLOCATED( us_av ) )  THEN
     
    812749                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    813750
    814                 CASE ( 'v_m' )
    815                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    816                    v_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    817                                 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    818 
    819751                CASE ( 'v_m_l' )
    820752                   IF ( k == 1 )  THEN
     
    870802                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    871803
    872                 CASE ( 'vpt_m' )
    873                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    874                    vpt_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    875                                 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    876 
    877804                CASE ( 'vsws' )
    878805                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     
    884811                   vswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    885812                        tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    886 
    887                 CASE ( 'vsws_m' )
    888                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    889                    vsws_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    890                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    891 
    892                 CASE ( 'vswst_m' )
    893                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    894                    vswst_m(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    895                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    896813
    897814                CASE ( 'w' )
     
    907824                   w_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    908825                               tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    909 
    910                 CASE ( 'w_m' )
    911                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    912                    w_m(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    913                                 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    914826
    915827                CASE ( 'w_m_l' )
Note: See TracChangeset for help on using the changeset viewer.