Changeset 3182 for palm/trunk/SOURCE/exchange_horiz_2d.f90

Timestamp:

Jul 27, 2018 1:36:03 PM (6 years ago)

Author:

suehring

Message:

New Inifor features: grid stretching, improved command-interface, support start dates in different formats in both YYYYMMDD and YYYYMMDDHH, Ability to manually control input file prefixes (--radiation-prefix, --soil-preifx, --flow-prefix, --soilmoisture-prefix) for compatiblity with DWD forcast naming scheme; GNU-style short and long option; Prepared output of large-scale forcing profiles (no computation yet); Added preprocessor flag netcdf4 to switch output format between netCDF 3 and 4; Updated netCDF variable names and attributes to comply with PIDS v1.9; Inifor bugfixes: Improved compatibility with older Intel Intel compilers by avoiding implicit array allocation; Added origin_lon/_lat values and correct reference time in dynamic driver global attributes; corresponding PALM changes: adjustments to revised Inifor; variables names in dynamic driver adjusted; enable geostrophic forcing also in offline nested mode; variable names in LES-LES and COSMO offline nesting changed; lateral boundary flags for nesting, in- and outflow conditions renamed

File:

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

palm/trunk/SOURCE/exchange_horiz_2d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in offline nesting mode and flags indicating lateral
+! boundary conditions
 ! 
 ! Former revisions:
@@ -84 +85 @@
 
     USE control_parameters,                                                    &
-        ONLY :  bc_lr_cyc, bc_ns_cyc, force_bound_l, force_bound_n,            &
-                force_bound_r, force_bound_s,                                  &
-                inflow_l, inflow_n, inflow_r, inflow_s,                        &
-                nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s,        &
-                outflow_l, outflow_n, outflow_r, outflow_s
+        ONLY :  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                &
+                bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,          &
+                bc_radiation_n, bc_radiation_r, bc_radiation_s 
                 
     USE cpulog,                                                                &
@@ -184 +183 @@
 !-- Neumann-conditions at inflow/outflow/nested boundaries
     IF ( nesting_mode /= 'vertical' )  THEN
-       IF ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.  force_bound_l )   &
-       THEN
+       IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
           DO  i = nbgp, 1, -1
              ar(:,nxl-i) = ar(:,nxl)
           ENDDO
        ENDIF
-       IF ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.  force_bound_r )   &
-       THEN
+       IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
           DO  i = 1, nbgp
              ar(:,nxr+i) = ar(:,nxr)
           ENDDO
        ENDIF
-       IF ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.  force_bound_s )   &
-       THEN
+       IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
           DO  i = nbgp, 1, -1
              ar(nys-i,:) = ar(nys,:)
           ENDDO
        ENDIF
-       IF ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.  force_bound_n )   &
-       THEN
+       IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
           DO  i = 1, nbgp
              ar(nyn+i,:) = ar(nyn,:)
@@ -227 +222 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, grid_level, force_bound_l, force_bound_n, &
-               force_bound_r, force_bound_s, nest_bound_l, nest_bound_n,       &
-               nest_bound_r, nest_bound_s
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,           &
+               bc_radiation_n, bc_radiation_r, bc_radiation_s, grid_level
         
     USE cpulog,                                                                &
@@ -326 +321 @@
 !
 !-- Neumann-conditions at inflow/outflow/nested boundaries
-    IF ( nest_bound_l  .OR.  force_bound_l )  THEN
+    IF ( bc_dirichlet_l )  THEN
        DO  i = nbgp_local, 1, -1
          ar(:,nxl_l-i) = ar(:,nxl_l)
        ENDDO
     ENDIF
-    IF ( nest_bound_r  .OR.  force_bound_r )  THEN
+    IF ( bc_dirichlet_r )  THEN
        DO  i = 1, nbgp_local
           ar(:,nxr_l+i) = ar(:,nxr_l)
        ENDDO
     ENDIF
-    IF ( nest_bound_s  .OR.  force_bound_s )  THEN
+    IF ( bc_dirichlet_s )  THEN
        DO  i = nbgp_local, 1, -1
          ar(nys_l-i,:) = ar(nys_l,:)
        ENDDO
     ENDIF
-    IF ( nest_bound_n  .OR.  force_bound_n )  THEN
+    IF ( bc_dirichlet_n )  THEN
        DO  i = 1, nbgp_local
          ar(nyn_l+i,:) = ar(nyn_l,:)