Ignore:
Timestamp:
Jul 27, 2018 1:36:03 PM (3 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:
1 edited

Legend:

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

    r2718 r3182  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Remove masking of geostrophic wind forcing in offline nesting case
    2323!
    2424! Former revisions:
     
    109109           
    110110       USE control_parameters,                                                 &
    111            ONLY:  f, forcing, fs, message_string
     111           ONLY:  f, fs, message_string, nesting_offline
    112112           
    113113       USE indices,                                                            &
     
    123123       INTEGER(iwp) ::  k          !< running index z direction
    124124
    125        REAL(wp)     ::  flag       !< flag to mask topography
    126        REAL(wp)     ::  flag_force !< flag to mask large-scale pressure gradient in case larger-scale forcing is applied
    127 
    128        flag_force = MERGE( 0.0_wp, 1.0_wp, forcing )
     125       REAL(wp)     ::  flag           !< flag to mask topography
     126
    129127!
    130128!--    Compute Coriolis terms for the three velocity components
     
    144142                      tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25_wp *          &
    145143                                   ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +    &
    146                                      v(k,j+1,i) ) - vg(k) * flag_force         &
    147                                                           ) * flag           &
     144                                     v(k,j+1,i) ) - vg(k) ) * flag             &
    148145                                                - fs *    ( 0.25_wp *          &
    149146                                   ( w(k-1,j,i-1) + w(k-1,j,i) + w(k,j,i-1) +  &
     
    167164                      tend(k,j,i) = tend(k,j,i) - f *     ( 0.25_wp *          &
    168165                                   ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +    &
    169                                      u(k,j,i+1) ) - ug(k) * flag_force         &
    170                                                           ) * flag
     166                                     u(k,j,i+1) ) - ug(k) ) * flag
    171167                   ENDDO
    172168                ENDDO
     
    212208           
    213209       USE control_parameters,                                                 &
    214            ONLY:  f, forcing, fs, message_string
     210           ONLY:  f, fs, message_string, nesting_offline
    215211           
    216212       USE indices,                                                            &
     
    227223
    228224       REAL(wp)     ::  flag       !< flag to mask topography
    229        REAL(wp)     ::  flag_force !< flag to mask large-scale pressure gradient in case larger-scale forcing is applied
    230 
    231        flag_force = MERGE( 0.0_wp, 1.0_wp, forcing )
     225
    232226!
    233227!--    Compute Coriolis terms for the three velocity components
     
    244238                tend(k,j,i) = tend(k,j,i) + f  *     ( 0.25_wp *               &
    245239                                ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +       &
    246                                   v(k,j+1,i) ) - vg(k) * flag_force            &
     240                                  v(k,j+1,i) ) - vg(k)                         &
    247241                                                     ) * flag                  &
    248242                                          - fs *     ( 0.25_wp *               &
     
    261255                tend(k,j,i) = tend(k,j,i) - f *        ( 0.25_wp *             &
    262256                                ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +       &
    263                                   u(k,j,i+1) ) - ug(k)   * flag_force          &
    264                                                        ) * flag
     257                                  u(k,j,i+1) ) - ug(k) ) * flag
    265258             ENDDO
    266259
Note: See TracChangeset for help on using the changeset viewer.