Changeset 3404 for palm/trunk/SOURCE/nesting_offl_mod.f90
- Timestamp:
- Oct 23, 2018 1:29:11 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/nesting_offl_mod.f90
r3364 r3404 25 25 ! ----------------- 26 26 ! 27 ! Consider time-dependent geostrophic wind components in offline nesting 28 ! 29 ! 27 30 ! Initial Revision: 28 31 ! - separate offline nesting from large_scale_nudging_mod … … 656 659 ! 657 660 !-- 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) 664 669 665 670 CALL cpu_log( log_point(58), 'offline nesting', 'stop' ) … … 959 964 960 965 !-- 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) ) 968 969 ! 969 970 !-- Allocate arrays for reading boundary values. Arrays will incorporate 2 … … 1008 1009 CALL netcdf_data_input_offline_nesting 1009 1010 ! 1010 !-- Initialize boundary data 1011 !-- Initialize boundary data. 1011 1012 IF ( bc_dirichlet_l ) THEN 1012 1013 u(nzb+1:nzt,nys:nyn,0) = nest_offl%u_left(0,nzb+1:nzt,nys:nyn) … … 1046 1047 ENDIF 1047 1048 ! 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 ! 1048 1062 !-- Initial calculation of the boundary layer depth from the prescribed 1049 1063 !-- boundary data. This is requiered for initialize the synthetic turbulence
Note: See TracChangeset
for help on using the changeset viewer.