Changeset 4768 for palm/trunk


Ignore:
Timestamp:
Nov 2, 2020 7:11:23 PM (4 years ago)
Author:
suehring
Message:

Enable 3D data output also with 64-bit precision

Location:
palm/trunk/SOURCE
Files:
17 edited

Legend:

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

    r4633 r4768  
    2626! -----------------
    2727! $Id$
     28! Enable 3D data output also with 64-bit precision
     29!
     30! 4633 2020-08-05 14:21:14Z suehring
    2831! Bugfix in check for humidity
    2932!
     
    11231126    LOGICAL, INTENT(OUT) ::  found   !< Output found?
    11241127
    1125     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< Temp. result grid to return
     1128    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< Temp. result grid to return
    11261129!
    11271130!-- Internal variables
  • palm/trunk/SOURCE/bulk_cloud_model_mod.f90

    r4755 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4755 2020-10-21 15:09:15Z schwenkel
    2629! Minor reformatting
    2730!
     
    45744577    REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
    45754578
    4576     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf   !< local
     4579    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf   !< local
    45774580                                                        !< array to which output data is resorted to
    45784581
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r4731 r4768  
    2626! -----------------
    2727! $Id$
     28! Enable 3D data output also with 64-bit precision
     29!
     30! 4731 2020-10-07 13:25:11Z schwenkel
    2831! Move exchange_horiz from time_integration to modules
    2932!
     
    15571560    REAL(wp)             ::  fill_value   !<
    15581561
    1559     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf
     1562    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf
    15601563!
    15611564!-- Local variables
  • palm/trunk/SOURCE/data_output_3d.f90

    r4559 r4768  
    1919! Current revisions:
    2020! ------------------
    21 ! 
    22 ! 
     21!
     22!
    2323! Former revisions:
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4559 2020-06-11 08:51:48Z raasch
    2629! file re-formatted to follow the PALM coding standard
    2730!
     
    170173    REAL(wp)     ::  s_r3      !< sum( particle-radius**3 )
    171174
    172     REAL(sp), DIMENSION(:,:,:), ALLOCATABLE ::  local_pf  !< output array
     175    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  local_pf  !< output array
    173176
    174177    REAL(wp), DIMENSION(:,:,:), POINTER ::  to_be_resorted  !< pointer to array which shall be
     
    261264!
    262265!--    Before each output, set array local_pf to fill value
    263        local_pf = fill_value
     266       local_pf = REAL( fill_value, KIND=wp )
    264267!
    265268!--    Set masking flag for topography for not resorted arrays
     
    674677!                count = (/ nxr-nxl+2, nyn-nys+2, nzt_do-nzb_do+1, 1 /) )
    675678!          ELSE
    676              nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                          &
    677                                      local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                      &
    678                                      start = (/ nxl+1, nys+1, nzb_do+1, do3d_time_count(av) /),    &
    679                                      count = (/ nxr-nxl+1, nyn-nys+1, nzt_do-nzb_do+1, 1 /)        &
    680                                    )
     679!
     680!--       Call for parallel runs
     681          nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                             &
     682                                   local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                        &
     683                                   start = (/ nxl+1, nys+1, nzb_do+1, do3d_time_count(av) /),      &
     684                                   count = (/ nxr-nxl+1, nyn-nys+1, nzt_do-nzb_do+1, 1 /)          &
     685                                 )
    681686!          ENDIF
    682687          CALL netcdf_handle_error( 'data_output_3d', 386 )
     
    685690#else
    686691#if defined( __netcdf )
     692!
     693!--    Call for non-parallel runs
    687694       nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                                &
    688695                               local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                            &
  • palm/trunk/SOURCE/diagnostic_output_quantities_mod.f90

    r4757 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4757 2020-10-26 10:23:38Z schwenkel
    2629! Implement relative humidity as diagnostic output quantity
    2730!
     
    12011204    REAL(wp) ::  fill_value       !< value for the _FillValue attribute
    12021205
    1203     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf        !<
     1206    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf        !<
    12041207    REAL(wp), DIMENSION(:,:,:), POINTER ::                 to_be_resorted  !< points to array which needs to be resorted for output
    12051208
  • palm/trunk/SOURCE/dynamics_mod.f90

    r4760 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4760 2020-10-26 13:26:47Z schwenkel
    2629! Implement relative humidity as diagnostic output quantity
    2730!
     
    18391842    REAL(wp), INTENT(IN) ::  fill_value    !< value for the _FillValue attribute
    18401843
    1841     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
     1844    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
    18421845
    18431846!
  • palm/trunk/SOURCE/gust_mod.f90

    r4646 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4646 2020-08-24 16:02:40Z raasch
    2629! file re-formatted to follow the PALM coding standard
    2730!
     
    543546       REAL(wp), INTENT(IN) ::  fill_value !< value for the _FillValue attribute
    544547
    545        REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf   !< local
     548       REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf   !< local
    546549                                                                                        !< array to which output data is resorted to
    547550
  • palm/trunk/SOURCE/indoor_model_mod.f90

    r4750 r4768  
    2525! -----------------
    2626! $Id$
     27! Enable 3D data output also with 64-bit precision
     28!
     29! 4750 2020-10-16 14:27:48Z suehring
    2730! - Namelist parameter added to switch-off/on the indoor model during wall/soil spinup
    2831! - Bugfix in window-wall treatment during spinup - in the urban-surface model the window fraction
     
    19961999    REAL(wp), INTENT(IN) ::  fill_value !< value for the _FillValue attribute
    19972000
    1998     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
     2001    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
    19992002
    20002003    local_pf = fill_value
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r4723 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4723 2020-10-06 15:16:11Z raasch
    2629! file re-formatted to follow the PALM coding standard
    2730!
     
    62236226    REAL(wp) ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
    62246227
    6225     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_soil:nzt_soil) ::  local_pf !<
     6228    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_soil:nzt_soil) ::  local_pf !<
    62266229
    62276230
  • palm/trunk/SOURCE/module_interface.f90

    r4757 r4768  
    2424! -----------------
    2525! $Id$
     26! Enable 3D data output also with 64-bit precision
     27!
     28! 4757 2020-10-26 10:23:38Z schwenkel
    2629! Implement doq calls
    2730!
     
    16151618    REAL(wp),          INTENT(IN)    ::  fill_value !< ToDo: refactor
    16161619
    1617     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf
     1620    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) ::  local_pf
    16181621
    16191622
  • palm/trunk/SOURCE/ocean_mod.f90

    r4731 r4768  
    2525! -----------------
    2626! $Id$
     27! Enable 3D data output also with 64-bit precision
     28!
     29! 4731 2020-10-07 13:25:11Z schwenkel
    2730! Move exchange_horiz from time_integration to modules
    2831!
     
    11211124    REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
    11221125
    1123     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< local
     1126    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< local
    11241127                                  !< array to which output data is resorted to
    11251128
  • palm/trunk/SOURCE/plant_canopy_model_mod.f90

    r4671 r4768  
    2727! -----------------
    2828! $Id$
     29! Enable 3D data output also with 64-bit precision
     30!
     31! 4671 2020-09-09 20:27:58Z pavelkrc
    2932! Implementation of downward facing USM and LSM surfaces
    3033!
     
    783786
    784787    REAL(wp)     ::  fill_value !< fill value
    785     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !< data output array
     788    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !< data output array
    786789
    787790
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r4741 r4768  
    2828! -----------------
    2929! $Id$
     30! Enable 3D data output also with 64-bit precision
     31!
     32! 4741 2020-10-14 14:32:50Z suehring
    3033! Add option to force calculation of horizontal mean profiles independent on data output
    3134!
     
    1163911642    REAL(wp)     ::  fill_value  !< value for the _FillValue attribute
    1164011643
    11641     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
     11644    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
    1164211645
    1164311646    CHARACTER (len=varnamelength)                   :: var, surfid
  • palm/trunk/SOURCE/salsa_mod.f90

    r4731 r4768  
    2626! -----------------
    2727! $Id$
     28! Enable 3D data output also with 64-bit precision
     29!
     30! 4731 2020-10-07 13:25:11Z schwenkel
    2831! Move exchange_horiz from time_integration to modules
    2932!
     
    1205912062    REAL(wp) ::  temp_bin                 !< temporary array for calculating output variables
    1206012063
    12061     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< local
     12064    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< local
    1206212065
    1206312066    found     = .TRUE.
  • palm/trunk/SOURCE/turbulence_closure_mod.f90

    r4717 r4768  
    2525! -----------------
    2626! $Id$
     27! Enable 3D data output also with 64-bit precision
     28!
     29! 4717 2020-09-30 22:27:40Z pavelkrc
    2730! Fixes and optimizations of OpenMP parallelization, formatting of OpenMP
    2831! directives (J. Resler)
     
    10641067    REAL(wp) ::  fill_value = -9999.0_wp  !< value for the _FillValue attribute
    10651068
    1066     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< local
     1069    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< local
    10671070       !< array to which output data is resorted to
    10681071
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r4750 r4768  
    2727! -----------------
    2828! $Id$
     29! Enable 3D data output also with 64-bit precision
     30!
     31! 4750 2020-10-16 14:27:48Z suehring
    2932! - bugfix in openmp directive
    3033! - make t_green_h and t_green_v public (required in indoor model)
     
    24552458    LOGICAL, INTENT(OUT)           ::  found     !<
    24562459
    2457     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< sp - it has to correspond to module data_output_3d
    2458     REAL(sp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf   !< temp array for urban surface output procedure
     2460    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< sp - it has to correspond to module data_output_3d
     2461    REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf   !< temp array for urban surface output procedure
    24592462
    24602463    found = .TRUE.
     
    25232526                j = surf_usm_h(l)%j(m)
    25242527                k = surf_usm_h(l)%k(m)
    2525                 temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = sp) )
     2528                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = wp) )
    25262529             ENDDO
    25272530          ELSE
     
    25302533                j = surf_usm_v(l)%j(m)
    25312534                k = surf_usm_v(l)%k(m)
    2532                 temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = sp) + 1.0_sp )
     2535                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = wp) + 1.0_sp )
    25332536             ENDDO
    25342537          ENDIF
  • palm/trunk/SOURCE/user_module.f90

    r4535 r4768  
    2020! Current revisions:
    2121! -----------------
    22 ! 
    23 ! 
     22!
     23!
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! Enable 3D data output also with 64-bit precision
     28!
     29! 4535 2020-05-15 12:07:23Z raasch
    2730! bugfix for restart data format query
    2831!
     
    970973!    REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
    971974
    972     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !<
     975    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !<
    973976
    974977!
Note: See TracChangeset for help on using the changeset viewer.