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/data_output_2d.f90

    r1552 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    139139               m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
    140140               qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
    141                qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
    142                t_soil_av, zs
     141               qsws_veg_eb, qsws_veg_eb_av, r_a, r_a_av, r_s, r_s_av, shf_eb,  &
     142               shf_eb_av, t_soil, t_soil_av, zs
    143143   
    144144    USE netcdf_control
     
    934934                   to_be_resorted => rho_av
    935935                ENDIF
     936
     937             CASE ( 'r_a*_xy' )        ! 2d-array
     938                IF ( av == 0 )  THEN
     939                   DO  i = nxlg, nxrg
     940                      DO  j = nysg, nyng
     941                         local_pf(i,j,nzb+1) = r_a(j,i)
     942                      ENDDO
     943                   ENDDO
     944                ELSE
     945                   DO  i = nxlg, nxrg
     946                      DO  j = nysg, nyng
     947                         local_pf(i,j,nzb+1) = r_a_av(j,i)
     948                      ENDDO
     949                   ENDDO
     950                ENDIF
     951                resorted = .TRUE.
     952                two_d = .TRUE.
     953                level_z(nzb+1) = zu(nzb+1)
     954
     955             CASE ( 'r_s*_xy' )        ! 2d-array
     956                IF ( av == 0 )  THEN
     957                   DO  i = nxlg, nxrg
     958                      DO  j = nysg, nyng
     959                         local_pf(i,j,nzb+1) = r_s(j,i)
     960                      ENDDO
     961                   ENDDO
     962                ELSE
     963                   DO  i = nxlg, nxrg
     964                      DO  j = nysg, nyng
     965                         local_pf(i,j,nzb+1) = r_s_av(j,i)
     966                      ENDDO
     967                   ENDDO
     968                ENDIF
     969                resorted = .TRUE.
     970                two_d = .TRUE.
     971                level_z(nzb+1) = zu(nzb+1)
    936972
    937973             CASE ( 's_xy', 's_xz', 's_yz' )
Note: See TracChangeset for help on using the changeset viewer.