Changeset 2798 for palm/trunk


Ignore:
Timestamp:
Feb 9, 2018 5:16:39 PM (4 years ago)
Author:
suehring
Message:

Bugfix initialization of %pt_surface array; Output of surface temperature also for default-type surfaces

Location:
palm/trunk/SOURCE
Files:
5 edited

Legend:

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

    r2797 r2798  
    2525! -----------------
    2626! $Id$
     27! Consider also default-type surfaces for surface temperature output.
     28!
     29! 2797 2018-02-08 13:24:35Z suehring
    2730! Enable output of ground-heat flux also at urban surfaces.
    2831!
     
    33683371          ENDIF
    33693372
    3370              IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'tsurf*'  .OR.  &
    3371                     TRIM( var ) == 'ghf*' )                                    &
     3373             IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'ghf*' )        &
    33723374                 .AND.  .NOT.  land_surface  .AND.  .NOT.  urban_surface )     &         
    33733375             THEN
  • palm/trunk/SOURCE/data_output_2d.f90

    r2797 r2798  
    2525! -----------------
    2626! $Id$
     27! Consider also default-type surfaces for surface temperature output.
     28!
     29! 2797 2018-02-08 13:24:35Z suehring
    2730! Enable output of ground-heat flux also at urban surfaces.
    2831!
     
    10371040             CASE ( 'tsurf*_xy' )        ! 2d-array
    10381041                IF ( av == 0 )  THEN
     1042                   DO  m = 1, surf_def_h(0)%ns
     1043                      i                   = surf_def_h(0)%i(m)           
     1044                      j                   = surf_def_h(0)%j(m)
     1045                      local_pf(i,j,nzb+1) = surf_def_h(0)%pt_surface(m)
     1046                   ENDDO
     1047
    10391048                   DO  m = 1, surf_lsm_h%ns
    10401049                      i                   = surf_lsm_h%i(m)           
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r2797 r2798  
    2525! -----------------
    2626! $Id$
     27! Minor bugfix for initialization of pt_surface
     28!
     29! 2797 2018-02-08 13:24:35Z suehring
    2730! Move output of ghf to general 2D output to output ghf also at urban-type
    2831! surfaces.
     
    40314034!--          layer
    40324035             IF ( surf_lsm_h%lambda_surface_s(m) == 0.0_wp )  THEN
    4033                 t_surface_h%var_1d(:)    = t_soil_h%var_2d(nzb_soil,m)
     4036                t_surface_h%var_1d(m)    = t_soil_h%var_2d(nzb_soil,m)
    40344037                surf_lsm_h%pt_surface(m) = t_soil_h%var_2d(nzb_soil,m) / exn
    40354038             ELSE
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r2797 r2798  
    2525! -----------------
    2626! $Id$
     27! Consider also default-type surfaces for surface temperature output.
     28!
     29! 2797 2018-02-08 13:24:35Z suehring
    2730! Enable output of ground-heat flux also at urban surfaces.
    2831!
     
    913916
    914917          CASE ( 'tsurf*' )
     918             DO  m = 1, surf_def_h(0)%ns
     919                i   = surf_def_h(0)%i(m)           
     920                j   = surf_def_h(0)%j(m)
     921                tsurf_av(j,i) = tsurf_av(j,i) + surf_def_h(0)%pt_surface(m)
     922             ENDDO
     923
    915924             DO  m = 1, surf_lsm_h%ns
    916925                i   = surf_lsm_h%i(m)           
  • palm/trunk/SOURCE/write_3d_binary.f90

    r2797 r2798  
    2525! -----------------
    2626! $Id$
     27! Consider also default-type surfaces for surface temperature output.
     28!
     29! 2797 2018-02-08 13:24:35Z suehring
    2730! Enable restarts for land- and urban-surface variables ghf_av, r_a_av and
    2831! tsurf_av.
     
    390393          WRITE ( 14 )  'r_a_av              ';  WRITE ( 14 )  r_a_av
    391394       ENDIF
    392        IF ( ALLOCATED( tsurf_av ) )  THEN
    393           WRITE ( 14 )  'tsurf_av            ';  WRITE ( 14 )  tsurf_av
    394        ENDIF
     395    ENDIF
     396    IF ( ALLOCATED( tsurf_av ) )  THEN
     397       WRITE ( 14 )  'tsurf_av            ';  WRITE ( 14 )  tsurf_av
    395398    ENDIF
    396399
Note: See TracChangeset for help on using the changeset viewer.