Changeset 2735 for palm/trunk/SOURCE/data_output_2d.f90
- Timestamp:
- Jan 11, 2018 12:01:27 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_2d.f90
r2718 r2735 25 25 ! ----------------- 26 26 ! $Id$ 27 ! output of r_a moved from land-surface to consider also urban-type surfaces 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 844 847 IF ( mode == 'xy' ) level_z = zu 845 848 846 849 CASE ( 'r_a*_xy' ) ! 2d-array 850 IF ( av == 0 ) THEN 851 DO m = 1, surf_lsm_h%ns 852 i = surf_lsm_h%i(m) 853 j = surf_lsm_h%j(m) 854 local_pf(i,j,nzb+1) = surf_lsm_h%r_a(m) 855 ENDDO 856 857 DO m = 1, surf_usm_h%ns 858 i = surf_usm_h%i(m) 859 j = surf_usm_h%j(m) 860 local_pf(i,j,nzb+1) = & 861 ( surf_usm_h%frac(0,m) * surf_usm_h%r_a(m) + & 862 surf_usm_h%frac(1,m) * surf_usm_h%r_a_green(m) + & 863 surf_usm_h%frac(2,m) * surf_usm_h%r_a_window(m) ) 864 ENDDO 865 866 ELSE 867 DO i = nxl, nxr 868 DO j = nys, nyn 869 local_pf(i,j,nzb+1) = r_a_av(j,i) 870 ENDDO 871 ENDDO 872 ENDIF 873 resorted = .TRUE. 874 two_d = .TRUE. 875 level_z(nzb+1) = zu(nzb+1) 847 876 848 877 CASE ( 'rho_ocean_xy', 'rho_ocean_xz', 'rho_ocean_yz' )
Note: See TracChangeset
for help on using the changeset viewer.