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

optimization of two-moments cloud physics

File:
1 edited

Legend:

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

    r1077 r1115  
    2020! Current revisions:
    2121! -----------------
     22! ql is calculated by calc_liquid_water_content
    2223!
    2324! Former revisions:
     
    483484             CASE ( 'qc_xy', 'qc_xz', 'qc_yz' )
    484485                IF ( av == 0 )  THEN
     486                   to_be_resorted => qc
     487                ELSE
     488                   to_be_resorted => qc_av
     489                ENDIF
     490                IF ( mode == 'xy' )  level_z = zu
     491
     492             CASE ( 'ql_xy', 'ql_xz', 'ql_yz' )
     493                IF ( av == 0 )  THEN
    485494                   to_be_resorted => ql
    486495                ELSE
    487496                   to_be_resorted => ql_av
    488                 ENDIF
    489                 IF ( mode == 'xy' )  level_z = zu
    490 
    491              CASE ( 'ql_xy', 'ql_xz', 'ql_yz' )
    492                 IF ( av == 0 )  THEN
    493                    IF ( cloud_physics  .AND.  icloud_scheme == 0 )  THEN
    494                       DO  i = nxlg, nxrg
    495                          DO  j = nysg, nyng
    496                             DO  k = nzb, nzt+1
    497                                local_pf(i,j,k) = ql(k,j,i) + qr(k,j,i)
    498                             ENDDO
    499                          ENDDO
    500                       ENDDO
    501                       resorted = .TRUE.
    502                    ELSE
    503                       to_be_resorted => ql
    504                    ENDIF
    505                 ELSE
    506                    IF ( cloud_physics  .AND.  icloud_scheme == 0 )  THEN
    507                       DO  i = nxlg, nxrg
    508                          DO  j = nysg, nyng
    509                             DO  k = nzb, nzt+1
    510                                local_pf(i,j,k) = ql_av(k,j,i) + qr_av(k,j,i)
    511                             ENDDO
    512                          ENDDO
    513                       ENDDO
    514                       resorted = .TRUE.
    515                    ELSE
    516                       to_be_resorted => ql_av
    517                    ENDIF
    518497                ENDIF
    519498                IF ( mode == 'xy' )  level_z = zu
Note: See TracChangeset for help on using the changeset viewer.