Ignore:
Timestamp:
Sep 24, 2018 3:42:55 PM (3 years ago)
Author:
knoop
Message:

Modularization of all bulk cloud physics code components

File:
1 edited

Legend:

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

    r3245 r3274  
    2525! -----------------
    2626! $Id$
     27! Modularization of all bulk cloud physics code components
     28!
     29! 3245 2018-09-13 14:08:16Z knoop
    2730! unused variables removed, shortest_distance has wp now
    2831!
     
    154157#endif
    155158
     159    USE basic_constants_and_equations_mod,                                     &
     160        ONLY:  g, kappa, lv_d_cp, lv_d_rd
     161
    156162    USE control_parameters,                                                    &
    157163        ONLY:  constant_diffusion, dt_3d, e_init, humidity,                    &
    158164               initializing_actions, intermediate_timestep_count,              &
    159                intermediate_timestep_count_max, kappa, km_constant,            &
     165               intermediate_timestep_count_max, km_constant,                   &
    160166               les_dynamic, les_mw, ocean, plant_canopy, prandtl_number,       &
    161167               prho_reference, pt_reference, rans_mode, rans_tke_e, rans_tke_l,&
     
    950956 SUBROUTINE tcm_init_arrays
    951957
    952     USE microphysics_mod,                                                      &
     958    USE bulk_cloud_model_mod,                                                  &
    953959        ONLY:  collision_turbulence
    954960
     
    25982604
    25992605    USE arrays_3d,                                                             &
    2600         ONLY:  ddzw, dd2zu, drho_air_zw, q, ql
    2601 
    2602     USE cloud_parameters,                                                      &
    2603         ONLY:  l_d_cp, l_d_r, pt_d_t, t_d_pt
     2606        ONLY:  ddzw, dd2zu, drho_air_zw, q, ql, d_exner, exner
    26042607
    26052608    USE control_parameters,                                                    &
    2606         ONLY:  cloud_droplets, cloud_physics, constant_flux_layer, g, neutral, &
     2609        ONLY:  cloud_droplets, constant_flux_layer, neutral,                  &
    26072610               rho_reference, use_single_reference_value, use_surface_fluxes,  &
    26082611               use_top_fluxes
     
    26102613    USE grid_variables,                                                        &
    26112614        ONLY:  ddx, dx, ddy, dy
     2615
     2616    USE bulk_cloud_model_mod,                                                  &
     2617        ONLY:  bulk_cloud_model
    26122618
    26132619    USE surface_mod,                                                           &
     
    30673073!--                Flag 9 is used to mask top fluxes, flag 30 to mask
    30683074!--                surface fluxes
    3069                    IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3075                   IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    30703076                      k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    30713077                      k2 = 0.61_wp * pt(k,j,i)
     
    30833089                                          BTEST( wall_flags_0(k,j,i), 9 )      &
    30843090                                        )
    3085                    ELSE IF ( cloud_physics )  THEN
     3091                   ELSE IF ( bulk_cloud_model )  THEN
    30863092                      IF ( ql(k,j,i) == 0.0_wp )  THEN
    30873093                         k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    30883094                         k2 = 0.61_wp * pt(k,j,i)
    30893095                      ELSE
    3090                          theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3091                          temp  = theta * t_d_pt(k)
     3096                         theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3097                         temp  = theta * exner(k)
    30923098                         k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                  &
    30933099                                       ( q(k,j,i) - ql(k,j,i) ) *              &
    3094                               ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /         &
    3095                               ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *           &
     3100                              ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
     3101                              ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
    30963102                              ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3097                          k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3103                         k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    30983104                      ENDIF
    30993105                      tend(k,j,i) = tend(k,j,i) - kh(k,j,i) *                  &
     
    31443150                         k = surf_def_h(l)%k(m)
    31453151
    3146                          IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3152                         IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    31473153                            k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    31483154                            k2 = 0.61_wp * pt(k,j,i)
    3149                          ELSE IF ( cloud_physics )  THEN
     3155                         ELSE IF ( bulk_cloud_model )  THEN
    31503156                            IF ( ql(k,j,i) == 0.0_wp )  THEN
    31513157                               k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    31523158                               k2 = 0.61_wp * pt(k,j,i)
    31533159                            ELSE
    3154                                theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3155                                temp  = theta * t_d_pt(k)
     3160                               theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3161                               temp  = theta * exner(k)
    31563162                               k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *            &
    31573163                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3158                                  ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /      &
    3159                                  ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *        &
     3164                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
     3165                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
    31603166                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3161                                k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3167                               k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    31623168                            ENDIF
    31633169                         ELSE IF ( cloud_droplets )  THEN
     
    31813187                      k = surf_lsm_h%k(m)
    31823188
    3183                       IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3189                      IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    31843190                         k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    31853191                         k2 = 0.61_wp * pt(k,j,i)
    3186                       ELSE IF ( cloud_physics )  THEN
     3192                      ELSE IF ( bulk_cloud_model )  THEN
    31873193                         IF ( ql(k,j,i) == 0.0_wp )  THEN
    31883194                            k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    31893195                            k2 = 0.61_wp * pt(k,j,i)
    31903196                         ELSE
    3191                             theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3192                             temp  = theta * t_d_pt(k)
     3197                            theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3198                            temp  = theta * exner(k)
    31933199                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    31943200                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3195                                  ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /      &
    3196                                  ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *        &
     3201                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
     3202                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
    31973203                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3198                             k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3204                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    31993205                         ENDIF
    32003206                      ELSE IF ( cloud_droplets )  THEN
     
    32173223                      k = surf_usm_h%k(m)
    32183224
    3219                       IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3225                      IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    32203226                         k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    32213227                         k2 = 0.61_wp * pt(k,j,i)
    3222                       ELSE IF ( cloud_physics )  THEN
     3228                      ELSE IF ( bulk_cloud_model )  THEN
    32233229                         IF ( ql(k,j,i) == 0.0_wp )  THEN
    32243230                            k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    32253231                            k2 = 0.61_wp * pt(k,j,i)
    32263232                         ELSE
    3227                             theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3228                             temp  = theta * t_d_pt(k)
     3233                            theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3234                            temp  = theta * exner(k)
    32293235                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    32303236                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3231                                  ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /      &
    3232                                  ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *        &
     3237                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
     3238                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
    32333239                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3234                             k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3240                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    32353241                         ENDIF
    32363242                      ELSE IF ( cloud_droplets )  THEN
     
    32603266                      k = surf_def_h(2)%k(m)
    32613267
    3262                       IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3268                      IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    32633269                         k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    32643270                         k2 = 0.61_wp * pt(k,j,i)
    3265                       ELSE IF ( cloud_physics )  THEN
     3271                      ELSE IF ( bulk_cloud_model )  THEN
    32663272                         IF ( ql(k,j,i) == 0.0_wp )  THEN
    32673273                            k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    32683274                            k2 = 0.61_wp * pt(k,j,i)
    32693275                         ELSE
    3270                             theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3271                             temp  = theta * t_d_pt(k)
     3276                            theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3277                            temp  = theta * exner(k)
    32723278                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    32733279                                       ( q(k,j,i) - ql(k,j,i) ) *              &
    3274                               ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /         &
    3275                               ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *           &
     3280                              ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
     3281                              ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
    32763282                              ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3277                             k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3283                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    32783284                         ENDIF
    32793285                      ELSE IF ( cloud_droplets )  THEN
     
    33163322
    33173323    USE arrays_3d,                                                             &
    3318         ONLY:  ddzw, dd2zu, drho_air_zw, q, ql
    3319 
    3320     USE cloud_parameters,                                                      &
    3321         ONLY:  l_d_cp, l_d_r, pt_d_t, t_d_pt
     3324        ONLY:  ddzw, dd2zu, drho_air_zw, q, ql, d_exner, exner
    33223325
    33233326    USE control_parameters,                                                    &
    3324         ONLY:  cloud_droplets, cloud_physics, constant_flux_layer, g, neutral, &
     3327        ONLY:  cloud_droplets, constant_flux_layer, neutral,                &
    33253328               rho_reference, use_single_reference_value, use_surface_fluxes,  &
    33263329               use_top_fluxes
     
    33283331    USE grid_variables,                                                        &
    33293332        ONLY:  ddx, dx, ddy, dy
     3333
     3334    USE bulk_cloud_model_mod,                                                  &
     3335        ONLY:  bulk_cloud_model
    33303336
    33313337    USE surface_mod,                                                           &
     
    38643870!
    38653871!--          Flag 9 is used to mask top fluxes, flag 30 to mask surface fluxes
    3866              IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets )  THEN
     3872             IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets )  THEN
    38673873                k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    38683874                k2 = 0.61_wp * pt(k,j,i)
     
    38803886                                        )
    38813887
    3882              ELSE IF ( cloud_physics )  THEN
     3888             ELSE IF ( bulk_cloud_model )  THEN
    38833889                IF ( ql(k,j,i) == 0.0_wp )  THEN
    38843890                   k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    38853891                   k2 = 0.61_wp * pt(k,j,i)
    38863892                ELSE
    3887                    theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3888                    temp  = theta * t_d_pt(k)
     3893                   theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3894                   temp  = theta * exner(k)
    38893895                   k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                        &
    38903896                                 ( q(k,j,i) - ql(k,j,i) ) *                    &
    3891                         ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /               &
    3892                         ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *                 &
     3897                        ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
     3898                        ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
    38933899                        ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3894                    k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3900                   k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    38953901                ENDIF
    38963902                tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g /                    &
     
    39343940                   k = surf_def_h(l)%k(m)
    39353941
    3936                    IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3942                   IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    39373943                      k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    39383944                      k2 = 0.61_wp * pt(k,j,i)
    3939                    ELSE IF ( cloud_physics )  THEN
     3945                   ELSE IF ( bulk_cloud_model )  THEN
    39403946                      IF ( ql(k,j,i) == 0.0_wp )  THEN
    39413947                         k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    39423948                         k2 = 0.61_wp * pt(k,j,i)
    39433949                      ELSE
    3944                         theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3945                         temp  = theta * t_d_pt(k)
     3950                        theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3951                        temp  = theta * exner(k)
    39463952                        k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                   &
    39473953                                   ( q(k,j,i) - ql(k,j,i) ) *                  &
    3948                           ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /             &
    3949                           ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *               &
     3954                          ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /           &
     3955                          ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *            &
    39503956                          ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3951                         k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3957                        k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    39523958                      ENDIF
    39533959                   ELSE IF ( cloud_droplets )  THEN
     
    39713977                k = surf_lsm_h%k(m)
    39723978
    3973                 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     3979                IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    39743980                    k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    39753981                    k2 = 0.61_wp * pt(k,j,i)
    3976                 ELSE IF ( cloud_physics )  THEN
     3982                ELSE IF ( bulk_cloud_model )  THEN
    39773983                    IF ( ql(k,j,i) == 0.0_wp )  THEN
    39783984                       k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    39793985                       k2 = 0.61_wp * pt(k,j,i)
    39803986                    ELSE
    3981                        theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    3982                        temp  = theta * t_d_pt(k)
     3987                       theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     3988                       temp  = theta * exner(k)
    39833989                       k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
    39843990                                     ( q(k,j,i) - ql(k,j,i) ) *                &
    3985                             ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /           &
    3986                             ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *             &
     3991                            ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
     3992                            ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
    39873993                            ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    3988                        k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     3994                       k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    39893995                   ENDIF
    39903996                ELSE IF ( cloud_droplets )  THEN
     
    40074013                k = surf_usm_h%k(m)
    40084014
    4009                 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
     4015                IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets ) THEN
    40104016                    k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    40114017                    k2 = 0.61_wp * pt(k,j,i)
    4012                 ELSE IF ( cloud_physics )  THEN
     4018                ELSE IF ( bulk_cloud_model )  THEN
    40134019                    IF ( ql(k,j,i) == 0.0_wp )  THEN
    40144020                       k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    40154021                       k2 = 0.61_wp * pt(k,j,i)
    40164022                    ELSE
    4017                        theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    4018                        temp  = theta * t_d_pt(k)
     4023                       theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     4024                       temp  = theta * exner(k)
    40194025                       k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
    40204026                                     ( q(k,j,i) - ql(k,j,i) ) *                &
    4021                             ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /           &
    4022                             ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *             &
     4027                            ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
     4028                            ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
    40234029                            ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    4024                        k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     4030                       k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    40254031                   ENDIF
    40264032                ELSE IF ( cloud_droplets )  THEN
     
    40474053
    40484054
    4049                 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets )  THEN
     4055                IF ( .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets )  THEN
    40504056                   k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    40514057                   k2 = 0.61_wp * pt(k,j,i)
    4052                 ELSE IF ( cloud_physics )  THEN
     4058                ELSE IF ( bulk_cloud_model )  THEN
    40534059                   IF ( ql(k,j,i) == 0.0_wp )  THEN
    40544060                      k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    40554061                      k2 = 0.61_wp * pt(k,j,i)
    40564062                   ELSE
    4057                       theta = pt(k,j,i) + pt_d_t(k) * l_d_cp * ql(k,j,i)
    4058                       temp  = theta * t_d_pt(k)
     4063                      theta = pt(k,j,i) + d_exner(k) * lv_d_cp * ql(k,j,i)
     4064                      temp  = theta * exner(k)
    40594065                      k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                     &
    40604066                                 ( q(k,j,i) - ql(k,j,i) ) *                    &
    4061                         ( 1.0_wp + 0.622_wp * l_d_r / temp ) ) /               &
    4062                         ( 1.0_wp + 0.622_wp * l_d_r * l_d_cp *                 &
     4067                        ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
     4068                        ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
    40634069                        ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    4064                       k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
     4070                      k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
    40654071                   ENDIF
    40664072                ELSE IF ( cloud_droplets )  THEN
     
    41004106        ONLY:  ddx2, ddy2
    41014107
    4102     USE microphysics_mod,                                                      &
     4108    USE bulk_cloud_model_mod,                                                  &
    41034109        ONLY:  collision_turbulence
    41044110
     
    42374243        ONLY:  ddx2, ddy2
    42384244       
    4239     USE microphysics_mod,                                                      &
     4245    USE bulk_cloud_model_mod,                                                  &
    42404246        ONLY:  collision_turbulence
    42414247
     
    44764482
    44774483    USE control_parameters,                                                    &
    4478         ONLY:  atmos_ocean_sign, g, use_single_reference_value,                &
     4484        ONLY:  atmos_ocean_sign, use_single_reference_value,                   &
    44794485               wall_adjustment, wall_adjustment_factor
    44804486
     
    45334539
    45344540    USE control_parameters,                                                    &
    4535         ONLY:  atmos_ocean_sign, g, use_single_reference_value
     4541        ONLY:  atmos_ocean_sign, use_single_reference_value
    45364542
    45374543    IMPLICIT NONE
Note: See TracChangeset for help on using the changeset viewer.