Ignore:
Timestamp:
Oct 23, 2018 1:29:11 PM (3 years ago)
Author:
suehring
Message:

Consider time-dependent geostrophic wind components in offline nesting

File:
1 edited

Legend:

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

    r3364 r3404  
    2525! -----------------
    2626!
     27! Consider time-dependent geostrophic wind components in offline nesting
     28!
     29!
    2730! Initial Revision:
    2831! - separate offline nesting from large_scale_nudging_mod
     
    656659!
    657660!-- Update geostrophic wind components from dynamic input file.
    658 !-- Note, at the moment this is not available. As a workaround, use reference
    659 !-- profiles derived from the boundary values.
    660 !-- This will be revised as soon as geostrophic wind components are provided by
    661 !-- dynamic input file.
    662     ug = u_init
    663     vg = v_init
     661    DO  k = nzb+1, nzt
     662       ug(k) = interpolate_in_time( nest_offl%ug(0,k), nest_offl%ug(1,k),      &
     663                                    fac_dt )
     664       vg(k) = interpolate_in_time( nest_offl%vg(0,k), nest_offl%vg(1,k),      &
     665                                    fac_dt )
     666    ENDDO
     667    ug(nzt+1) = ug(nzt)
     668    vg(nzt+1) = vg(nzt)
    664669   
    665670    CALL  cpu_log( log_point(58), 'offline nesting', 'stop' )
     
    959964
    960965!--    Allocate arrays for geostrophic wind components. Arrays will
    961 !--    incorporate 2 time levels in order to interpolate in between. Please
    962 !--    note, forcing using geostrophic wind components is only required in
    963 !--    case of cyclic boundary conditions.
    964 !        IF ( bc_lr_cyc  .AND.  bc_ns_cyc )  THEN
    965 !           ALLOCATE( nest_offl%ug(0:1,nzb:nzt+1) )
    966 !           ALLOCATE( nest_offl%vg(0:1,nzb:nzt+1) )
    967 !        ENDIF
     966!--    incorporate 2 time levels in order to interpolate in between.
     967       ALLOCATE( nest_offl%ug(0:1,1:nzt) )
     968       ALLOCATE( nest_offl%vg(0:1,1:nzt) )
    968969!
    969970!--    Allocate arrays for reading boundary values. Arrays will incorporate 2
     
    10081009       CALL netcdf_data_input_offline_nesting
    10091010!
    1010 !--    Initialize boundary data 
     1011!--    Initialize boundary data.
    10111012       IF ( bc_dirichlet_l )  THEN
    10121013          u(nzb+1:nzt,nys:nyn,0)    = nest_offl%u_left(0,nzb+1:nzt,nys:nyn)
     
    10461047       ENDIF
    10471048!
     1049!--    Initialize geostrophic wind components. Actually this is already done in
     1050!--    init_3d_model when initializing_action = 'inifor', however, in speical
     1051!--    case of user-defined initialization this will be done here again, in
     1052!--    order to have a consistent initialization.
     1053       ug(nzb+1:nzt) = nest_offl%ug(0,nzb+1:nzt)
     1054       vg(nzb+1:nzt) = nest_offl%vg(0,nzb+1:nzt)
     1055!
     1056!--    Set bottom and top boundary condition for geostrophic wind components
     1057       ug(nzt+1) = ug(nzt)
     1058       vg(nzt+1) = vg(nzt)
     1059       ug(nzb)   = ug(nzb+1)
     1060       vg(nzb)   = vg(nzb+1)
     1061!
    10481062!--    Initial calculation of the boundary layer depth from the prescribed
    10491063!--    boundary data. This is requiered for initialize the synthetic turbulence
Note: See TracChangeset for help on using the changeset viewer.