Changeset 329
- Timestamp:
- Jun 2, 2009 2:16:11 PM (15 years ago)
- Location:
- palm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/cross_sections.ncl
r286 r329 1019 1019 1020 1020 if (xyc .EQ. 1) then 1021 d= (ye-ys+1)/(major_ticks_y-1)1021 d= (ye-ys+1)/(major_ticks_y-1) 1022 1022 e=(xe-xs+1)/(major_ticks_x-1) 1023 1023 array_yl =new(major_ticks_y,integer) … … 1036 1036 if (ar .GT. 0) 1037 1037 do min_ar=0,3 1038 array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)1038 array_minor_yl(4*(ar-1)+min_ar)= y_d(array_yl(ar-1))+(y_d(array_yl(ar))-y_d(array_yl(ar-1)))/5*(min_ar+1) 1039 1039 end do 1040 1040 end if … … 1045 1045 if (br .GT. 0) 1046 1046 do min_br=0,3 1047 array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)1047 array_minor_xb(4*(br-1)+min_br)= x_d(array_xb(br-1))+(x_d(array_xb(br))-x_d(array_xb(br-1)))/5*(min_br+1) 1048 1048 end do 1049 1049 end if … … 1069 1069 if (ar .GT. 0) 1070 1070 do min_ar=0,3 1071 array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)1071 array_minor_yl(4*(ar-1)+min_ar)= z_d(array_yl(ar-1))+(z_d(array_yl(ar))-z_d(array_yl(ar-1)))/5*(min_ar+1) 1072 1072 end do 1073 1073 end if … … 1078 1078 if (br .GT. 0) 1079 1079 do min_br=0,3 1080 array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)1080 array_minor_xb(4*(br-1)+min_br)= x_d(array_xb(br-1))+(x_d(array_xb(br))-x_d(array_xb(br-1)))/5*(min_br+1) 1081 1081 end do 1082 1082 end if … … 1099 1099 do ar=1,major_ticks_y-1 1100 1100 array_yl(ar)=d*(ar-1)+d-1 1101 array_yl_labels(ar) = y_d(array_yl(ar))1101 array_yl_labels(ar) = z_d(array_yl(ar)) 1102 1102 if (ar .GT. 0) 1103 1103 do min_ar=0,3 1104 array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)1104 array_minor_yl(4*(ar-1)+min_ar)= z_d(array_yl(ar-1))+(z_d(array_yl(ar))-z_d(array_yl(ar-1)))/5*(min_ar+1) 1105 1105 end do 1106 1106 end if … … 1108 1108 do br=1,major_ticks_x-1 1109 1109 array_xb(br)=e*(br-1)+e-1 1110 array_xb_labels(br) = x_d(array_xb(br))1110 array_xb_labels(br) = y_d(array_xb(br)) 1111 1111 if (br .GT. 0) 1112 1112 do min_br=0,3 1113 array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)1113 array_minor_xb(4*(br-1)+min_br)= y_d(array_xb(br-1))+(y_d(array_xb(br))-y_d(array_xb(br-1)))/5*(min_br+1) 1114 1114 end do 1115 1115 end if 1116 1116 end do 1117 1117 end if 1118 1118 1119 1119 if (axes_explicit .EQ. 1)then 1120 1120 cs_res@tmYLMode = "Explicit" 1121 1121 cs_res@tmXBMode = "Explicit" 1122 cs_res@tmYLValues = array_yl 1123 if (xyc .EQ. 1)then 1122 if (xyc .EQ. 1)then 1123 cs_res@tmYLValues = y_d(array_yl) 1124 cs_res@tmXBValues = x_d(array_xb) 1124 1125 cs_res@tmYLLabels = array_yl_labels/norm_y 1125 1126 cs_res@tmXBLabels = array_xb_labels/norm_x … … 1136 1137 end if 1137 1138 if (xzc .EQ. 1)then 1139 cs_res@tmYLValues = z_d(array_yl) 1140 cs_res@tmXBValues = x_d(array_xb) 1138 1141 cs_res@tmYLLabels = array_yl_labels/norm_z 1139 1142 cs_res@tmXBLabels = array_xb_labels/norm_x … … 1150 1153 end if 1151 1154 if (yzc .EQ. 1)then 1155 cs_res@tmYLValues = z_d(array_yl) 1156 cs_res@tmXBValues = y_d(array_xb) 1152 1157 cs_res@tmYLLabels = array_yl_labels/norm_z 1153 1158 cs_res@tmXBLabels = array_xb_labels/norm_y … … 1163 1168 end if 1164 1169 end if 1165 cs_res@tmXBValues = array_xb1166 1170 cs_res@tmYLMinorValues = array_minor_yl 1167 1171 cs_res@tmXBMinorValues = array_minor_xb … … 1301 1305 temp = f[:]->$vNam(varn)$ 1302 1306 data_att = f_att->$vNam(varn)$ 1303 data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe) 1307 if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") 1308 ;these variables depend von zu1_xy and that's why they have only one z-layer 1309 data(varn,:,0,:,:)=temp(:,0,ys:ye,xs:xe) 1310 else 1311 data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe) 1312 end if 1313 delete(temp) 1304 1314 end if 1305 1315 if ( xzc .eq. 1 ) then … … 1330 1340 1331 1341 unit(varn) = data_att@units 1342 delete(data_att) 1332 1343 1333 1344 end if -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r328 r329 43 43 Changed: 44 44 ------- 45 Small change in formatting of the message handling routine concering the output in the 46 job protocoll. (message) 47 45 48 initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill', now 46 49 independent of turbulent_inflow (check_parameters, header, init_3d_model) … … 85 88 Default value of grid_matching changed to strict. 86 89 87 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list90 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list 88 91 89 92 -
palm/trunk/SOURCE/message.f90
r311 r329 138 138 WRITE( file_id, '(4X,A)' ) TRIM( information_string_1 ) 139 139 WRITE( file_id, '(4X,A)' ) TRIM( information_string_2 ) 140 WRITE( *, '(4X,A)' ) '' 140 141 ! 141 142 !-- Flush buffer, if requested
Note: See TracChangeset
for help on using the changeset viewer.