Changeset 1239 for palm/trunk/SOURCE/microphysics.f90

Timestamp:

Oct 29, 2013 10:11:53 AM (10 years ago)

Author:

heinze

Message:

routines for nudging and large scale forcing from external file added

File:

• palm/trunk/SOURCE/microphysics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/microphysics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! hyp and rho have to be calculated at each time step if data from external
+! file LSF_DATA are used
 !
 ! Former revisions:
@@ -103 +104 @@
 
        USE arrays_3d
-       USE control_parameters
+       USE cloud_parameters
+       USE control_parameters
+       USE grid_variables
        USE indices
        USE statistics
@@ -296 +299 @@
        USE cloud_parameters
        USE control_parameters
+       USE grid_variables
+       USE indices
        USE statistics
 
        IMPLICIT NONE
 
-       INTEGER ::  i, j
+       INTEGER ::  i, j, k 
+       REAL    ::  t_surface
+
+       IF ( large_scale_forcing .AND. lsf_surf ) THEN
+!
+!--       Calculate:
+!--       pt / t : ratio of potential and actual temperature (pt_d_t)
+!--       t / pt : ratio of actual and potential temperature (t_d_pt)
+!--       p_0(z) : vertical profile of the hydrostatic pressure (hyp)
+          t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
+          DO  k = nzb, nzt+1
+             hyp(k)    = surface_pressure * 100.0 * &
+                         ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
+             pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
+             t_d_pt(k) = 1.0 / pt_d_t(k)
+             hyrho(k)  = hyp(k) / ( r_d * t_d_pt(k) * pt_init(k) )       
+          ENDDO
+!
+!--       Compute reference density
+          rho_surface = surface_pressure * 100.0 / ( r_d * t_surface )
+       ENDIF
+
 
        dt_micro = dt_3d * weight_pres(intermediate_timestep_count)