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

Timestamp:

Mar 4, 2015 5:44:27 PM (9 years ago)

Author:

maronga

Message:

LSM output of r_a and r_s added

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added output of r_a and r_s.
 ! 
 ! Former revisions:
@@ -160 +160 @@
     USE land_surface_model_mod,                                                &
         ONLY:   dots_soil, ghf_eb, land_surface, m_soil, nzb_soil, nzt_soil,   &
-                qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, shf_eb,       &
-                t_soil
+                qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, r_a, r_s,     &
+                shf_eb, t_soil
 
     USE pegrid
@@ -717 +717 @@
 
              IF ( land_surface )  THEN
-                sums_l(nzb,93,tn) = sums_l(nzb,93,tn) + ghf_eb(j,i)
-                sums_l(nzb,94,tn) = sums_l(nzb,94,tn) + shf_eb(j,i)
-                sums_l(nzb,95,tn) = sums_l(nzb,95,tn) + qsws_eb(j,i)
-                sums_l(nzb,96,tn) = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
-                sums_l(nzb,97,tn) = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
-                sums_l(nzb,98,tn) = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
+                sums_l(nzb,93,tn)  = sums_l(nzb,93,tn) + ghf_eb(j,i)
+                sums_l(nzb,94,tn)  = sums_l(nzb,94,tn) + shf_eb(j,i)
+                sums_l(nzb,95,tn)  = sums_l(nzb,95,tn) + qsws_eb(j,i)
+                sums_l(nzb,96,tn)  = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
+                sums_l(nzb,97,tn)  = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
+                sums_l(nzb,98,tn)  = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
+                sums_l(nzb,99,tn)  = sums_l(nzb,99,tn) + r_a(j,i)
+                sums_l(nzb,100,tn) = sums_l(nzb,100,tn)+ r_s(j,i)
              ENDIF
 
              IF ( radiation )  THEN
-                sums_l(nzb,99,tn)  = sums_l(nzb,99,tn)  + rad_net(j,i)
-                sums_l(nzb,100,tn) = sums_l(nzb,100,tn) + rad_sw_in(j,i)
+                sums_l(nzb,101,tn)  = sums_l(nzb,101,tn)  + rad_net(j,i)
+                sums_l(nzb,102,tn) = sums_l(nzb,102,tn) + rad_sw_in(j,i)
              ENDIF
 
@@ -1293 +1295 @@
           hom(:,1,91,sr) = sums(:,91)              ! m_soil
                                                    ! 92 is initial m_soil profile
-          hom(:,1,93,sr) = sums(:,93)              ! ghf_eb
-          hom(:,1,94,sr) = sums(:,94)              ! shf_eb
-          hom(:,1,95,sr) = sums(:,95)              ! qsws_eb
-          hom(:,1,96,sr) = sums(:,96)              ! qsws_liq_eb
-          hom(:,1,97,sr) = sums(:,97)              ! qsws_soil_eb
-          hom(:,1,98,sr) = sums(:,98)              ! qsws_veg_eb
+          hom(:,1,93,sr)  = sums(:,93)             ! ghf_eb
+          hom(:,1,94,sr)  = sums(:,94)             ! shf_eb
+          hom(:,1,95,sr)  = sums(:,95)             ! qsws_eb
+          hom(:,1,96,sr)  = sums(:,96)             ! qsws_liq_eb
+          hom(:,1,97,sr)  = sums(:,97)             ! qsws_soil_eb
+          hom(:,1,98,sr)  = sums(:,98)             ! qsws_veg_eb
+          hom(:,1,99,sr)  = sums(:,99)             ! r_a
+          hom(:,1,100,sr) = sums(:,100)            ! r_s
+
        ENDIF
 
        IF ( radiation )  THEN
-          hom(:,1,99 ,sr) = sums(:,99)             ! rad_net
-          hom(:,1,100,sr) = sums(:,100)            ! rad_sw_in
+          hom(:,1,101 ,sr) = sums(:,101)            ! rad_net
+          hom(:,1,102,sr)  = sums(:,102)            ! rad_sw_in
        ENDIF
 
@@ -1460 +1465 @@
           ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
           ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
+          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
+          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
        ENDIF
 !
@@ -3263 +3270 @@
           ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
           ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
+          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
+          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
        ENDIF
 !