Changeset 3773 for palm/trunk/SOURCE/sum_up_3d_data.f90
- Timestamp:
- Mar 1, 2019 8:56:57 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/sum_up_3d_data.f90
r3761 r3773 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Added output of theta_2m*_xy_av 28 ! 29 ! 3761 2019-02-25 15:31:42Z raasch 27 30 ! unused variables removed 28 31 ! … … 920 923 ENDIF 921 924 925 CASE ( 'theta_2m*' ) 926 IF ( ALLOCATED( pt_2m_av ) ) THEN 927 DO i = nxl, nxr 928 DO j = nys, nyn 929 match_def = surf_def_h(0)%start_index(j,i) <= & 930 surf_def_h(0)%end_index(j,i) 931 match_lsm = surf_lsm_h%start_index(j,i) <= & 932 surf_lsm_h%end_index(j,i) 933 match_usm = surf_usm_h%start_index(j,i) <= & 934 surf_usm_h%end_index(j,i) 935 936 IF ( match_def ) THEN 937 m = surf_def_h(0)%end_index(j,i) 938 pt_2m_av(j,i) = pt_2m_av(j,i) + & 939 surf_def_h(0)%pt_2m(m) 940 ELSEIF ( match_lsm .AND. .NOT. match_usm ) THEN 941 m = surf_lsm_h%end_index(j,i) 942 pt_2m_av(j,i) = pt_2m_av(j,i) + & 943 surf_lsm_h%pt_2m(m) 944 ELSEIF ( match_usm ) THEN 945 m = surf_usm_h%end_index(j,i) 946 pt_2m_av(j,i) = pt_2m_av(j,i) + & 947 surf_usm_h%pt_2m(m) 948 ENDIF 949 ENDDO 950 ENDDO 951 ENDIF 952 953 922 954 CASE ( 't*' ) 923 955 IF ( ALLOCATED( ts_av ) ) THEN
Note: See TracChangeset
for help on using the changeset viewer.