Changeset 3404 for palm/trunk/SOURCE/nesting_offl_mod.f90

Timestamp:

Oct 23, 2018 1:29:11 PM (4 years ago)

Author:

suehring

Message:

Consider time-dependent geostrophic wind components in offline nesting

File:

• palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 !
+! Consider time-dependent geostrophic wind components in offline nesting
+! 
+! 
 ! Initial Revision: 
 ! - separate offline nesting from large_scale_nudging_mod
@@ -656 +659 @@
 !
 !-- Update geostrophic wind components from dynamic input file. 
-!-- Note, at the moment this is not available. As a workaround, use reference
-!-- profiles derived from the boundary values. 
-!-- This will be revised as soon as geostrophic wind components are provided by
-!-- dynamic input file.
-    ug = u_init
-    vg = v_init
+    DO  k = nzb+1, nzt
+       ug(k) = interpolate_in_time( nest_offl%ug(0,k), nest_offl%ug(1,k),      &
+                                    fac_dt )
+       vg(k) = interpolate_in_time( nest_offl%vg(0,k), nest_offl%vg(1,k),      &
+                                    fac_dt )
+    ENDDO
+    ug(nzt+1) = ug(nzt)
+    vg(nzt+1) = vg(nzt)
     
     CALL  cpu_log( log_point(58), 'offline nesting', 'stop' )
@@ -959 +964 @@
 
 !--    Allocate arrays for geostrophic wind components. Arrays will 
-!--    incorporate 2 time levels in order to interpolate in between. Please
-!--    note, forcing using geostrophic wind components is only required in 
-!--    case of cyclic boundary conditions.
-!        IF ( bc_lr_cyc  .AND.  bc_ns_cyc )  THEN
-!           ALLOCATE( nest_offl%ug(0:1,nzb:nzt+1) )
-!           ALLOCATE( nest_offl%vg(0:1,nzb:nzt+1) )
-!        ENDIF
+!--    incorporate 2 time levels in order to interpolate in between. 
+       ALLOCATE( nest_offl%ug(0:1,1:nzt) )
+       ALLOCATE( nest_offl%vg(0:1,1:nzt) )
 !
 !--    Allocate arrays for reading boundary values. Arrays will incorporate 2 
@@ -1008 +1009 @@
        CALL netcdf_data_input_offline_nesting
 !
-!--    Initialize boundary data 
+!--    Initialize boundary data.
        IF ( bc_dirichlet_l )  THEN
           u(nzb+1:nzt,nys:nyn,0)    = nest_offl%u_left(0,nzb+1:nzt,nys:nyn)
@@ -1046 +1047 @@
        ENDIF
 !
+!--    Initialize geostrophic wind components. Actually this is already done in
+!--    init_3d_model when initializing_action = 'inifor', however, in speical 
+!--    case of user-defined initialization this will be done here again, in 
+!--    order to have a consistent initialization. 
+       ug(nzb+1:nzt) = nest_offl%ug(0,nzb+1:nzt)
+       vg(nzb+1:nzt) = nest_offl%vg(0,nzb+1:nzt)
+!
+!--    Set bottom and top boundary condition for geostrophic wind components
+       ug(nzt+1) = ug(nzt)
+       vg(nzt+1) = vg(nzt)
+       ug(nzb)   = ug(nzb+1)
+       vg(nzb)   = vg(nzb+1)
+!
 !--    Initial calculation of the boundary layer depth from the prescribed 
 !--    boundary data. This is requiered for initialize the synthetic turbulence