Ignore:
Timestamp:
Mar 4, 2015 5:44:27 PM (7 years ago)
Author:
maronga
Message:

LSM output of r_a and r_s added

File:
1 edited

Legend:

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

    r1552 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    102102               m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
    103103               qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
    104                qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
    105                t_soil_av
     104               qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, r_a, r_a_av,    &
     105               r_s, r_s_av, t_soil, t_soil_av
    106106
    107107    USE particle_attributes,                                                   &
     
    334334                ENDIF
    335335                rho_av = 0.0_wp
     336
     337             CASE ( 'r_a*' )
     338                IF ( .NOT. ALLOCATED( r_a_av ) )  THEN
     339                   ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
     340                ENDIF
     341                r_a_av = 0.0_wp
     342
     343             CASE ( 'r_s*' )
     344                IF ( .NOT. ALLOCATED( r_s_av ) )  THEN
     345                   ALLOCATE( r_s_av(nysg:nyng,nxlg:nxrg) )
     346                ENDIF
     347                r_s_av = 0.0_wp
    336348
    337349             CASE ( 's' )
     
    724736             ENDDO
    725737
     738          CASE ( 'r_a*' )
     739             DO  i = nxlg, nxrg
     740                DO  j = nysg, nyng
     741                   r_a_av(j,i) = r_a_av(j,i) + r_a(j,i)
     742                ENDDO
     743             ENDDO
     744
     745          CASE ( 'r_s*' )
     746             DO  i = nxlg, nxrg
     747                DO  j = nysg, nyng
     748                   r_s_av(j,i) = r_s_av(j,i) + r_s(j,i)
     749                ENDDO
     750             ENDDO
     751
    726752          CASE ( 'rho' )
    727753             DO  i = nxlg, nxrg
Note: See TracChangeset for help on using the changeset viewer.