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

Timestamp:

Oct 29, 2013 1:21:31 PM (10 years ago)

Author:

heinze

Message:

Undoing commit 1239

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:
@@ -104 +103 @@
 
        USE arrays_3d
-       USE cloud_parameters
-       USE control_parameters
-       USE grid_variables
+       USE control_parameters
        USE indices
        USE statistics
@@ -299 +296 @@
        USE cloud_parameters
        USE control_parameters
-       USE grid_variables
-       USE indices
        USE statistics
 
        IMPLICIT NONE
 
-       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
-
+       INTEGER ::  i, j
 
        dt_micro = dt_3d * weight_pres(intermediate_timestep_count)