Ignore:
Timestamp:
Nov 20, 2018 10:53:14 AM (4 years ago)
Author:
eckhard
Message:

inifor: COSMO-D2 support

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/UTIL/inifor/src/inifor_grid.f90

    r3456 r3537  
    2626! -----------------
    2727! $Id$
     28! Read COSMO domain extents and soil depths from input files
     29! Report averaging mode and debugging mode in log
     30!
     31!
     32! 3456 2018-10-30 14:29:54Z eckhard
    2833! Remove surface forcing from netCDF output (radiation + precipitation)
    2934! NetCDf output of internal arrays only with --debug option
     
    8388               RHO_L, OMEGA, HECTO
    8489    USE io,                                                                    &
    85         ONLY:  get_netcdf_variable, get_netcdf_attribute,                      &
    86                parse_command_line_arguments, validate_config
     90        ONLY:  get_netcdf_attribute, get_netcdf_dim_vector,                    &
     91               get_netcdf_variable, parse_command_line_arguments,              &
     92               validate_config
    8793    USE netcdf,                                                                &
    8894        ONLY:  NF90_MAX_NAME, NF90_MAX_VAR_DIMS
     
    154160    REAL(dp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  hfl             !< heights of full layers (cell centres) above sea level in COSMO-DE, computed as arithmetic average of hhl
    155161    REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  depths          !< COSMO-DE's TERRA-ML soil layer depths
    156     REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  d_depth_rho_inv !< COSMO-DE's TERRA-ML soil layer thicknesses
     162    REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  d_depth         !< COSMO-DE's TERRA-ML soil layer thicknesses
     163    REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  d_depth_rho_inv !< inverted soil water mass
    157164    REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  rlon            !< longitudes of COSMO-DE's rotated-pole grid
    158165    REAL(dp), DIMENSION(:), ALLOCATABLE, TARGET     ::  rlat            !< latitudes of COSMO-DE's rotated-pole grid
     
    177184    INTEGER ::  nlat  !< number of latitudal points in target grid (COSMO-DE)
    178185    INTEGER ::  nlev  !< number of levels in target grid (COSMO-DE)
    179     INTEGER ::  layers !< number of COSMO-DE soil layers
     186    INTEGER ::  ndepths !< number of COSMO-DE soil layers
    180187    INTEGER ::  start_hour_flow         !< start of flow forcing in number of hours relative to start_date
    181188    INTEGER ::  start_hour_soil         !< start of soil forcing in number of hours relative to start_date, typically equals start_hour_flow
     
    293300       start_hour_soilmoisture = - (4 * 7 * 24) - 2
    294301
    295        lonmin =  -5.0_dp * TO_RADIANS
    296        lonmax =   5.5_dp * TO_RADIANS
    297        latmin =  -5.0_dp * TO_RADIANS
    298        latmax =   6.5_dp * TO_RADIANS
    299 
    300302       ! COSMO-DE default rotated pole
    301303       phi_n     =   40.0_dp * TO_RADIANS
    302304       phi_equat =   50.0_dp * TO_RADIANS
    303305       lambda_n  = -170.0_dp * TO_RADIANS
    304 
    305        ! COMSMO-DE soil layers
    306        layers = 8   
    307        ALLOCATE( depths(layers), d_depth_rho_inv(layers) )
    308        depths = (/0.005_dp, 0.02_dp, 0.06_dp, 0.18_dp, 0.54_dp, 1.62_dp, 4.86_dp, 14.58_dp/)
    309        d_depth_rho_inv = 1.0_dp / &
    310           ( (/0.01_dp, 0.02_dp, 0.06_dp, 0.18_dp, 0.54_dp, 1.62_dp, 4.86_dp, 14.58_dp/) * RHO_L )
    311306
    312307       ! Defaultmain centre (_c) of the PALM-4U grid in the geographical system (_g)
     
    395390       CALL report('setup_parameters', "initialization mode: " // TRIM(cfg % ic_mode))
    396391       CALL report('setup_parameters', "       forcing mode: " // TRIM(cfg % bc_mode))
     392       CALL report('setup_parameters', "     averaging mode: " // TRIM(cfg % averaging_mode))
    397393       CALL report('setup_parameters', "          data path: " // TRIM(cfg % input_path))
    398394       CALL report('setup_parameters', "           hhl file: " // TRIM(cfg % hhl_file))
     
    400396       CALL report('setup_parameters', "      namelist file: " // TRIM(cfg % namelist_file))
    401397       CALL report('setup_parameters', "   output data file: " // TRIM(output_file % name))
     398       IF (cfg % debug )  CALL report('setup_parameters', "     debugging mode: enabled")
    402399
    403400 CALL run_control('time', 'init')
     
    472469       cosmo_var % name = 'HHL'
    473470       CALL get_netcdf_variable(cfg % hhl_file, cosmo_var, hhl)
     471       CALL get_netcdf_dim_vector(cfg % hhl_file, 'rlon', rlon)
     472       CALL get_netcdf_dim_vector(cfg % hhl_file, 'rlat', rlat)
     473       CALL get_netcdf_dim_vector(soil_files(1), 'depth_2', depths)
    474474 CALL run_control('time', 'read')
    475475
     
    478478       nlat = SIZE(hhl, 2)
    479479       nlev = SIZE(hhl, 3)
    480 
     480       ndepths = SIZE(depths)
     481
     482       lonmin = MINVAL(rlon) * TO_RADIANS
     483       lonmax = MAXVAL(rlon) * TO_RADIANS
     484       latmin = MINVAL(rlat) * TO_RADIANS
     485       latmax = MAXVAL(rlat) * TO_RADIANS
    481486 CALL run_control('time', 'comp')
    482487
    483488       ! Appoximate COSMO-DE heights of full layers (cell centres)
    484489       ALLOCATE( hfl(nlon, nlat, nlev-1) )
     490       ALLOCATE( d_depth(ndepths), d_depth_rho_inv(ndepths) )
     491
    485492 CALL run_control('time', 'alloc')
     493       CALL get_soil_layer_thickness( depths, d_depth )
     494       d_depth_rho_inv = 1.0_dp / ( d_depth * RHO_L )
     495
     496       ! Appoximate COSMO-DE heights of full layers (cell centres)
    486497       DO k = 1, nlev-1
    487498          hfl(:,:,k) = 0.5_dp * ( hhl(:,:,k) +                                 &
    488499                                  hhl(:,:,k+1) )
    489500       END DO
     501 CALL run_control('time', 'comp')
     502
     503
    490504
    491505!------------------------------------------------------------------------------
     
    14891503       END IF
    14901504
    1491        ! Copy averaged grid to all COSMO columnts, leads to computing the same
     1505       ! Copy averaged grid to all COSMO columns, leads to computing the same
    14921506       ! vertical interpolation weights for all columns and to COSMO grid level
    14931507       ! based averaging onto the averaged COSMO heights.
     
    15421556
    15431557       imin = CEILING( (avg_grid % lon(1) - cosmo_lon(0)) / dlon )
    1544        imax = FLOOR (avg_grid % lon(2) - cosmo_lon(0)) / dlon )
     1558       imax = FLOOR  ( (avg_grid % lon(2) - cosmo_lon(0)) / dlon )
    15451559
    15461560       jmin = CEILING( (avg_grid % lat(1) - cosmo_lat(0)) / dlat )
    1547        jmax = FLOOR (avg_grid % lat(2) - cosmo_lat(0)) / dlat )
     1561       jmax = FLOOR  ( (avg_grid % lat(2) - cosmo_lat(0)) / dlat )
    15481562       
    1549        message = "Averaging over "//                                           &
    1550                  TRIM(str(imin)) // " < i < " // TRIM(str(imax)) //            &
     1563       message = "Averaging '" // TRIM(avg_grid % kind) // "' over "//         &
     1564                 TRIM(str(imin)) // " <= i <= " // TRIM(str(imax)) //          &
    15511565                 " and " //                                                    &
    1552                  TRIM(str(jmin)) // " < j < " // TRIM(str(jmax))
     1566                 TRIM(str(jmin)) // " <= j <= " // TRIM(str(jmax))
    15531567       CALL report( 'get_cosmo_averaging_region', message )
    15541568
     
    15571571       avg_grid % n_columns = nx * ny
    15581572
    1559        ALLOCATE( avg_grid % iii(avg_grid % n_columns),                                 &
     1573       ALLOCATE( avg_grid % iii(avg_grid % n_columns),                         &
    15601574                 avg_grid % jjj(avg_grid % n_columns) )
    15611575
     
    39333947    END SUBROUTINE fill_water_cells
    39343948
     3949
     3950!------------------------------------------------------------------------------!
     3951! Description:
     3952! ------------
     3953!> Fills the thickness array of the COSMO soil layers. Since COSMO's (i.e.
     3954!> TERRA_ML's [1]) soil layer boundaries follow the rule
     3955!>
     3956!>    depth(0) = 0.0, and
     3957!>    depth(k) = 0.01 * 3**(k-1), k in [1,2,3,...,7]
     3958!>
     3959!> and full levels are defined as the midpoints between two layer boundaries,
     3960!> all except the first layer thicknesses equal the depth of the midpoint.
     3961!>
     3962!> [1] A Description of the Nonhydrostatic Regional COSMO Model Part II :
     3963!>     Physical Parameterization*, Sect. 11 TERRA_ML.
     3964!>     http://www.cosmo-model.org/content/model/documentation/core/cosmoPhysParamtr.pdf)
     3965!>
     3966!> Input parameters:
     3967!> -----------------
     3968!>
     3969!> depths: array of full soil layer depths (cell midpoints)
     3970!>
     3971!>
     3972!> Output parameters:
     3973!> ------------------
     3974!>
     3975!> d_depth: array of soil layer thicknesses
     3976!>
     3977!------------------------------------------------------------------------------!
     3978    SUBROUTINE get_soil_layer_thickness(depths, d_depth)
     3979       
     3980       REAL(dp), INTENT(IN)  ::  depths(:)
     3981       REAL(dp), INTENT(OUT) ::  d_depth(:)
     3982
     3983       d_depth(:) = depths(:)
     3984       d_depth(1) = 2.0_dp * depths(1)
     3985
     3986    END SUBROUTINE get_soil_layer_thickness
     3987
    39353988 END MODULE grid
Note: See TracChangeset for help on using the changeset viewer.