Ignore:
Timestamp:
Nov 13, 2012 5:11:03 PM (9 years ago)
Author:
hoffmann
Message:

two-moment cloud physics implemented

Location:
palm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk

  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/read_3d_binary.f90

    r1037 r1053  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! necessary expansions according to the two new prognostic equations (nr, qr)
     23! of the two-moment cloud physics scheme:
     24! +prr, prr_av, *, *_av, *s, *sws, *swst
    2325!
    2426! Former revisions:
     
    279281!--    First compare the version numbers
    280282       READ ( 13 )  version_on_file
    281        binary_version = '3.2'
     283       binary_version = '3.3'
    282284       IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
    283285          WRITE( message_string, * ) 'version mismatch concerning data ',      &
     
    457459                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    458460
     461                CASE ( 'prr' )
     462                   IF ( .NOT. ALLOCATED( prr ) )  THEN
     463                      ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     464                   ENDIF
     465                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     466                   prr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     467                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     468
     469                CASE ( 'prr_av' )
     470                   IF ( .NOT. ALLOCATED( prr_av ) )  THEN
     471                      ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     472                   ENDIF
     473                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     474                   prr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     475                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     476
    459477                CASE ( 'precipitation_amount' )
    460478                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     
    508526                   ql_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    509527                                       tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     528
     529                CASE ( 'qr' )
     530                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     531                   qr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     532                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     533
     534                CASE ( 'nr' )
     535                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     536                   nr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     537                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     538
     539                CASE ( 'qr_av' )
     540                   IF ( .NOT. ALLOCATED( qr_av ) )  THEN
     541                      ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     542                   ENDIF
     543                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     544                   qr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     545                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     546
     547                CASE ( 'nr_av' )
     548                   IF ( .NOT. ALLOCATED( nr_av ) )  THEN
     549                      ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     550                   ENDIF
     551                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     552                   nr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     553                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     554
     555                CASE ( 'qrs' )
     556                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     557                   qrs(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     558                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     559
     560                CASE ( 'qrsws' )
     561                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     562                   qrsws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     563                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     564
     565                CASE ( 'qrswst' )
     566                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     567                   qrswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     568                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     569
     570                CASE ( 'nrs' )
     571                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     572                   nrs(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     573                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     574
     575                CASE ( 'nrsws' )
     576                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     577                   nrsws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     578                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     579
     580                CASE ( 'nrswst' )
     581                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     582                   nrswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     583                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    510584
    511585                CASE ( 'ql_c_av' )
Note: See TracChangeset for help on using the changeset viewer.