Changeset 87 for palm/trunk/SOURCE
- Timestamp:
- May 22, 2007 3:46:47 PM (18 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 13 edited
-
CURRENT_MODIFICATIONS (modified) (2 diffs)
-
check_parameters.f90 (modified) (17 diffs)
-
data_output_profiles.f90 (modified) (14 diffs)
-
flow_statistics.f90 (modified) (11 diffs)
-
header.f90 (modified) (3 diffs)
-
init_3d_model.f90 (modified) (3 diffs)
-
modules.f90 (modified) (6 diffs)
-
parin.f90 (modified) (3 diffs)
-
read_var_list.f90 (modified) (7 diffs)
-
run_control.f90 (modified) (2 diffs)
-
time_integration.f90 (modified) (2 diffs)
-
user_interface.f90 (modified) (7 diffs)
-
write_var_list.f90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r86 r87 1 1 New: 2 2 --- 3 Calculation and output of user-defined profiles. New &userpar parameters data_output_pr_user and max_pr_user. 4 5 check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list 3 6 4 7 5 8 Changed: 6 9 ------- 7 Division through dt_3d replaced by multiplication of the inverse. For performance optimisation, this is done in the loop calculating the divergence instead of using a seperate loop. (pres.f90) 10 Division through dt_3d replaced by multiplication of the inverse. For performance optimisation, this is done in the loop calculating the divergence instead of using a seperate loop. (pres.f90) var_hom and var_sum renamed pr_palm. 8 11 9 pres 12 data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration 10 13 11 14 … … 16 19 Bugfix: field_chr renamed field_char (user_interface). 17 20 18 poisfft, user_interface 21 Bugfix: output of use_upstream_for_tke (header). 22 23 header, poisfft, user_interface -
palm/trunk/SOURCE/check_parameters.f90
r77 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Check for user-defined profiles 7 7 ! 8 8 ! Former revisions: … … 1214 1214 CASE ( 'u', '#u' ) 1215 1215 dopr_index(i) = 1 1216 dopr_unit(i) = 'm/s' 1216 1217 hom(:,2,1,:) = SPREAD( zu, 2, statistic_regions+1 ) 1217 1218 IF ( data_output_pr(i)(1:1) == '#' ) THEN … … 1223 1224 CASE ( 'v', '#v' ) 1224 1225 dopr_index(i) = 2 1225 hom(:,2,2,:) = SPREAD( zu, 2, statistic_regions+1 ) 1226 dopr_unit(i) = 'm/s' 1227 hom(:,2,2,:) = SPREAD( zu, 2, statistic_regions+1 ) 1226 1228 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1227 1229 dopr_initial_index(i) = 6 … … 1232 1234 CASE ( 'w' ) 1233 1235 dopr_index(i) = 3 1234 hom(:,2,3,:) = SPREAD( zw, 2, statistic_regions+1 ) 1236 dopr_unit(i) = 'm/s' 1237 hom(:,2,3,:) = SPREAD( zw, 2, statistic_regions+1 ) 1235 1238 1236 1239 CASE ( 'pt', '#pt' ) 1237 1240 IF ( .NOT. cloud_physics ) THEN 1238 1241 dopr_index(i) = 4 1242 dopr_unit(i) = 'K' 1243 hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 ) 1244 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1245 dopr_initial_index(i) = 7 1246 hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 ) 1247 hom(nzb,2,7,:) = 0.0 ! because zu(nzb) is negative 1248 data_output_pr(i) = data_output_pr(i)(2:) 1249 ENDIF 1250 ELSE 1251 dopr_index(i) = 43 1252 dopr_unit(i) = 'K' 1253 hom(:,2,43,:) = SPREAD( zu, 2, statistic_regions+1 ) 1254 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1255 dopr_initial_index(i) = 28 1256 hom(:,2,28,:) = SPREAD( zu, 2, statistic_regions+1 ) 1257 hom(nzb,2,28,:) = 0.0 ! because zu(nzb) is negative 1258 data_output_pr(i) = data_output_pr(i)(2:) 1259 ENDIF 1260 ENDIF 1261 1262 CASE ( 'e' ) 1263 dopr_index(i) = 8 1264 dopr_unit(i) = 'm2/s2' 1265 hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 ) 1266 hom(nzb,2,8,:) = 0.0 1267 1268 CASE ( 'km', '#km' ) 1269 dopr_index(i) = 9 1270 dopr_unit(i) = 'm2/s' 1271 hom(:,2,9,:) = SPREAD( zu, 2, statistic_regions+1 ) 1272 hom(nzb,2,9,:) = 0.0 1273 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1274 dopr_initial_index(i) = 23 1275 hom(:,2,23,:) = hom(:,2,9,:) 1276 data_output_pr(i) = data_output_pr(i)(2:) 1277 ENDIF 1278 1279 CASE ( 'kh', '#kh' ) 1280 dopr_index(i) = 10 1281 dopr_unit(i) = 'm2/s' 1282 hom(:,2,10,:) = SPREAD( zu, 2, statistic_regions+1 ) 1283 hom(nzb,2,10,:) = 0.0 1284 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1285 dopr_initial_index(i) = 24 1286 hom(:,2,24,:) = hom(:,2,10,:) 1287 data_output_pr(i) = data_output_pr(i)(2:) 1288 ENDIF 1289 1290 CASE ( 'l', '#l' ) 1291 dopr_index(i) = 11 1292 dopr_unit(i) = 'm' 1293 hom(:,2,11,:) = SPREAD( zu, 2, statistic_regions+1 ) 1294 hom(nzb,2,11,:) = 0.0 1295 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1296 dopr_initial_index(i) = 25 1297 hom(:,2,25,:) = hom(:,2,11,:) 1298 data_output_pr(i) = data_output_pr(i)(2:) 1299 ENDIF 1300 1301 CASE ( 'w"u"' ) 1302 dopr_index(i) = 12 1303 dopr_unit(i) = 'm2/s2' 1304 hom(:,2,12,:) = SPREAD( zw, 2, statistic_regions+1 ) 1305 IF ( prandtl_layer ) hom(nzb,2,12,:) = zu(1) 1306 1307 CASE ( 'w*u*' ) 1308 dopr_index(i) = 13 1309 dopr_unit(i) = 'm2/s2' 1310 hom(:,2,13,:) = SPREAD( zw, 2, statistic_regions+1 ) 1311 1312 CASE ( 'w"v"' ) 1313 dopr_index(i) = 14 1314 dopr_unit(i) = 'm2/s2' 1315 hom(:,2,14,:) = SPREAD( zw, 2, statistic_regions+1 ) 1316 IF ( prandtl_layer ) hom(nzb,2,14,:) = zu(1) 1317 1318 CASE ( 'w*v*' ) 1319 dopr_index(i) = 15 1320 dopr_unit(i) = 'm2/s2' 1321 hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 ) 1322 1323 CASE ( 'w"pt"' ) 1324 dopr_index(i) = 16 1325 dopr_unit(i) = 'K m/s' 1326 hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 ) 1327 1328 CASE ( 'w*pt*' ) 1329 dopr_index(i) = 17 1330 dopr_unit(i) = 'K m/s' 1331 hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 ) 1332 1333 CASE ( 'wpt' ) 1334 dopr_index(i) = 18 1335 dopr_unit(i) = 'K m/s' 1336 hom(:,2,18,:) = SPREAD( zw, 2, statistic_regions+1 ) 1337 1338 CASE ( 'wu' ) 1339 dopr_index(i) = 19 1340 dopr_unit(i) = 'm2/s2' 1341 hom(:,2,19,:) = SPREAD( zw, 2, statistic_regions+1 ) 1342 IF ( prandtl_layer ) hom(nzb,2,19,:) = zu(1) 1343 1344 CASE ( 'wv' ) 1345 dopr_index(i) = 20 1346 dopr_unit(i) = 'm2/s2' 1347 hom(:,2,20,:) = SPREAD( zw, 2, statistic_regions+1 ) 1348 IF ( prandtl_layer ) hom(nzb,2,20,:) = zu(1) 1349 1350 CASE ( 'w*pt*BC' ) 1351 dopr_index(i) = 21 1352 dopr_unit(i) = 'K m/s' 1353 hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 ) 1354 1355 CASE ( 'wptBC' ) 1356 dopr_index(i) = 22 1357 dopr_unit(i) = 'K m/s' 1358 hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 ) 1359 1360 CASE ( 'u*2' ) 1361 dopr_index(i) = 30 1362 dopr_unit(i) = 'm2/s2' 1363 hom(:,2,30,:) = SPREAD( zu, 2, statistic_regions+1 ) 1364 1365 CASE ( 'v*2' ) 1366 dopr_index(i) = 31 1367 dopr_unit(i) = 'm2/s2' 1368 hom(:,2,31,:) = SPREAD( zu, 2, statistic_regions+1 ) 1369 1370 CASE ( 'w*2' ) 1371 dopr_index(i) = 32 1372 dopr_unit(i) = 'm2/s2' 1373 hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 ) 1374 1375 CASE ( 'pt*2' ) 1376 dopr_index(i) = 33 1377 dopr_unit(i) = 'K2' 1378 hom(:,2,33,:) = SPREAD( zu, 2, statistic_regions+1 ) 1379 1380 CASE ( 'e*' ) 1381 dopr_index(i) = 34 1382 dopr_unit(i) = 'm2/s2' 1383 hom(:,2,34,:) = SPREAD( zu, 2, statistic_regions+1 ) 1384 1385 CASE ( 'w*2pt*' ) 1386 dopr_index(i) = 35 1387 dopr_unit(i) = 'K m2/s2' 1388 hom(:,2,35,:) = SPREAD( zw, 2, statistic_regions+1 ) 1389 1390 CASE ( 'w*pt*2' ) 1391 dopr_index(i) = 36 1392 dopr_unit(i) = 'K2 m/s' 1393 hom(:,2,36,:) = SPREAD( zw, 2, statistic_regions+1 ) 1394 1395 CASE ( 'w*e*' ) 1396 dopr_index(i) = 37 1397 dopr_unit(i) = 'm3/s3' 1398 hom(:,2,37,:) = SPREAD( zw, 2, statistic_regions+1 ) 1399 1400 CASE ( 'w*3' ) 1401 dopr_index(i) = 38 1402 dopr_unit(i) = 'm3/s3' 1403 hom(:,2,38,:) = SPREAD( zw, 2, statistic_regions+1 ) 1404 1405 CASE ( 'Sw' ) 1406 dopr_index(i) = 39 1407 hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 ) 1408 1409 CASE ( 'q', '#q' ) 1410 IF ( .NOT. cloud_physics ) THEN 1411 IF ( myid == 0 ) THEN 1412 PRINT*, '+++ check_parameters: data_output_pr = ', & 1413 data_output_pr(i), & 1414 ' is not implemented for cloud_physics = FALSE' 1415 ENDIF 1416 CALL local_stop 1417 ELSE 1418 dopr_index(i) = 41 1419 dopr_unit(i) = 'kg/kg' 1420 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 ) 1421 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1422 dopr_initial_index(i) = 26 1423 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 ) 1424 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist 1425 data_output_pr(i) = data_output_pr(i)(2:) 1426 ENDIF 1427 ENDIF 1428 1429 CASE ( 's', '#s' ) 1430 IF ( .NOT. passive_scalar ) THEN 1431 IF ( myid == 0 ) THEN 1432 PRINT*, '+++ check_parameters: data_output_pr = ', & 1433 data_output_pr(i), & 1434 ' is not implemented for passive_scalar = FALSE' 1435 ENDIF 1436 CALL local_stop 1437 ELSE 1438 dopr_index(i) = 41 1439 dopr_unit(i) = 'kg/m3' 1440 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 ) 1441 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1442 dopr_initial_index(i) = 26 1443 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 ) 1444 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist 1445 data_output_pr(i) = data_output_pr(i)(2:) 1446 ENDIF 1447 ENDIF 1448 1449 CASE ( 'qv', '#qv' ) 1450 IF ( .NOT. cloud_physics ) THEN 1451 dopr_index(i) = 41 1452 dopr_unit(i) = 'kg/kg' 1453 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 ) 1454 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1455 dopr_initial_index(i) = 26 1456 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 ) 1457 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist 1458 data_output_pr(i) = data_output_pr(i)(2:) 1459 ENDIF 1460 ELSE 1461 dopr_index(i) = 42 1462 dopr_unit(i) = 'kg/kg' 1463 hom(:,2,42,:) = SPREAD( zu, 2, statistic_regions+1 ) 1464 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1465 dopr_initial_index(i) = 27 1466 hom(:,2,27,:) = SPREAD( zu, 2, statistic_regions+1 ) 1467 hom(nzb,2,27,:) = 0.0 ! weil zu(nzb) negativ ist 1468 data_output_pr(i) = data_output_pr(i)(2:) 1469 ENDIF 1470 ENDIF 1471 1472 CASE ( 'lpt', '#lpt' ) 1473 IF ( .NOT. cloud_physics ) THEN 1474 IF ( myid == 0 ) THEN 1475 PRINT*, '+++ check_parameters: data_output_pr = ', & 1476 data_output_pr(i), & 1477 ' is not implemented for cloud_physics = FALSE' 1478 ENDIF 1479 CALL local_stop 1480 ELSE 1481 dopr_index(i) = 4 1482 dopr_unit(i) = 'K' 1239 1483 hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 ) 1240 1484 IF ( data_output_pr(i)(1:1) == '#' ) THEN … … 1244 1488 data_output_pr(i) = data_output_pr(i)(2:) 1245 1489 ENDIF 1246 ELSE1247 dopr_index(i) = 431248 hom(:,2,43,:) = SPREAD( zu, 2, statistic_regions+1 )1249 IF ( data_output_pr(i)(1:1) == '#' ) THEN1250 dopr_initial_index(i) = 281251 hom(:,2,28,:) = SPREAD( zu, 2, statistic_regions+1 )1252 hom(nzb,2,28,:) = 0.0 ! weil zu(nzb) negativ ist1253 data_output_pr(i) = data_output_pr(i)(2:)1254 ENDIF1255 ENDIF1256 1257 CASE ( 'e' )1258 dopr_index(i) = 81259 hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 )1260 hom(nzb,2,8,:) = 0.01261 1262 CASE ( 'km', '#km' )1263 dopr_index(i) = 91264 hom(:,2,9,:) = SPREAD( zu, 2, statistic_regions+1 )1265 hom(nzb,2,9,:) = 0.01266 IF ( data_output_pr(i)(1:1) == '#' ) THEN1267 dopr_initial_index(i) = 231268 hom(:,2,23,:) = hom(:,2,9,:)1269 data_output_pr(i) = data_output_pr(i)(2:)1270 ENDIF1271 1272 CASE ( 'kh', '#kh' )1273 dopr_index(i) = 101274 hom(:,2,10,:) = SPREAD( zu, 2, statistic_regions+1 )1275 hom(nzb,2,10,:) = 0.01276 IF ( data_output_pr(i)(1:1) == '#' ) THEN1277 dopr_initial_index(i) = 241278 hom(:,2,24,:) = hom(:,2,10,:)1279 data_output_pr(i) = data_output_pr(i)(2:)1280 ENDIF1281 1282 CASE ( 'l', '#l' )1283 dopr_index(i) = 111284 hom(:,2,11,:) = SPREAD( zu, 2, statistic_regions+1 )1285 hom(nzb,2,11,:) = 0.01286 IF ( data_output_pr(i)(1:1) == '#' ) THEN1287 dopr_initial_index(i) = 251288 hom(:,2,25,:) = hom(:,2,11,:)1289 data_output_pr(i) = data_output_pr(i)(2:)1290 ENDIF1291 1292 CASE ( 'w"u"' )1293 dopr_index(i) = 121294 hom(:,2,12,:) = SPREAD( zw, 2, statistic_regions+1 )1295 IF ( prandtl_layer ) hom(nzb,2,12,:) = zu(1)1296 1297 CASE ( 'w*u*' )1298 dopr_index(i) = 131299 hom(:,2,13,:) = SPREAD( zw, 2, statistic_regions+1 )1300 1301 CASE ( 'w"v"' )1302 dopr_index(i) = 141303 hom(:,2,14,:) = SPREAD( zw, 2, statistic_regions+1 )1304 IF ( prandtl_layer ) hom(nzb,2,14,:) = zu(1)1305 1306 CASE ( 'w*v*' )1307 dopr_index(i) = 151308 hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 )1309 1310 CASE ( 'w"pt"' )1311 dopr_index(i) = 161312 hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 )1313 1314 CASE ( 'w*pt*' )1315 dopr_index(i) = 171316 hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 )1317 1318 CASE ( 'wpt' )1319 dopr_index(i) = 181320 hom(:,2,18,:) = SPREAD( zw, 2, statistic_regions+1 )1321 1322 CASE ( 'wu' )1323 dopr_index(i) = 191324 hom(:,2,19,:) = SPREAD( zw, 2, statistic_regions+1 )1325 IF ( prandtl_layer ) hom(nzb,2,19,:) = zu(1)1326 1327 CASE ( 'wv' )1328 dopr_index(i) = 201329 hom(:,2,20,:) = SPREAD( zw, 2, statistic_regions+1 )1330 IF ( prandtl_layer ) hom(nzb,2,20,:) = zu(1)1331 1332 CASE ( 'w*pt*BC' )1333 dopr_index(i) = 211334 hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 )1335 1336 CASE ( 'wptBC' )1337 dopr_index(i) = 221338 hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 )1339 1340 CASE ( 'u*2' )1341 dopr_index(i) = 301342 hom(:,2,30,:) = SPREAD( zu, 2, statistic_regions+1 )1343 1344 CASE ( 'v*2' )1345 dopr_index(i) = 311346 hom(:,2,31,:) = SPREAD( zu, 2, statistic_regions+1 )1347 1348 CASE ( 'w*2' )1349 dopr_index(i) = 321350 hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 )1351 1352 CASE ( 'pt*2' )1353 dopr_index(i) = 331354 hom(:,2,33,:) = SPREAD( zu, 2, statistic_regions+1 )1355 1356 CASE ( 'e*' )1357 dopr_index(i) = 341358 hom(:,2,34,:) = SPREAD( zu, 2, statistic_regions+1 )1359 1360 CASE ( 'w*2pt*' )1361 dopr_index(i) = 351362 hom(:,2,35,:) = SPREAD( zw, 2, statistic_regions+1 )1363 1364 CASE ( 'w*pt*2' )1365 dopr_index(i) = 361366 hom(:,2,36,:) = SPREAD( zw, 2, statistic_regions+1 )1367 1368 CASE ( 'w*e*' )1369 dopr_index(i) = 371370 hom(:,2,37,:) = SPREAD( zw, 2, statistic_regions+1 )1371 1372 CASE ( 'w*3' )1373 dopr_index(i) = 381374 hom(:,2,38,:) = SPREAD( zw, 2, statistic_regions+1 )1375 1376 CASE ( 'Sw' )1377 dopr_index(i) = 391378 hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )1379 1380 CASE ( 'q', '#q' )1381 IF ( .NOT. cloud_physics ) THEN1382 IF ( myid == 0 ) THEN1383 PRINT*, '+++ check_parameters: data_output_pr = ', &1384 data_output_pr(i), &1385 ' is not implemented for cloud_physics = FALSE'1386 ENDIF1387 CALL local_stop1388 ELSE1389 dopr_index(i) = 411390 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )1391 IF ( data_output_pr(i)(1:1) == '#' ) THEN1392 dopr_initial_index(i) = 261393 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )1394 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist1395 data_output_pr(i) = data_output_pr(i)(2:)1396 ENDIF1397 ENDIF1398 1399 CASE ( 's', '#s' )1400 IF ( .NOT. passive_scalar ) THEN1401 IF ( myid == 0 ) THEN1402 PRINT*, '+++ check_parameters: data_output_pr = ', &1403 data_output_pr(i), &1404 ' is not implemented for passive_scalar = FALSE'1405 ENDIF1406 CALL local_stop1407 ELSE1408 dopr_index(i) = 411409 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )1410 IF ( data_output_pr(i)(1:1) == '#' ) THEN1411 dopr_initial_index(i) = 261412 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )1413 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist1414 data_output_pr(i) = data_output_pr(i)(2:)1415 ENDIF1416 ENDIF1417 1418 CASE ( 'qv', '#qv' )1419 IF ( .NOT. cloud_physics ) THEN1420 dopr_index(i) = 411421 hom(:,2,41,:) = SPREAD( zu, 2, statistic_regions+1 )1422 IF ( data_output_pr(i)(1:1) == '#' ) THEN1423 dopr_initial_index(i) = 261424 hom(:,2,26,:) = SPREAD( zu, 2, statistic_regions+1 )1425 hom(nzb,2,26,:) = 0.0 ! weil zu(nzb) negativ ist1426 data_output_pr(i) = data_output_pr(i)(2:)1427 ENDIF1428 ELSE1429 dopr_index(i) = 421430 hom(:,2,42,:) = SPREAD( zu, 2, statistic_regions+1 )1431 IF ( data_output_pr(i)(1:1) == '#' ) THEN1432 dopr_initial_index(i) = 271433 hom(:,2,27,:) = SPREAD( zu, 2, statistic_regions+1 )1434 hom(nzb,2,27,:) = 0.0 ! weil zu(nzb) negativ ist1435 data_output_pr(i) = data_output_pr(i)(2:)1436 ENDIF1437 ENDIF1438 1439 CASE ( 'lpt', '#lpt' )1440 IF ( .NOT. cloud_physics ) THEN1441 IF ( myid == 0 ) THEN1442 PRINT*, '+++ check_parameters: data_output_pr = ', &1443 data_output_pr(i), &1444 ' is not implemented for cloud_physics = FALSE'1445 ENDIF1446 CALL local_stop1447 ELSE1448 dopr_index(i) = 41449 hom(:,2,4,:) = SPREAD( zu, 2, statistic_regions+1 )1450 IF ( data_output_pr(i)(1:1) == '#' ) THEN1451 dopr_initial_index(i) = 71452 hom(:,2,7,:) = SPREAD( zu, 2, statistic_regions+1 )1453 hom(nzb,2,7,:) = 0.0 ! weil zu(nzb) negativ ist1454 data_output_pr(i) = data_output_pr(i)(2:)1455 ENDIF1456 1490 ENDIF 1457 1491 1458 1492 CASE ( 'vpt', '#vpt' ) 1459 1493 dopr_index(i) = 44 1460 hom(:,2,44,:) = SPREAD( zu, 2, statistic_regions+1 ) 1494 dopr_unit(i) = 'K' 1495 hom(:,2,44,:) = SPREAD( zu, 2, statistic_regions+1 ) 1461 1496 IF ( data_output_pr(i)(1:1) == '#' ) THEN 1462 1497 dopr_initial_index(i) = 29 … … 1468 1503 CASE ( 'w"vpt"' ) 1469 1504 dopr_index(i) = 45 1505 dopr_unit(i) = 'K m/s' 1470 1506 hom(:,2,45,:) = SPREAD( zw, 2, statistic_regions+1 ) 1471 1507 1472 1508 CASE ( 'w*vpt*' ) 1473 1509 dopr_index(i) = 46 1510 dopr_unit(i) = 'K m/s' 1474 1511 hom(:,2,46,:) = SPREAD( zw, 2, statistic_regions+1 ) 1475 1512 1476 1513 CASE ( 'wvpt' ) 1477 1514 dopr_index(i) = 47 1515 dopr_unit(i) = 'K m/s' 1478 1516 hom(:,2,47,:) = SPREAD( zw, 2, statistic_regions+1 ) 1479 1517 … … 1488 1526 ELSE 1489 1527 dopr_index(i) = 48 1528 dopr_unit(i) = 'kg/kg m/s' 1490 1529 hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 ) 1491 1530 ENDIF … … 1501 1540 ELSE 1502 1541 dopr_index(i) = 49 1542 dopr_unit(i) = 'kg/kg m/s' 1503 1543 hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 ) 1504 1544 ENDIF … … 1514 1554 ELSE 1515 1555 dopr_index(i) = 50 1556 dopr_unit(i) = 'kg/kg m/s' 1516 1557 hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 ) 1517 1558 ENDIF … … 1527 1568 ELSE 1528 1569 dopr_index(i) = 48 1570 dopr_unit(i) = 'kg/m3 m/s' 1529 1571 hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 ) 1530 1572 ENDIF … … 1540 1582 ELSE 1541 1583 dopr_index(i) = 49 1584 dopr_unit(i) = 'kg/m3 m/s' 1542 1585 hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 ) 1543 1586 ENDIF … … 1553 1596 ELSE 1554 1597 dopr_index(i) = 50 1598 dopr_unit(i) = 'kg/m3 m/s' 1555 1599 hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 ) 1556 1600 ENDIF … … 1560 1604 THEN 1561 1605 dopr_index(i) = 48 1606 dopr_unit(i) = 'kg/kg m/s' 1562 1607 hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 ) 1563 1608 ELSEIF( humidity .AND. cloud_physics ) THEN 1564 1609 dopr_index(i) = 51 1610 dopr_unit(i) = 'kg/kg m/s' 1565 1611 hom(:,2,51,:) = SPREAD( zw, 2, statistic_regions+1 ) 1566 1612 ELSE … … 1578 1624 THEN 1579 1625 dopr_index(i) = 49 1626 dopr_unit(i) = 'kg/kg m/s' 1580 1627 hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 ) 1581 1628 ELSEIF( humidity .AND. cloud_physics ) THEN 1582 1629 dopr_index(i) = 52 1630 dopr_unit(i) = 'kg/kg m/s' 1583 1631 hom(:,2,52,:) = SPREAD( zw, 2, statistic_regions+1 ) 1584 1632 ELSE … … 1596 1644 THEN 1597 1645 dopr_index(i) = 50 1646 dopr_unit(i) = 'kg/kg m/s' 1598 1647 hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 ) 1599 1648 ELSEIF( humidity .AND. cloud_physics ) THEN 1600 1649 dopr_index(i) = 53 1650 dopr_unit(i) = 'kg/kg m/s' 1601 1651 hom(:,2,53,:) = SPREAD( zw, 2, statistic_regions+1 ) 1602 1652 ELSE … … 1620 1670 ELSE 1621 1671 dopr_index(i) = 54 1672 dopr_unit(i) = 'kg/kg' 1622 1673 hom(:,2,54,:) = SPREAD( zu, 2, statistic_regions+1 ) 1623 1674 ENDIF … … 1625 1676 CASE ( 'w*u*u*/dz' ) 1626 1677 dopr_index(i) = 55 1678 dopr_unit(i) = 'm2/s3' 1627 1679 hom(:,2,55,:) = SPREAD( zu, 2, statistic_regions+1 ) 1628 1680 1629 1681 CASE ( 'w*p*/dz' ) 1630 1682 dopr_index(i) = 56 1683 dopr_unit(i) = 'm2/s3' 1631 1684 hom(:,2,56,:) = SPREAD( zu, 2, statistic_regions+1 ) 1632 1685 1633 1686 CASE ( 'w"e/dz' ) 1634 1687 dopr_index(i) = 57 1688 dopr_unit(i) = 'm2/s3' 1635 1689 hom(:,2,57,:) = SPREAD( zu, 2, statistic_regions+1 ) 1636 1690 1637 1691 CASE ( 'u"pt"' ) 1638 1692 dopr_index(i) = 58 1693 dopr_unit(i) = 'K m/s' 1639 1694 hom(:,2,58,:) = SPREAD( zu, 2, statistic_regions+1 ) 1640 1695 1641 1696 CASE ( 'u*pt*' ) 1642 1697 dopr_index(i) = 59 1698 dopr_unit(i) = 'K m/s' 1643 1699 hom(:,2,59,:) = SPREAD( zu, 2, statistic_regions+1 ) 1644 1700 1645 1701 CASE ( 'upt_t' ) 1646 1702 dopr_index(i) = 60 1703 dopr_unit(i) = 'K m/s' 1647 1704 hom(:,2,60,:) = SPREAD( zu, 2, statistic_regions+1 ) 1648 1705 1649 1706 CASE ( 'v"pt"' ) 1650 1707 dopr_index(i) = 61 1708 dopr_unit(i) = 'K m/s' 1651 1709 hom(:,2,61,:) = SPREAD( zu, 2, statistic_regions+1 ) 1652 1710 1653 1711 CASE ( 'v*pt*' ) 1654 1712 dopr_index(i) = 62 1713 dopr_unit(i) = 'K m/s' 1655 1714 hom(:,2,62,:) = SPREAD( zu, 2, statistic_regions+1 ) 1656 1715 1657 1716 CASE ( 'vpt_t' ) 1658 1717 dopr_index(i) = 63 1718 dopr_unit(i) = 'K m/s' 1659 1719 hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 ) 1660 1720 1661 1721 1662 1722 CASE DEFAULT 1663 IF ( myid == 0 ) THEN 1664 PRINT*, '+++ check_parameters: unknown output profile: ', & 1665 'data_output_pr = ', data_output_pr(i) 1666 ENDIF 1667 CALL local_stop 1723 1724 CALL user_check_data_output_pr( data_output_pr(i), i, unit ) 1725 1726 IF ( unit == 'illegal' ) THEN 1727 IF ( myid == 0 ) THEN 1728 IF ( data_output_pr_user(1) /= ' ' ) THEN 1729 PRINT*, '+++ check_parameters: illegal value for data_',& 1730 'output_pr or data_output_pr_user: "', & 1731 TRIM( data_output_pr(i) ), '"' 1732 ELSE 1733 PRINT*, '+++ check_parameters: illegal value for data_',& 1734 'output_pr: "', TRIM( data_output_pr(i) ),'"' 1735 ENDIF 1736 ENDIF 1737 CALL local_stop 1738 ENDIF 1668 1739 1669 1740 END SELECT -
palm/trunk/SOURCE/data_output_profiles.f90
r4 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! var_hom renamed pr_palm 7 7 ! 8 8 ! Former revisions: … … 111 111 112 112 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2 113 (/ hom_sum(nzb+8, var_hom,normalizing_region)**2 /), &113 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), & 114 114 start = (/ 1 /), count = (/ 1 /) ) 115 115 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 331 ) 116 116 117 117 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2 118 (/ hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &118 (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 119 119 start = (/ 1 /), count = (/ 1 /) ) 120 120 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 332 ) 121 121 122 122 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3 123 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 /), &123 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), & 124 124 start = (/ 1 /), count = (/ 1 /) ) 125 125 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 333 ) 126 126 127 127 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw 128 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 * &129 hom_sum(nzb+3, var_hom,normalizing_region) /), &128 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 * & 129 hom_sum(nzb+3,pr_palm,normalizing_region) /), & 130 130 start = (/ 1 /), count = (/ 1 /) ) 131 131 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 334 ) 132 132 133 133 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2 134 (/ hom_sum(nzb+8, var_hom,normalizing_region) * &135 hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &134 (/ hom_sum(nzb+8,pr_palm,normalizing_region) * & 135 hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 136 136 start = (/ 1 /), count = (/ 1 /) ) 137 137 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 335 ) 138 138 139 139 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i 140 (/ hom_sum(nzb+6, var_hom,normalizing_region) /), &140 (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), & 141 141 start = (/ 1 /), count = (/ 1 /) ) 142 142 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 336 ) … … 237 237 CASE ( 'ts2' ) 238 238 cross_normx_factor(k,j) = & 239 ( hom_sum(nzb+3, var_hom,normalizing_region) )**2239 ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2 240 240 CASE ( 'wpt0' ) 241 241 cross_normx_factor(k,j) = & … … 243 243 CASE ( 'wsts2' ) 244 244 cross_normx_factor(k,j) = & 245 hom_sum(nzb+8, var_hom,normalizing_region) &246 * ( hom_sum(nzb+3, var_hom,normalizing_region) )**2245 hom_sum(nzb+8,pr_palm,normalizing_region) & 246 * ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2 247 247 CASE ( 'ws2' ) 248 248 cross_normx_factor(k,j) = & 249 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2249 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 250 250 CASE ( 'ws2ts' ) 251 251 cross_normx_factor(k,j) = & 252 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2 &253 * hom_sum(nzb+3, var_hom,normalizing_region)252 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 & 253 * hom_sum(nzb+3,pr_palm,normalizing_region) 254 254 CASE ( 'ws3' ) 255 255 cross_normx_factor(k,j) = & 256 ( hom_sum(nzb+8, var_hom,normalizing_region) )**3256 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3 257 257 258 258 END SELECT … … 262 262 CASE ( 'z_i' ) 263 263 cross_normy_factor(k,j) = & 264 hom_sum(nzb+6, var_hom,normalizing_region)264 hom_sum(nzb+6,pr_palm,normalizing_region) 265 265 266 266 END SELECT … … 351 351 352 352 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2 353 (/ hom_sum(nzb+8, var_hom,normalizing_region)**2 /), &353 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), & 354 354 start = (/ dopr_time_count /), & 355 355 count = (/ 1 /) ) … … 357 357 358 358 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2 359 (/ hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &359 (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 360 360 start = (/ dopr_time_count /), & 361 361 count = (/ 1 /) ) … … 363 363 364 364 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3 365 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 /), &365 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), & 366 366 start = (/ dopr_time_count /), & 367 367 count = (/ 1 /) ) … … 369 369 370 370 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw 371 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 * &372 hom_sum(nzb+3, var_hom,normalizing_region) /), &371 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 * & 372 hom_sum(nzb+3,pr_palm,normalizing_region) /), & 373 373 start = (/ dopr_time_count /), & 374 374 count = (/ 1 /) ) … … 376 376 377 377 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2 378 (/ hom_sum(nzb+8, var_hom,normalizing_region) * &379 hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &378 (/ hom_sum(nzb+8,pr_palm,normalizing_region) * & 379 hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 380 380 start = (/ dopr_time_count /), & 381 381 count = (/ 1 /) ) … … 383 383 384 384 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i 385 (/ hom_sum(nzb+6, var_hom,normalizing_region) /), &385 (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), & 386 386 start = (/ dopr_time_count /), & 387 387 count = (/ 1 /) ) … … 482 482 CASE ( 'tsw2' ) 483 483 cross_normx_factor(k,j) = & 484 ( hom_sum(nzb+11, var_hom,normalizing_region) )**2484 ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2 485 485 CASE ( 'wpt0' ) 486 486 cross_normx_factor(k,j) = & … … 488 488 CASE ( 'wstsw2' ) 489 489 cross_normx_factor(k,j) = & 490 hom_sum(nzb+8, var_hom,normalizing_region) &491 * ( hom_sum(nzb+11, var_hom,normalizing_region) )**2490 hom_sum(nzb+8,pr_palm,normalizing_region) & 491 * ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2 492 492 CASE ( 'ws2' ) 493 493 cross_normx_factor(k,j) = & 494 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2494 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 495 495 CASE ( 'ws2tsw' ) 496 496 cross_normx_factor(k,j) = & 497 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2&498 * hom_sum(nzb+11, var_hom,normalizing_region)497 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2& 498 * hom_sum(nzb+11,pr_palm,normalizing_region) 499 499 CASE ( 'ws3' ) 500 500 cross_normx_factor(k,j) = & 501 ( hom_sum(nzb+8, var_hom,normalizing_region) )**3501 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3 502 502 503 503 END SELECT … … 506 506 CASE ( 'z_i' ) 507 507 cross_normy_factor(k,j) = & 508 hom_sum(nzb+6, var_hom,normalizing_region)508 hom_sum(nzb+6,pr_palm,normalizing_region) 509 509 510 510 END SELECT -
palm/trunk/SOURCE/flow_statistics.f90
r83 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Two more arguments added to user_statistics, which is now also called for 7 ! user-defined profiles, 8 ! var_hom and var_sum renamed pr_palm 7 9 ! 8 10 ! Former revisions: … … 87 89 !-- WARNING: next line still has to be adjusted for OpenMP 88 90 sums_l(:,21,0) = sums_wsts_bc_l(:,sr) ! heat flux from advec_s_bc 89 sums_l(nzb+9, var_sum,0) = sums_divold_l(sr) ! old divergence from pres90 sums_l(nzb+10, var_sum,0) = sums_divnew_l(sr) ! new divergence from pres91 sums_l(nzb+9,pr_palm,0) = sums_divold_l(sr) ! old divergence from pres 92 sums_l(nzb+10,pr_palm,0) = sums_divnew_l(sr) ! new divergence from pres 91 93 !-- WARNING: next four lines still may have to be adjusted for OpenMP 92 sums_l(nzb:nzb+2, var_sum-1,0) = sums_up_fraction_l(1,1:3,sr)! upstream93 sums_l(nzb+3:nzb+5, var_sum-1,0) = sums_up_fraction_l(2,1:3,sr)! parts94 sums_l(nzb+6:nzb+8, var_sum-1,0) = sums_up_fraction_l(3,1:3,sr)! from95 sums_l(nzb+9:nzb+11, var_sum-1,0) = sums_up_fraction_l(4,1:3,sr)! spline94 sums_l(nzb:nzb+2,pr_palm-1,0) = sums_up_fraction_l(1,1:3,sr)! upstream 95 sums_l(nzb+3:nzb+5,pr_palm-1,0) = sums_up_fraction_l(2,1:3,sr)! parts 96 sums_l(nzb+6:nzb+8,pr_palm-1,0) = sums_up_fraction_l(3,1:3,sr)! from 97 sums_l(nzb+9:nzb+11,pr_palm-1,0) = sums_up_fraction_l(4,1:3,sr)! spline 96 98 97 99 ! … … 316 318 !-- quantities is seperated from the previous loop in order to 317 319 !-- allow vectorization of that loop. 318 sums_l(nzb+4, var_sum,tn) = sums_l(nzb+4,var_sum,tn) + sums_l_etot319 sums_l(nzb+5, var_sum,tn) = sums_l(nzb+5,var_sum,tn) + sums_l_eper320 sums_l(nzb+4,pr_palm,tn) = sums_l(nzb+4,pr_palm,tn) + sums_l_etot 321 sums_l(nzb+5,pr_palm,tn) = sums_l(nzb+5,pr_palm,tn) + sums_l_eper 320 322 ! 321 323 !-- 2D-arrays (being collected in the last column of sums_l) 322 sums_l(nzb, var_sum,tn) = sums_l(nzb,var_sum,tn) + &324 sums_l(nzb,pr_palm,tn) = sums_l(nzb,pr_palm,tn) + & 323 325 us(j,i) * rmask(j,i,sr) 324 sums_l(nzb+1, var_sum,tn) = sums_l(nzb+1,var_sum,tn) + &326 sums_l(nzb+1,pr_palm,tn) = sums_l(nzb+1,pr_palm,tn) + & 325 327 usws(j,i) * rmask(j,i,sr) 326 sums_l(nzb+2, var_sum,tn) = sums_l(nzb+2,var_sum,tn) + &328 sums_l(nzb+2,pr_palm,tn) = sums_l(nzb+2,pr_palm,tn) + & 327 329 vsws(j,i) * rmask(j,i,sr) 328 sums_l(nzb+3, var_sum,tn) = sums_l(nzb+3,var_sum,tn) + &330 sums_l(nzb+3,pr_palm,tn) = sums_l(nzb+3,pr_palm,tn) + & 329 331 ts(j,i) * rmask(j,i,sr) 330 332 ENDDO … … 652 654 653 655 ENDIF 656 657 ! 658 !-- Calculate the user-defined profiles 659 CALL user_statistics( 'profiles', sr, tn ) 654 660 !$OMP END PARALLEL 655 661 … … 660 666 sums_l(:,3,0) = sums_l(:,3,0) + sums_l(:,3,i) 661 667 sums_l(:,4:40,0) = sums_l(:,4:40,0) + sums_l(:,4:40,i) 662 sums_l(:,45:var_sum,0) = sums_l(:,45:var_sum,0) + & 663 sums_l(:,45:var_sum,i) 668 sums_l(:,45:pr_palm,0) = sums_l(:,45:pr_palm,0) + & 669 sums_l(:,45:pr_palm,i) 670 IF ( max_pr_user > 0 ) THEN 671 sums_l(:,pr_palm+1:pr_palm+max_pr_user,0) = & 672 sums_l(:,pr_palm+1:pr_palm+max_pr_user,0) + & 673 sums_l(:,pr_palm+1:pr_palm+max_pr_user,i) 674 ENDIF 664 675 ENDDO 665 676 ENDIF … … 679 690 !-- Profiles: 680 691 DO k = nzb, nzt+1 681 sums(k,: var_sum-2) = sums(k,:var_sum-2) / ngp_2dh_outer(k,sr)692 sums(k,:pr_palm-2) = sums(k,:pr_palm-2) / ngp_2dh_outer(k,sr) 682 693 ENDDO 683 694 !-- Upstream-parts 684 sums(nzb:nzb+11, var_sum-1) = sums(nzb:nzb+11,var_sum-1) / ngp_3d(sr)695 sums(nzb:nzb+11,pr_palm-1) = sums(nzb:nzb+11,pr_palm-1) / ngp_3d(sr) 685 696 !-- u* and so on 686 !-- As sums(nzb:nzb+3, var_sum) are full 2D arrays (us, usws, vsws, ts) whose697 !-- As sums(nzb:nzb+3,pr_palm) are full 2D arrays (us, usws, vsws, ts) whose 687 698 !-- size is always ( nx + 1 ) * ( ny + 1 ), defined at the first grid layer 688 699 !-- above the topography, they are being divided by ngp_2dh(sr) 689 sums(nzb:nzb+3, var_sum) = sums(nzb:nzb+3,var_sum) / &700 sums(nzb:nzb+3,pr_palm) = sums(nzb:nzb+3,pr_palm) / & 690 701 ngp_2dh(sr) 691 702 !-- eges, e* 692 sums(nzb+4:nzb+5, var_sum) = sums(nzb+4:nzb+5,var_sum) / &703 sums(nzb+4:nzb+5,pr_palm) = sums(nzb+4:nzb+5,pr_palm) / & 693 704 ngp_3d_inner(sr) 694 705 !-- Old and new divergence 695 sums(nzb+9:nzb+10, var_sum) = sums(nzb+9:nzb+10,var_sum) / &706 sums(nzb+9:nzb+10,pr_palm) = sums(nzb+9:nzb+10,pr_palm) / & 696 707 ngp_3d_inner(sr) 708 709 !-- User-defined profiles 710 IF ( max_pr_user > 0 ) THEN 711 DO k = nzb, nzt+1 712 sums(k,pr_palm+1:pr_palm+max_pr_user) = & 713 sums(k,pr_palm+1:pr_palm+max_pr_user) / & 714 ngp_2dh_outer(k,sr) 715 ENDDO 716 ENDIF 697 717 698 718 ! … … 748 768 hom(:,1,63,sr) = sums(:,61) + sums(:,62) ! vpt_t 749 769 750 hom(:,1, var_hom-1,sr) = sums(:,var_sum-1)770 hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1) 751 771 ! upstream-parts u_x, u_y, u_z, v_x, 752 772 ! v_y, usw. (in last but one profile) 753 hom(:,1, var_hom,sr) = sums(:,var_sum)773 hom(:,1,pr_palm,sr) = sums(:,pr_palm) 754 774 ! u*, w'u', w'v', t* (in last profile) 775 776 IF ( max_pr_user > 0 ) THEN ! user-defined profiles 777 hom(:,1,pr_palm+1:pr_palm+max_pr_user,sr) = & 778 sums(:,pr_palm+1:pr_palm+max_pr_user) 779 ENDIF 755 780 756 781 ! … … 792 817 ENDDO 793 818 794 hom(nzb+6,1, var_hom,sr) = z_i(1)795 hom(nzb+7,1, var_hom,sr) = z_i(2)819 hom(nzb+6,1,pr_palm,sr) = z_i(1) 820 hom(nzb+7,1,pr_palm,sr) = z_i(2) 796 821 797 822 ! … … 800 825 !-- The horizontal average at nzb+1 is input for the average temperature. 801 826 IF ( hom(nzb,1,18,sr) > 0.0 .AND. z_i(1) /= 0.0 ) THEN 802 hom(nzb+8,1, var_hom,sr) = ( g / hom(nzb+1,1,4,sr) * &827 hom(nzb+8,1,pr_palm,sr) = ( g / hom(nzb+1,1,4,sr) * & 803 828 hom(nzb,1,18,sr) * z_i(1) )**0.333333333 804 829 !-- so far this only works if Prandtl layer is used 805 hom(nzb+11,1, var_hom,sr) = hom(nzb,1,16,sr) / hom(nzb+8,1,var_hom,sr)830 hom(nzb+11,1,pr_palm,sr) = hom(nzb,1,16,sr) / hom(nzb+8,1,pr_palm,sr) 806 831 ELSE 807 hom(nzb+8,1, var_hom,sr) = 0.0808 hom(nzb+11,1, var_hom,sr) = 0.0832 hom(nzb+8,1,pr_palm,sr) = 0.0 833 hom(nzb+11,1,pr_palm,sr) = 0.0 809 834 ENDIF 810 835 811 836 ! 812 837 !-- Collect the time series quantities 813 ts_value(1,sr) = hom(nzb+4,1, var_hom,sr) ! E814 ts_value(2,sr) = hom(nzb+5,1, var_hom,sr) ! E*838 ts_value(1,sr) = hom(nzb+4,1,pr_palm,sr) ! E 839 ts_value(2,sr) = hom(nzb+5,1,pr_palm,sr) ! E* 815 840 ts_value(3,sr) = dt_3d 816 ts_value(4,sr) = hom(nzb,1, var_hom,sr) ! u*817 ts_value(5,sr) = hom(nzb+3,1, var_hom,sr) ! th*841 ts_value(4,sr) = hom(nzb,1,pr_palm,sr) ! u* 842 ts_value(5,sr) = hom(nzb+3,1,pr_palm,sr) ! th* 818 843 ts_value(6,sr) = u_max 819 844 ts_value(7,sr) = v_max 820 845 ts_value(8,sr) = w_max 821 ts_value(9,sr) = hom(nzb+10,1, var_sum,sr) ! new divergence822 ts_value(10,sr) = hom(nzb+9,1, var_hom,sr) ! old Divergence823 ts_value(11,sr) = hom(nzb+6,1, var_hom,sr) ! z_i(1)824 ts_value(12,sr) = hom(nzb+7,1, var_hom,sr) ! z_i(2)825 ts_value(13,sr) = hom(nzb+8,1, var_hom,sr) ! w*846 ts_value(9,sr) = hom(nzb+10,1,pr_palm,sr) ! new divergence 847 ts_value(10,sr) = hom(nzb+9,1,pr_palm,sr) ! old Divergence 848 ts_value(11,sr) = hom(nzb+6,1,pr_palm,sr) ! z_i(1) 849 ts_value(12,sr) = hom(nzb+7,1,pr_palm,sr) ! z_i(2) 850 ts_value(13,sr) = hom(nzb+8,1,pr_palm,sr) ! w* 826 851 ts_value(14,sr) = hom(nzb,1,16,sr) ! w'pt' at k=0 827 852 ts_value(15,sr) = hom(nzb+1,1,16,sr) ! w'pt' at k=1 … … 829 854 ts_value(17,sr) = hom(nzb,1,4,sr) ! pt(0) 830 855 ts_value(18,sr) = hom(nzb+1,1,4,sr) ! pt(zp) 831 ts_value(19,sr) = hom(nzb+9,1, var_hom-1,sr) ! splptx832 ts_value(20,sr) = hom(nzb+10,1, var_hom-1,sr) ! splpty833 ts_value(21,sr) = hom(nzb+11,1, var_hom-1,sr) ! splptz856 ts_value(19,sr) = hom(nzb+9,1,pr_palm-1,sr) ! splptx 857 ts_value(20,sr) = hom(nzb+10,1,pr_palm-1,sr) ! splpty 858 ts_value(21,sr) = hom(nzb+11,1,pr_palm-1,sr) ! splptz 834 859 IF ( ts_value(5,sr) /= 0.0 ) THEN 835 860 ts_value(22,sr) = ts_value(4,sr)**2 / & … … 841 866 ! 842 867 !-- Calculate additional statistics provided by the user interface 843 CALL user_statistics( sr)868 CALL user_statistics( 'time_series', sr, 0 ) 844 869 845 870 ENDDO ! loop of the subregions -
palm/trunk/SOURCE/header.f90
r83 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Bugfix: output of use_upstream_for_tke 7 7 ! 8 8 ! Former revisions: … … 224 224 WRITE ( io, 122 ) timestep_scheme 225 225 ENDIF 226 IF ( use_upstream_for_tke ) WRITE ( io, 143 ) 226 227 IF ( rayleigh_damping_factor /= 0.0 ) THEN 227 228 WRITE ( io, 123 ) rayleigh_damping_height, rayleigh_damping_factor … … 1195 1196 142 FORMAT (' perturbation pressure is calculated at every Runge-Kutta ', & 1196 1197 'step') 1198 143 FORMAT (' Euler/upstream scheme is used for the SGS turbulent ', & 1199 'kinetic energy') 1197 1200 150 FORMAT (' --> Volume flow at the right and north boundary will be ', & 1198 1201 'conserved') -
palm/trunk/SOURCE/init_3d_model.f90
r77 r87 7 7 ! Actual revisions: 8 8 ! ----------------- 9 ! 9 ! var_hom and var_sum renamed pr_palm 10 10 ! 11 11 ! Former revisions: … … 81 81 sums_divold_l(0:statistic_regions) ) 82 82 ALLOCATE( rdf(nzb+1:nzt) ) 83 ALLOCATE( hom_sum(nzb:nzt+1, var_hom,0:statistic_regions),&83 ALLOCATE( hom_sum(nzb:nzt+1,pr_palm+max_pr_user,0:statistic_regions), & 84 84 ngp_2dh_outer(nzb:nzt+1,0:statistic_regions), & 85 85 ngp_2dh_outer_l(nzb:nzt+1,0:statistic_regions), & 86 86 rmask(nys-1:nyn+1,nxl-1:nxr+1,0:statistic_regions), & 87 sums(nzb:nzt+1, var_sum),&88 sums_l(nzb:nzt+1, var_sum,0:threads_per_task-1),&87 sums(nzb:nzt+1,pr_palm+max_pr_user), & 88 sums_l(nzb:nzt+1,pr_palm+max_pr_user,0:threads_per_task-1), & 89 89 sums_l_l(nzb:nzt+1,0:statistic_regions,0:threads_per_task-1), & 90 90 sums_up_fraction_l(10,3,0:statistic_regions), & … … 931 931 ngp_3d_inner = 0 932 932 ngp_3d = 0 933 ngp_sums = ( nz + 2 ) * var_sum933 ngp_sums = ( nz + 2 ) * ( pr_palm + max_pr_user ) 934 934 935 935 DO sr = 0, statistic_regions -
palm/trunk/SOURCE/modules.f90
r83 r87 5 5 ! Actual revisions: 6 6 ! ----------------- 7 ! 7 ! +data_output_pr_user, max_pr_user, num_user_profiles, size of data_output_pr, 8 ! dopr_index, dopr_initial_index and dopr_unit enlarged, 9 ! var_hom and var_sum renamed pr_palm 8 10 ! 9 11 ! Former revisions: … … 250 252 CHARACTER (LEN=10), DIMENSION(10) :: data_output_format = ' ' 251 253 CHARACTER (LEN=10), DIMENSION(100) :: data_output = ' ', & 252 data_output_pr = ' ', &253 254 data_output_user = ' ', doav = ' ' 255 CHARACTER (LEN=10), DIMENSION(300) :: data_output_pr = ' ' 256 CHARACTER (LEN=10), DIMENSION(200) :: data_output_pr_user = ' ' 254 257 CHARACTER (LEN=20), DIMENSION(10) :: netcdf_precision = ' ' 255 258 … … 269 272 intermediate_timestep_count, intermediate_timestep_count_max, & 270 273 iran = -1234567, last_dt_change = 0, maximum_grid_level, & 271 mgcycles = 0, mg_cycles = -1, mg_switch_to_pe0_level = 0, & 272 ngsrb = 2, nsor = 20, nsor_ini = 100, n_sor, & 273 normalizing_region = 0, nz_do1d, nz_do3d = -9999, & 274 max_pr_user = 0, mgcycles = 0, mg_cycles = -1, & 275 mg_switch_to_pe0_level = 0, ngsrb = 2, nsor = 20, & 276 nsor_ini = 100, n_sor, normalizing_region = 0, & 277 num_user_profiles = 0, nz_do1d, nz_do3d = -9999, & 274 278 outflow_damping_width = -1, prt_time_count = 0, runnr = 0, & 275 279 skip_do_avs = 0, timestep_count = 0 … … 719 723 ( 'unknown', i9 = 1, 78 ) /) 720 724 721 CHARACTER (LEN=7), DIMENSION( 100) :: dopr_unit = 'unknown'725 CHARACTER (LEN=7), DIMENSION(300) :: dopr_unit = 'unknown' 722 726 723 727 CHARACTER (LEN=7), DIMENSION(0:1,100) :: do2d_unit, do3d_unit … … 966 970 cross_ts_numbers(crmax,crmax) = 0, & 967 971 cross_ts_number_count(crmax) = 0, dopr_crossindex(100) = 0, & 968 dopr_index( 100) = 0, dopr_initial_index(100) = 0, &972 dopr_index(300) = 0, dopr_initial_index(300) = 0, & 969 973 dots_crossindex(100) = 0, dots_index(100) = 0, & 970 974 linecolors(10) = (/ 2, 3, 4, 5, 7, 8, 12, 15, 16, 23 /), & … … 1050 1054 1051 1055 CHARACTER (LEN=40) :: region(0:9) 1052 INTEGER :: statistic_regions = 0, var_hom = 80, var_sum = 80, var_ts = 1001056 INTEGER :: pr_palm = 80, statistic_regions = 0, var_ts = 100 1053 1057 INTEGER :: u_max_ijk(3), v_max_ijk(3), w_max_ijk(3) 1054 1058 LOGICAL :: flow_statistics_called = .FALSE. -
palm/trunk/SOURCE/parin.f90
r83 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Size of hom increased by the maximum number of user-defined profiles, 7 ! var_hom renamed pr_palm 7 8 ! 8 9 ! Former revisions: … … 179 180 runnr = runnr + 1 180 181 181 ELSE 182 ! 183 !-- This is not a restart job. 184 !-- Check, if the grid point numbers are well defined. 182 ENDIF 183 184 ! 185 !-- Definition of names of areas used for computing statistics. They must 186 !-- be defined at this place, because they are allowed to be redefined by 187 !-- the user in user_parin. 188 region = 'total domain' 189 190 ! 191 !-- Read runtime parameters given by the user for this run (namelist "d3par"). 192 !-- The namelist "d3par" can be omitted. In that case, default values are 193 !-- used for the parameters. 194 READ ( 11, d3par, END=20 ) 195 196 ! 197 !-- Read control parameters for optionally used model software packages 198 20 CALL package_parin 199 200 ! 201 !-- Read user-defined variables 202 CALL user_parin 203 204 ! 205 !-- Check in case of initial run, if the grid point numbers are well defined 206 !-- and allocate some arrays which are already needed in init_pegrid or 207 !-- check_parameters. During restart jobs, these arrays will be allocated 208 !-- in read_var_list. All other arrays are allocated in init_3d_model. 209 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 210 185 211 IF ( nx <= 0 ) THEN 186 212 IF ( myid == 0 ) THEN … … 202 228 ENDIF 203 229 204 !205 !-- Allocate arrays which will be already initialized in init_pegrid or206 !-- check_parameters. During restart jobs, these arrays will be allocated207 !-- in read_var_list. All other arrays are allocated in init_3d_model.208 230 ALLOCATE( ug(0:nz+1), vg(0:nz+1), & 209 231 pt_init(0:nz+1), q_init(0:nz+1), u_init(0:nz+1), & 210 232 v_init(0:nz+1), & 211 hom(0:nz+1,2, var_hom,0:statistic_regions) )233 hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions) ) 212 234 hom = 0.0 213 235 214 236 ENDIF 215 216 !217 !-- Definition of names of areas used for computing statistics. They must218 !-- be defined at this place, because they are allowed to be redefined by219 !-- the user in user_parin.220 region = 'total domain'221 222 !223 !-- Read runtime parameters given by the user for this run (namelist "d3par").224 !-- The namelist "d3par" can be omitted. In that case, default values are225 !-- used for the parameters.226 READ ( 11, d3par, END=20 )227 228 !229 !-- Read control parameters for optionally used model software packages230 20 CALL package_parin231 232 !233 !-- Read user-defined variables234 CALL user_parin235 237 236 238 ! -
palm/trunk/SOURCE/read_var_list.f90
r77 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! +max_pr_user (version 3.1), var_hom renamed pr_palm 7 7 ! 8 8 ! Former revisions: … … 32 32 USE arrays_3d 33 33 USE averaging 34 USE control_parameters 34 35 USE grid_variables 35 36 USE indices … … 38 39 USE profil_parameter 39 40 USE statistics 40 USE control_parameters41 41 42 42 IMPLICIT NONE … … 44 44 CHARACTER (LEN=10) :: binary_version, version_on_file 45 45 CHARACTER (LEN=30) :: variable_chr 46 47 INTEGER :: max_pr_user_on_file 46 48 47 49 CALL check_open( 13 ) … … 50 52 !-- Make version number check first 51 53 READ ( 13 ) version_on_file 52 binary_version = '3. 0'54 binary_version = '3.1' 53 55 IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN 54 56 IF ( myid == 0 ) THEN … … 74 76 ENDIF 75 77 READ ( 13 ) nz 78 79 READ ( 13 ) variable_chr 80 IF ( TRIM( variable_chr ) /= 'max_pr_user' ) THEN 81 PRINT*, '+++ read_var_list: max_pr_user not found in data from ', & 82 'prior run on PE ', myid 83 CALL local_stop 84 ENDIF 85 READ ( 13 ) max_pr_user_on_file 86 IF ( max_pr_user_on_file /= max_pr_user ) THEN 87 IF ( myid == 0 ) THEN 88 PRINT*, '+++ read_var_list: version mismatch concerning maximum', & 89 ' number of user profiles' 90 PRINT*, ' max_pr_user on file = "', & 91 max_pr_user_on_file, '"' 92 PRINT*, ' max_pr_user from run = "', & 93 max_pr_user, '"' 94 ENDIF 95 CALL local_stop 96 ENDIF 97 76 98 READ ( 13 ) variable_chr 77 99 IF ( TRIM( variable_chr ) /= 'statistic_regions' ) THEN … … 83 105 ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1), v_init(0:nz+1), & 84 106 pt_init(0:nz+1), q_init(0:nz+1), & 85 hom(0:nz+1,2, var_hom,0:statistic_regions) )107 hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions) ) 86 108 87 109 ! -
palm/trunk/SOURCE/run_control.f90
r83 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! var_hom renamed pr_palm 7 7 ! 8 8 ! Former revisions: … … 86 86 simulated_time-INT( simulated_time ), dt_3d, & 87 87 timestep_reason, change_chr, u_max, disturb_chr, & 88 v_max, disturb_chr, w_max, hom(nzb,1, var_hom,0), &89 hom(nzb+8,1, var_hom,0), hom(nzb+3,1,var_hom,0), &90 hom(nzb+6,1, var_hom,0), hom(nzb+4,1,var_hom,0), &91 hom(nzb+5,1, var_hom,0), hom(nzb+9,1,var_hom,0), &92 hom(nzb+10,1, var_hom,0), u_max_ijk(1:3), &88 v_max, disturb_chr, w_max, hom(nzb,1,pr_palm,0), & 89 hom(nzb+8,1,pr_palm,0), hom(nzb+3,1,pr_palm,0), & 90 hom(nzb+6,1,pr_palm,0), hom(nzb+4,1,pr_palm,0), & 91 hom(nzb+5,1,pr_palm,0), hom(nzb+9,1,pr_palm,0), & 92 hom(nzb+10,1,pr_palm,0), u_max_ijk(1:3), & 93 93 v_max_ijk(1:3), w_max_ijk(1:3), & 94 94 advected_distance_x/1000.0, & -
palm/trunk/SOURCE/time_integration.f90
r77 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! var_hom renamed pr_palm 7 7 ! 8 8 ! Former revisions: … … 187 187 time_disturb = time_disturb + dt_3d 188 188 IF ( time_disturb >= dt_disturb ) THEN 189 IF ( hom(nzb+5,1, var_hom,0) < disturbance_energy_limit ) THEN189 IF ( hom(nzb+5,1,pr_palm,0) < disturbance_energy_limit ) THEN 190 190 CALL disturb_field( nzb_u_inner, tend, u ) 191 191 CALL disturb_field( nzb_v_inner, tend, v ) -
palm/trunk/SOURCE/user_interface.f90
r86 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Calculation and output of user-defined profiles: new routine 7 ! user_check_data_output_pr, +data_output_pr_user, max_pr_user in userpar, 8 ! routine user_statistics has got two more arguments 6 9 ! Bugfix: field_chr renamed field_char 7 10 ! … … 56 59 57 60 USE control_parameters 61 USE pegrid 58 62 USE statistics 59 63 USE user … … 63 67 CHARACTER (LEN=80) :: zeile 64 68 65 66 NAMELIST /userpar/ data_output_user, region 69 INTEGER :: i, j 70 71 72 NAMELIST /userpar/ data_output_pr_user, data_output_user, max_pr_user, & 73 region 67 74 68 75 ! … … 82 89 READ ( 11, userpar ) 83 90 user_defined_namelist_found = .TRUE. 91 92 ! 93 !-- Determine the maximum number of user-profiles allowed to be calculated. 94 !-- This number may be (much) larger than num_user_profiles which only gives 95 !-- the number of user-defined profiles to be output. 96 IF ( data_output_pr_user(1) /= ' ' .AND. max_pr_user == 0 ) THEN 97 max_pr_user = 100 98 ENDIF 99 100 ! 101 !-- Determine the number of user-defined profiles and append them to the 102 !-- standard data output (data_output_pr) 103 IF ( data_output_pr_user(1) /= ' ' ) THEN 104 i = 1 105 DO WHILE ( data_output_pr(i) /= ' ' .AND. i <= 100 ) 106 i = i + 1 107 ENDDO 108 j = 1 109 DO WHILE ( data_output_pr_user(j) /= ' ' .AND. j <= 100 ) 110 data_output_pr(i) = data_output_pr_user(j) 111 num_user_profiles = num_user_profiles + 1 112 i = i + 1 113 j = j + 1 114 ENDDO 115 ENDIF 116 117 IF ( num_user_profiles > max_pr_user ) THEN 118 IF ( myid == 0 ) THEN 119 PRINT*, '+++ user_parin: number of user-defined profiles given by ' 120 PRINT*, ' data_output_pr_user (', num_user_profiles, & 121 ') is larger than the allowed maximum' 122 PRINT*, ' number of profiles: max_pr_user(', & 123 max_pr_user, ')' 124 ENDIF 125 CALL local_stop 126 ENDIF 84 127 85 128 100 RETURN … … 447 490 448 491 449 SUBROUTINE user_statistics( sr ) 450 451 !------------------------------------------------------------------------------! 452 ! 453 ! Description: 454 ! ------------ 455 ! Calculation of user-defined statistics 492 SUBROUTINE user_statistics( mode, sr, tn ) 493 494 !------------------------------------------------------------------------------! 495 ! 496 ! Description: 497 ! ------------ 498 ! Calculation of user-defined statistics, i.e. horizontally averaged profiles 499 ! and time series. 456 500 ! This routine is called for every statistic region sr defined by the user, 457 501 ! but at least for the region "total domain" (sr=0). 458 !------------------------------------------------------------------------------! 459 502 ! See section 3.5.4 on how to define, calculate, and output user defined 503 ! quantities. 504 !------------------------------------------------------------------------------! 505 506 USE arrays_3d 507 USE indices 460 508 USE statistics 461 509 USE user … … 463 511 IMPLICIT NONE 464 512 465 INTEGER :: sr 466 467 ! 468 !-- Sample on how to add values for the user-defined time series quantities. 469 !-- These have to be defined before in routine user_init. This sample 470 !-- creates two time series for the absolut values of the horizontal 471 !-- velocities u and v. 472 ! ts_value(dots_num_palm+1,sr) = ABS( u_max ) 473 ! ts_value(dots_num_palm+2,sr) = ABS( v_max ) 513 CHARACTER (LEN=*) :: mode 514 515 INTEGER :: i, j, k, sr, tn 516 517 518 IF ( mode == 'profiles' ) THEN 519 520 ! 521 !-- Sample on how to calculate horizontally averaged profiles of user- 522 !-- defined quantities. Each quantity is identified by the index 523 !-- "pr_palm+#" where "#" is an integer starting from 1. These 524 !-- user-profile-numbers must also be assigned to the respective strings 525 !-- given by data_output_pr_user in routine user_check_data_output_pr. 526 ! !$OMP DO 527 ! DO i = nxl, nxr 528 ! DO j = nys, nyn 529 ! DO k = nzb_s_outer(j,i)+1, nzt 530 !! 531 !!-- Sample on how to calculate the profile for vertical velocity 532 ! sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) + & 533 ! w(k,j,i) & 534 ! * rmask(j,i,sr) 535 !! 536 !!-- Further profiles can be defined and calculated by increasing 537 !!-- the second index of array sums_l (replace ... appropriately) 538 ! sums_l(k,pr_palm+2,tn) = sums_l(k,pr_palm+2,tn) + ... & 539 ! * rmask(j,i,sr) 540 ! ENDDO 541 ! ENDDO 542 ! ENDDO 543 544 ELSEIF ( mode == 'time_series' ) THEN 545 546 ! 547 !-- Sample on how to add values for the user-defined time series quantities. 548 !-- These have to be defined before in routine user_init. This sample 549 !-- creates two time series for the absolut values of the horizontal 550 !-- velocities u and v. 551 ! ts_value(dots_num_palm+1,sr) = ABS( u_max ) 552 ! ts_value(dots_num_palm+2,sr) = ABS( v_max ) 553 554 ENDIF 474 555 475 556 END SUBROUTINE user_statistics … … 649 730 650 731 END SUBROUTINE user_check_data_output 732 733 734 735 SUBROUTINE user_check_data_output_pr( variable, var_count, unit ) 736 737 !------------------------------------------------------------------------------! 738 ! 739 ! Description: 740 ! ------------ 741 ! Set the unit of user defined profile output quantities. For those variables 742 ! not recognized by the user, the parameter unit is set to "illegal", which 743 ! tells the calling routine that the output variable is not defined and leads 744 ! to a program abort. 745 !------------------------------------------------------------------------------! 746 747 USE arrays_3d 748 USE indices 749 USE netcdf_control 750 USE profil_parameter 751 USE statistics 752 USE user 753 754 IMPLICIT NONE 755 756 CHARACTER (LEN=*) :: unit, variable 757 758 INTEGER :: index, var_count 759 760 761 SELECT CASE ( TRIM( variable ) ) 762 763 ! 764 !-- Uncomment and extend the following lines, if necessary. 765 !-- Add additional CASE statements depending on the number of quantities 766 !-- for which profiles are to be calculated. The respective calculations 767 !-- to be performed have to be added in routine user_statistics. 768 !-- The quantities are (internally) identified by a user-profile-number 769 !-- (see variable "index" below). The first user-profile must be assigned 770 !-- the number "pr_palm+1", the second one "pr_palm+2", etc. The respective 771 !-- user-profile-numbers have also to be used in routine user_statistics! 772 ! CASE ( 'w_user' ) ! quantity string as given in 773 ! ! data_output_pr_user 774 ! index = pr_palm + 1 775 ! dopr_index(var_count) = index ! the quantities' user-profile-number 776 ! dopr_unit(var_count) = 'm/s' ! quantity unit 777 ! hom(:,2,index,:) = SPREAD( zw, 2, statistic_regions+1 ) 778 ! ! grid on which the quantity is 779 ! ! defined (use zu or zw) 780 781 CASE DEFAULT 782 unit = 'illegal' 783 784 END SELECT 785 786 787 END SUBROUTINE user_check_data_output_pr 651 788 652 789 -
palm/trunk/SOURCE/write_var_list.f90
r77 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! +max_pr_user (version 3.1) 7 7 ! 8 8 ! Former revisions: … … 32 32 USE arrays_3d 33 33 USE averaging 34 USE control_parameters 34 35 USE grid_variables 35 36 USE indices … … 38 39 USE profil_parameter 39 40 USE statistics 40 USE control_parameters41 41 42 42 IMPLICIT NONE … … 46 46 47 47 48 binary_version = '3. 0'48 binary_version = '3.1' 49 49 50 50 WRITE ( 14 ) binary_version ! opened in write_3d_binary … … 52 52 WRITE ( 14 ) 'nz ' 53 53 WRITE ( 14 ) nz 54 WRITE ( 14 ) 'max_pr_user ' 55 WRITE ( 14 ) max_pr_user 54 56 WRITE ( 14 ) 'statistic_regions ' 55 57 WRITE ( 14 ) statistic_regions
Note: See TracChangeset
for help on using the changeset viewer.
