Changeset 1380 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Apr 28, 2014 12:40:45 PM (10 years ago)

Author:

heinze

Message:

Upper boundary conditions for pt and q in case of nudging adjusted

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (7 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! CALL of nudge_ref added 
+! bc_pt_t_val and bc_q_t_val are updated in case nudging is used
 ! 
 ! Former revisions:
@@ -131 +132 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  diss, e_p, nr_p, prho, pt, pt_p, q, ql, ql_c, ql_v, ql_vp, qr_p,&
-               q_p, ref_state, rho, sa_p, tend, u, u_p, v, vpt, v_p, w_p
+        ONLY:  diss, dzu, e_p, nr_p, prho, pt, pt_p, pt_init, q_init, q, ql,   &
+               ql_c, ql_v, ql_vp, qr_p, q_p, ref_state, rho, sa_p, tend, u,    &
+               u_p, v, vpt, v_p, w_p
 
     USE calc_mean_profile_mod,                                                 &
@@ -140 +142 @@
         ONLY:  advected_distance_x, advected_distance_y, average_count_3d,     &
                average_count_sp, averaging_interval, averaging_interval_pr,    &
-               averaging_interval_sp, bc_lr_cyc, bc_ns_cyc,                    &
-               call_psolver_at_all_substeps, cloud_droplets, cloud_physics,    &
-               constant_heatflux, create_disturbances, dopr_n,                 &
+               averaging_interval_sp, bc_lr_cyc, bc_ns_cyc, bc_pt_t_val,       &
+               bc_q_t_val, call_psolver_at_all_substeps, cloud_droplets,       &
+               cloud_physics, constant_heatflux, create_disturbances, dopr_n,  &
                constant_diffusion, coupling_mode, coupling_start_time,         &
                current_timestep_number, disturbance_created,                   &
@@ -176 +178 @@
     USE indices,                                                               &
         ONLY:  i_left, i_right, j_north, j_south, nbgp, nx, nxl, nxlg, nxr,    &
-               nxrg, nyn, nys, nzb, nzb_u_inner, nzb_v_inner 
+               nxrg, nyn, nys, nzb, nzt, nzb_u_inner, nzb_v_inner 
 
     USE interaction_droplets_ptq_mod,                                          &
@@ -187 +189 @@
 
     USE nudge_mod,                                                             &
-        ONLY:  calc_tnudge
+        ONLY:  calc_tnudge, nudge_ref
 
     USE particle_attributes,                                                   &
@@ -267 +269 @@
 
 !
+!--    Set pt_init and q_init to the current profiles taken from 
+!--    NUDGING_DATA 
+       IF ( nudging )  THEN
+           CALL nudge_ref ( simulated_time )
+!
+!--        Store temperature gradient at the top boundary for possible Neumann
+!--        boundary condition
+           bc_pt_t_val = ( pt_init(nzt+1) - pt_init(nzt) ) / dzu(nzt+1)
+           bc_q_t_val  = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
+       ENDIF
+
+!
 !--    Execute the user-defined actions
        CALL user_actions( 'before_timestep' )
@@ -312 +326 @@
 !
 !--       In case of nudging calculate current nudging time scale and horizontal
-!--       means of u,v,pt and q
+!--       means of u, v, pt and q
           IF ( nudging )  THEN
              CALL calc_tnudge( simulated_time )