Changeset 1949 for palm/trunk/SOURCE

Timestamp:

Jun 17, 2016 7:19:16 AM (8 years ago)

Author:

maronga

Message:

bugfixes in output of land surface variables

Location:

palm/trunk/SOURCE

Files:

• land_surface_model_mod.f90 (modified) (3 diffs)
• sum_up_3d_data.f90 (modified) (5 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: calculation of qsws_soil_eb with precipitation = .TRUE. gave
+! qsws_soil_eb = 0 due to a typo
 ! 
 ! Former revisions:
@@ -1415 +1416 @@
 !--                Add precipitation to bare soil according to the bare soil
 !--                coverage.
-                   qsws_soil_eb(j,i) = qsws_soil_eb(j,i) * (1.0_wp             &
+                   qsws_soil_eb(j,i) = qsws_soil_eb(j,i) + (1.0_wp             &
                                        - c_veg(j,i)) * prr(k,j,i) * hyrho(k)   &
                                        * 0.001_wp * rho_l * l_v
@@ -1432 +1433 @@
                       qsws_soil_eb(j,i) = qsws_soil_eb(j,i)                    &
                                            + qsws_liq_eb(j,i)
+
                       qsws_liq_eb(j,i)  = 0.0_wp
                    ENDIF

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: calculation of lai_av, c_veg_av and c_liq_av.
 ! 
 ! Former revisions:
@@ -530 +530 @@
              DO  i = nxlg, nxrg
                 DO  j = nysg, nyng
-                   c_liq_av(j,i) = c_liq_av(j,i)
+                   c_liq_av(j,i) = c_liq_av(j,i) + c_liq(j,i)
                 ENDDO
              ENDDO
@@ -537 +537 @@
              DO  i = nxlg, nxrg
                 DO  j = nysg, nyng
-                   c_soil_av(j,i) = c_soil_av(j,i) + (1.0_wp - c_veg(j,i))
+                   c_soil_av(j,i) = c_soil_av(j,i) + (1.0 - c_veg(j,i))
                 ENDDO
              ENDDO
@@ -544 +544 @@
              DO  i = nxlg, nxrg
                 DO  j = nysg, nyng
-                   c_veg_av(j,i) = c_veg_av(j,i)
+                   c_veg_av(j,i) = c_veg_av(j,i) + c_veg(j,i)
                 ENDDO
              ENDDO
@@ -567 +567 @@
              DO  i = nxlg, nxrg
                 DO  j = nysg, nyng
-                   lai_av(j,i) = lai_av(j,i)
+                   lai_av(j,i) = lai_av(j,i) + lai(j,i)
                 ENDDO
              ENDDO