Changeset 1701 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Nov 2, 2015 7:43:04 AM (8 years ago)

Author:

maronga

Message:

minor bugfixes for radiation model. bugfix in subjob

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (5 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -1196 +1196 @@
                    sums_l(k,108,tn)  = sums_l(k,108,tn)  + rad_sw_cs_hr(k,j,i) &
                                        * rmask(j,i,sr)
-                   sums_l(k,109,tn)  = sums_l(k,108,tn)  + rad_sw_hr(k,j,i)    &
+                   sums_l(k,109,tn)  = sums_l(k,109,tn)  + rad_sw_hr(k,j,i)    &
                                        * rmask(j,i,sr)
                 ENDDO
@@ -1401 +1401 @@
 
           IF ( radiation_scheme == 'rrtmg' )  THEN
+#if defined ( __rrtmg )
              hom(:,1,106,sr) = sums(:,106)            ! rad_lw_cs_hr
              hom(:,1,107,sr) = sums(:,107)            ! rad_lw_hr
@@ -1406 +1407 @@
              hom(:,1,109,sr) = sums(:,109)            ! rad_sw_hr
 
-#if defined ( __rrtmg )
              hom(:,1,110,sr) = sums(:,110)            ! rrtm_aldif
              hom(:,1,111,sr) = sums(:,111)            ! rrtm_aldir
@@ -1588 +1588 @@
           ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
           ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
-          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
-          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
+          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
+          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
 
 #if defined ( __rrtmg )
@@ -3529 +3529 @@
           ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
           ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
-          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
-          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
+          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
+          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
 
 #if defined ( __rrtmg )