Changeset 1353 for palm/trunk/SOURCE/calc_radiation.f90

Timestamp:

Apr 8, 2014 3:21:23 PM (10 years ago)

Author:

heinze

Message:

REAL constants provided with KIND-attribute

File:

• palm/trunk/SOURCE/calc_radiation.f90 (modified) (8 diffs)

palm/trunk/SOURCE/calc_radiation.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! REAL constants provided with KIND-attribute
 !
 ! Former revisions:
@@ -55 +55 @@
     
     LOGICAL, SAVE ::  first_call = .TRUE. !:
-    REAL(wp), SAVE ::  sigma = 5.67E-08   !:
+    REAL(wp), SAVE ::  sigma = 5.67E-08_wp   !:
 
     REAL(wp), DIMENSION(:), ALLOCATABLE, SAVE ::  lwp_ground         !:
@@ -126 +126 @@
 !--          Compute the liquid water path (LWP) and blackbody_emission
 !--          at all vertical levels
-             lwp_ground(nzb) = 0.0
+             lwp_ground(nzb) = 0.0_wp
              lwp_top(nzt+1)  = rho_surface * ql(nzt+1,j,i) * dzw(nzt+1)
 
@@ -156 +156 @@
 !
 !--          See Chlond '92, this is just a first guess
-             impinging_flux_at_top = blackbody_emission(nzb) - 100.0
+             impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
 
              DO  k = nzb_2d(j,i)+1, nzt
@@ -162 +162 @@
 !--             Save some computational time, but this may cause load
 !--             imbalances if ql is not distributed uniformly
-                IF ( ql(k,j,i) /= 0.0 ) THEN
+                IF ( ql(k,j,i) /= 0.0_wp ) THEN
 !
 !--                Compute effective emissivities
-                   effective_emission_up_p   = 1.0 -                           &
-                                               EXP( -130.0 * lwp_ground(k+1) )
-                   effective_emission_up_m   = 1.0 -                           &
-                                               EXP( -130.0 * lwp_ground(k-1) )
-                   effective_emission_down_p = 1.0 -                           &
-                                               EXP( -158.0 * lwp_top(k+1) )
-                   effective_emission_down_m = 1.0 -                           &
-                                               EXP( -158.0 * lwp_top(k-1) )  
+                   effective_emission_up_p   = 1.0_wp -                        &
+                                               EXP( -130.0_wp * lwp_ground(k+1) )
+                   effective_emission_up_m   = 1.0_wp -                        &
+                                               EXP( -130.0_wp * lwp_ground(k-1) )
+                   effective_emission_down_p = 1.0_wp -                        &
+                                               EXP( -158.0_wp * lwp_top(k+1) )
+                   effective_emission_down_m = 1.0_wp -                        &
+                                               EXP( -158.0_wp * lwp_top(k-1) )  
 
 !
@@ -266 +266 @@
 !--    Compute the liquid water path (LWP) and blackbody_emission
 !--    at all vertical levels
-       lwp_ground(nzb) = 0.0
+       lwp_ground(nzb) = 0.0_wp
        lwp_top(nzt+1)  = rho_surface * ql(nzt+1,j,i) * dzw(nzt+1)
 
@@ -293 +293 @@
 !
 !--    See Chlond '92, this is just a first guess
-       impinging_flux_at_top = blackbody_emission(nzb) - 100.0
+       impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
 
        DO  k = nzb_2d(j,i)+1, nzt
@@ -299 +299 @@
 !--       Store some computational time,
 !--       this may cause load imbalances if ql is not distributed uniformly
-          IF ( ql(k,j,i) /= 0.0 ) THEN
+          IF ( ql(k,j,i) /= 0.0_wp ) THEN
 !
 !--          Compute effective emissivities
-             effective_emission_up_p   = 1.0 -                                 &
-                                         EXP( -130.0 * lwp_ground(k+1) )
-             effective_emission_up_m   = 1.0 -                                 &
-                                         EXP( -130.0 * lwp_ground(k-1) )
-             effective_emission_down_p = 1.0 -                                 &
-                                         EXP( -158.0 * lwp_top(k+1) )
-             effective_emission_down_m = 1.0 -                                 &
-                                         EXP( -158.0 * lwp_top(k-1) )  
+             effective_emission_up_p   = 1.0_wp -                              &
+                                         EXP( -130.0_wp * lwp_ground(k+1) )
+             effective_emission_up_m   = 1.0_wp -                              &
+                                         EXP( -130.0_wp * lwp_ground(k-1) )
+             effective_emission_down_p = 1.0_wp -                              &
+                                         EXP( -158.0_wp * lwp_top(k+1) )
+             effective_emission_down_m = 1.0_wp -                              &
+                                         EXP( -158.0_wp * lwp_top(k-1) )  
              
 !