Changeset 48 for palm/trunk/SOURCE/data_output_tseries.f90
 Timestamp:
 Mar 6, 2007 12:28:36 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/data_output_tseries.f90
r4 r48 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! Collection of time series quantities moved to routine flow_statistics, 7 ! output for "profil" removed 7 8 ! 8 9 ! Former revisions: … … 37 38 38 39 INTEGER :: file_id, i, j, sr 39 REAL :: ts_value(30)40 40 41 41 … … 70 70 ! applicable) 71 71 DO sr = 0, statistic_regions 72 !73 ! Open file for time series output.74 IF ( profil_output ) THEN75 file_id = 50 + sr76 CALL check_open( file_id )77 ENDIF78 79 !80 ! Collect and printout all time series quantities in a single line.81 ts_value(1) = hom(nzb+4,1,var_hom,sr) ! E82 ts_value(2) = hom(nzb+5,1,var_hom,sr) ! E*83 ts_value(3) = dt_3d84 ts_value(4) = hom(nzb,1,var_hom,sr) ! u*85 ts_value(5) = hom(nzb+3,1,var_hom,sr) ! th*86 ts_value(6) = u_max87 ts_value(7) = v_max88 ts_value(8) = w_max89 ts_value(9) = hom(nzb+10,1,var_sum,sr) ! new divergence90 ts_value(10) = hom(nzb+9,1,var_hom,sr) ! old Divergence91 ts_value(11) = hom(nzb+6,1,var_hom,sr) ! z_i(1)92 ts_value(12) = hom(nzb+7,1,var_hom,sr) ! z_i(2)93 ts_value(13) = hom(nzb+8,1,var_hom,sr) ! w*94 ts_value(14) = hom(nzb,1,16,sr) ! w'pt' at k=095 ts_value(15) = hom(nzb+1,1,16,sr) ! w'pt' at k=196 ts_value(16) = hom(nzb+1,1,18,sr) ! wpt at k=197 ts_value(17) = hom(nzb,1,4,sr) ! pt(0)98 ts_value(18) = hom(nzb+1,1,4,sr) ! pt(zp)99 ts_value(19) = hom(nzb+9,1,var_hom1,sr) ! splptx100 ts_value(20) = hom(nzb+10,1,var_hom1,sr) ! splpty101 ts_value(21) = hom(nzb+11,1,var_hom1,sr) ! splptz102 IF ( ts_value(5) /= 0.0 ) THEN103 ts_value(22) = ts_value(4)**2 / &104 ( kappa * g * ts_value(5) / ts_value(18) ) ! L105 ELSE106 ts_value(22) = 10000.0107 ENDIF108 72 109 73 #if defined( __netcdf ) … … 111 75 DO i = 1, dots_num 112 76 nc_stat = NF90_PUT_VAR( id_set_ts, id_var_dots(i,sr), & 113 (/ ts_value(i ) /),&77 (/ ts_value(i,sr) /), & 114 78 start = (/ dots_time_count /), & 115 79 count = (/ 1 /) ) … … 119 83 #endif 120 84 121 IF ( profil_output ) THEN 122 WRITE ( file_id, 500 ) simulated_time, ts_value(1:22) 123 ! 124 ! yvalue range of the crosses to be drawn by PROFIL 125 ! If required, enlarge them, provided they have not yet been 126 ! specified in 127 ! check_parameters 128 DO i = 1, dots_n 129 130 j = dots_crossindex(i) 131 132 IF ( cross_ts_uymin(j) == 999.999 ) THEN 133 ! 134 ! When the value range of the first line in the corresponding 135 ! cross is determined, its value range is simply adopted. 136 IF ( cross_ts_uymin_computed(j) == 999.999 ) & 137 THEN 138 cross_ts_uymin_computed(j) = ts_value(dots_index(i)) 139 ELSE 140 cross_ts_uymin_computed(j) = & 141 MIN(cross_ts_uymin_computed(j),ts_value(dots_index(i))) 142 ENDIF 143 ENDIF 144 145 IF ( cross_ts_uymax(j) == 999.999 ) THEN 146 ! 147 ! When the value range of the first line in the corresponding 148 ! cross is determined, its value range is simply adopted. 149 IF ( cross_ts_uymax_computed(j) == 999.999 ) & 150 THEN 151 cross_ts_uymax_computed(j) = ts_value(dots_index(i)) 152 ELSE 153 cross_ts_uymax_computed(j) = & 154 MAX(cross_ts_uymax_computed(j),ts_value(dots_index(i))) 155 ENDIF 156 ENDIF 157 158 ENDDO 159 160 ENDIF 161 162 ENDDO ! Loop of subregions. 85 ENDDO 163 86 164 87 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.