Ignore:
Timestamp:
Mar 26, 2013 6:16:16 PM (8 years ago)
Author:
hoffmann
Message:

optimization of two-moments cloud physics

File:
1 edited

Legend:

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

    r1054 r1115  
    2020! Current revisions:
    2121! -----------------
     22! unused variables removed
    2223!
    2324! Former revisions:
     
    283284!--    First compare the version numbers
    284285       READ ( 13 )  version_on_file
    285        binary_version = '3.3'
     286       binary_version = '3.4'
    286287       IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
    287288          WRITE( message_string, * ) 'version mismatch concerning data ',      &
     
    432433                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    433434
     435                CASE ( 'nr' )
     436                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     437                   nr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     438                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     439
     440                CASE ( 'nr_av' )
     441                   IF ( .NOT. ALLOCATED( nr_av ) )  THEN
     442                      ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     443                   ENDIF
     444                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     445                   nr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     446                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     447
     448                CASE ( 'nrs' )
     449                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     450                   nrs(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     451                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     452
     453                CASE ( 'nrsws' )
     454                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     455                   nrsws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     456                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     457
     458                CASE ( 'nrswst' )
     459                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     460                   nrswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     461                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     462
    434463                CASE ( 'p' )
    435464                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    516545                                     tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    517546
     547                CASE ( 'qc' )
     548                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     549                   qc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     550                                       tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     551
     552                CASE ( 'qc_av' )
     553                   IF ( .NOT. ALLOCATED( qc_av ) )  THEN
     554                      ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     555                   ENDIF
     556                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     557                   qc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     558                                       tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     559
    518560                CASE ( 'ql' )
    519561                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    534576                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    535577
    536                 CASE ( 'nr' )
    537                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    538                    nr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    539                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    540 
    541578                CASE ( 'qr_av' )
    542579                   IF ( .NOT. ALLOCATED( qr_av ) )  THEN
     
    547584                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    548585
    549                 CASE ( 'nr_av' )
    550                    IF ( .NOT. ALLOCATED( nr_av ) )  THEN
    551                       ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    552                    ENDIF
    553                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    554                    nr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    555                                     tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    556 
    557586                CASE ( 'qrs' )
    558587                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     
    568597                   IF ( k == 1 )  READ ( 13 )  tmp_2d
    569598                   qrswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    570                                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    571 
    572                 CASE ( 'nrs' )
    573                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    574                    nrs(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    575                                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    576 
    577                 CASE ( 'nrsws' )
    578                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    579                    nrsws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    580                                           tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    581 
    582                 CASE ( 'nrswst' )
    583                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    584                    nrswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    585599                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    586600
Note: See TracChangeset for help on using the changeset viewer.