Changeset 1701


Ignore:
Timestamp:
Nov 2, 2015 7:43:04 AM (9 years ago)
Author:
maronga
Message:

minor bugfixes for radiation model. bugfix in subjob

Location:
palm/trunk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/subjob

    r1623 r1701  
    2323# Current revisions:
    2424# ------------------
    25 #
     25# Bugfix: added missing init_cmds for lccrayh/lccrayb
    2626#
    2727# Former revisions:
     
    841841$featuredir
    842842
     843$init_cmds
    843844$module_calls
    844845
  • palm/trunk/SOURCE/check_parameters.f90

    r1692 r1701  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfix: definition of rad_net timeseries was missing
    2222!
    2323! Former revisions:
     
    31163116             ENDIF
    31173117
     3118          CASE ( 'rad_net' )
     3119             IF ( (.NOT. radiation) .OR. radiation_scheme == 'constant' )  THEN
     3120                message_string = 'data_output_pr = ' //                        &
     3121                                 TRIM( data_output_pr(i) ) // ' is not ava' // &
     3122                                 'lable for radiation = .FALSE. or ' //        &
     3123                                 'radiation_scheme = "constant"'
     3124                CALL message( 'check_parameters', 'PA0408', 1, 2, 0, 6, 0 )
     3125             ELSE
     3126                dopr_index(i) = 101
     3127                dopr_unit(i)  = 'W/m2'
     3128                hom(:,2,101,:)  = SPREAD( zw, 2, statistic_regions+1 )
     3129             ENDIF
     3130
    31183131          CASE ( 'rad_lw_in' )
    31193132             IF ( (.NOT. radiation) .OR. radiation_scheme == 'constant' )  THEN
  • palm/trunk/SOURCE/data_output_2d.f90

    r1692 r1701  
    1515!
    1616! Copyright 1997-2015 Leibniz Universitaet Hannover
    17 !--------------------------------------------------------------------------------!
     17 !--------------------------------------------------------------------------------!
    1818!
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfix in output of RRTGM data
    2222!
    2323! Former revisions:
     
    953953                   to_be_resorted => rad_lw_in_av
    954954                ENDIF
     955                IF ( mode == 'xy' )  level_z = zu
    955956
    956957             CASE ( 'rad_lw_out_xy', 'rad_lw_out_xz', 'rad_lw_out_yz' )
     
    960961                   to_be_resorted => rad_lw_out_av
    961962                ENDIF
     963                IF ( mode == 'xy' )  level_z = zu
    962964
    963965             CASE ( 'rad_lw_cs_hr_xy', 'rad_lw_cs_hr_xz', 'rad_lw_cs_hr_yz' )
     
    967969                   to_be_resorted => rad_lw_cs_hr_av
    968970                ENDIF
     971                IF ( mode == 'xy' )  level_z = zw
    969972
    970973             CASE ( 'rad_lw_hr_xy', 'rad_lw_hr_xz', 'rad_lw_hr_yz' )
     
    974977                   to_be_resorted => rad_lw_hr_av
    975978                ENDIF
     979                IF ( mode == 'xy' )  level_z = zw
    976980
    977981             CASE ( 'rad_sw_in_xy', 'rad_sw_in_xz', 'rad_sw_in_yz' )
     
    981985                   to_be_resorted => rad_sw_in_av
    982986                ENDIF
     987                IF ( mode == 'xy' )  level_z = zu
    983988
    984989             CASE ( 'rad_sw_out_xy', 'rad_sw_out_xz', 'rad_sw_out_yz' )
     
    988993                   to_be_resorted => rad_sw_out_av
    989994                ENDIF
     995                IF ( mode == 'xy' )  level_z = zu
    990996
    991997             CASE ( 'rad_sw_cs_hr_xy', 'rad_sw_cs_hr_xz', 'rad_sw_cs_hr_yz' )
     
    9951001                   to_be_resorted => rad_sw_cs_hr_av
    9961002                ENDIF
     1003                IF ( mode == 'xy' )  level_z = zw
    9971004
    9981005             CASE ( 'rad_sw_hr_xy', 'rad_sw_hr_xz', 'rad_sw_hr_yz' )
     
    10021009                   to_be_resorted => rad_sw_hr_av
    10031010                ENDIF
     1011                IF ( mode == 'xy' )  level_z = zw
    10041012
    10051013             CASE ( 'rho_xy', 'rho_xz', 'rho_yz' )
  • palm/trunk/SOURCE/flow_statistics.f90

    r1692 r1701  
    11961196                   sums_l(k,108,tn)  = sums_l(k,108,tn)  + rad_sw_cs_hr(k,j,i) &
    11971197                                       * rmask(j,i,sr)
    1198                    sums_l(k,109,tn)  = sums_l(k,108,tn)  + rad_sw_hr(k,j,i)    &
     1198                   sums_l(k,109,tn)  = sums_l(k,109,tn)  + rad_sw_hr(k,j,i)    &
    11991199                                       * rmask(j,i,sr)
    12001200                ENDDO
     
    14011401
    14021402          IF ( radiation_scheme == 'rrtmg' )  THEN
     1403#if defined ( __rrtmg )
    14031404             hom(:,1,106,sr) = sums(:,106)            ! rad_lw_cs_hr
    14041405             hom(:,1,107,sr) = sums(:,107)            ! rad_lw_hr
     
    14061407             hom(:,1,109,sr) = sums(:,109)            ! rad_sw_hr
    14071408
    1408 #if defined ( __rrtmg )
    14091409             hom(:,1,110,sr) = sums(:,110)            ! rrtm_aldif
    14101410             hom(:,1,111,sr) = sums(:,111)            ! rrtm_aldir
     
    15881588          ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
    15891589          ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
    1590           ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
    1591           ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
     1590          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
     1591          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
    15921592
    15931593#if defined ( __rrtmg )
     
    35293529          ts_value(dots_rad+1,sr) = hom(nzb,1,102,sr)          ! rad_lw_in
    35303530          ts_value(dots_rad+2,sr) = hom(nzb,1,103,sr)          ! rad_lw_out
    3531           ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_lw_in
    3532           ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_lw_out
     3531          ts_value(dots_rad+3,sr) = hom(nzb,1,104,sr)          ! rad_sw_in
     3532          ts_value(dots_rad+4,sr) = hom(nzb,1,105,sr)          ! rad_sw_out
    35333533
    35343534#if defined ( __rrtmg )
  • palm/trunk/SOURCE/radiation_model.f90

    r1692 r1701  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfixes: wrong index for output of timeseries, setting of nz_snd_end
    2222!
    2323! Former revisions:
     
    628628       dots_num = dots_num + 5
    629629
    630        dots_label(dots_rad+1) = "rad_net"
    631        dots_label(dots_rad+2) = "rad_lw_in"
    632        dots_label(dots_rad+3) = "rad_lw_out"
    633        dots_label(dots_rad+4) = "rad_sw_in"
    634        dots_label(dots_rad+5) = "rad_sw_out"
     630       dots_label(dots_rad) = "rad_net"
     631       dots_label(dots_rad+1) = "rad_lw_in"
     632       dots_label(dots_rad+2) = "rad_lw_out"
     633       dots_label(dots_rad+3) = "rad_sw_in"
     634       dots_label(dots_rad+4) = "rad_sw_out"
    635635       dots_unit(dots_rad:dots_rad+4) = "W/m2"
    636636
     
    11611161!
    11621162! !--    Allocate temporary array for storing pressure data
    1163        ALLOCATE( hyp_snd_tmp(nzb+1:nz_snd) )
     1163       ALLOCATE( hyp_snd_tmp(1:nz_snd) )
    11641164       hyp_snd_tmp = 0.0_wp
    11651165
     
    11731173!
    11741174!--    Allocate temporary array for storing temperature data
    1175        ALLOCATE( t_snd_tmp(nzb+1:nz_snd) )
     1175       ALLOCATE( t_snd_tmp(1:nz_snd) )
    11761176       t_snd_tmp = 0.0_wp
    11771177
     
    11861186!--    Calculate start of sounding data
    11871187       nz_snd_start = nz_snd + 1
    1188        nz_snd_end   = nz_snd_end
     1188       nz_snd_end   = nz_snd + 1
    11891189
    11901190!
     
    11991199
    12001200       IF ( nz_snd_start <= nz_snd )  THEN
    1201           nz_snd_end = nz_snd - 1
     1201          nz_snd_end = nz_snd
    12021202       END IF
    12031203
     
    12051205!
    12061206!--    Calculate of total grid points for RRTMG calculations
    1207        nzt_rad = nzt + nz_snd_end - nz_snd_start + 2
     1207       nzt_rad = nzt + nz_snd_end - nz_snd_start + 1
    12081208
    12091209!
Note: See TracChangeset for help on using the changeset viewer.