Changeset 1949 for palm/trunk/SOURCE


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r1857 r1949  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfix: calculation of qsws_soil_eb with precipitation = .TRUE. gave
     22! qsws_soil_eb = 0 due to a typo
    2223!
    2324! Former revisions:
     
    14151416!--                Add precipitation to bare soil according to the bare soil
    14161417!--                coverage.
    1417                    qsws_soil_eb(j,i) = qsws_soil_eb(j,i) * (1.0_wp             &
     1418                   qsws_soil_eb(j,i) = qsws_soil_eb(j,i) + (1.0_wp             &
    14181419                                       - c_veg(j,i)) * prr(k,j,i) * hyrho(k)   &
    14191420                                       * 0.001_wp * rho_l * l_v
     
    14321433                      qsws_soil_eb(j,i) = qsws_soil_eb(j,i)                    &
    14331434                                           + qsws_liq_eb(j,i)
     1435
    14341436                      qsws_liq_eb(j,i)  = 0.0_wp
    14351437                   ENDIF
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r1852 r1949  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfix: calculation of lai_av, c_veg_av and c_liq_av.
    2222!
    2323! Former revisions:
     
    530530             DO  i = nxlg, nxrg
    531531                DO  j = nysg, nyng
    532                    c_liq_av(j,i) = c_liq_av(j,i)
     532                   c_liq_av(j,i) = c_liq_av(j,i) + c_liq(j,i)
    533533                ENDDO
    534534             ENDDO
     
    537537             DO  i = nxlg, nxrg
    538538                DO  j = nysg, nyng
    539                    c_soil_av(j,i) = c_soil_av(j,i) + (1.0_wp - c_veg(j,i))
     539                   c_soil_av(j,i) = c_soil_av(j,i) + (1.0 - c_veg(j,i))
    540540                ENDDO
    541541             ENDDO
     
    544544             DO  i = nxlg, nxrg
    545545                DO  j = nysg, nyng
    546                    c_veg_av(j,i) = c_veg_av(j,i)
     546                   c_veg_av(j,i) = c_veg_av(j,i) + c_veg(j,i)
    547547                ENDDO
    548548             ENDDO
     
    567567             DO  i = nxlg, nxrg
    568568                DO  j = nysg, nyng
    569                    lai_av(j,i) = lai_av(j,i)
     569                   lai_av(j,i) = lai_av(j,i) + lai(j,i)
    570570                ENDDO
    571571             ENDDO
Note: See TracChangeset for help on using the changeset viewer.