Changeset 4521 for palm


Ignore:
Timestamp:
May 6, 2020 11:39:49 AM (7 months ago)
Author:
schwenkel
Message:

add error number

Location:
palm/trunk/SOURCE
Files:
6 edited

Legend:

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

    r4517 r4521  
    2525! -----------------
    2626! $Id$
     27! Rename variable
     28!
     29! 4517 2020-05-03 14:29:30Z raasch
    2730! added restart with MPI-IO for reading local arrays
    2831!
     
    197200                surf_bulk_cloud_model,                                         &
    198201                surf_microphysics_morrison, surf_microphysics_seifert,         &
    199                 surf_microphysics_ice_extension,                               &
     202                surf_microphysics_ice_phase,                               &
    200203                surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h,    &
    201204                surf_usm_v
     
    220223
    221224    LOGICAL ::  call_microphysics_at_all_substeps = .FALSE.      !< namelist parameter
    222     LOGICAL ::  microphysics_ice_extension = .FALSE.              !< use ice microphysics scheme
     225    LOGICAL ::  microphysics_ice_phase = .FALSE.              !< use ice microphysics scheme
    223226    LOGICAL ::  microphysics_sat_adjust = .FALSE.                !< use saturation adjust bulk scheme?
    224227    LOGICAL ::  microphysics_kessler = .FALSE.                   !< use kessler bulk scheme?
     
    325328           microphysics_sat_adjust, &
    326329           microphysics_seifert, &
    327            microphysics_ice_extension, &
     330           microphysics_ice_phase, &
    328331           na_init, &
    329332           nc_const, &
     
    454457          ventilation_effect, &
    455458          ice_crystal_sedimentation, &
    456           microphysics_ice_extension, &
     459          microphysics_ice_phase, &
    457460          start_ice_microphysics, &
    458461          in_init
     
    582585       surf_microphysics_morrison = microphysics_morrison
    583586       surf_microphysics_seifert = microphysics_seifert
    584        surf_microphysics_ice_extension = microphysics_ice_extension
     587       surf_microphysics_ice_phase = microphysics_ice_phase
    585588!
    586589!--    Check aerosol
     
    629632
    630633          CASE ( 'ni' )
    631              IF ( .NOT.  microphysics_ice_extension )  THEN
     634             IF ( .NOT.  microphysics_ice_phase )  THEN
    632635                message_string = 'output of "' // TRIM( var ) // '" ' //       &
    633636                                 'requires ' //                                &
    634                                  'microphysics_ice_extension = ".TRUE."'
     637                                 'microphysics_ice_phase = ".TRUE."'
    635638                CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    636639             ENDIF
     
    659662
    660663          CASE ( 'qi' )
    661              IF ( .NOT.  microphysics_ice_extension ) THEN
     664             IF ( .NOT.  microphysics_ice_phase ) THEN
    662665                message_string = 'output of "' // TRIM( var ) // '" ' //       &
    663666                                 'requires ' //                                &
    664                                  'microphysics_ice_extension = ".TRUE."'
     667                                 'microphysics_ice_phase = ".TRUE."'
    665668                CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    666669             ENDIF
     
    758761
    759762          CASE ( 'ni' )
    760              IF ( .NOT.  microphysics_ice_extension )  THEN
     763             IF ( .NOT.  microphysics_ice_phase )  THEN
    761764                message_string = 'data_output_pr = ' //                        &
    762765                                 TRIM( data_output_pr(var_count) ) //          &
    763766                                 ' is not implemented for' //                  &
    764                                  ' microphysics_ice_extension = ".F."'
     767                                 ' microphysics_ice_phase = ".F."'
    765768                CALL message( 'check_parameters', 'PA0358', 1, 2, 0, 6, 0 )
    766769             ENDIF
     
    807810
    808811          CASE ( 'qi' )
    809              IF ( .NOT.  microphysics_ice_extension )  THEN
     812             IF ( .NOT.  microphysics_ice_phase )  THEN
    810813                message_string = 'data_output_pr = ' //                        &
    811814                                 TRIM( data_output_pr(var_count) ) //          &
    812815                                 ' is not implemented for' //                  &
    813                                  ' microphysics_ice_extension = ".F."'
     816                                 ' microphysics_ice_phase = ".F."'
    814817                CALL message( 'check_parameters', 'PA0358', 1, 2, 0, 6, 0 )
    815818             ENDIF
     
    894897       ENDIF
    895898
    896        IF ( microphysics_ice_extension )  THEN
     899       IF ( microphysics_ice_phase )  THEN
    897900!
    898901!--       3D-cloud drop water content, cloud drop concentration arrays
     
    918921             sums_wsnrs_ws_l = 0.0_wp
    919922          ENDIF
    920           IF ( microphysics_ice_extension )  THEN
     923          IF ( microphysics_ice_phase )  THEN
    921924             ALLOCATE( sums_wsqis_ws_l(nzb:nzt+1,0:threads_per_task-1) )
    922925             ALLOCATE( sums_wsnis_ws_l(nzb:nzt+1,0:threads_per_task-1) )
     
    953956                          diss_l_nr(nzb+1:nzt,nys:nyn,0:threads_per_task-1) )
    954957             ENDIF
    955              IF ( microphysics_ice_extension )  THEN
     958             IF ( microphysics_ice_phase )  THEN
    956959                ALLOCATE( flux_s_qi(nzb+1:nzt,0:threads_per_task-1),           &
    957960                          diss_s_qi(nzb+1:nzt,0:threads_per_task-1),           &
     
    980983          nr => nr_1;  nr_p  => nr_2;  tnr_m  => nr_3
    981984       ENDIF
    982        IF ( microphysics_ice_extension )  THEN
     985       IF ( microphysics_ice_phase )  THEN
    983986          qi => qi_1;  qi_p  => qi_2;  tqi_m  => qi_3
    984987          ni => ni_1;  ni_p  => ni_2;  tni_m  => ni_3
     
    10271030!
    10281031!--          Initialize the remaining quantities
    1029              IF ( microphysics_ice_extension )  THEN
     1032             IF ( microphysics_ice_phase )  THEN
    10301033                DO  i = nxlg, nxrg
    10311034                   DO  j = nysg, nyng
     
    10561059                nr_p  = nr
    10571060             ENDIF
    1058              IF ( microphysics_ice_extension )  THEN
     1061             IF ( microphysics_ice_phase )  THEN
    10591062                tqi_m = 0.0_wp
    10601063                tni_m = 0.0_wp
     
    12031206                sums_wsnrs_ws_l = 0.0_wp
    12041207             ENDIF
    1205              IF ( microphysics_ice_extension )  THEN
     1208             IF ( microphysics_ice_phase )  THEN
    12061209                sums_wsqis_ws_l = 0.0_wp
    12071210                sums_wsnis_ws_l = 0.0_wp
     
    12471250                sums_wsnrs_ws_l = 0.0_wp
    12481251             ENDIF
    1249              IF ( microphysics_ice_extension )  THEN
     1252             IF ( microphysics_ice_phase )  THEN
    12501253                sums_wsqis_ws_l = 0.0_wp
    12511254                sums_wsnis_ws_l = 0.0_wp
     
    13191322          ELSEIF ( microphysics_seifert  .AND.  .NOT. microphysics_morrison )  THEN
    13201323             CALL adjust_cloud
    1321              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1324             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13221325             CALL ice_nucleation
    1323              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1326             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13241327             CALL ice_deposition
    13251328             CALL autoconversion
     
    13291332             CALL sedimentation_rain
    13301333             IF ( cloud_water_sedimentation )  CALL sedimentation_cloud
    1331              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1334             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13321335             CALL adjust_ice
    1333              IF ( ice_crystal_sedimentation .AND. microphysics_ice_extension &
     1336             IF ( ice_crystal_sedimentation .AND. microphysics_ice_phase &
    13341337                  .AND. simulated_time > start_ice_microphysics ) CALL sedimentation_ice
    13351338
     
    13401343             CALL activation
    13411344             CALL condensation
    1342              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1345             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13431346             CALL adjust_ice
    1344              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1347             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13451348             CALL ice_nucleation
    1346              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1349             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13471350             CALL ice_deposition
    13481351             IF ( cloud_water_sedimentation )  CALL sedimentation_cloud
     
    13551358             CALL activation
    13561359             CALL condensation
    1357              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1360             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13581361             CALL adjust_ice
    1359              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1362             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13601363             CALL ice_nucleation
    1361              IF ( microphysics_ice_extension  .AND.  simulated_time > start_ice_microphysics )     &
     1364             IF ( microphysics_ice_phase  .AND.  simulated_time > start_ice_microphysics )     &
    13621365             CALL ice_deposition
    13631366             CALL autoconversion
     
    14341437          ELSEIF ( microphysics_seifert  .AND.  .NOT. microphysics_morrison )  THEN
    14351438             CALL adjust_cloud_ij( i,j )
    1436              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1439             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14371440                 CALL ice_nucleation_ij( i,j )
    1438              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1441             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14391442                 CALL ice_deposition_ij( i,j )
    14401443             CALL autoconversion_ij( i,j )
     
    14441447             CALL sedimentation_rain_ij( i,j )
    14451448             IF ( cloud_water_sedimentation )  CALL sedimentation_cloud_ij( i,j )
    1446              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1449             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14471450                 CALL adjust_ice_ij ( i,j )
    1448              IF ( ice_crystal_sedimentation .AND. microphysics_ice_extension &
     1451             IF ( ice_crystal_sedimentation .AND. microphysics_ice_phase &
    14491452                  .AND. simulated_time > start_ice_microphysics )  CALL sedimentation_ice_ij ( i,j )
    14501453!
     
    14541457             CALL activation_ij( i,j )
    14551458             CALL condensation_ij( i,j )
    1456              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1459             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14571460                 CALL adjust_ice_ij ( i,j )
    1458              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1461             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14591462                 CALL ice_nucleation_ij( i,j )
    1460              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1463             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14611464                 CALL ice_deposition_ij( i,j )
    14621465             IF ( cloud_water_sedimentation )  CALL sedimentation_cloud_ij( i,j )
     
    14691472             CALL activation_ij( i,j )
    14701473             CALL condensation_ij( i,j )
    1471              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1474             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14721475                 CALL adjust_ice_ij ( i,j )
    1473              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1476             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14741477                 CALL ice_nucleation_ij( i,j )
    1475              IF ( microphysics_ice_extension .AND. simulated_time > start_ice_microphysics )       &
     1478             IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics )       &
    14761479                 CALL ice_deposition_ij( i,j )
    14771480             CALL autoconversion_ij( i,j )
     
    15141517             CALL exchange_horiz( nr, nbgp )
    15151518          ENDIF
    1516           IF ( microphysics_ice_extension ) THEN
     1519          IF ( microphysics_ice_phase ) THEN
    15171520             CALL exchange_horiz( qi, nbgp )
    15181521             CALL exchange_horiz( ni, nbgp )
     
    17361739!--    If required, calculate prognostic equations for ice crystal content
    17371740!--    and ice crystal concentration
    1738        IF ( microphysics_ice_extension )  THEN
     1741       IF ( microphysics_ice_phase )  THEN
    17391742
    17401743          CALL cpu_log( log_point(70), 'qi-equation', 'start' )
     
    22622265!--    If required, calculate prognostic equations for ice crystal mixing ratio
    22632266!--    and ice crystal concentration
    2264        IF ( microphysics_ice_extension )  THEN
     2267       IF ( microphysics_ice_phase )  THEN
    22652268!
    22662269!--       Calculate prognostic equation for ice crystal mixing ratio
     
    25372540                   nr => nr_1;    nr_p => nr_2
    25382541                ENDIF
    2539                 IF ( microphysics_ice_extension )  THEN
     2542                IF ( microphysics_ice_phase )  THEN
    25402543                   qi => qi_1;    qi_p => qi_2
    25412544                   ni => ni_1;    ni_p => ni_2
     
    25522555                   nr => nr_2;    nr_p => nr_1
    25532556                ENDIF
    2554                 IF ( microphysics_ice_extension )  THEN
     2557                IF ( microphysics_ice_phase )  THEN
    25552558                   qi => qi_2;    qi_p => qi_1
    25562559                   ni => ni_2;    ni_p => ni_1
     
    26012604       ENDIF
    26022605
    2603        IF ( microphysics_ice_extension )  THEN
     2606       IF ( microphysics_ice_phase )  THEN
    26042607!
    26052608!--       Surface conditions ice crysral (Dirichlet)
     
    26612664             nr_p(:,nys-1,:) = nr_p(:,nys,:)
    26622665          ENDIF
    2663           IF ( microphysics_ice_extension )  THEN
     2666          IF ( microphysics_ice_phase )  THEN
    26642667             qi_p(:,nys-1,:) = qi_p(:,nys,:)
    26652668             ni_p(:,nys-1,:) = ni_p(:,nys,:)
     
    26742677             nr_p(:,nyn+1,:) = nr_p(:,nyn,:)
    26752678          ENDIF
    2676           IF ( microphysics_ice_extension )  THEN
     2679          IF ( microphysics_ice_phase )  THEN
    26772680             qi_p(:,nyn+1,:) = qi_p(:,nyn,:)
    26782681             ni_p(:,nyn+1,:) = ni_p(:,nyn,:)
     
    26872690             nr_p(:,:,nxl-1) = nr_p(:,:,nxl)
    26882691          ENDIF
    2689           IF ( microphysics_ice_extension )  THEN
     2692          IF ( microphysics_ice_phase )  THEN
    26902693             qi_p(:,:,nxl-1) = qi_p(:,:,nxl)
    26912694             ni_p(:,:,nxl-1) = ni_p(:,:,nxl)
     
    27002703             nr_p(:,:,nxr+1) = nr_p(:,:,nxr)
    27012704          ENDIF
    2702           IF ( microphysics_ice_extension )  THEN
     2705          IF ( microphysics_ice_phase )  THEN
    27032706             qi_p(:,:,nxr+1) = qi_p(:,:,nxr)
    27042707             ni_p(:,:,nxr+1) = ni_p(:,:,nxr)
     
    33453348             READ ( 13 )  curvature_solution_effects_bulk
    33463349
    3347           CASE ( 'microphysics_ice_extension' )
    3348              READ ( 13 )  microphysics_ice_extension
     3350          CASE ( 'microphysics_ice_phase' )
     3351             READ ( 13 )  microphysics_ice_phase
    33493352
    33503353          CASE ( 'ice_crystal_sedimentation' )
     
    33893392       CALL rrd_mpi_io( 'curvature_solution_effects_bulk', curvature_solution_effects_bulk )
    33903393       CALL rrd_mpi_io( 'start_ice_microphysics', start_ice_microphysics )
    3391        CALL rrd_mpi_io( 'microphysics_ice_extension', microphysics_ice_extension )
     3394       CALL rrd_mpi_io( 'microphysics_ice_phase', microphysics_ice_phase )
    33923395       CALL rrd_mpi_io( 'in_init', in_init )
    33933396       CALL rrd_mpi_io( 'ice_crystal_sedimentation', ice_crystal_sedimentation )
     
    36423645       ENDIF
    36433646
    3644        IF ( microphysics_ice_extension )  THEN
     3647       IF ( microphysics_ice_phase )  THEN
    36453648
    36463649          CALL rrd_mpi_io( 'ni', ni )
     
    37243727          WRITE ( 14 )  start_ice_microphysics
    37253728
    3726           CALL wrd_write_string( 'microphysics_ice_extension' )
    3727           WRITE ( 14 )  microphysics_ice_extension
     3729          CALL wrd_write_string( 'microphysics_ice_phase' )
     3730          WRITE ( 14 )  microphysics_ice_phase
    37283731
    37293732          CALL wrd_write_string( 'in_init' )
     
    37503753          CALL wrd_mpi_io( 'curvature_solution_effects_bulk', curvature_solution_effects_bulk )
    37513754          CALL wrd_mpi_io( 'start_ice_microphysics', start_ice_microphysics )
    3752           CALL wrd_mpi_io( 'microphysics_ice_extension', microphysics_ice_extension )
     3755          CALL wrd_mpi_io( 'microphysics_ice_phase', microphysics_ice_phase )
    37533756          CALL wrd_mpi_io( 'in_init', in_init )
    37543757          CALL wrd_mpi_io( 'ice_crystal_sedimentation', ice_crystal_sedimentation )
     
    38153818          ENDIF
    38163819
    3817           IF ( microphysics_ice_extension )  THEN
     3820          IF ( microphysics_ice_phase )  THEN
    38183821
    38193822             CALL wrd_write_string( 'ni' )
     
    38773880             IF ( ALLOCATED( qr_av ) )  CALL wrd_mpi_io( 'qr_av', qr_av )
    38783881          ENDIF
    3879           IF ( microphysics_ice_extension )  THEN
     3882          IF ( microphysics_ice_phase )  THEN
    38803883             CALL wrd_mpi_io( 'ni', ni )
    38813884             IF ( ALLOCATED( ni_av ) )  CALL wrd_mpi_io( 'ni_av', ni_av )
     
    44524455!--             Actual temperature, t_l is calculated directly before
    44534456!--             in supersaturation
    4454                 IF ( microphysics_ice_extension ) THEN
     4457                IF ( microphysics_ice_phase ) THEN
    44554458                   temp = t_l + lv_d_cp * ql(k,j,i) + ls_d_cp * qi(k,j,i)
    44564459                ELSE
     
    45394542!--       Actual temperature, t_l is calculated directly before
    45404543!--       in supersaturation
    4541           IF ( microphysics_ice_extension ) THEN
     4544          IF ( microphysics_ice_phase ) THEN
    45424545             temp = t_l + lv_d_cp * ql(k,j,i) + ls_d_cp * qi(k,j,i)
    45434546          ELSE
     
    65516554       q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
    65526555
    6553        IF ( .NOT. microphysics_ice_extension ) THEN
     6556       IF ( .NOT. microphysics_ice_phase ) THEN
    65546557          IF ( microphysics_seifert ) THEN
    65556558             sat = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
     
    66336636!
    66346637!--             Compute the liquid water content
    6635                 IF ( .NOT. microphysics_ice_extension ) THEN
     6638                IF ( .NOT. microphysics_ice_phase ) THEN
    66366639                   IF ( microphysics_seifert  .AND. .NOT.                      &
    66376640                        microphysics_morrison )  THEN
  • palm/trunk/SOURCE/compute_vpt.f90

    r4502 r4521  
    2525! -----------------
    2626! $Id$
     27! Rename variable
     28!
     29! 4502 2020-04-17 16:14:16Z schwenkel
    2730! Implementation of ice microphysics
    2831!
     
    5962
    6063    USE bulk_cloud_model_mod,                                                  &
    61         ONLY:  bulk_cloud_model, microphysics_ice_extension
     64        ONLY:  bulk_cloud_model, microphysics_ice_phase
    6265
    6366    IMPLICIT NONE
     
    6770    IF ( .NOT. bulk_cloud_model  .AND.  .NOT. cloud_droplets )  THEN
    6871       vpt = pt * ( 1.0_wp + 0.61_wp * q )
    69     ELSEIF ( bulk_cloud_model  .AND.  .NOT. microphysics_ice_extension )  THEN
     72    ELSEIF ( bulk_cloud_model  .AND.  .NOT. microphysics_ice_phase )  THEN
    7073       DO  k = nzb, nzt+1
    7174              vpt(k,:,:) = ( pt(k,:,:) + d_exner(k) * lv_d_cp * ql(k,:,:) ) *  &
    7275                       ( 1.0_wp + 0.61_wp * q(k,:,:) - 1.61_wp *  ql(k,:,:)  )
    7376       ENDDO
    74     ELSEIF ( bulk_cloud_model  .AND.  microphysics_ice_extension ) THEN
     77    ELSEIF ( bulk_cloud_model  .AND.  microphysics_ice_phase ) THEN
    7578       DO  k = nzb, nzt+1
    7679          vpt(k,:,:) = ( pt(k,:,:) + d_exner(k) * lv_d_cp * ql(k,:,:)   +      &
  • palm/trunk/SOURCE/flow_statistics.f90

    r4502 r4521  
    2525! -----------------
    2626! $Id$
     27! Rename variable
     28!
     29! 4502 2020-04-17 16:14:16Z schwenkel
    2730! Implementation of ice microphysics
    2831!
     
    97100    USE bulk_cloud_model_mod,                                                  &
    98101        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert,   &
    99               microphysics_ice_extension
     102              microphysics_ice_phase
    100103
    101104    USE chem_modules,                                                          &
     
    14181421                                                                flag
    14191422                         ENDIF
    1420                          IF ( microphysics_ice_extension )  THEN
     1423                         IF ( microphysics_ice_phase )  THEN
    14211424                            sums_l(k,124,tn) = sums_l(k,124,tn) + ni(k,j,i) *  &
    14221425                                                                rmask(j,i,sr) *&
  • palm/trunk/SOURCE/init_masks.f90

    r4502 r4521  
    2525! -----------------
    2626! $Id$
     27! Rename variable
     28!
     29! 4502 2020-04-17 16:14:16Z schwenkel
    2730! Implementation of ice microphysics
    2831!
     
    6265    USE bulk_cloud_model_mod,                                                  &
    6366        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert,   &
    64               microphysics_ice_extension
     67              microphysics_ice_phase
    6568
    6669    USE control_parameters,                                                    &
     
    277280                        '" requires bulk_cloud_model = .TRUE.'
    278281                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
    279                  ELSEIF ( .NOT. microphysics_ice_extension ) THEN
     282                 ELSEIF ( .NOT. microphysics_ice_phase ) THEN
    280283                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
    281                          'requires  microphysics_ice_extension = .TRUE.'
     284                         'requires  microphysics_ice_phase = .TRUE.'
    282285                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    283286                ENDIF
     
    355358                            'requires bulk_cloud_model = .TRUE.'
    356359                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
    357                 ELSEIF ( .NOT. microphysics_ice_extension ) THEN
     360                ELSEIF ( .NOT. microphysics_ice_phase ) THEN
    358361                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
    359                             'requires microphysics_ice_extension = .TRUE.'
     362                            'requires microphysics_ice_phase = .TRUE.'
    360363                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    361364                ENDIF
  • palm/trunk/SOURCE/surface_mod.f90

    r4517 r4521  
    2626! -----------------
    2727! $Id$
     28! Rename variable
     29!
     30! 4517 2020-05-03 14:29:30Z raasch
    2831! added restart with MPI-IO for reading local arrays
    2932!
     
    479482    LOGICAL ::  surf_microphysics_morrison = .FALSE.   !< use 2-moment Morrison (add. prog. eq. for nc and qc)
    480483    LOGICAL ::  surf_microphysics_seifert = .FALSE.    !< use 2-moment Seifert and Beheng scheme
    481     LOGICAL ::  surf_microphysics_ice_extension = .FALSE. !< use 2-moment Seifert and Beheng scheme
     484    LOGICAL ::  surf_microphysics_ice_phase = .FALSE. !< use 2-moment Seifert and Beheng scheme
    482485
    483486
     
    535538           surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, surf_usm_v, surf_type,      &
    536539           vertical_surfaces_exist, surf_bulk_cloud_model, surf_microphysics_morrison,             &
    537            surf_microphysics_seifert, surf_microphysics_ice_extension
     540           surf_microphysics_seifert, surf_microphysics_ice_phase
    538541!
    539542!-- Public subroutines and functions
     
    13151318!
    13161319!--
    1317        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     1320       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    13181321          DEALLOCATE ( surfaces%qis )
    13191322          DEALLOCATE ( surfaces%nis )
     
    14521455!
    14531456!--
    1454        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     1457       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    14551458          ALLOCATE ( surfaces%qis(1:surfaces%ns)   )
    14561459          ALLOCATE ( surfaces%nis(1:surfaces%ns)   )
     
    16111614!
    16121615!--
    1613        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     1616       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    16141617          DEALLOCATE ( surfaces%qisws )
    16151618          DEALLOCATE ( surfaces%nisws )
     
    16901693!
    16911694!--
    1692        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     1695       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    16931696          ALLOCATE ( surfaces%qisws(1:surfaces%ns) )
    16941697          ALLOCATE ( surfaces%nisws(1:surfaces%ns) )
     
    18751878       ENDIF
    18761879
    1877        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     1880       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    18781881          DEALLOCATE ( surfaces%qis )
    18791882          DEALLOCATE ( surfaces%nis )
     
    20082011       ENDIF
    20092012
    2010        IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     2013       IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    20112014          ALLOCATE ( surfaces%qis(1:surfaces%ns)   )
    20122015          ALLOCATE ( surfaces%nis(1:surfaces%ns)   )
     
    26572660                ENDIF
    26582661
    2659                 IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     2662                IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    26602663                   surf%qis(num_h) = 0.0_wp
    26612664                   surf%nis(num_h) = 0.0_wp
     
    28382841                   surf%qrsws(num_h) = 0.0_wp
    28392842                ENDIF
    2840                 IF ( surf_bulk_cloud_model  .AND.  surf_microphysics_ice_extension ) THEN
     2843                IF ( surf_bulk_cloud_model  .AND.  surf_microphysics_ice_phase ) THEN
    28412844                   surf%nisws(num_h) = 0.0_wp
    28422845                   surf%qisws(num_h) = 0.0_wp
     
    29822985                   surf%nrsws(num_v) = 0.0_wp
    29832986                ENDIF
    2984                 IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_extension)  THEN
     2987                IF ( surf_bulk_cloud_model .AND. surf_microphysics_ice_phase)  THEN
    29852988                   surf%qis(num_v) = 0.0_wp
    29862989                   surf%nis(num_v) = 0.0_wp
  • palm/trunk/SOURCE/time_integration.f90

    r4511 r4521  
    2525! -----------------
    2626! $Id$
     27! Rename variable
     28!
     29! 4511 2020-04-30 12:20:40Z raasch
    2730! chemistry decycling replaced by explicit setting of lateral boundary conditions
    2831!
     
    233236    USE bulk_cloud_model_mod,                                                                      &
    234237        ONLY: bulk_cloud_model, calc_liquid_water_content, collision_turbulence,                   &
    235               microphysics_ice_extension, microphysics_morrison, microphysics_seifert
     238              microphysics_ice_phase, microphysics_morrison, microphysics_seifert
    236239
    237240    USE calc_mean_profile_mod,                                                                     &
     
    808811                CALL exchange_horiz( nr_p, nbgp )
    809812             ENDIF
    810              IF ( bulk_cloud_model .AND. microphysics_ice_extension )  THEN
     813             IF ( bulk_cloud_model .AND. microphysics_ice_phase )  THEN
    811814                CALL exchange_horiz( qi_p, nbgp )
    812815                CALL exchange_horiz( ni_p, nbgp )
     
    907910                       CALL exchange_horiz( nr, nbgp )
    908911                   ENDIF
    909                    IF ( bulk_cloud_model  .AND.  microphysics_ice_extension )  THEN
     912                   IF ( bulk_cloud_model  .AND.  microphysics_ice_phase )  THEN
    910913                      CALL exchange_horiz( qi, nbgp )
    911914                      CALL exchange_horiz( ni, nbgp )
Note: See TracChangeset for help on using the changeset viewer.