Changeset 3011 for palm


Ignore:
Timestamp:
May 7, 2018 2:38:42 PM (7 years ago)
Author:
schwenkel
Message:

removed redundant if statement

File:
1 edited

Legend:

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

    r3004 r3011  
    2525! -----------------
    2626! $Id$
     27! removed redundant if statement
     28!
     29! 3004 2018-04-27 12:33:25Z Giersch
    2730! precipitation_rate removed
    2831!
     
    691694       ALLOCATE( q(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                             &
    692695                 q_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
    693                  tq_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     696                 tq_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
     697                 vpt(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    694698#else
    695699       ALLOCATE( q_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
    696700                 q_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
    697                  q_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     701                 q_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
     702                 vpt_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    698703#endif
    699704
    700 !
    701 !--    3D-arrays needed for humidity
    702        IF ( humidity )  THEN
    703 #if defined( __nopointer )
    704           ALLOCATE( vpt(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    705 #else
    706           ALLOCATE( vpt_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    707 #endif
    708 
    709           IF ( cloud_physics )  THEN
     705       IF ( cloud_physics )  THEN
    710706!
    711707!--          Liquid water content
    712708#if defined( __nopointer )
    713              ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     709          ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    714710#else
    715              ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     711          ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    716712#endif
    717713
    718714!
    719 !--          3D-cloud water content
    720              IF ( .NOT. microphysics_morrison )  THEN
     715!--       3D-cloud water content
     716          IF ( .NOT. microphysics_morrison )  THEN
    721717#if defined( __nopointer )
    722                 ALLOCATE( qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     718             ALLOCATE( qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    723719#else
    724                 ALLOCATE( qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     720             ALLOCATE( qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    725721#endif
    726              ENDIF
    727 !
    728 !--          Precipitation amount and rate (only needed if output is switched)
    729              ALLOCATE( precipitation_amount(nysg:nyng,nxlg:nxrg) )
    730 
    731 !
    732 !--          3d-precipitation rate
    733              ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    734 
    735              IF ( microphysics_morrison )  THEN
    736 !
    737 !--             3D-cloud drop water content, cloud drop concentration arrays
     722          ENDIF
     723!
     724!--       Precipitation amount and rate (only needed if output is switched)
     725          ALLOCATE( precipitation_amount(nysg:nyng,nxlg:nxrg) )
     726
     727!
     728!--       3d-precipitation rate
     729          ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     730
     731          IF ( microphysics_morrison )  THEN
     732!
     733!--          3D-cloud drop water content, cloud drop concentration arrays
    738734#if defined( __nopointer )
    739                 ALLOCATE( nc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                &
    740                           nc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    741                           qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                &
    742                           qc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    743                           tnc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),             &
    744                           tqc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     735             ALLOCATE( nc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     736                       nc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     737                       qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     738                       qc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     739                       tnc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                   &
     740                       tqc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    745741#else
    746                 ALLOCATE( nc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    747                           nc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    748                           nc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    749                           qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    750                           qc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    751                           qc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     742             ALLOCATE( nc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     743                       nc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     744                       nc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     745                       qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     746                       qc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     747                       qc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    752748#endif
    753              ENDIF
    754 
    755              IF ( microphysics_seifert )  THEN
    756 !
    757 !--             3D-rain water content, rain drop concentration arrays
     749          ENDIF
     750
     751          IF ( microphysics_seifert )  THEN
     752!
     753!--          3D-rain water content, rain drop concentration arrays
    758754#if defined( __nopointer )
    759                 ALLOCATE( nr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                &
    760                           nr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    761                           qr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                &
    762                           qr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    763                           tnr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),             &
    764                           tqr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     755             ALLOCATE( nr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     756                       nr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     757                       qr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     758                       qr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     759                       tnr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                   &
     760                       tqr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    765761#else
    766                 ALLOCATE( nr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    767                           nr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    768                           nr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    769                           qr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    770                           qr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),              &
    771                           qr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     762             ALLOCATE( nr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     763                       nr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     764                       nr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     765                       qr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     766                       qr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                    &
     767                       qr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    772768#endif
    773              ENDIF
    774 
    775           ENDIF
    776 
    777           IF ( cloud_droplets )  THEN
    778 !
    779 !--          Liquid water content, change in liquid water content
     769          ENDIF
     770
     771       ENDIF
     772
     773       IF ( cloud_droplets )  THEN
     774!
     775!--       Liquid water content, change in liquid water content
    780776#if defined( __nopointer )
    781              ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                     &
    782                         ql_c(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     777          ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                        &
     778                     ql_c(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    783779#else
    784              ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                   &
    785                         ql_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     780          ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     781                     ql_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    786782#endif
    787783!
    788 !--          Real volume of particles (with weighting), volume of particles
    789              ALLOCATE ( ql_v(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                   &
    790                         ql_vp(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    791           ENDIF
    792 
    793        ENDIF
    794 
    795     ENDIF
    796    
     784!--       Real volume of particles (with weighting), volume of particles
     785          ALLOCATE ( ql_v(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                      &
     786                     ql_vp(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     787       ENDIF
     788
     789    ENDIF   
    797790   
    798791    IF ( passive_scalar )  THEN
Note: See TracChangeset for help on using the changeset viewer.