Changeset 4824 for palm/trunk


Ignore:
Timestamp:
Dec 12, 2020 8:45:23 PM (4 years ago)
Author:
moh.hefny
Message:

bugfix in surface_data_output_mod and radiation_model_mod

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4801 r4824  
    2727! -----------------
    2828! $Id$
     29! Bugfix: change the unit of radiation fluxes, few formating changes
     30!
     31! 4801 2020-11-27 16:06:40Z raasch
    2932! bugfix for r4800
    3033!
     
    15391542    ELSE
    15401543       SELECT CASE ( var )
    1541           CASE ( 'rad_lw_cs_hr', 'rad_lw_hr', 'rad_lw_in', 'rad_lw_out', 'rad_sw_cs_hr',           &
    1542                  'rad_sw_hr', 'rad_sw_in', 'rad_sw_out' )
     1544          CASE ( 'rad_lw_cs_hr', 'rad_lw_hr', 'rad_sw_cs_hr', 'rad_sw_hr' )
    15431545             IF ( .NOT.  radiation  .OR.  radiation_scheme /= 'rrtmg' )  THEN
    15441546                message_string = '"output of "' // var // '" requires radiation = .TRUE. and ' //  &
     
    15471549             ENDIF
    15481550             unit = 'K/h'
     1551
     1552          CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_sw_in', 'rad_sw_out' )
     1553             IF ( .NOT.  radiation  .OR.  radiation_scheme /= 'rrtmg' )  THEN
     1554                message_string = '"output of "' // var // '" requires radiation = .TRUE. and ' //  &
     1555                                 'radiation_scheme = "rrtmg"'
     1556                CALL message( 'check_parameters', 'PA0406', 1, 2, 0, 6, 0 )
     1557             ENDIF
     1558             unit = 'W/m2'
    15491559
    15501560          CASE ( 'rad_net*', 'rad_lw_in*', 'rad_lw_out*', 'rad_sw_in*', 'rad_sw_out*' )
     
    62966306!--     Reset for the input from next reflective pass
    62976307        surfins = 0.0_wp
    6298         surfinl = 0.00_wp
     6308        surfinl = 0.0_wp
    62996309!
    63006310!--     Reflected radiation
     
    66376647                                           uvec        !< z, y, x unit vector of incoming flux
    66386648
    6639 
    6640 
    6641 
    6642 
    6643 
    6644 
    6645 
    6646 
    6647 
    6648 
    6649 
     6649   
    66506650    xshift = uvec(3) / uvec(1) * boxsize(1)
    66516651    xmin = MIN( 0._wp, - xshift )
  • palm/trunk/SOURCE/surface_data_output_mod.f90

    r4805 r4824  
    2525! -----------------
    2626! $Id$
     27! Bugfix: correct the down-ward surf_lsm index in the routine surface_data_output_collect_2d
     28!
     29! 4805 2020-12-02 12:05:34Z gronemeier
    2730! Bugfix: add time_dosurf to restart data
    2831!
     
    40404043       ENDDO
    40414044    ELSE
    4042        n_surf = n_surf + surf_lsm_h(0)%ns
     4045       n_surf = n_surf + surf_lsm_h(1)%ns
    40434046    ENDIF
    40444047    IF ( ALLOCATED( var_usm_h0 ) )  THEN
Note: See TracChangeset for help on using the changeset viewer.