Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

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

    r1958 r1960  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Separate humidity and passive scalar
     22! +bc_s_t_val, diss_s_s, diss_l_s, flux_s_s, flux_l_s, s, sp, s1, s2, s3, ssws_av,
     23!  s_init, s_surf, sums_wsss_ws_l, ss, ssws, sswst, ts_m, wall_sflux
     24! +constant_scalarflux, ibc_s_b, ibc_s_t, s_vertical_gradient_level_ind
     25!
     26! Unused variables removed
     27! -gamma_x, gamma_y, gamma_z, var_x, var_y, var_z
     28!
     29! Change default values (in order to unify gradient calculation):
     30! pt_vertical_gradient_level, sa_vertical_gradient_level
    2231!
    2332! Former revisions:
     
    411420          dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, l_grid,           &
    412421          ptdf_x, ptdf_y, p_surf, pt_surf, pt_init, qsws_surf, q_init, q_surf, &
    413           rdf, rdf_sc, ref_state, sa_init, shf_surf, timenudge, time_surf,     &
    414           time_vert, tmp_tnudge, ug, u_init, u_nzb_p1_for_vfc, vg, v_init,     &
    415           v_nzb_p1_for_vfc, w_subs, zu, zw
     422          rdf, rdf_sc, ref_state, s_init, s_surf, sa_init, shf_surf,           &
     423          timenudge, time_surf, time_vert, tmp_tnudge, ug, u_init,             &
     424          u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
    416425
    417426    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::                                   &
    418427          c_u, c_v, c_w, diss_s_e, diss_s_nr, diss_s_pt, diss_s_q,             &
    419           diss_s_qr, diss_s_sa, diss_s_u, diss_s_v, diss_s_w, dzu_mg, dzw_mg,  &
    420           flux_s_e, flux_s_nr, flux_s_pt, flux_s_q, flux_s_qr, flux_s_sa,      &
    421           flux_s_u, flux_s_v, flux_s_w, f1_mg, f2_mg, f3_mg,                   &
    422           mean_inflow_profiles, nrs, nrsws, nrswst,                            &
     428          diss_s_qr, diss_s_s, diss_s_sa, diss_s_u, diss_s_v, diss_s_w, dzu_mg,&
     429          dzw_mg, flux_s_e, flux_s_nr, flux_s_pt, flux_s_q, flux_s_qr,         &
     430          flux_s_s, flux_s_sa, flux_s_u, flux_s_v, flux_s_w, f1_mg, f2_mg,     &
     431          f3_mg, mean_inflow_profiles, nrs, nrsws, nrswst,                     &
    423432          ol,                                                                  & !< Obukhov length
    424433          precipitation_amount, precipitation_rate, ptnudge, pt_slope_ref,     &
    425434          qnudge, qs, qsws, qswst, qswst_remote, qrs, qrsws, qrswst,           &
    426           saswsb, saswst, shf, tnudge, td_lsa_lpt, td_lsa_q, td_sub_lpt,       &
     435          saswsb, saswst, shf, ss, ssws, sswst, tnudge, td_lsa_lpt, td_lsa_q,  &
     436          td_sub_lpt,                                                          &
    427437          td_sub_q, total_2d_a, total_2d_o, ts, tswst, ug_vert, unudge, us,    &
    428438          usws, uswst, vnudge, vg_vert, vsws, vswst, wnudge, wsubs_vert,       &
     
    433443    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::                                 &
    434444          d, de_dx, de_dy, de_dz, diss, diss_l_e,                              &
    435           diss_l_nr, diss_l_pt, diss_l_q, diss_l_qr, diss_l_sa, diss_l_u,      &
    436           diss_l_v, diss_l_w, flux_l_e, flux_l_nr, flux_l_pt, flux_l_q,        &
    437           flux_l_qr, flux_l_sa, flux_l_u, flux_l_v, flux_l_w, kh, km,          &
    438           l_wall, prr, p_loc, tend, tric,                                      &
     445          diss_l_nr, diss_l_pt, diss_l_q, diss_l_qr, diss_l_s, diss_l_sa,      &
     446          diss_l_u, diss_l_v, diss_l_w, flux_l_e, flux_l_nr, flux_l_pt,        &
     447          flux_l_q, flux_l_qr, flux_l_s, flux_l_sa, flux_l_u, flux_l_v,        &
     448          flux_l_w, kh, km, l_wall, prr, p_loc, tend, tric,                    &
    439449          u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n, v_m_r, v_m_s, w_m_l,       &
    440450          w_m_n, w_m_r, w_m_s
     
    443453    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                         &
    444454          e, e_p, nr, nr_p, p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, ql_v,     &
    445           ql_vp, qr, qr_p, rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m,     &
    446           tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
     455          ql_vp, qr, qr_p, rho, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m,    &
     456          tqr_m, ts_m, tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
    447457#else
    448458    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                         &
    449459          e_1, e_2, e_3, p, prho_1, nr_1, nr_2, nr_3, pt_1, pt_2, pt_3, q_1,   &
    450460          q_2, q_3, qc_1, ql_v, ql_vp, ql_1, ql_2, qr_1, qr_2, qr_3, rho_1,    &
    451           sa_1, sa_2, sa_3, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1, w_1, w_2, w_3
     461          s_1, s_2, s_3, sa_1, sa_2, sa_3, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1,&
     462          w_1, w_2, w_3
    452463
    453464    REAL(wp), DIMENSION(:,:,:), POINTER ::                                     &
    454465          e, e_p, nr, nr_p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, qr, qr_p,    &
    455           rho, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m, tu_m, tv_m,   &
    456           tw_m, u, u_p, v, v_p, vpt, w, w_p
     466          rho, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, ts_m,        &
     467          tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
    457468#endif
    458469
    459470    REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::  rif_wall, tri
    460 
    461     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: var_x, var_y, var_z, gamma_x,   &
    462                                            gamma_y, gamma_z
    463471
    464472
     
    481489                                              ol_av,                 & !< Avg. of Obukhov length
    482490                                              qsws_av,               & !< Avg. of surface moisture flux
     491                                              ssws_av,               & !< Avg. of surface scalar flux
    483492                                              shf_av,                & !< Avg. of surface heat flux
    484493                                              ts_av,                 & !< Avg. of characteristic temperature scale
     
    631640                     dz_stretch_level_index, ensemble_member_nr = 0, gamma_mg, gathered_size, &
    632641                     grid_level, ibc_e_b, ibc_p_b, ibc_p_t, &
    633                      ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, &
     642                     ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b, ibc_s_t, &
    634643                     ibc_sa_t, ibc_uv_b, ibc_uv_t, &
    635644                     inflow_disturbance_begin = -1, inflow_disturbance_end = -1, &
     
    657666                     pt_vertical_gradient_level_ind(10) = -9999, &
    658667                     q_vertical_gradient_level_ind(10) = -9999, &
     668                     s_vertical_gradient_level_ind(10) = -9999, &                     
    659669                     sa_vertical_gradient_level_ind(10) = -9999, &
    660670                     section(100,3), section_xy(100) = -9999, &
     
    685695                constant_top_momentumflux = .FALSE., &
    686696                constant_top_salinityflux = .TRUE., &
     697                constant_scalarflux = .TRUE., &             
    687698                constant_waterflux = .TRUE., create_disturbances = .TRUE., &
    688699                data_output_2d_on_each_pe = .TRUE., &
     
    739750                 alpha_surface = 0.0_wp, atmos_ocean_sign = 1.0_wp, &
    740751                 averaging_interval = 0.0_wp, averaging_interval_pr = 9999999.9_wp, &
    741                  bc_pt_t_val, bc_q_t_val, bottom_salinityflux = 0.0_wp, &
     752                 bc_pt_t_val, bc_q_t_val, bc_s_t_val, bottom_salinityflux = 0.0_wp, &
    742753                 building_height = 50.0_wp, building_length_x = 50.0_wp, &
    743754                 building_length_y = 50.0_wp, building_wall_left = 9999999.9_wp, &
     
    813824                 mask_scale(3), &
    814825                 pt_vertical_gradient(10) = 0.0_wp, &
    815                  pt_vertical_gradient_level(10) = -9999999.9_wp, &
     826                 pt_vertical_gradient_level(10) = -1.0_wp, &
    816827                 q_vertical_gradient(10) = 0.0_wp, &
    817828                 q_vertical_gradient_level(10) = -1.0_wp, &
     
    819830                 s_vertical_gradient_level(10) = -1.0_wp, &
    820831                 sa_vertical_gradient(10) = 0.0_wp, &
    821                  sa_vertical_gradient_level(10) = -9999999.9_wp, &
     832                 sa_vertical_gradient_level(10) = -1.0_wp, &
    822833                 skip_time_domask(max_masks) = 9999999.9_wp, threshold(20) = 0.0_wp, &
    823834                 time_domask(max_masks) = 0.0_wp, &
     
    833844                 wall_humidityflux(0:4) = 0.0_wp, wall_nrflux(0:4) = 0.0_wp, &
    834845                 wall_qflux(0:4) = 0.0_wp, wall_qrflux(0:4) = 0.0_wp, &
    835                  wall_salinityflux(0:4) = 0.0_wp, wall_scalarflux(0:4) = 0.0_wp, &
     846                 wall_salinityflux(0:4) = 0.0_wp, wall_sflux(0:4) = 0.0_wp, wall_scalarflux(0:4) = 0.0_wp, &
    836847                 subs_vertical_gradient(10) = 0.0_wp, &
    837848                 subs_vertical_gradient_level(10) = -9999999.9_wp
     
    12551266                                                  sums_wsnrs_ws_l,                &
    12561267                                                  sums_wspts_ws_l,                &
     1268                                                  sums_wsqs_ws_l,                 &   
     1269                                                  sums_wsqrs_ws_l,                &
    12571270                                                  sums_wssas_ws_l,                &
    1258                                                   sums_wsqs_ws_l,                 &
    1259                                                   sums_wsqrs_ws_l,                &
     1271                                                  sums_wsss_ws_l,                 &
    12601272                                                  sums_ls_l
    12611273                                             
Note: See TracChangeset for help on using the changeset viewer.