Changeset 87


Ignore:
Timestamp:
May 22, 2007 3:46:47 PM (14 years ago)
Author:
raasch
Message:

Preliminary update for user defined profiles

Location:
palm/trunk
Files:
14 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_3.5.4.html

    r62 r87  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    22<html><head>
    3 <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM
    4 chapter 3.5.5</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;14050943"> <meta name="CHANGED" content="20041117;12180008"> <meta name="KEYWORDS" content="parallel LES model"> <style>
     3<meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 3.5.5</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;14050943"> <meta name="CHANGED" content="20041117;12180008"> <meta name="KEYWORDS" content="parallel LES model"> <style>
    54<!--
    65@page { size: 21cm 29.7cm }
     
    1413quantities which are not part of PALM's standard output. The basic user
    1514interface includes a number of subroutines which allow the calculation
    16 of user-defined quantities and output of these quantities as 2d cross
     15of user-defined quantities and output of these quantities as time series, (horizontally averaged) vertical profile, 2d cross
    1716section or 3d volume data. The respective subroutines
    18 contain&nbsp;as
    19 an example code lines (written as comment lines) for calculating and
    20 output the square of the u-component velocity (note: this quantity
     17contain sample code lines (written as comment lines) for defining, calculating and
     18output of such quantities.<br><br>The
     19rest of this chapter explains step-by-step how to modify/extend the
     20default file user_interface.f90 in order to generate the respective output.<br><br><br><br><h4><span style="font-weight: bold;">Output of 2d cross sections or 3d volume data</span></h4><br>This example shows the output of the
     21quantity "square of the u-component" (<span style="font-weight: bold;">Note:</span> this quantity
    2122could of course easily be calculated from the u-component by
    2223postprocessing the PALM output so that calculation within PALM is not
    23 necessarily required).<br><br>The
    24 rest of this chapter explains step-by-step how to modify/extend the
    25 default file user_interface.f90 in order to generate output of the
    26 quantity "square of the u-component". If more than one user-defined
     24necessarily required). If more than one user-defined
    2725quantity shall be output, these steps have to be carried out in the
    2826same way for each of the quantities.<br><br><ol><li>The
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r86 r87  
    11New:
    22---
     3Calculation and output of user-defined profiles. New &userpar parameters data_output_pr_user and max_pr_user.
     4
     5check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list
    36
    47
    58Changed:
    69-------
    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)
     10Division 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.
    811
    9 pres
     12data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration
    1013
    1114
     
    1619Bugfix: field_chr renamed field_char (user_interface).
    1720
    18 poisfft, user_interface
     21Bugfix: output of use_upstream_for_tke (header).
     22
     23header, poisfft, user_interface
  • palm/trunk/SOURCE/check_parameters.f90

    r77 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Check for user-defined profiles
    77!
    88! Former revisions:
     
    12141214          CASE ( 'u', '#u' )
    12151215             dopr_index(i) = 1
     1216             dopr_unit(i)  = 'm/s'
    12161217             hom(:,2,1,:)  = SPREAD( zu, 2, statistic_regions+1 )
    12171218             IF ( data_output_pr(i)(1:1) == '#' )  THEN
     
    12231224          CASE ( 'v', '#v' )
    12241225             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 )
    12261228             IF ( data_output_pr(i)(1:1) == '#' )  THEN
    12271229                dopr_initial_index(i) = 6
     
    12321234          CASE ( 'w' )
    12331235             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 )
    12351238
    12361239          CASE ( 'pt', '#pt' )
    12371240             IF ( .NOT. cloud_physics ) THEN
    12381241                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'
    12391483                hom(:,2,4,:)  = SPREAD( zu, 2, statistic_regions+1 )
    12401484                IF ( data_output_pr(i)(1:1) == '#' )  THEN
     
    12441488                   data_output_pr(i)     = data_output_pr(i)(2:)
    12451489                ENDIF
    1246              ELSE
    1247                 dopr_index(i) = 43
    1248                 hom(:,2,43,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1249                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1250                    dopr_initial_index(i) = 28
    1251                    hom(:,2,28,:)         = SPREAD( zu, 2, statistic_regions+1 )
    1252                    hom(nzb,2,28,:)       = 0.0    ! weil zu(nzb) negativ ist
    1253                    data_output_pr(i)     = data_output_pr(i)(2:)
    1254                 ENDIF
    1255              ENDIF
    1256 
    1257           CASE ( 'e' )
    1258              dopr_index(i)  = 8
    1259              hom(:,2,8,:)   = SPREAD( zu, 2, statistic_regions+1 )
    1260              hom(nzb,2,8,:) = 0.0
    1261 
    1262           CASE ( 'km', '#km' )
    1263              dopr_index(i)  = 9
    1264              hom(:,2,9,:)   = SPREAD( zu, 2, statistic_regions+1 )
    1265              hom(nzb,2,9,:) = 0.0
    1266              IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1267                 dopr_initial_index(i) = 23
    1268                 hom(:,2,23,:)         = hom(:,2,9,:)
    1269                 data_output_pr(i)     = data_output_pr(i)(2:)
    1270              ENDIF
    1271 
    1272           CASE ( 'kh', '#kh' )
    1273              dopr_index(i)   = 10
    1274              hom(:,2,10,:)   = SPREAD( zu, 2, statistic_regions+1 )
    1275              hom(nzb,2,10,:) = 0.0
    1276              IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1277                 dopr_initial_index(i) = 24
    1278                 hom(:,2,24,:)         = hom(:,2,10,:)
    1279                 data_output_pr(i)     = data_output_pr(i)(2:)
    1280              ENDIF
    1281 
    1282           CASE ( 'l', '#l' )
    1283              dopr_index(i)   = 11
    1284              hom(:,2,11,:)   = SPREAD( zu, 2, statistic_regions+1 )
    1285              hom(nzb,2,11,:) = 0.0
    1286              IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1287                 dopr_initial_index(i) = 25
    1288                 hom(:,2,25,:)         = hom(:,2,11,:)
    1289                 data_output_pr(i)     = data_output_pr(i)(2:)
    1290              ENDIF
    1291 
    1292           CASE ( 'w"u"' )
    1293              dopr_index(i) = 12
    1294              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) = 13
    1299              hom(:,2,13,:) = SPREAD( zw, 2, statistic_regions+1 )
    1300 
    1301           CASE ( 'w"v"' )
    1302              dopr_index(i) = 14
    1303              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) = 15
    1308              hom(:,2,15,:) = SPREAD( zw, 2, statistic_regions+1 )
    1309 
    1310           CASE ( 'w"pt"' )
    1311              dopr_index(i) = 16
    1312              hom(:,2,16,:) = SPREAD( zw, 2, statistic_regions+1 )
    1313 
    1314           CASE ( 'w*pt*' )
    1315              dopr_index(i) = 17
    1316              hom(:,2,17,:) = SPREAD( zw, 2, statistic_regions+1 )
    1317 
    1318           CASE ( 'wpt' )
    1319              dopr_index(i) = 18
    1320              hom(:,2,18,:) = SPREAD( zw, 2, statistic_regions+1 )
    1321 
    1322           CASE ( 'wu' )
    1323              dopr_index(i) = 19
    1324              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) = 20
    1329              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) = 21
    1334              hom(:,2,21,:) = SPREAD( zw, 2, statistic_regions+1 )
    1335 
    1336           CASE ( 'wptBC' )
    1337              dopr_index(i) = 22
    1338              hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 )
    1339 
    1340           CASE ( 'u*2' )
    1341              dopr_index(i) = 30
    1342              hom(:,2,30,:) = SPREAD( zu, 2, statistic_regions+1 )
    1343 
    1344           CASE ( 'v*2' )
    1345              dopr_index(i) = 31
    1346              hom(:,2,31,:) = SPREAD( zu, 2, statistic_regions+1 )
    1347 
    1348           CASE ( 'w*2' )
    1349              dopr_index(i) = 32
    1350              hom(:,2,32,:) = SPREAD( zw, 2, statistic_regions+1 )
    1351 
    1352           CASE ( 'pt*2' )
    1353              dopr_index(i) = 33
    1354              hom(:,2,33,:) = SPREAD( zu, 2, statistic_regions+1 )
    1355 
    1356           CASE ( 'e*' )
    1357              dopr_index(i) = 34
    1358              hom(:,2,34,:) = SPREAD( zu, 2, statistic_regions+1 )
    1359 
    1360           CASE ( 'w*2pt*' )
    1361              dopr_index(i) = 35
    1362              hom(:,2,35,:) = SPREAD( zw, 2, statistic_regions+1 )
    1363 
    1364           CASE ( 'w*pt*2' )
    1365              dopr_index(i) = 36
    1366              hom(:,2,36,:) = SPREAD( zw, 2, statistic_regions+1 )
    1367 
    1368           CASE ( 'w*e*' )
    1369              dopr_index(i) = 37
    1370              hom(:,2,37,:) = SPREAD( zw, 2, statistic_regions+1 )
    1371 
    1372           CASE ( 'w*3' )
    1373              dopr_index(i) = 38
    1374              hom(:,2,38,:) = SPREAD( zw, 2, statistic_regions+1 )
    1375 
    1376           CASE ( 'Sw' )
    1377              dopr_index(i) = 39
    1378              hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 )
    1379 
    1380           CASE ( 'q', '#q' )
    1381              IF ( .NOT. cloud_physics )  THEN
    1382                 IF ( myid == 0 )  THEN
    1383                    PRINT*, '+++ check_parameters:  data_output_pr = ', &
    1384                            data_output_pr(i),                          &
    1385                            '    is not implemented for cloud_physics = FALSE'
    1386                 ENDIF
    1387                 CALL local_stop
    1388              ELSE
    1389                 dopr_index(i) = 41
    1390                 hom(:,2,41,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1391                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1392                    dopr_initial_index(i) = 26
    1393                    hom(:,2,26,:)         = SPREAD( zu, 2, statistic_regions+1 )
    1394                    hom(nzb,2,26,:)       = 0.0    ! weil zu(nzb) negativ ist
    1395                    data_output_pr(i)     = data_output_pr(i)(2:)
    1396                 ENDIF
    1397              ENDIF
    1398 
    1399           CASE ( 's', '#s' )
    1400              IF ( .NOT. passive_scalar )  THEN
    1401                 IF ( myid == 0 )  THEN
    1402                    PRINT*, '+++ check_parameters:  data_output_pr = ', &
    1403                            data_output_pr(i),                          &
    1404                            '    is not implemented for passive_scalar = FALSE'
    1405                 ENDIF
    1406                 CALL local_stop
    1407              ELSE
    1408                 dopr_index(i) = 41
    1409                 hom(:,2,41,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1410                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1411                    dopr_initial_index(i) = 26
    1412                    hom(:,2,26,:)         = SPREAD( zu, 2, statistic_regions+1 )
    1413                    hom(nzb,2,26,:)       = 0.0    ! weil zu(nzb) negativ ist
    1414                    data_output_pr(i)     = data_output_pr(i)(2:)
    1415                 ENDIF
    1416              ENDIF
    1417 
    1418           CASE ( 'qv', '#qv' )
    1419              IF ( .NOT. cloud_physics ) THEN
    1420                 dopr_index(i) = 41
    1421                 hom(:,2,41,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1422                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1423                    dopr_initial_index(i) = 26
    1424                    hom(:,2,26,:)         = SPREAD( zu, 2, statistic_regions+1 )
    1425                    hom(nzb,2,26,:)       = 0.0    ! weil zu(nzb) negativ ist
    1426                    data_output_pr(i)     = data_output_pr(i)(2:)
    1427                 ENDIF
    1428              ELSE
    1429                 dopr_index(i) = 42
    1430                 hom(:,2,42,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1431                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1432                    dopr_initial_index(i) = 27
    1433                    hom(:,2,27,:)         = SPREAD( zu, 2, statistic_regions+1 )
    1434                    hom(nzb,2,27,:)       = 0.0    ! weil zu(nzb) negativ ist
    1435                    data_output_pr(i)     = data_output_pr(i)(2:)
    1436                 ENDIF
    1437              ENDIF
    1438 
    1439           CASE ( 'lpt', '#lpt' )
    1440              IF ( .NOT. cloud_physics ) THEN
    1441                 IF ( myid == 0 )  THEN
    1442                    PRINT*, '+++ check_parameters:  data_output_pr = ', &
    1443                            data_output_pr(i),                          &
    1444                            '    is not implemented for cloud_physics = FALSE'
    1445                 ENDIF
    1446                 CALL local_stop
    1447              ELSE
    1448                 dopr_index(i) = 4
    1449                 hom(:,2,4,:)  = SPREAD( zu, 2, statistic_regions+1 )
    1450                 IF ( data_output_pr(i)(1:1) == '#' )  THEN
    1451                    dopr_initial_index(i) = 7
    1452                    hom(:,2,7,:)          = SPREAD( zu, 2, statistic_regions+1 )
    1453                    hom(nzb,2,7,:)        = 0.0    ! weil zu(nzb) negativ ist
    1454                    data_output_pr(i)     = data_output_pr(i)(2:)
    1455                 ENDIF
    14561490             ENDIF
    14571491
    14581492          CASE ( 'vpt', '#vpt' )
    14591493             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 )
    14611496             IF ( data_output_pr(i)(1:1) == '#' )  THEN
    14621497                dopr_initial_index(i) = 29
     
    14681503          CASE ( 'w"vpt"' )
    14691504             dopr_index(i) = 45
     1505             dopr_unit(i)  = 'K m/s'
    14701506             hom(:,2,45,:) = SPREAD( zw, 2, statistic_regions+1 )
    14711507
    14721508          CASE ( 'w*vpt*' )
    14731509             dopr_index(i) = 46
     1510             dopr_unit(i)  = 'K m/s'
    14741511             hom(:,2,46,:) = SPREAD( zw, 2, statistic_regions+1 )
    14751512
    14761513          CASE ( 'wvpt' )
    14771514             dopr_index(i) = 47
     1515             dopr_unit(i)  = 'K m/s'
    14781516             hom(:,2,47,:) = SPREAD( zw, 2, statistic_regions+1 )
    14791517
     
    14881526             ELSE
    14891527                dopr_index(i) = 48
     1528                dopr_unit(i)  = 'kg/kg m/s'
    14901529                hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
    14911530             ENDIF
     
    15011540             ELSE
    15021541                dopr_index(i) = 49
     1542                dopr_unit(i)  = 'kg/kg m/s'
    15031543                hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
    15041544             ENDIF
     
    15141554             ELSE
    15151555                dopr_index(i) = 50
     1556                dopr_unit(i)  = 'kg/kg m/s'
    15161557                hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
    15171558             ENDIF
     
    15271568             ELSE
    15281569                dopr_index(i) = 48
     1570                dopr_unit(i)  = 'kg/m3 m/s'
    15291571                hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
    15301572             ENDIF
     
    15401582             ELSE
    15411583                dopr_index(i) = 49
     1584                dopr_unit(i)  = 'kg/m3 m/s'
    15421585                hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
    15431586             ENDIF
     
    15531596             ELSE
    15541597                dopr_index(i) = 50
     1598                dopr_unit(i)  = 'kg/m3 m/s'
    15551599                hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
    15561600             ENDIF
     
    15601604             THEN
    15611605                dopr_index(i) = 48
     1606                dopr_unit(i)  = 'kg/kg m/s'
    15621607                hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
    15631608             ELSEIF( humidity .AND. cloud_physics ) THEN
    15641609                dopr_index(i) = 51
     1610                dopr_unit(i)  = 'kg/kg m/s'
    15651611                hom(:,2,51,:) = SPREAD( zw, 2, statistic_regions+1 )
    15661612             ELSE
     
    15781624             THEN
    15791625                dopr_index(i) = 49
     1626                dopr_unit(i)  = 'kg/kg m/s'
    15801627                hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
    15811628             ELSEIF( humidity .AND. cloud_physics ) THEN
    15821629                dopr_index(i) = 52
     1630                dopr_unit(i)  = 'kg/kg m/s'
    15831631                hom(:,2,52,:) = SPREAD( zw, 2, statistic_regions+1 )
    15841632             ELSE
     
    15961644             THEN
    15971645                dopr_index(i) = 50
     1646                dopr_unit(i)  = 'kg/kg m/s'
    15981647                hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
    15991648             ELSEIF( humidity .AND. cloud_physics ) THEN
    16001649                dopr_index(i) = 53
     1650                dopr_unit(i)  = 'kg/kg m/s'
    16011651                hom(:,2,53,:) = SPREAD( zw, 2, statistic_regions+1 )
    16021652             ELSE
     
    16201670             ELSE
    16211671                dopr_index(i) = 54
     1672                dopr_unit(i)  = 'kg/kg'
    16221673                hom(:,2,54,:)  = SPREAD( zu, 2, statistic_regions+1 )
    16231674             ENDIF
     
    16251676          CASE ( 'w*u*u*/dz' )
    16261677             dopr_index(i) = 55
     1678             dopr_unit(i)  = 'm2/s3'
    16271679             hom(:,2,55,:) = SPREAD( zu, 2, statistic_regions+1 )
    16281680
    16291681          CASE ( 'w*p*/dz' )
    16301682             dopr_index(i) = 56
     1683             dopr_unit(i)  = 'm2/s3'
    16311684             hom(:,2,56,:) = SPREAD( zu, 2, statistic_regions+1 )
    16321685
    16331686          CASE ( 'w"e/dz' )
    16341687             dopr_index(i) = 57
     1688             dopr_unit(i)  = 'm2/s3'
    16351689             hom(:,2,57,:) = SPREAD( zu, 2, statistic_regions+1 )
    16361690
    16371691          CASE ( 'u"pt"' )
    16381692             dopr_index(i) = 58
     1693             dopr_unit(i)  = 'K m/s'
    16391694             hom(:,2,58,:) = SPREAD( zu, 2, statistic_regions+1 )
    16401695
    16411696          CASE ( 'u*pt*' )
    16421697             dopr_index(i) = 59
     1698             dopr_unit(i)  = 'K m/s'
    16431699             hom(:,2,59,:) = SPREAD( zu, 2, statistic_regions+1 )
    16441700
    16451701          CASE ( 'upt_t' )
    16461702             dopr_index(i) = 60
     1703             dopr_unit(i)  = 'K m/s'
    16471704             hom(:,2,60,:) = SPREAD( zu, 2, statistic_regions+1 )
    16481705
    16491706          CASE ( 'v"pt"' )
    16501707             dopr_index(i) = 61
     1708             dopr_unit(i)  = 'K m/s'
    16511709             hom(:,2,61,:) = SPREAD( zu, 2, statistic_regions+1 )
    16521710             
    16531711          CASE ( 'v*pt*' )
    16541712             dopr_index(i) = 62
     1713             dopr_unit(i)  = 'K m/s'
    16551714             hom(:,2,62,:) = SPREAD( zu, 2, statistic_regions+1 )
    16561715
    16571716          CASE ( 'vpt_t' )
    16581717             dopr_index(i) = 63
     1718             dopr_unit(i)  = 'K m/s'
    16591719             hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 )
    16601720
    16611721
    16621722          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
    16681739
    16691740       END SELECT
  • palm/trunk/SOURCE/data_output_profiles.f90

    r4 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! var_hom renamed pr_palm
    77!
    88! Former revisions:
     
    111111
    112112             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 /), &
    114114                                     start = (/ 1 /), count = (/ 1 /) )
    115115             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 331 )
    116116
    117117             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 /), &
    119119                                     start = (/ 1 /), count = (/ 1 /) )
    120120             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 332 )
    121121
    122122             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 /), &
    124124                                     start = (/ 1 /), count = (/ 1 /) )
    125125             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 333 )
    126126
    127127             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)    /), &
    130130                                     start = (/ 1 /), count = (/ 1 /) )
    131131             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 334 )
    132132
    133133             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 /), &
    136136                                     start = (/ 1 /), count = (/ 1 /) )
    137137             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 335 )
    138138
    139139             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) /), &
    141141                                     start = (/ 1 /), count = (/ 1 /) )
    142142             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 336 )
     
    237237                         CASE ( 'ts2' )
    238238                            cross_normx_factor(k,j) = &
    239                              ( hom_sum(nzb+3,var_hom,normalizing_region) )**2
     239                             ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
    240240                         CASE ( 'wpt0' )
    241241                            cross_normx_factor(k,j) = &
     
    243243                         CASE ( 'wsts2' )
    244244                            cross_normx_factor(k,j) = &
    245                              hom_sum(nzb+8,var_hom,normalizing_region)  &
    246                            * ( hom_sum(nzb+3,var_hom,normalizing_region) )**2
     245                             hom_sum(nzb+8,pr_palm,normalizing_region)  &
     246                           * ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
    247247                         CASE ( 'ws2' )
    248248                            cross_normx_factor(k,j) = &
    249                              ( hom_sum(nzb+8,var_hom,normalizing_region) )**2
     249                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
    250250                         CASE ( 'ws2ts' )
    251251                            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)
    254254                         CASE ( 'ws3' )
    255255                            cross_normx_factor(k,j) = &
    256                              ( hom_sum(nzb+8,var_hom,normalizing_region) )**3
     256                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
    257257
    258258                      END SELECT
     
    262262                         CASE ( 'z_i' )
    263263                            cross_normy_factor(k,j) = &
    264                                     hom_sum(nzb+6,var_hom,normalizing_region)
     264                                    hom_sum(nzb+6,pr_palm,normalizing_region)
    265265
    266266                      END SELECT
     
    351351
    352352          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 /), &
    354354                                  start = (/ dopr_time_count /),               &
    355355                                  count = (/ 1 /) )
     
    357357
    358358          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 /), &
    360360                                  start = (/ dopr_time_count /),               &
    361361                                  count = (/ 1 /) )
     
    363363
    364364          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 /), &
    366366                                  start = (/ dopr_time_count /),               &
    367367                                  count = (/ 1 /) )
     
    369369
    370370          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)    /), &
    373373                                  start = (/ dopr_time_count /),               &
    374374                                  count = (/ 1 /) )
     
    376376
    377377          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 /), &
    380380                                  start = (/ dopr_time_count /),               &
    381381                                  count = (/ 1 /) )
     
    383383
    384384          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) /), &
    386386                                  start = (/ dopr_time_count /),               &
    387387                                  count = (/ 1 /) )
     
    482482                   CASE ( 'tsw2' )
    483483                      cross_normx_factor(k,j) = &
    484                             ( hom_sum(nzb+11,var_hom,normalizing_region) )**2
     484                            ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
    485485                   CASE ( 'wpt0' )
    486486                      cross_normx_factor(k,j) = &
     
    488488                   CASE ( 'wstsw2' )
    489489                      cross_normx_factor(k,j) = &
    490                               hom_sum(nzb+8,var_hom,normalizing_region)  &
    491                           * ( hom_sum(nzb+11,var_hom,normalizing_region) )**2
     490                              hom_sum(nzb+8,pr_palm,normalizing_region)  &
     491                          * ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
    492492                   CASE ( 'ws2' )
    493493                      cross_normx_factor(k,j) = &
    494                             ( hom_sum(nzb+8,var_hom,normalizing_region) )**2
     494                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
    495495                   CASE ( 'ws2tsw' )
    496496                      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)
    499499                   CASE ( 'ws3' )
    500500                      cross_normx_factor(k,j) = &
    501                             ( hom_sum(nzb+8,var_hom,normalizing_region) )**3
     501                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
    502502
    503503                END SELECT
     
    506506                   CASE ( 'z_i' )
    507507                      cross_normy_factor(k,j) = &
    508                                    hom_sum(nzb+6,var_hom,normalizing_region)
     508                                   hom_sum(nzb+6,pr_palm,normalizing_region)
    509509
    510510                END SELECT
  • palm/trunk/SOURCE/flow_statistics.f90

    r83 r87  
    44! Actual revisions:
    55! -----------------
    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
    79!
    810! Former revisions:
     
    8789!--    WARNING: next line still has to be adjusted for OpenMP
    8890       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 pres
    90        sums_l(nzb+10,var_sum,0) = sums_divnew_l(sr)  ! new divergence from pres
     91       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
    9193!--    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)! upstream
    93        sums_l(nzb+3:nzb+5,var_sum-1,0)  = sums_up_fraction_l(2,1:3,sr)! parts
    94        sums_l(nzb+6:nzb+8,var_sum-1,0)  = sums_up_fraction_l(3,1:3,sr)! from
    95        sums_l(nzb+9:nzb+11,var_sum-1,0) = sums_up_fraction_l(4,1:3,sr)! spline
     94       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
    9698
    9799!
     
    316318!--          quantities is seperated from the previous loop in order to
    317319!--          allow vectorization of that loop.
    318              sums_l(nzb+4,var_sum,tn) = sums_l(nzb+4,var_sum,tn) + sums_l_etot
    319              sums_l(nzb+5,var_sum,tn) = sums_l(nzb+5,var_sum,tn) + sums_l_eper
     320             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
    320322!
    321323!--          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) +   &
    323325                                        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) + &
    325327                                        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) + &
    327329                                        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) + &
    329331                                        ts(j,i)   * rmask(j,i,sr)
    330332          ENDDO
     
    652654
    653655       ENDIF
     656
     657!
     658!--    Calculate the user-defined profiles
     659       CALL user_statistics( 'profiles', sr, tn )
    654660       !$OMP END PARALLEL
    655661
     
    660666             sums_l(:,3,0)          = sums_l(:,3,0) + sums_l(:,3,i)
    661667             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
    664675          ENDDO
    665676       ENDIF
     
    679690!--    Profiles:
    680691       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)
    682693       ENDDO
    683694!--    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)
    685696!--    u* and so on
    686 !--    As sums(nzb:nzb+3,var_sum) are full 2D arrays (us, usws, vsws, ts) whose
     697!--    As sums(nzb:nzb+3,pr_palm) are full 2D arrays (us, usws, vsws, ts) whose
    687698!--    size is always ( nx + 1 ) * ( ny + 1 ), defined at the first grid layer
    688699!--    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)    / &
    690701                                    ngp_2dh(sr)
    691702!--    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)  / &
    693704                                    ngp_3d_inner(sr)
    694705!--    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) / &
    696707                                    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
    697717
    698718!
     
    748768       hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
    749769
    750        hom(:,1,var_hom-1,sr) = sums(:,var_sum-1)
     770       hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)
    751771                                       ! upstream-parts u_x, u_y, u_z, v_x,
    752772                                       ! 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) 
    754774                                       ! 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
    755780
    756781!
     
    792817       ENDDO
    793818
    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)
    796821
    797822!
     
    800825!--    The horizontal average at nzb+1 is input for the average temperature.
    801826       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) * &
    803828                                       hom(nzb,1,18,sr) * z_i(1) )**0.333333333
    804829!--       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)
    806831       ELSE
    807           hom(nzb+8,1,var_hom,sr)  = 0.0
    808           hom(nzb+11,1,var_hom,sr) = 0.0
     832          hom(nzb+8,1,pr_palm,sr)  = 0.0
     833          hom(nzb+11,1,pr_palm,sr) = 0.0
    809834       ENDIF
    810835
    811836!
    812837!--    Collect the time series quantities
    813        ts_value(1,sr) = hom(nzb+4,1,var_hom,sr)     ! E
    814        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*
    815840       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*
    818843       ts_value(6,sr) = u_max
    819844       ts_value(7,sr) = v_max
    820845       ts_value(8,sr) = w_max
    821        ts_value(9,sr) = hom(nzb+10,1,var_sum,sr)    ! new divergence
    822        ts_value(10,sr) = hom(nzb+9,1,var_hom,sr)    ! old Divergence
    823        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*
    826851       ts_value(14,sr) = hom(nzb,1,16,sr)           ! w'pt'   at k=0
    827852       ts_value(15,sr) = hom(nzb+1,1,16,sr)         ! w'pt'   at k=1
     
    829854       ts_value(17,sr) = hom(nzb,1,4,sr)            ! pt(0)
    830855       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)  ! splptx
    832        ts_value(20,sr) = hom(nzb+10,1,var_hom-1,sr) ! splpty
    833        ts_value(21,sr) = hom(nzb+11,1,var_hom-1,sr) ! splptz
     856       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
    834859       IF ( ts_value(5,sr) /= 0.0 )  THEN
    835860          ts_value(22,sr) = ts_value(4,sr)**2 / &
     
    841866!
    842867!--    Calculate additional statistics provided by the user interface
    843        CALL user_statistics( sr )
     868       CALL user_statistics( 'time_series', sr, 0 )
    844869
    845870    ENDDO    ! loop of the subregions
  • palm/trunk/SOURCE/header.f90

    r83 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: output of use_upstream_for_tke
    77!
    88! Former revisions:
     
    224224       WRITE ( io, 122 )  timestep_scheme
    225225    ENDIF
     226    IF ( use_upstream_for_tke )  WRITE ( io, 143 )
    226227    IF ( rayleigh_damping_factor /= 0.0 )  THEN
    227228       WRITE ( io, 123 )  rayleigh_damping_height, rayleigh_damping_factor
     
    11951196142 FORMAT ('     perturbation pressure is calculated at every Runge-Kutta ', &
    11961197                  'step')
     1198143 FORMAT ('     Euler/upstream scheme is used for the SGS turbulent ', &
     1199                  'kinetic energy')
    11971200150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
    11981201                  'conserved')
  • palm/trunk/SOURCE/init_3d_model.f90

    r77 r87  
    77! Actual revisions:
    88! -----------------
    9 !
     9! var_hom and var_sum renamed pr_palm
    1010!
    1111! Former revisions:
     
    8181              sums_divold_l(0:statistic_regions) )
    8282    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),   &
    8484              ngp_2dh_outer(nzb:nzt+1,0:statistic_regions),                 &
    8585              ngp_2dh_outer_l(nzb:nzt+1,0:statistic_regions),               &
    8686              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),   &
    8989              sums_l_l(nzb:nzt+1,0:statistic_regions,0:threads_per_task-1), &
    9090              sums_up_fraction_l(10,3,0:statistic_regions),                 &
     
    931931    ngp_3d_inner    = 0
    932932    ngp_3d          = 0
    933     ngp_sums        = ( nz + 2 ) * var_sum
     933    ngp_sums        = ( nz + 2 ) * ( pr_palm + max_pr_user )
    934934
    935935    DO  sr = 0, statistic_regions
  • palm/trunk/SOURCE/modules.f90

    r83 r87  
    55! Actual revisions:
    66! -----------------
    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
    810!
    911! Former revisions:
     
    250252    CHARACTER (LEN=10), DIMENSION(10)  ::  data_output_format = ' '
    251253    CHARACTER (LEN=10), DIMENSION(100) ::  data_output = ' ',    &
    252                                            data_output_pr = ' ', &
    253254                                           data_output_user = ' ', doav = ' '
     255    CHARACTER (LEN=10), DIMENSION(300) ::  data_output_pr = ' '
     256    CHARACTER (LEN=10), DIMENSION(200) ::  data_output_pr_user = ' '
    254257    CHARACTER (LEN=20), DIMENSION(10)  ::  netcdf_precision = ' '
    255258
     
    269272                intermediate_timestep_count, intermediate_timestep_count_max, &
    270273                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, &
    274278                outflow_damping_width = -1, prt_time_count = 0, runnr = 0, &
    275279                skip_do_avs = 0, timestep_count = 0
     
    719723             ( 'unknown', i9 = 1, 78 ) /)
    720724
    721     CHARACTER (LEN=7), DIMENSION(100) ::  dopr_unit = 'unknown'
     725    CHARACTER (LEN=7), DIMENSION(300) ::  dopr_unit = 'unknown'
    722726
    723727    CHARACTER (LEN=7), DIMENSION(0:1,100) ::  do2d_unit, do3d_unit
     
    966970                cross_ts_numbers(crmax,crmax) = 0, &
    967971                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, &
    969973                dots_crossindex(100) = 0, dots_index(100) = 0, &
    970974                linecolors(10) = (/ 2, 3, 4,  5, 7, 8, 12, 15, 16, 23 /), &
     
    10501054
    10511055    CHARACTER (LEN=40) ::  region(0:9)
    1052     INTEGER ::  statistic_regions = 0, var_hom = 80, var_sum = 80, var_ts = 100
     1056    INTEGER ::  pr_palm = 80, statistic_regions = 0, var_ts = 100
    10531057    INTEGER ::  u_max_ijk(3), v_max_ijk(3), w_max_ijk(3)
    10541058    LOGICAL ::  flow_statistics_called = .FALSE.
  • palm/trunk/SOURCE/parin.f90

    r83 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Size of hom increased by the maximum number of user-defined profiles,
     7! var_hom renamed pr_palm
    78!
    89! Former revisions:
     
    179180       runnr = runnr + 1
    180181
    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
    185211       IF ( nx <= 0 )  THEN
    186212          IF ( myid == 0 )  THEN
     
    202228       ENDIF
    203229
    204 !
    205 !--    Allocate arrays which will be already initialized in init_pegrid or
    206 !--    check_parameters. During restart jobs, these arrays will be allocated
    207 !--    in read_var_list. All other arrays are allocated in init_3d_model.
    208230       ALLOCATE( ug(0:nz+1), vg(0:nz+1), &
    209231                 pt_init(0:nz+1), q_init(0:nz+1), u_init(0:nz+1), &
    210232                 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) )
    212234       hom = 0.0
    213235
    214236    ENDIF
    215 
    216 !
    217 !-- Definition of names of areas used for computing statistics. They must
    218 !-- be defined at this place, because they are allowed to be redefined by
    219 !-- 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 are
    225 !-- used for the parameters.
    226     READ ( 11, d3par, END=20 )
    227 
    228 !
    229 !-- Read control parameters for optionally used model software packages
    230  20 CALL package_parin
    231 
    232 !
    233 !-- Read user-defined variables
    234     CALL user_parin
    235237
    236238!
  • palm/trunk/SOURCE/read_var_list.f90

    r77 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +max_pr_user (version 3.1), var_hom renamed pr_palm
    77!
    88! Former revisions:
     
    3232    USE arrays_3d
    3333    USE averaging
     34    USE control_parameters
    3435    USE grid_variables
    3536    USE indices
     
    3839    USE profil_parameter
    3940    USE statistics
    40     USE control_parameters
    4141
    4242    IMPLICIT NONE
     
    4444    CHARACTER (LEN=10) ::  binary_version, version_on_file
    4545    CHARACTER (LEN=30) ::  variable_chr
     46
     47    INTEGER            ::  max_pr_user_on_file
    4648   
    4749    CALL check_open( 13 )
     
    5052!-- Make version number check first
    5153    READ ( 13 )  version_on_file
    52     binary_version = '3.0'
     54    binary_version = '3.1'
    5355    IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
    5456       IF ( myid == 0 )  THEN
     
    7476    ENDIF
    7577    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
    7698    READ ( 13 )  variable_chr
    7799    IF ( TRIM( variable_chr ) /= 'statistic_regions' )  THEN
     
    83105    ALLOCATE( ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1), v_init(0:nz+1), &
    84106              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) )
    86108
    87109!
  • palm/trunk/SOURCE/run_control.f90

    r83 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! var_hom renamed pr_palm
    77!
    88! Former revisions:
     
    8686                          simulated_time-INT( simulated_time ), dt_3d,         &
    8787                          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),             &
    9393                          v_max_ijk(1:3), w_max_ijk(1:3),                      &
    9494                          advected_distance_x/1000.0,                          &
  • palm/trunk/SOURCE/time_integration.f90

    r77 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! var_hom renamed pr_palm
    77!
    88! Former revisions:
     
    187187             time_disturb = time_disturb + dt_3d
    188188             IF ( time_disturb >= dt_disturb )  THEN
    189                 IF ( hom(nzb+5,1,var_hom,0) < disturbance_energy_limit )  THEN
     189                IF ( hom(nzb+5,1,pr_palm,0) < disturbance_energy_limit )  THEN
    190190                   CALL disturb_field( nzb_u_inner, tend, u )
    191191                   CALL disturb_field( nzb_v_inner, tend, v )
  • palm/trunk/SOURCE/user_interface.f90

    r86 r87  
    44! Actual revisions:
    55! -----------------
     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
    69! Bugfix: field_chr renamed field_char
    710!
     
    5659
    5760    USE control_parameters
     61    USE pegrid
    5862    USE statistics
    5963    USE user
     
    6367    CHARACTER (LEN=80) ::  zeile
    6468
    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
    6774
    6875!
     
    8289    READ ( 11, userpar )
    8390    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
    84127
    85128100 RETURN
     
    447490
    448491
    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.
    456500! This routine is called for every statistic region sr defined by the user,
    457501! 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
    460508    USE statistics
    461509    USE user
     
    463511    IMPLICIT NONE
    464512
    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
    474555
    475556 END SUBROUTINE user_statistics
     
    649730
    650731 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
    651788
    652789
  • palm/trunk/SOURCE/write_var_list.f90

    r77 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +max_pr_user (version 3.1)
    77!
    88! Former revisions:
     
    3232    USE arrays_3d
    3333    USE averaging
     34    USE control_parameters
    3435    USE grid_variables
    3536    USE indices
     
    3839    USE profil_parameter
    3940    USE statistics
    40     USE control_parameters
    4141
    4242    IMPLICIT NONE
     
    4646
    4747
    48     binary_version = '3.0'
     48    binary_version = '3.1'
    4949
    5050    WRITE ( 14 )  binary_version    ! opened in write_3d_binary
     
    5252    WRITE ( 14 )  'nz                            '
    5353    WRITE ( 14 )  nz
     54    WRITE ( 14 )  'max_pr_user                   '
     55    WRITE ( 14 )  max_pr_user
    5456    WRITE ( 14 )  'statistic_regions             '
    5557    WRITE ( 14 )  statistic_regions
Note: See TracChangeset for help on using the changeset viewer.