Ignore:
Timestamp:
Jul 30, 2019 2:47:10 PM (23 months ago)
Author:
suehring
Message:

Merge with branch resler: biomet- output of bio_mrt added; plant_canopy - separate vertical dimension for 3D output (to save disk space); radiation - remove unused plant canopy variables; urban-surface model - do not add anthropogenic heat during wall spin-up

File:
1 edited

Legend:

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

    r4069 r4127  
    2525! -----------------
    2626! $Id$
     27! -Introduce new vertical dimension for plant-canopy output.
     28! -Temporarlily disable masked output for soil (merge from branch resler)
     29!
     30! 4069 2019-07-01 14:05:51Z Giersch
    2731! Masked output running index mid has been introduced as a local variable to
    2832! avoid runtime error (Loop variable has been modified) in time_integration
     
    538542                    id_dim_y_xz, id_dim_yv_xz, id_dim_y_yz, id_dim_yv_yz, &
    539543                    id_dim_y_3d, id_dim_yv_3d, id_dim_zs_xy, id_dim_zs_xz, &
    540                     id_dim_zs_yz, id_dim_zs_3d, id_dim_zu_xy, id_dim_zu1_xy, &
     544                    id_dim_zs_yz, id_dim_zs_3d, id_dim_zpc_3d, &
     545                    id_dim_zu_xy, id_dim_zu1_xy, &
    541546                    id_dim_zu_xz, id_dim_zu_yz, id_dim_zu_3d, id_dim_zw_xy, &
    542547                    id_dim_zw_xz, id_dim_zw_yz, id_dim_zw_3d, id_set_xy, &
     
    548553                    id_var_yv_xy, id_var_y_xz, id_var_yv_xz, id_var_y_yz, &
    549554                    id_var_yv_yz, id_var_y_3d, id_var_yv_3d, id_var_zs_xy, &
    550                     id_var_zs_xz, id_var_zs_yz, id_var_zs_3d, id_var_zusi_xy, &
    551                     id_var_zusi_3d, id_var_zu_xy, id_var_zu1_xy, id_var_zu_xz, &
     555                    id_var_zs_xz, id_var_zs_yz, id_var_zs_3d, id_var_zpc_3d, &
     556                    id_var_zusi_xy, id_var_zusi_3d, id_var_zu_xy, id_var_zu1_xy, id_var_zu_xz, &
    552557                    id_var_zu_yz, id_var_zu_3d, id_var_zwwi_xy, id_var_zwwi_3d, &
    553558                    id_var_zw_xy, id_var_zw_xz, id_var_zw_yz, id_var_zw_3d
     
    738743
    739744    USE plant_canopy_model_mod,                                                &
    740         ONLY:  pcm_define_netcdf_grid
     745        ONLY:  pch_index, pcm_define_netcdf_grid
    741746
    742747    USE profil_parameter,                                                      &
     
    15521557
    15531558          ENDIF
    1554 
    1555           IF ( land_surface )  THEN
     1559!
     1560!--       soil is not in masked output for now - disable temporary this block
     1561!          IF ( land_surface )  THEN
    15561562!
    15571563!--          Write zs data (vertical axes for soil model), use negative values
    15581564!--          to indicate soil depth
    1559              ALLOCATE( netcdf_data(mask_size(mid,3)) )
    1560 
    1561              netcdf_data = zs( mask_k_global(mid,:mask_size(mid,3)) )
    1562 
    1563              nc_stat = NF90_PUT_VAR( id_set_mask(mid,av), id_var_zs_mask(mid,av), &
    1564                                      netcdf_data, start = (/ 1 /), &
    1565                                      count = (/ mask_size(mid,3) /) )
    1566              CALL netcdf_handle_error( 'netcdf_define_header', 538 )
    1567 
    1568              DEALLOCATE( netcdf_data )
    1569 
    1570           ENDIF
     1565!             ALLOCATE( netcdf_data(mask_size(mid,3)) )
     1566!
     1567!             netcdf_data = zs( mask_k_global(mid,:mask_size(mid,3)) )
     1568!
     1569!             nc_stat = NF90_PUT_VAR( id_set_mask(mid,av), id_var_zs_mask(mid,av), &
     1570!                                     netcdf_data, start = (/ 1 /), &
     1571!                                     count = (/ mask_size(mid,3) /) )
     1572!             CALL netcdf_handle_error( 'netcdf_define_header', 538 )
     1573!
     1574!             DEALLOCATE( netcdf_data )
     1575!
     1576!          ENDIF
    15711577
    15721578!
     
    18471853          ENDIF
    18481854
     1855          IF ( plant_canopy )  THEN
     1856!
     1857!--          Define vertical coordinate grid (zpc grid)
     1858             CALL netcdf_create_dim( id_set_3d(av), 'zpc_3d',                  &
     1859                                     pch_index+1, id_dim_zpc_3d(av), 70 )
     1860             !netcdf_create_dim(ncid, dim_name, ncdim_type, ncdim_id, error_no)
     1861             CALL netcdf_create_var( id_set_3d(av), (/ id_dim_zpc_3d(av) /),    &
     1862                                     'zpc_3d', NF90_DOUBLE, id_var_zpc_3d(av),   &
     1863                                     'meters', '', 71, 72, 00 )
     1864
     1865          ENDIF
     1866
    18491867!
    18501868!--       Define the variables
     
    20062024             ELSEIF ( grid_z == 'zs' )  THEN
    20072025                id_z = id_dim_zs_3d(av)
     2026             ELSEIF ( grid_z == 'zpc' )  THEN
     2027                id_z = id_dim_zpc_3d(av)
    20082028             ENDIF
    20092029
     
    22482268                                        - zs(nzb_soil:nzt_soil), start = (/ 1 /), &
    22492269                                        count = (/ nzt_soil-nzb_soil+1 /) )
     2270                CALL netcdf_handle_error( 'netcdf_define_header', 86 )
     2271             ENDIF
     2272
     2273             IF ( plant_canopy )  THEN
     2274!
     2275!--             Write zpc grid
     2276                nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zpc_3d(av),  &
     2277                                        zu(nzb:nzb+pch_index), start = (/ 1 /), &
     2278                                        count = (/ pch_index+1 /) )
    22502279                CALL netcdf_handle_error( 'netcdf_define_header', 86 )
    22512280             ENDIF
Note: See TracChangeset for help on using the changeset viewer.