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

Timestamp:

Mar 20, 2014 4:38:49 PM (10 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

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

palm/trunk/SOURCE/calc_radiation.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! exponent 4.0 changed to integer
 !
 ! Former revisions:
@@ -127 +127 @@
 
              temperature     = pt(nzb,j,i) * t_d_pt(nzb) + l_d_cp * ql(nzb,j,i)
-             blackbody_emission(nzb) = sigma * temperature**4.0
+             blackbody_emission(nzb) = sigma * temperature**4
 
              DO  k = nzb_2d(j,i)+1, nzt
@@ -139 +139 @@
 
                 temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
-                blackbody_emission(k) = sigma * temperature**4.0
+                blackbody_emission(k) = sigma * temperature**4
 
              ENDDO
@@ -149 +149 @@
              temperature       = pt(nzt+1,j,i) * t_d_pt(nzt+1) + l_d_cp *      &
                                  ql(nzt+1,j,i)
-             blackbody_emission(nzt+1) = sigma * temperature**4.0
+             blackbody_emission(nzt+1) = sigma * temperature**4
 
 !
@@ -267 +267 @@
 
        temperature     = pt(nzb,j,i) * t_d_pt(nzb) + l_d_cp * ql(nzb,j,i)
-       blackbody_emission(nzb) = sigma * temperature**4.0
+       blackbody_emission(nzb) = sigma * temperature**4
 
        DO  k = nzb_2d(j,i)+1, nzt
@@ -277 +277 @@
 
           temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
-          blackbody_emission(k) = sigma * temperature**4.0
+          blackbody_emission(k) = sigma * temperature**4
 
        ENDDO
@@ -286 +286 @@
        temperature       = pt(nzt+1,j,i) * t_d_pt(nzt+1) + l_d_cp *            &
                            ql(nzt+1,j,i)
-       blackbody_emission(nzt+1) = sigma * temperature**4.0
+       blackbody_emission(nzt+1) = sigma * temperature**4
 
 !