Changeset 4504 for palm


Ignore:
Timestamp:
Apr 20, 2020 12:11:24 PM (7 months ago)
Author:
raasch
Message:

file re-formatted to follow the PALM coding standard, hint for setting rmask arrays added

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/user_module.f90

    r4497 r4504  
    2121! -----------------
    2222!
    23 ! 
     23!
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! hint for setting rmask arrays added
     28!
     29! 4497 2020-04-15 10:20:51Z raasch
    2730! file re-formatted to follow the PALM coding standard
    28 !
    2931!
    3032! 4495 2020-04-13 20:11:20Z raasch
     
    464466!
    465467!-- Example for defining a statistic region:
     468!-- ATTENTION: rmask = 0 is required at the ghost boundaries to guarantee correct statistic
     469!--            evaluations (otherwise ghost points would be counted twice). This setting has
     470!--            already been cared for in routine init_3d_model. Please don't set the ghost points
     471!--            /= 0. i.e. run the following loop only over nxl,nxr and nys,nyn.
    466472!     IF ( statistic_regions >= 1 )  THEN
    467473!        region = 1
  • palm/trunk/SOURCE/virtual_measurement_mod.f90

    r4498 r4504  
    2121! -----------------
    2222!
    23 ! 
     23!
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
    2727! file re-formatted to follow the PALM coding standard
    28 !
    2928!
    3029! 4481 2020-03-31 18:55:54Z maronga
     
    456455                        'If virtual measurements are taken, the number ' //                        &
    457456                        'of surrounding grid points must not be larger ' //                        &
    458                         'than the number of ghost points - 1, which is: ',                         &
    459                         nbgp - 1
     457                        'than the number of ghost points - 1, which is: ', nbgp - 1
    460458       CALL message( 'vm_check_parameters', 'PA0705', 1, 2, 0, 6, 0 )
    461459    ENDIF
     
    13601358!
    13611359!-- Determine the number of output timesteps.
    1362     ntimesteps = CEILING(                                                                          &
    1363                   ( end_time - MAX( vm_time_start, time_since_reference_point )                    &
    1364                   ) / dt_virtual_measurement )
     1360    ntimesteps = CEILING( ( end_time - MAX( vm_time_start, time_since_reference_point )            &
     1361                          ) / dt_virtual_measurement )
    13651362!
    13661363!-- Create directory where output files will be stored.
     
    13841381                                       dum_lon, dum_lat )
    13851382
    1386        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1387                                    attribute_name = 'site',                                        &
     1383       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'site',                   &
    13881384                                   value = TRIM( vmea(l)%site ) )
    1389        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1390                                    attribute_name = 'title',                                       &
     1385       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'title',                  &
    13911386                                   value = 'Virtual measurement output')
    1392        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1393                                    attribute_name = 'source',                                      &
     1387       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'source',                 &
    13941388                                   value = 'PALM-4U')
    1395        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1396                                    attribute_name = 'institution',                                 &
     1389       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'institution',            &
    13971390                                   value = input_file_atts%institution )
    1398        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1399                                    attribute_name = 'acronym',                                     &
     1391       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'acronym',                &
    14001392                                   value = input_file_atts%acronym )
    1401        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1402                                    attribute_name = 'author',                                      &
     1393       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'author',                 &
    14031394                                   value = input_file_atts%author )
    1404        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1405                                    attribute_name = 'contact_person',                              &
     1395       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'contact_person',         &
    14061396                                   value = input_file_atts%contact_person )
    1407        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1408                                    attribute_name = 'iop',                                         &
     1397       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'iop',                    &
    14091398                                   value = input_file_atts%campaign )
    1410        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1411                                    attribute_name = 'campaign',                                    &
     1399       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'campaign',               &
    14121400                                   value = 'PALM-4U' )
    1413        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1414                                    attribute_name = 'origin_time ',                                &
     1401       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_time ',           &
    14151402                                   value = origin_date_time)
    1416        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1417                                    attribute_name = 'location',                                    &
     1403       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'location',               &
    14181404                                   value = input_file_atts%location )
    1419        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1420                                    attribute_name = 'origin_x',                                    &
     1405       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_x',               &
    14211406                                   value = vmea(l)%origin_x_obs )
    1422        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1423                                    attribute_name = 'origin_y',                                    &
     1407       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_y',               &
    14241408                                   value = vmea(l)%origin_y_obs )
    1425        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1426                                    attribute_name = 'origin_lon',                                  &
     1409       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_lon',             &
    14271410                                   value = dum_lon )
    1428        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1429                                    attribute_name = 'origin_lat',                                  &
     1411       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_lat',             &
    14301412                                   value = dum_lat )
    1431        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1432                                    attribute_name = 'origin_z',                                    &
    1433                                    value = 0.0 )
    1434        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1435                                    attribute_name = 'rotation_angle',                              &
     1413       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'origin_z', value = 0.0 )
     1414       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'rotation_angle',         &
    14361415                                   value = input_file_atts%rotation_angle )
    1437        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1438                                    attribute_name = 'featureType',                                 &
     1416       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'featureType',            &
    14391417                                   value = TRIM( vmea(l)%feature_type_out ) )
    1440        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1441                                    attribute_name = 'data_content',                                &
     1418       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'data_content',           &
    14421419                                   value = TRIM( vmea(l)%data_content ) )
    1443        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1444                                    attribute_name = 'creation_time',                               &
     1420       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'creation_time',          &
    14451421                                   value = input_file_atts%creation_time )
    1446        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1447                                    attribute_name = 'version',                                     &
    1448                                    value = 1 ) !input_file_atts%version )
    1449        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1450                                    attribute_name = 'creation_time',                               &
     1422       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'version', value = 1 ) !input_file_atts%version
     1423       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'creation_time',          &
    14511424                                   value = TRIM( vmea(l)%site ) )
    1452        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1453                                    attribute_name = 'Conventions',                                 &
     1425       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'Conventions',            &
    14541426                                   value = input_file_atts%conventions )
    1455        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1456                                    attribute_name = 'dependencies',                                &
     1427       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'dependencies',           &
    14571428                                   value = input_file_atts%dependencies )
    1458        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1459                                    attribute_name = 'history',                                     &
     1429       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'history',                &
    14601430                                   value = input_file_atts%history )
    1461        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1462                                    attribute_name = 'references',                                  &
     1431       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'references',             &
    14631432                                   value = input_file_atts%references )
    1464        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1465                                    attribute_name = 'comment',                                     &
     1433       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'comment',                &
    14661434                                   value = input_file_atts%comment )
    1467        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1468                                    attribute_name = 'keywords',                                    &
     1435       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'keywords',               &
    14691436                                   value = input_file_atts%keywords )
    1470        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1471                                    attribute_name = 'licence',                                     &
     1437       return_value = dom_def_att( vmea(l)%nc_filename, attribute_name = 'licence',                &
    14721438                                   value = '[UC]2 Open Licence; see [UC]2 ' //                     &
    14731439                                           'data policy available at ' //                          &
     
    15841550!--    of attributes.
    15851551!--    Units
    1586        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1587                                    variable_name = 'time',                                         &
    1588                                    attribute_name = char_unit,                                     &
    1589                                    value = 'seconds since ' // origin_date_time )
    1590        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1591                                    variable_name = 'z',                                            &
    1592                                    attribute_name = char_unit,                                     &
    1593                                    value = 'm' )
    1594        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1595                                    variable_name = 'station_h',                                    &
    1596                                    attribute_name = char_unit,                                     &
    1597                                    value = 'm' )
    1598        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1599                                    variable_name = 'x',                                            &
    1600                                    attribute_name = char_unit,                                     &
    1601                                    value = 'm' )
    1602        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1603                                    variable_name = 'y',                                            &
    1604                                    attribute_name = char_unit,                                     &
    1605                                    value = 'm' )
    1606        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1607                                    variable_name = 'E_UTM',                                        &
    1608                                    attribute_name = char_unit,                                     &
    1609                                    value = 'm' )
    1610        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1611                                    variable_name = 'N_UTM',                                        &
    1612                                    attribute_name = char_unit,                                     &
    1613                                    value = 'm' )
    1614        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1615                                    variable_name = 'lat',                                          &
    1616                                    attribute_name = char_unit,                                     &
    1617                                    value = 'degrees_north' )
    1618        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1619                                    variable_name = 'lon',                                          &
    1620                                    attribute_name = char_unit,                                     &
    1621                                    value = 'degrees_east' )
     1552       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1553                                   attribute_name = char_unit, value = 'seconds since ' //         &
     1554                                   origin_date_time )
     1555       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z',                       &
     1556                                   attribute_name = char_unit, value = 'm' )
     1557       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h',               &
     1558                                   attribute_name = char_unit, value = 'm' )
     1559       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x',                       &
     1560                                   attribute_name = char_unit, value = 'm' )
     1561       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y',                       &
     1562                                   attribute_name = char_unit, value = 'm' )
     1563       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM',                   &
     1564                                   attribute_name = char_unit, value = 'm' )
     1565       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM',                   &
     1566                                   attribute_name = char_unit, value = 'm' )
     1567       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat',                     &
     1568                                   attribute_name = char_unit, value = 'degrees_north' )
     1569       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon',                     &
     1570                                   attribute_name = char_unit, value = 'degrees_east' )
    16221571!
    16231572!--    Long name
    1624        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1625                                    variable_name = 'station_name',                                 &
    1626                                    attribute_name = char_long,                                     &
    1627                                    value = 'station name')
    1628        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1629                                    variable_name = 'time',                                         &
    1630                                    attribute_name = char_long,                                     &
    1631                                    value = 'time')
    1632        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1633                                    variable_name = 'z',                                            &
    1634                                    attribute_name = char_long,                                     &
    1635                                    value = 'height above origin' )
    1636        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1637                                    variable_name = 'station_h',                                    &
    1638                                    attribute_name = char_long,                                     &
    1639                                    value = 'surface altitude' )
    1640        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1641                                    variable_name = 'x',                                            &
     1573       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name',            &
     1574                                   attribute_name = char_long, value = 'station name')
     1575       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1576                                   attribute_name = char_long, value = 'time')
     1577       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z',                       &
     1578                                   attribute_name = char_long, value = 'height above origin' )
     1579       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h',               &
     1580                                   attribute_name = char_long, value = 'surface altitude' )
     1581       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x',                       &
    16421582                                   attribute_name = char_long,                                     &
    16431583                                   value = 'distance to origin in x-direction')
    1644        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1645                                    variable_name = 'y',                                            &
     1584       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y',                       &
    16461585                                   attribute_name = char_long,                                     &
    16471586                                   value = 'distance to origin in y-direction')
    1648        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1649                                    variable_name = 'E_UTM',                                        &
    1650                                    attribute_name = char_long,                                     &
    1651                                    value = 'easting' )
    1652        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1653                                    variable_name = 'N_UTM',                                        &
    1654                                    attribute_name = char_long,                                     &
    1655                                    value = 'northing' )
    1656        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1657                                    variable_name = 'lat',                                          &
    1658                                    attribute_name = char_long,                                     &
    1659                                    value = 'latitude' )
    1660        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1661                                    variable_name = 'lon',                                          &
    1662                                    attribute_name = char_long,                                     &
    1663                                    value = 'longitude' )
     1587       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM',                   &
     1588                                   attribute_name = char_long, value = 'easting' )
     1589       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM',                   &
     1590                                   attribute_name = char_long, value = 'northing' )
     1591       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat',                     &
     1592                                   attribute_name = char_long, value = 'latitude' )
     1593       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon',                     &
     1594                                   attribute_name = char_long, value = 'longitude' )
    16641595!
    16651596!--    Standard name
    1666        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1667                                    variable_name = 'station_name',                                 &
    1668                                    attribute_name = char_standard,                                 &
    1669                                    value = 'platform_name')
    1670        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1671                                    variable_name = 'time',                                         &
    1672                                    attribute_name = char_standard,                                 &
    1673                                    value = 'time')
    1674        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1675                                    variable_name = 'z',                                            &
     1597       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name',            &
     1598                                   attribute_name = char_standard, value = 'platform_name')
     1599       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1600                                   attribute_name = char_standard, value = 'time')
     1601       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z',                       &
    16761602                                   attribute_name = char_standard,                                 &
    16771603                                   value = 'height_above_mean_sea_level' )
    1678        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1679                                    variable_name = 'station_h',                                    &
    1680                                    attribute_name = char_standard,                                 &
    1681                                    value = 'surface_altitude' )
    1682        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1683                                    variable_name = 'E_UTM',                                        &
     1604       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h',               &
     1605                                   attribute_name = char_standard, value = 'surface_altitude' )
     1606       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM',                   &
    16841607                                   attribute_name = char_standard,                                 &
    16851608                                   value = 'projection_x_coordinate' )
    1686        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1687                                    variable_name = 'N_UTM',                                        &
     1609       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM',                   &
    16881610                                   attribute_name = char_standard,                                 &
    16891611                                   value = 'projection_y_coordinate' )
    1690        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1691                                    variable_name = 'lat',                                          &
    1692                                    attribute_name = char_standard,                                 &
    1693                                    value = 'latitude' )
    1694        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1695                                    variable_name = 'lon',                                          &
    1696                                    attribute_name = char_standard,                                 &
    1697                                    value = 'longitude' )
     1612       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat',                     &
     1613                                   attribute_name = char_standard, value = 'latitude' )
     1614       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon',                     &
     1615                                   attribute_name = char_standard, value = 'longitude' )
    16981616!
    16991617!--    Axis
    1700        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1701                                    variable_name = 'time',                                         &
    1702                                    attribute_name = 'axis',                                        &
    1703                                    value = 'T')
    1704        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1705                                    variable_name = 'z',                                            &
    1706                                    attribute_name = 'axis',                                        &
    1707                                    value = 'Z' )
    1708        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1709                                    variable_name = 'x',                                            &
    1710                                    attribute_name = 'axis',                                        &
    1711                                    value = 'X' )
    1712        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1713                                    variable_name = 'y',                                            &
    1714                                    attribute_name = 'axis',                                        &
    1715                                    value = 'Y' )
     1618       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1619                                   attribute_name = 'axis', value = 'T')
     1620       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z',                       &
     1621                                   attribute_name = 'axis', value = 'Z' )
     1622       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x',                       &
     1623                                   attribute_name = 'axis', value = 'X' )
     1624       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y',                       &
     1625                                   attribute_name = 'axis', value = 'Y' )
    17161626!
    17171627!--    Set further individual attributes for the coordinate variables.
    17181628!--    For station name
    1719        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1720                                    variable_name = 'station_name',                                 &
    1721                                    attribute_name = 'cf_role',                                     &
    1722                                    value = 'timeseries_id' )
     1629       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name',            &
     1630                                   attribute_name = 'cf_role', value = 'timeseries_id' )
    17231631!
    17241632!--    For time
    1725        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1726                                    variable_name = 'time',                                         &
    1727                                    attribute_name = 'calendar',                                    &
    1728                                    value = 'proleptic_gregorian' )
    1729        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1730                                    variable_name = 'time',                                         &
    1731                                    attribute_name = 'bounds',                                      &
    1732                                    value = 'time_bounds' )
     1633       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1634                                   attribute_name = 'calendar', value = 'proleptic_gregorian' )
     1635       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time',                    &
     1636                                   attribute_name = 'bounds', value = 'time_bounds' )
    17331637!
    17341638!--    For vertical reference system
    1735        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1736                                    variable_name = 'vrs',                                          &
    1737                                    attribute_name = char_long,                                     &
    1738                                    value = 'vertical reference system' )
    1739        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1740                                    variable_name = 'vrs',                                          &
    1741                                    attribute_name = 'system_name',                                 &
    1742                                    value = 'DHHN2016' )
     1639       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'vrs',                     &
     1640                                   attribute_name = char_long, value = 'vertical reference system' )
     1641       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'vrs',                     &
     1642                                   attribute_name = 'system_name', value = 'DHHN2016' )
    17431643!
    17441644!--    For z
    1745        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1746                                    variable_name = 'z',                                            &
    1747                                    attribute_name = 'positive',                                    &
    1748                                    value = 'up' )
     1645       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z',                       &
     1646                                   attribute_name = 'positive', value = 'up' )
    17491647!
    17501648!--    For coordinate reference system
    1751        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1752                                    variable_name = 'crs',                                          &
    1753                                    attribute_name = 'epsg_code',                                   &
    1754                                    value = coord_ref_sys%epsg_code )
    1755        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1756                                    variable_name = 'crs',                                          &
     1649       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
     1650                                   attribute_name = 'epsg_code', value = coord_ref_sys%epsg_code )
     1651       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17571652                                   attribute_name = 'false_easting',                               &
    17581653                                   value = coord_ref_sys%false_easting )
    1759        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1760                                    variable_name = 'crs',                                          &
     1654       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17611655                                   attribute_name = 'false_northing',                              &
    17621656                                   value = coord_ref_sys%false_northing )
    1763        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1764                                    variable_name = 'crs',                                          &
     1657       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17651658                                   attribute_name = 'grid_mapping_name',                           &
    17661659                                   value = coord_ref_sys%grid_mapping_name )
    1767        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1768                                    variable_name = 'crs',                                          &
     1660       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17691661                                   attribute_name = 'inverse_flattening',                          &
    17701662                                   value = coord_ref_sys%inverse_flattening )
    1771        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1772                                    variable_name = 'crs',                                          &
     1663       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17731664                                   attribute_name = 'latitude_of_projection_origin',&
    17741665                                   value = coord_ref_sys%latitude_of_projection_origin )
    1775        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1776                                    variable_name = 'crs',                                          &
    1777                                    attribute_name = char_long,                                     &
    1778                                    value = coord_ref_sys%long_name )
    1779        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1780                                    variable_name = 'crs',                                          &
     1666       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
     1667                                   attribute_name = char_long, value = coord_ref_sys%long_name )
     1668       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17811669                                   attribute_name = 'longitude_of_central_meridian',               &
    17821670                                   value = coord_ref_sys%longitude_of_central_meridian )
    1783        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1784                                    variable_name = 'crs',                                          &
     1671       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17851672                                   attribute_name = 'longitude_of_prime_meridian',                 &
    17861673                                   value = coord_ref_sys%longitude_of_prime_meridian )
    1787        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1788                                    variable_name = 'crs',                                          &
     1674       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17891675                                   attribute_name = 'scale_factor_at_central_meridian',            &
    17901676                                   value = coord_ref_sys%scale_factor_at_central_meridian )
    1791        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1792                                    variable_name = 'crs',                                          &
     1677       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
    17931678                                   attribute_name = 'semi_major_axis',                             &
    17941679                                   value = coord_ref_sys%semi_major_axis )
    1795        return_value = dom_def_att( vmea(l)%nc_filename,                                            &
    1796                                    variable_name = 'crs',                                          &
    1797                                    attribute_name = char_unit,                                     &
    1798                                    value = coord_ref_sys%units )
     1680       return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'crs',                     &
     1681                                   attribute_name = char_unit, value = coord_ref_sys%units )
    17991682!
    18001683!--    In case of sampled soil quantities, define further dimensions and coordinates.
     
    18071690          ENDDO
    18081691
    1809           return_value = dom_def_dim( vmea(l)%nc_filename,                                         &
    1810                                       dimension_name = 'station_soil',                             &
     1692          return_value = dom_def_dim( vmea(l)%nc_filename, dimension_name = 'station_soil',        &
    18111693                                      output_type = 'int32',                                       &
    1812                                       bounds = (/1_iwp,vmea(l)%ns_soil_tot/),                      &
    1813                                       values_int32 = ndim )
     1694                                      bounds = (/1_iwp,vmea(l)%ns_soil_tot/), values_int32 = ndim )
    18141695          DEALLOCATE( ndim )
    18151696!
     
    18201701          ENDDO
    18211702
    1822           return_value = dom_def_dim( vmea(l)%nc_filename,                                         &
    1823                                       dimension_name = 'ntime_soil',                               &
    1824                                       output_type = 'int32',                                       &
    1825                                       bounds = (/1_iwp,ntimesteps/),                               &
     1703          return_value = dom_def_dim( vmea(l)%nc_filename, dimension_name = 'ntime_soil',          &
     1704                                      output_type = 'int32', bounds = (/1_iwp,ntimesteps/),        &
    18261705                                      values_int32 = ndim )
    18271706          DEALLOCATE( ndim )
     
    18291708!--       time for soil
    18301709          variable_name = 'time_soil'
    1831           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1832                                       variable_name = variable_name,                               &
    1833                                       dimension_names = (/'station_soil',                          &
    1834                                                           'ntime_soil  '/),                        &
     1710          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1711                                      dimension_names = (/'station_soil', 'ntime_soil  '/),        &
    18351712                                      output_type = 'real32' )
    18361713!
    18371714!--       station_name for soil
    18381715          variable_name = 'station_name_soil'
    1839           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1840                                       variable_name = variable_name,                               &
    1841                                       dimension_names = (/ 'max_name_len',                         &
    1842                                                            'station_soil' /),                      &
     1716          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1717                                      dimension_names = (/ 'max_name_len', 'station_soil' /),      &
    18431718                                      output_type = 'char' )
    18441719!
    18451720!--       z
    18461721          variable_name = 'z_soil'
    1847           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1848                                       variable_name = variable_name,                               &
    1849                                       dimension_names = (/'station_soil'/),                        &
    1850                                       output_type = 'real32' )
     1722          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1723                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18511724!
    18521725!--       station_h for soil
    18531726          variable_name = 'station_h_soil'
    1854           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1855                                       variable_name = variable_name,                               &
    1856                                       dimension_names = (/'station_soil'/),                        &
    1857                                       output_type = 'real32' )
     1727          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1728                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18581729!
    18591730!--       x soil
    18601731          variable_name = 'x_soil'
    1861           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1862                                       variable_name = variable_name,                               &
    1863                                       dimension_names = (/'station_soil'/),                        &
    1864                                       output_type = 'real32' )
     1732          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1733                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18651734!
    18661735!-        y soil
    18671736          variable_name = 'y_soil'
    1868           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1869                                       variable_name = variable_name,                               &
    1870                                       dimension_names = (/'station_soil'/),                        &
    1871                                       output_type = 'real32' )
     1737          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1738                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18721739!
    18731740!--       E-UTM soil
    18741741          variable_name = 'E_UTM_soil'
    1875           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1876                                       variable_name = variable_name,                               &
    1877                                       dimension_names = (/'station_soil'/),                        &
    1878                                       output_type = 'real32' )
     1742          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1743                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18791744!
    18801745!--       N-UTM soil
    18811746          variable_name = 'N_UTM_soil'
    1882           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1883                                       variable_name = variable_name,                               &
    1884                                       dimension_names = (/'station_soil'/),                        &
    1885                                       output_type = 'real32' )
     1747          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1748                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18861749!
    18871750!--       latitude soil
    18881751          variable_name = 'lat_soil'
    1889           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1890                                       variable_name = variable_name,                               &
    1891                                       dimension_names = (/'station_soil'/),                        &
    1892                                       output_type = 'real32' )
     1752          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1753                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    18931754!
    18941755!--       longitude soil
    18951756          variable_name = 'lon_soil'
    1896           return_value = dom_def_var( vmea(l)%nc_filename,                                         &
    1897                                       variable_name = variable_name,                               &
    1898                                       dimension_names = (/'station_soil'/),                        &
    1899                                       output_type = 'real32' )
     1757          return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,          &
     1758                                      dimension_names = (/'station_soil'/), output_type = 'real32' )
    19001759!
    19011760!--       Set attributes for the coordinate variables. Note, not all coordinates have the same
    19021761!--       number of attributes.
    19031762!--       Units
    1904           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1905                                       variable_name = 'time_soil',                                 &
    1906                                       attribute_name = char_unit,                                  &
    1907                                       value = 'seconds since ' // origin_date_time )
    1908           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1909                                       variable_name = 'z_soil',                                    &
    1910                                       attribute_name = char_unit,                                  &
    1911                                       value = 'm' )
    1912           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1913                                       variable_name = 'station_h_soil',                            &
    1914                                       attribute_name = char_unit,                                  &
    1915                                       value = 'm' )
    1916           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1917                                       variable_name = 'x_soil',                                    &
    1918                                       attribute_name = char_unit,                                  &
    1919                                       value = 'm' )
    1920           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1921                                       variable_name = 'y_soil',                                    &
    1922                                       attribute_name = char_unit,                                  &
    1923                                       value = 'm' )
    1924           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1925                                       variable_name = 'E_UTM_soil',                                &
    1926                                       attribute_name = char_unit,                                  &
    1927                                       value = 'm' )
    1928           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1929                                       variable_name = 'N_UTM_soil',                                &
    1930                                       attribute_name = char_unit,                                  &
    1931                                       value = 'm' )
    1932           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1933                                       variable_name = 'lat_soil',                                  &
    1934                                       attribute_name = char_unit,                                  &
    1935                                       value = 'degrees_north' )
    1936           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1937                                       variable_name = 'lon_soil',                                  &
    1938                                       attribute_name = char_unit,                                  &
    1939                                       value = 'degrees_east' )
     1763          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1764                                      attribute_name = char_unit, value = 'seconds since ' //      &
     1765                                      origin_date_time )
     1766          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z_soil',               &
     1767                                      attribute_name = char_unit, value = 'm' )
     1768          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h_soil',       &
     1769                                      attribute_name = char_unit, value = 'm' )
     1770          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x_soil',               &
     1771                                      attribute_name = char_unit, value = 'm' )
     1772          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y_soil',               &
     1773                                      attribute_name = char_unit, value = 'm' )
     1774          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM_soil',           &
     1775                                      attribute_name = char_unit, value = 'm' )
     1776          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM_soil',           &
     1777                                      attribute_name = char_unit, value = 'm' )
     1778          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat_soil',             &
     1779                                      attribute_name = char_unit, value = 'degrees_north' )
     1780          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon_soil',             &
     1781                                      attribute_name = char_unit, value = 'degrees_east' )
    19401782!
    19411783!--       Long name
    1942           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1943                                       variable_name = 'station_name_soil',                         &
    1944                                       attribute_name = char_long,                                  &
    1945                                       value = 'station name')
    1946           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1947                                       variable_name = 'time_soil',                                 &
    1948                                       attribute_name = char_long,                                  &
    1949                                       value = 'time')
    1950           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1951                                       variable_name = 'z_soil',                                    &
    1952                                       attribute_name = char_long,                                  &
    1953                                       value = 'height above origin' )
    1954           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1955                                       variable_name = 'station_h_soil',                            &
    1956                                       attribute_name = char_long,                                  &
    1957                                       value = 'surface altitude' )
    1958           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1959                                       variable_name = 'x_soil',                                    &
     1784          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name_soil',    &
     1785                                      attribute_name = char_long, value = 'station name')
     1786          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1787                                      attribute_name = char_long, value = 'time')
     1788          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z_soil',               &
     1789                                      attribute_name = char_long, value = 'height above origin' )
     1790          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h_soil',       &
     1791                                      attribute_name = char_long, value = 'surface altitude' )
     1792          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x_soil',               &
    19601793                                      attribute_name = char_long,                                  &
    19611794                                      value = 'distance to origin in x-direction' )
    1962           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1963                                       variable_name = 'y_soil',                                    &
     1795          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y_soil',               &
    19641796                                      attribute_name = char_long,                                  &
    19651797                                      value = 'distance to origin in y-direction' )
    1966           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1967                                       variable_name = 'E_UTM_soil',                                &
    1968                                       attribute_name = char_long,                                  &
    1969                                       value = 'easting' )
    1970           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1971                                       variable_name = 'N_UTM_soil',                                &
    1972                                       attribute_name = char_long,                                  &
    1973                                       value = 'northing' )
    1974           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1975                                       variable_name = 'lat_soil',                                  &
    1976                                       attribute_name = char_long,                                  &
    1977                                       value = 'latitude' )
    1978           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1979                                       variable_name = 'lon_soil',                                  &
    1980                                       attribute_name = char_long,                                  &
    1981                                       value = 'longitude' )
     1798          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM_soil',           &
     1799                                      attribute_name = char_long, value = 'easting' )
     1800          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM_soil',           &
     1801                                      attribute_name = char_long, value = 'northing' )
     1802          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat_soil',             &
     1803                                      attribute_name = char_long, value = 'latitude' )
     1804          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon_soil',             &
     1805                                      attribute_name = char_long, value = 'longitude' )
    19821806!
    19831807!--       Standard name
    1984           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1985                                       variable_name = 'station_name_soil',                         &
    1986                                       attribute_name = char_standard,                              &
    1987                                       value = 'platform_name')
    1988           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1989                                       variable_name = 'time_soil',                                 &
    1990                                       attribute_name = char_standard,                              &
    1991                                       value = 'time')
    1992           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1993                                       variable_name = 'z_soil',                                    &
     1808          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name_soil',    &
     1809                                      attribute_name = char_standard, value = 'platform_name')
     1810          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1811                                      attribute_name = char_standard, value = 'time')
     1812          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z_soil',               &
    19941813                                      attribute_name = char_standard,                              &
    19951814                                      value = 'height_above_mean_sea_level' )
    1996           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    1997                                       variable_name = 'station_h_soil',                            &
    1998                                       attribute_name = char_standard,                              &
    1999                                       value = 'surface_altitude' )
    2000           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2001                                       variable_name = 'E_UTM_soil',                                &
     1815          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_h_soil',       &
     1816                                      attribute_name = char_standard, value = 'surface_altitude' )
     1817          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'E_UTM_soil',           &
    20021818                                      attribute_name = char_standard,                              &
    20031819                                      value = 'projection_x_coordinate' )
    2004           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2005                                       variable_name = 'N_UTM_soil',                                &
     1820          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'N_UTM_soil',           &
    20061821                                      attribute_name = char_standard,                              &
    20071822                                      value = 'projection_y_coordinate' )
    2008           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2009                                       variable_name = 'lat_soil',                                  &
    2010                                       attribute_name = char_standard,                              &
    2011                                       value = 'latitude' )
    2012           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2013                                       variable_name = 'lon_soil',                                  &
    2014                                       attribute_name = char_standard,                              &
    2015                                       value = 'longitude' )
     1823          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lat_soil',             &
     1824                                      attribute_name = char_standard, value = 'latitude' )
     1825          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'lon_soil',             &
     1826                                      attribute_name = char_standard, value = 'longitude' )
    20161827!
    20171828!--       Axis
    2018           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2019                                       variable_name = 'time_soil',                                 &
    2020                                       attribute_name = 'axis',                                     &
    2021                                       value = 'T')
    2022           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2023                                       variable_name = 'z_soil',                                    &
    2024                                       attribute_name = 'axis',                                     &
    2025                                       value = 'Z' )
    2026           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2027                                       variable_name = 'x_soil',                                    &
    2028                                       attribute_name = 'axis',                                     &
    2029                                       value = 'X' )
    2030           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2031                                       variable_name = 'y_soil',                                    &
    2032                                       attribute_name = 'axis',                                     &
    2033                                       value = 'Y' )
     1829          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1830                                      attribute_name = 'axis', value = 'T')
     1831          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z_soil',               &
     1832                                      attribute_name = 'axis', value = 'Z' )
     1833          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'x_soil',               &
     1834                                      attribute_name = 'axis', value = 'X' )
     1835          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'y_soil',               &
     1836                                      attribute_name = 'axis', value = 'Y' )
    20341837!
    20351838!--       Set further individual attributes for the coordinate variables.
    20361839!--       For station name soil
    2037           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2038                                       variable_name = 'station_name_soil',                         &
    2039                                       attribute_name = 'cf_role',                                  &
    2040                                       value = 'timeseries_id' )
     1840          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'station_name_soil',    &
     1841                                      attribute_name = 'cf_role', value = 'timeseries_id' )
    20411842!
    20421843!--       For time soil
    2043           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2044                                       variable_name = 'time_soil',                                 &
    2045                                       attribute_name = 'calendar',                                 &
    2046                                       value = 'proleptic_gregorian' )
    2047           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2048                                       variable_name = 'time_soil',                                 &
    2049                                       attribute_name = 'bounds',                                   &
    2050                                       value = 'time_bounds' )
     1844          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1845                                      attribute_name = 'calendar', value = 'proleptic_gregorian' )
     1846          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'time_soil',            &
     1847                                      attribute_name = 'bounds', value = 'time_bounds' )
    20511848!
    20521849!--       For z soil
    2053           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2054                                       variable_name = 'z_soil',                                    &
    2055                                       attribute_name = 'positive',                                 &
    2056                                       value = 'up' )
     1850          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = 'z_soil',               &
     1851                                      attribute_name = 'positive', value = 'up' )
    20571852       ENDIF
    20581853!
     
    20661861               ANY( TRIM( vmea(l)%var_atts(n)%name) == soil_vars ) )  THEN
    20671862
    2068              return_value = dom_def_var( vmea(l)%nc_filename,                                      &
    2069                                          variable_name = variable_name,                            &
    2070                                          dimension_names = (/'station_soil',                       &
    2071                                                              'ntime_soil  '/),                     &
     1863             return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,       &
     1864                                         dimension_names = (/'station_soil', 'ntime_soil  '/),     &
    20721865                                         output_type = 'real32' )
    20731866          ELSE
    20741867
    2075              return_value = dom_def_var( vmea(l)%nc_filename,                                      &
    2076                                          variable_name = variable_name,                            &
    2077                                          dimension_names = (/'station',                            &
    2078                                                              'ntime  '/),                          &
     1868             return_value = dom_def_var( vmea(l)%nc_filename, variable_name = variable_name,       &
     1869                                         dimension_names = (/'station', 'ntime  '/),               &
    20791870                                         output_type = 'real32' )
    20801871          ENDIF
     
    20851876
    20861877          IF ( vmea(l)%var_atts(n)%long_name /= 'none' )  THEN
    2087              return_value = dom_def_att( vmea(l)%nc_filename,                                      &
    2088                                          variable_name = variable_name,                            &
     1878             return_value = dom_def_att( vmea(l)%nc_filename,  variable_name = variable_name,      &
    20891879                                         attribute_name = char_long,                               &
    20901880                                         value = TRIM( vmea(l)%var_atts(n)%long_name ) )
    20911881          ENDIF
    20921882          IF ( vmea(l)%var_atts(n)%standard_name /= 'none' )  THEN
    2093              return_value = dom_def_att( vmea(l)%nc_filename,                                      &
    2094                                          variable_name = variable_name,                            &
     1883             return_value = dom_def_att( vmea(l)%nc_filename, variable_name = variable_name,       &
    20951884                                         attribute_name = char_standard,                           &
    20961885                                         value = TRIM( vmea(l)%var_atts(n)%standard_name ) )
    20971886          ENDIF
    20981887          IF ( vmea(l)%var_atts(n)%units /= 'none' )  THEN
    2099              return_value = dom_def_att( vmea(l)%nc_filename,                                      &
    2100                                          variable_name = variable_name,                            &
     1888             return_value = dom_def_att( vmea(l)%nc_filename, variable_name = variable_name,       &
    21011889                                         attribute_name = char_unit,                               &
    21021890                                         value = TRIM( vmea(l)%var_atts(n)%units ) )
    21031891          ENDIF
    21041892
    2105           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2106                                       variable_name = variable_name,                               &
     1893          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = variable_name,          &
    21071894                                      attribute_name = 'grid_mapping',                             &
    21081895                                      value = TRIM( vmea(l)%var_atts(n)%grid_mapping ) )
    21091896
    2110           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2111                                       variable_name = variable_name,                               &
     1897          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = variable_name,          &
    21121898                                      attribute_name = 'coordinates',                              &
    21131899                                      value = TRIM( vmea(l)%var_atts(n)%coordinates ) )
    21141900
    2115           return_value = dom_def_att( vmea(l)%nc_filename,                                         &
    2116                                       variable_name = variable_name,                               &
     1901          return_value = dom_def_att( vmea(l)%nc_filename, variable_name = variable_name,          &
    21171902                                      attribute_name = char_fill,                                  &
    21181903                                      value = REAL( vmea(l)%var_atts(n)%fill_value, KIND=4 ) )
     
    21721957          output_values_1d_pointer => output_values_1d_target
    21731958
    2174           return_value =                                                                           &
    2175                   dom_write_var( vmea(l)%nc_filename,                                              &
    2176                                  'E_UTM',                                                          &
    2177                                  values_realwp_1d = output_values_1d_pointer,                      &
    2178                                  bounds_start = (/vmea(l)%start_coord_a/),                         &
    2179                                  bounds_end   = (/vmea(l)%end_coord_a  /) )
     1959          return_value = dom_write_var( vmea(l)%nc_filename, 'E_UTM',                              &
     1960                                        values_realwp_1d = output_values_1d_pointer,               &
     1961                                        bounds_start = (/vmea(l)%start_coord_a/),                  &
     1962                                        bounds_end   = (/vmea(l)%end_coord_a  /) )
    21801963!
    21811964!--       Output of Northing coordinate. Before output, recalculate NUTM.
     
    21871970
    21881971          output_values_1d_pointer => output_values_1d_target
    2189           return_value =                                                                           &
    2190                   dom_write_var( vmea(l)%nc_filename,                                              &
    2191                                  'N_UTM',                                                          &
    2192                                  values_realwp_1d = output_values_1d_pointer,                      &
    2193                                  bounds_start = (/vmea(l)%start_coord_a/),                         &
    2194                                  bounds_end   = (/vmea(l)%end_coord_a  /) )
     1972          return_value = dom_write_var( vmea(l)%nc_filename, 'N_UTM',                              &
     1973                                        values_realwp_1d = output_values_1d_pointer,               &
     1974                                        bounds_start = (/vmea(l)%start_coord_a/),                  &
     1975                                        bounds_end   = (/vmea(l)%end_coord_a  /) )
    21951976!
    21961977!--       Output of relative height coordinate.
     
    22031984          output_values_1d_target = vmea(l)%zar(1:vmea(l)%ns) + oro_rel(:)
    22041985          output_values_1d_pointer => output_values_1d_target
    2205           return_value =                                                                           &
    2206                   dom_write_var( vmea(l)%nc_filename,                                              &
    2207                                  'z',                                                              &
    2208                                  values_realwp_1d = output_values_1d_pointer,                      &
    2209                                  bounds_start = (/vmea(l)%start_coord_a/),                         &
    2210                                  bounds_end   = (/vmea(l)%end_coord_a  /) )
     1986          return_value = dom_write_var( vmea(l)%nc_filename, 'z',                                  &
     1987                                        values_realwp_1d = output_values_1d_pointer,               &
     1988                                        bounds_start = (/vmea(l)%start_coord_a/),                  &
     1989                                        bounds_end   = (/vmea(l)%end_coord_a  /) )
    22111990!
    22121991!--       Write surface altitude for the station. Note, since z is already a relative observation
     
    22141993          output_values_1d_target = oro_rel(:)
    22151994          output_values_1d_pointer => output_values_1d_target
    2216           return_value =                                                                           &
    2217                   dom_write_var( vmea(l)%nc_filename,                                              &
    2218                                  'station_h',                                                      &
    2219                                  values_realwp_1d = output_values_1d_pointer,                      &
    2220                                  bounds_start = (/vmea(l)%start_coord_a/),                         &
    2221                                  bounds_end   = (/vmea(l)%end_coord_a  /) )
     1995          return_value = dom_write_var( vmea(l)%nc_filename, 'station_h',                          &
     1996                                        values_realwp_1d = output_values_1d_pointer,               &
     1997                                        bounds_start = (/vmea(l)%start_coord_a/),                  &
     1998                                        bounds_end   = (/vmea(l)%end_coord_a  /) )
    22221999
    22232000          DEALLOCATE( oro_rel )
     
    22392016          output_values_2d_char_pointer => output_values_2d_char_target
    22402017
    2241           return_value =                                                                           &
    2242                   dom_write_var( vmea(l)%nc_filename,                                              &
    2243                                  'station_name',                                                   &
    2244                                  values_char_2d = output_values_2d_char_pointer,                   &
    2245                                  bounds_start = (/ 1,                   vmea(l)%start_coord_a /),  &
    2246                                  bounds_end   = (/ maximum_name_length, vmea(l)%end_coord_a /) )
     2018          return_value = dom_write_var( vmea(l)%nc_filename, 'station_name',                       &
     2019                                        values_char_2d = output_values_2d_char_pointer,            &
     2020                                        bounds_start = (/ 1, vmea(l)%start_coord_a /),             &
     2021                                        bounds_end   = (/ maximum_name_length,                     &
     2022                                        vmea(l)%end_coord_a /) )
    22472023
    22482024          DEALLOCATE( station_name )
     
    22602036               * SIN( init_model%rotation_angle * pi / 180.0_wp )
    22612037             output_values_1d_pointer => output_values_1d_target
    2262              return_value =                                                                        &
    2263                   dom_write_var( vmea(l)%nc_filename,                                              &
    2264                                  'E_UTM_soil',                                                     &
    2265                                  values_realwp_1d = output_values_1d_pointer,                      &
    2266                                  bounds_start = (/vmea(l)%start_coord_s/),                         &
    2267                                  bounds_end   = (/vmea(l)%end_coord_s  /) )
     2038             return_value = dom_write_var( vmea(l)%nc_filename, 'E_UTM_soil',                      &
     2039                                           values_realwp_1d = output_values_1d_pointer,            &
     2040                                           bounds_start = (/vmea(l)%start_coord_s/),               &
     2041                                           bounds_end   = (/vmea(l)%end_coord_s  /) )
    22682042!
    22692043!--          Output of Northing coordinate. Before output, recalculate NUTM.
     
    22752049
    22762050             output_values_1d_pointer => output_values_1d_target
    2277              return_value =                                                                        &
    2278                   dom_write_var( vmea(l)%nc_filename,                                              &
    2279                                  'N_UTM_soil',                                                     &
    2280                                  values_realwp_1d = output_values_1d_pointer,                      &
    2281                                  bounds_start = (/vmea(l)%start_coord_s/),                         &
    2282                                  bounds_end   = (/vmea(l)%end_coord_s  /) )
     2051             return_value = dom_write_var( vmea(l)%nc_filename, 'N_UTM_soil',                      &
     2052                                           values_realwp_1d = output_values_1d_pointer,            &
     2053                                           bounds_start = (/vmea(l)%start_coord_s/),               &
     2054                                           bounds_end   = (/vmea(l)%end_coord_s  /) )
    22832055!
    22842056!--          Output of relative height coordinate.
     
    22912063             output_values_1d_target = vmea(l)%depth(1:vmea(l)%ns_soil) + oro_rel(:)
    22922064             output_values_1d_pointer => output_values_1d_target
    2293              return_value =                                                                        &
    2294                   dom_write_var( vmea(l)%nc_filename,                                              &
    2295                                  'z_soil',                                                         &
    2296                                  values_realwp_1d = output_values_1d_pointer,                      &
    2297                                  bounds_start = (/vmea(l)%start_coord_s/),                         &
    2298                                  bounds_end   = (/vmea(l)%end_coord_s  /) )
     2065             return_value = dom_write_var( vmea(l)%nc_filename, 'z_soil',                          &
     2066                                           values_realwp_1d = output_values_1d_pointer,            &
     2067                                           bounds_start = (/vmea(l)%start_coord_s/),               &
     2068                                           bounds_end   = (/vmea(l)%end_coord_s  /) )
    22992069!
    23002070!--          Write surface altitude for the station. Note, since z is already a relative observation
     
    23022072             output_values_1d_target = oro_rel(:)
    23032073             output_values_1d_pointer => output_values_1d_target
    2304              return_value =                                                                        &
    2305                   dom_write_var( vmea(l)%nc_filename,                                              &
    2306                                  'station_h_soil',                                                 &
    2307                                  values_realwp_1d = output_values_1d_pointer,                      &
    2308                                  bounds_start = (/vmea(l)%start_coord_s/),                         &
    2309                                  bounds_end   = (/vmea(l)%end_coord_s  /) )
     2074             return_value = dom_write_var( vmea(l)%nc_filename, 'station_h_soil',                  &
     2075                                           values_realwp_1d = output_values_1d_pointer,            &
     2076                                           bounds_start = (/vmea(l)%start_coord_s/),               &
     2077                                           bounds_end   = (/vmea(l)%end_coord_s  /) )
    23102078
    23112079             DEALLOCATE( oro_rel )
     
    23262094             output_values_2d_char_pointer => output_values_2d_char_target
    23272095
    2328              return_value =                                                                        &
    2329                   dom_write_var( vmea(l)%nc_filename,                                              &
    2330                                  'station_name_soil',                                              &
    2331                                  values_char_2d = output_values_2d_char_pointer,                   &
    2332                                  bounds_start = (/ 1,                   vmea(l)%start_coord_s /),  &
    2333                                  bounds_end   = (/ maximum_name_length, vmea(l)%end_coord_s   /) )
     2096             return_value = dom_write_var( vmea(l)%nc_filename, 'station_name_soil',               &
     2097                                           values_char_2d = output_values_2d_char_pointer,         &
     2098                                           bounds_start = (/ 1, vmea(l)%start_coord_s /),          &
     2099                                           bounds_end   = (/ maximum_name_length,                  &
     2100                                           vmea(l)%end_coord_s   /) )
    23342101
    23352102             DEALLOCATE( station_name )
     
    23632130             output_values_2d_pointer => output_values_2d_target
    23642131
    2365              return_value = dom_write_var( vmea(l)%nc_filename,                                    &
    2366                                            variable_name,                                          &
    2367                                            values_realwp_2d = output_values_2d_pointer, &
     2132             return_value = dom_write_var( vmea(l)%nc_filename, variable_name,                     &
     2133                                           values_realwp_2d = output_values_2d_pointer,            &
    23682134                                           bounds_start = (/vmea(l)%start_coord_s, t_ind/),        &
    23692135                                           bounds_end   = (/vmea(l)%end_coord_s, t_ind /) )
     
    23722138             output_values_2d_target(t_ind,:) = vmea(l)%measured_vars_soil(:,n)
    23732139             output_values_2d_pointer => output_values_2d_target
    2374              return_value =                                                                        &
    2375                       dom_write_var( vmea(l)%nc_filename,                                          &
    2376                                      variable_name,                                                &
    2377                                      values_realwp_2d = output_values_2d_pointer,                  &
    2378                                      bounds_start = (/vmea(l)%start_coord_s, t_ind/),              &
    2379                                      bounds_end   = (/vmea(l)%end_coord_s, t_ind  /) )
     2140             return_value = dom_write_var( vmea(l)%nc_filename, variable_name,                     &
     2141                                           values_realwp_2d = output_values_2d_pointer,            &
     2142                                           bounds_start = (/vmea(l)%start_coord_s, t_ind/),        &
     2143                                           bounds_end   = (/vmea(l)%end_coord_s, t_ind  /) )
    23802144             DEALLOCATE( output_values_2d_target )
    23812145          ELSE
     
    23872151             output_values_2d_pointer => output_values_2d_target
    23882152
    2389              return_value = dom_write_var( vmea(l)%nc_filename,                                    &
    2390                                            variable_name,                                          &
     2153             return_value = dom_write_var( vmea(l)%nc_filename, variable_name,                     &
    23912154                                           values_realwp_2d = output_values_2d_pointer,            &
    23922155                                           bounds_start = (/vmea(l)%start_coord_a, t_ind/),        &
     
    23972160             output_values_2d_target(t_ind,:) = vmea(l)%measured_vars(:,n)
    23982161             output_values_2d_pointer => output_values_2d_target
    2399              return_value =                                                                        &
    2400                       dom_write_var( vmea(l)%nc_filename,                                          &
    2401                                      variable_name,                                                &
    2402                                      values_realwp_2d = output_values_2d_pointer,                  &
    2403                                      bounds_start = (/ vmea(l)%start_coord_a, t_ind /),            &
    2404                                      bounds_end   = (/ vmea(l)%end_coord_a, t_ind /) )
     2162             return_value = dom_write_var( vmea(l)%nc_filename, variable_name,                     &
     2163                                           values_realwp_2d = output_values_2d_pointer,            &
     2164                                           bounds_start = (/ vmea(l)%start_coord_a, t_ind /),      &
     2165                                           bounds_end   = (/ vmea(l)%end_coord_a, t_ind /) )
    24052166
    24062167             DEALLOCATE( output_values_2d_target )
     
    25632324                    j = vmea(l)%j(m)
    25642325                    i = vmea(l)%i(m)
    2565                     vmea(l)%measured_vars(m,n) = SQRT(                                             &
    2566                                                  ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) )**2 +       &
    2567                                                  ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) )**2         &
     2326                    vmea(l)%measured_vars(m,n) = SQRT(   ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) )**2 &
     2327                                                       + ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) )**2 &
    25682328                                                     )
    25692329                 ENDDO
     
    25762336
    25772337                    vmea(l)%measured_vars(m,n) = 180.0_wp + 180.0_wp / pi * ATAN2(                 &
    2578                                                  0.5_wp * ( v(k,j,i) + v(k,j+1,i) ),              &
    2579                                                  0.5_wp * ( u(k,j,i) + u(k,j,i+1) )                &
    2580                                                                                   )
     2338                                                               0.5_wp * ( v(k,j,i) + v(k,j+1,i) ), &
     2339                                                               0.5_wp * ( u(k,j,i) + u(k,j,i+1) )  &
     2340                                                                                 )
    25812341                 ENDDO
    25822342
Note: See TracChangeset for help on using the changeset viewer.