Ignore:
Timestamp:
Jun 1, 2018 5:15:08 PM (3 years ago)
Author:
suehring
Message:

Bugfix, initialization of surface elements also in case of restart runs and cyclic fill

File:
1 edited

Legend:

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

    r3026 r3055  
    2626! -----------------
    2727! $Id$
     28! Initialize Obukhov length, friction velocity and momentum fluxes also
     29! in case of restarts and cyclic fill
     30!
     31! 3026 2018-05-22 10:30:53Z schwenkel
    2832! Changed the name specific humidity to mixing ratio, since we are computing
    2933! mixing ratios.
     
    20532057                ENDIF
    20542058!
    2055 !--          Initialization in case of constant profiles
    2056              ELSEIF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0&
    2057                 .OR.  INDEX(initializing_actions, 'inifor' ) /= 0 )  THEN
     2059!--          Initialization in all other cases
     2060             ELSE
    20582061
    20592062                surf%ol(num_h)   = surf%z_mo(num_h) / zeta_min
     
    20672070             ENDIF
    20682071
    2069              surf%rib(num_h)   = 0.0_wp
     2072             surf%rib(num_h)     = 0.0_wp
    20702073             surf%uvw_abs(num_h) = 0.0_wp
    20712074
     
    33213324
    33223325       found              = .TRUE.
    3323 
    33243326
    33253327       SELECT CASE ( restart_string(1:length) )
     
    39463948
    39473949       END SELECT
    3948 
    39493950!
    39503951!--    Redistribute surface elements on its respective type.
     
    39523953            .NOT. INDEX( restart_string(1:length), '%start_index' ) /= 0 )     &
    39533954       THEN
     3955       
    39543956          ic = nxlc
    39553957          DO  i = nxlf, nxrf
     
    39653967                                    surf_usm_h%start_index(jc,ic) )            &
    39663968                            .AND.  l == 0
    3967 
     3969                           
    39683970                IF ( surf_match_def )  THEN
    39693971                   mm = surf_def_h(l)%start_index(jc,ic)
     
    40174019                surf_match_usm  = surf_usm_v(l)%end_index(jc,ic) >=            &
    40184020                                  surf_usm_v(l)%start_index(jc,ic)
    4019 
     4021                                 
    40204022                IF ( surf_match_def )  THEN
    40214023                   mm = surf_def_v(l)%start_index(jc,ic)
Note: See TracChangeset for help on using the changeset viewer.