Changeset 2797


Ignore:
Timestamp:
Feb 8, 2018 1:24:35 PM (4 years ago)
Author:
suehring
Message:

Output of ground-heat flux at natural- and urban-type surfaces in one output variable; enable restart data of _av variables that belong to both land- and urban-surface model

Location:
palm/trunk/SOURCE
Files:
9 edited

Legend:

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

    r2766 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable output of ground-heat flux also at urban surfaces.
     28!
     29! 2766 2018-01-22 17:17:47Z kanani
    2730! Removed preprocessor directive __chem
    2831!
     
    200203             ENDDO
    201204
     205          CASE ( 'ghf*' )
     206             DO  i = nxl, nxr
     207                DO  j = nys, nyn
     208                   ghf_av(j,i) = ghf_av(j,i)                                   &
     209                                 / REAL( average_count_3d, KIND=wp )
     210                ENDDO
     211             ENDDO
     212             CALL exchange_horiz_2d( ghf_av, nbgp )
     213
    202214          CASE ( 'qsws*' )
    203215             DO  i = nxlg, nxrg
  • palm/trunk/SOURCE/check_parameters.f90

    r2776 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable output of ground-heat flux also at urban surfaces.
     28!
     29! 2776 2018-01-31 10:44:42Z Giersch
    2730! Variable synthetic_turbulence_generator has been abbreviated
    2831!
     
    33223325             CONTINUE
    33233326
    3324           CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'r_a*', 'shf*',       &
    3325                  'ssws*', 't*', 'tsurf*', 'u*', 'z0*', 'z0h*', 'z0q*' )
     3327          CASE ( 'ghf*', 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'r_a*',       &
     3328                 'shf*', 'ssws*', 't*', 'tsurf*', 'u*', 'z0*', 'z0h*', 'z0q*' )
    33263329             IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
    33273330                message_string = 'illegal value for data_output: "' //         &
     
    33583361                CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
    33593362             ENDIF
    3360              IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'tsurf*' )      &
     3363
     3364          IF ( TRIM( var ) == 'ghf*'  .AND.  .NOT.  land_surface )  THEN
     3365             message_string = 'output of "' // TRIM( var ) // '" requi' //     &
     3366                              'res land_surface = .TRUE.'
     3367             CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
     3368          ENDIF
     3369
     3370             IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'tsurf*'  .OR.  &
     3371                    TRIM( var ) == 'ghf*' )                                    &
    33613372                 .AND.  .NOT.  land_surface  .AND.  .NOT.  urban_surface )     &         
    33623373             THEN
     
    33723383             ENDIF
    33733384
     3385             IF ( TRIM( var ) == 'ghf*'   )  unit = 'W/m2'
    33743386             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/m2'
    33753387             IF ( TRIM( var ) == 'ol*'    )  unit = 'm'
  • palm/trunk/SOURCE/data_output_2d.f90

    r2743 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable output of ground-heat flux also at urban surfaces.
     28!
     29! 2743 2018-01-12 16:03:39Z suehring
    2730! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
    2831!
     
    545548                ENDIF
    546549                IF ( mode == 'xy' )  level_z = zu
     550
     551             CASE ( 'ghf*_xy' )        ! 2d-array
     552                IF ( av == 0 )  THEN
     553                   DO  m = 1, surf_lsm_h%ns
     554                      i                   = surf_lsm_h%i(m)           
     555                      j                   = surf_lsm_h%j(m)
     556                      local_pf(i,j,nzb+1) = surf_lsm_h%ghf(m)
     557                   ENDDO
     558                   DO  m = 1, surf_usm_h%ns
     559                      i                   = surf_usm_h%i(m)           
     560                      j                   = surf_usm_h%j(m)
     561                      local_pf(i,j,nzb+1) = surf_usm_h%frac(0,m)     *          &
     562                                            surf_usm_h%wghf_eb(m)        +      &
     563                                            surf_usm_h%frac(1,m)     *          &
     564                                            surf_usm_h%wghf_eb_green(m)  +      &
     565                                            surf_usm_h%frac(2,m)     *          &
     566                                            surf_usm_h%wghf_eb_window(m)
     567                   ENDDO
     568                ELSE
     569                   DO  i = nxl, nxr
     570                      DO  j = nys, nyn
     571                         local_pf(i,j,nzb+1) = ghf_av(j,i)
     572                      ENDDO
     573                   ENDDO
     574                ENDIF
     575
     576                resorted = .TRUE.
     577                two_d = .TRUE.
     578                level_z(nzb+1) = zu(nzb+1)
    547579
    548580             CASE ( 'ol*_xy' )        ! 2d-array
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r2765 r2797  
    2525! -----------------
    2626! $Id$
     27! Move output of ghf to general 2D output to output ghf also at urban-type
     28! surfaces.
     29! Move restart data of ghf_av to read/write_3d_binary, as this is not a
     30! exclusively LSM variable anymore.   
     31!
     32! 2765 2018-01-22 11:34:58Z maronga
    2733! Major bugfix in calculation of f_shf for vertical surfaces
    2834!
     
    598604              c_soil_av,        & !< average of c_soil
    599605              c_veg_av,         & !< average of c_veg
    600               ghf_av,           & !< average of ghf
    601606              lai_av,           & !< average of lai
    602607              qsws_liq_av,      & !< average of qsws_liq
     
    10441049          unit = 'K'   
    10451050             
    1046        CASE ( 'lai*', 'c_liq*', 'c_soil*', 'c_veg*', 'ghf*', 'm_liq*',         &
     1051       CASE ( 'lai*', 'c_liq*', 'c_soil*', 'c_veg*', 'm_liq*',                 &
    10471052              'qsws_liq*', 'qsws_soil*', 'qsws_veg*', 'r_s*' )
    10481053          IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
     
    10721077             CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    10731078          ENDIF
    1074           IF ( TRIM( var ) == 'ghf*'  .AND.  .NOT.  land_surface )  THEN
    1075              message_string = 'output of "' // TRIM( var ) // '" requi' //     &
    1076                               'res land_surface = .TRUE.'
    1077              CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    1078           ENDIF
    10791079          IF ( TRIM( var ) == 'm_liq*'  .AND.  .NOT.  land_surface )  THEN
    10801080             message_string = 'output of "' // TRIM( var ) // '" requi' //     &
     
    11111111          IF ( TRIM( var ) == 'c_soil*')      unit = 'none'
    11121112          IF ( TRIM( var ) == 'c_veg*' )      unit = 'none'
    1113           IF ( TRIM( var ) == 'ghf*')         unit = 'W/m2'
    11141113          IF ( TRIM( var ) == 'm_liq*'     )  unit = 'm'
    11151114          IF ( TRIM( var ) == 'qsws_liq*'  )  unit = 'W/m2'
     
    50425041                c_veg_av = 0.0_wp
    50435042
    5044              CASE ( 'ghf*' )
    5045                 IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
    5046                    ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
    5047                 ENDIF
    5048                 ghf_av = 0.0_wp
    5049 
    50505043             CASE ( 'lai*' )
    50515044                IF ( .NOT. ALLOCATED( lai_av ) )  THEN
     
    51245117                j   = surf_lsm_h%j(m)
    51255118                c_veg_av(j,i) = c_veg_av(j,i) + surf_lsm_h%c_veg(m)
    5126              ENDDO
    5127 
    5128           CASE ( 'ghf*' )
    5129              DO  m = 1, surf_lsm_h%ns
    5130                 i   = surf_lsm_h%i(m)           
    5131                 j   = surf_lsm_h%j(m)
    5132                 ghf_av(j,i) = ghf_av(j,i) + surf_lsm_h%ghf(m)
    51335119             ENDDO
    51345120
     
    52265212                   c_veg_av(j,i) = c_veg_av(j,i)                               &
    52275213                                   / REAL( average_count_3d, KIND=wp )
    5228                 ENDDO
    5229              ENDDO
    5230 
    5231           CASE ( 'ghf*' )
    5232              DO  i = nxl, nxr
    5233                 DO  j = nys, nyn
    5234                    ghf_av(j,i) = ghf_av(j,i)                                   &
    5235                                  / REAL( average_count_3d, KIND=wp )
    52365214                ENDDO
    52375215             ENDDO
     
    54395417                DO  j = nys, nyn
    54405418                   local_pf(i,j,nzb+1) = c_veg_av(j,i)
    5441                 ENDDO
    5442              ENDDO
    5443           ENDIF
    5444 
    5445           two_d = .TRUE.
    5446           grid = 'zu1'
    5447 
    5448        CASE ( 'ghf*_xy' )        ! 2d-array
    5449           IF ( av == 0 )  THEN
    5450              DO  m = 1, surf_lsm_h%ns
    5451                 i                   = surf_lsm_h%i(m)           
    5452                 j                   = surf_lsm_h%j(m)
    5453                 local_pf(i,j,nzb+1) = surf_lsm_h%ghf(m)
    5454              ENDDO
    5455           ELSE
    5456              DO  i = nxl, nxr
    5457                 DO  j = nys, nyn
    5458                    local_pf(i,j,nzb+1) = ghf_av(j,i)
    54595419                ENDDO
    54605420             ENDDO
     
    57505710       IF ( ALLOCATED( c_veg_av ) )  THEN
    57515711          WRITE ( 14 )  'c_veg_av            ';  WRITE ( 14 ) c_veg_av
    5752        ENDIF
    5753        IF ( ALLOCATED( ghf_av ) )  THEN
    5754           WRITE ( 14 )  'ghf_av              ';  WRITE ( 14 )  ghf_av
    57555712       ENDIF
    57565713       IF ( ALLOCATED( lai_av ) )  THEN
     
    59935950                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    59945951
    5995                 CASE ( 'ghf_av' )
    5996                    IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
    5997                       ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
    5998                    ENDIF
    5999                    IF ( k == 1 )  READ ( 13 )  tmp_2d
    6000                    ghf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =           &
    6001                                   tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    6002 
    60035952                CASE ( 'lai_av' )
    60045953                   IF ( .NOT. ALLOCATED( lai_av ) )  THEN
  • palm/trunk/SOURCE/modules.f90

    r2776 r2797  
    2525! -----------------
    2626! $Id$
     27! +ghf_av
     28!
     29! 2776 2018-01-31 10:44:42Z Giersch
    2730! Variable synthetic_turbulence_generator has been abbreviated and _prerun flags
    2831! for skipping module related restart data has beed introduced
     
    896899    USE kinds
    897900
     901    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  ghf_av                 !< avg. ground heat flux
    898902    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  lwp_av                 !< avg. liquid water path
    899903    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  precipitation_rate_av  !< avg. precipitation rate
  • palm/trunk/SOURCE/read_3d_binary.f90

    r2766 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable restarts for land- and urban-surface variables ghf_av, r_a_av and
     28! tsurf_av.
     29!
     30! 2766 2018-01-22 17:17:47Z kanani
    2731! Removed preprocessor directive __chem
    2832!
     
    481485             SELECT CASE ( TRIM( field_chr ) )
    482486
     487                CASE ( 'ghf_av' )
     488                   IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
     489                      ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
     490                   ENDIF
     491                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     492                   ghf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =           &
     493                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     494
    483495                CASE ( 'e' )
    484496                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     
    729741                   qv_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    730742                             tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     743
     744                CASE ( 'r_a_av' )
     745                   IF ( .NOT. ALLOCATED( r_a_av ) )  THEN
     746                      ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
     747                   ENDIF
     748                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     749                   r_a_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =           &
     750                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    731751
    732752                CASE ( 'random_iv' )  ! still unresolved issue
     
    840860                   ts_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
    841861                        tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     862
     863                CASE ( 'tsurf_av' )
     864                   IF ( .NOT. ALLOCATED( tsurf_av ) )  THEN
     865                      ALLOCATE( tsurf_av(nysg:nyng,nxlg:nxrg) )
     866                   ENDIF
     867                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     868                   tsurf_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =         &
     869                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    842870
    843871                CASE ( 'u' )
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r2790 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable output of ground-heat flux also at urban surfaces.
     28!
     29! 2790 2018-02-06 11:57:19Z suehring
    2730! Bugfix in summation of surface sensible and latent heat flux
    2831!
     
    172175
    173176    USE averaging,                                                             &
    174         ONLY:  diss_av, e_av, kh_av, km_av, lpt_av, lwp_av, nc_av, nr_av,      &
     177        ONLY:  diss_av, e_av, ghf_av, kh_av, km_av, lpt_av, lwp_av, nc_av,     &
     178               nr_av,                                                          &
    175179               ol_av, p_av, pc_av, pr_av, prr_av, precipitation_rate_av, pt_av,&
    176180               q_av, qc_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, &
     
    255259          SELECT CASE ( trimvar )
    256260
     261             CASE ( 'ghf*' )
     262                IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
     263                   ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
     264                ENDIF
     265                ghf_av = 0.0_wp
     266
    257267             CASE ( 'e' )
    258268                IF ( .NOT. ALLOCATED( e_av ) )  THEN
     
    538548!--    Store the array chosen on the temporary array.
    539549       SELECT CASE ( trimvar )
     550
     551          CASE ( 'ghf*' )
     552             DO  m = 1, surf_lsm_h%ns
     553                i   = surf_lsm_h%i(m)           
     554                j   = surf_lsm_h%j(m)
     555                ghf_av(j,i) = ghf_av(j,i) + surf_lsm_h%ghf(m)
     556             ENDDO
     557
     558             DO  m = 1, surf_usm_h%ns
     559                i   = surf_usm_h%i(m)           
     560                j   = surf_usm_h%j(m)
     561                ghf_av(j,i) = ghf_av(j,i) + surf_usm_h%frac(0,m)     *          &
     562                                            surf_usm_h%wghf_eb(m)        +      &
     563                                            surf_usm_h%frac(1,m)     *          &
     564                                            surf_usm_h%wghf_eb_green(m)  +      &
     565                                            surf_usm_h%frac(2,m)     *          &
     566                                            surf_usm_h%wghf_eb_window(m)
     567             ENDDO
    540568
    541569          CASE ( 'e' )
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r2766 r2797  
    2626! -----------------
    2727! $Id$
     28! Comment concerning output of ground-heat flux added.
     29!
     30! 2766 2018-01-22 17:17:47Z kanani
    2831! Removed redundant commas, added some blanks
    2932!
     
    13321335                       ENDIF
    13331336                    ENDDO
    1334 
     1337!
     1338!--             Please note, the following output quantities belongs to the
     1339!--             individual tile fractions - ground heat flux at wall-, window-,
     1340!--             and green fraction. Aggregated ground-heat flux is treated
     1341!--             accordingly in average_3d_data, sum_up_3d_data, etc..
    13351342                CASE ( 'usm_wghf' )
    13361343!--                 array of heat flux from ground (wall, roof, land)
  • palm/trunk/SOURCE/write_3d_binary.f90

    r2718 r2797  
    2525! -----------------
    2626! $Id$
     27! Enable restarts for land- and urban-surface variables ghf_av, r_a_av and
     28! tsurf_av.
     29!
     30! 2718 2018-01-02 08:49:38Z maronga
    2731! Corrected "Former revisions" section
    2832!
     
    135139
    136140    USE control_parameters,                                                    &
    137         ONLY:  iran, humidity, passive_scalar, cloud_physics, cloud_droplets,  &
    138                microphysics_morrison, microphysics_seifert, ocean, topography
     141        ONLY:  cloud_physics, cloud_droplets, humidity, iran,                  &
     142               land_surface, microphysics_morrison, microphysics_seifert,      &
     143               ocean, passive_scalar, topography, urban_surface
    139144               
    140145    USE indices,                                                               &
     
    378383       WRITE ( 14 )  'z0q_av              ';  WRITE ( 14 )  z0q_av
    379384    ENDIF
     385    IF ( land_surface  .OR.  urban_surface )  THEN
     386       IF ( ALLOCATED( ghf_av ) )  THEN
     387          WRITE ( 14 )  'ghf_av              ';  WRITE ( 14 )  ghf_av
     388       ENDIF
     389       IF ( ALLOCATED( r_a_av ) )  THEN
     390          WRITE ( 14 )  'r_a_av              ';  WRITE ( 14 )  r_a_av
     391       ENDIF
     392       IF ( ALLOCATED( tsurf_av ) )  THEN
     393          WRITE ( 14 )  'tsurf_av            ';  WRITE ( 14 )  tsurf_av
     394       ENDIF
     395    ENDIF
    380396
    381397!
Note: See TracChangeset for help on using the changeset viewer.