Ignore:
Timestamp:
Mar 26, 2012 2:18:34 PM (12 years ago)
Author:
suehring
Message:

WS5 is available in combination with topography. Version number changed from 3.8 to 3.8a. Modification in init_pt_anomaly.

File:
1 edited

Legend:

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

    r844 r861  
    44! Current revisions:
    55! -----------------
     6! Set wall_flags_0. The array is needed for degradation in ws-scheme near walls,
     7! inflow and outflow boundaries as well as near the bottom and the top of the
     8! model domain.
     9!
     10! Initialization of nzb_s_inner and nzb_w_inner.
     11! gls has to be at least nbgp to do not exceed the array bounds of nzb_local
     12! while setting wall_flags_0
    613!
    714! Former revisions:
     
    104111                                             corner_sr, wall_l, wall_n, wall_r,&
    105112                                             wall_s, nzb_local, nzb_tmp
     113
     114    LOGICAL :: flag_set = .FALSE.
    106115
    107116    REAL    ::  dx_l, dy_l, dz_stretched
     
    311320!-- the flag arrays needed for the multigrid method
    312321    gls = 2**( maximum_grid_level )
     322    IF ( gls < nbgp )  gls = nbgp
    313323
    314324    ALLOCATE( corner_nl(nys:nyn,nxl:nxr), corner_nr(nys:nyn,nxl:nxr),       &
     
    589599
    590600    END SELECT
    591 
     601!
     602!-- Determine the maximum level of topography. Furthermore it is used for
     603!-- steering the degradation of order of the applied advection scheme.
     604    nzb_max = MAXVAL( nzb_local )
    592605!
    593606!-- Consistency checks and index array initialization are only required for
     
    653666!
    654667!--    Initialize index arrays nzb_s_inner and nzb_w_inner
    655        nzb_s_inner = nzb_local(nys-1:nyn+1,nxl-1:nxr+1)
    656        nzb_w_inner = nzb_local(nys-1:nyn+1,nxl-1:nxr+1)
     668       nzb_s_inner = nzb_local(nysg:nyng,nxlg:nxrg)
     669       nzb_w_inner = nzb_local(nysg:nyng,nxlg:nxrg)
    657670
    658671!
     
    10181031                ENDDO
    10191032             ENDDO
    1020           ENDDO 
     1033          ENDDO
    10211034
    10221035!
     
    10361049
    10371050    ENDIF
     1051!
     1052!-- Allocate flags needed for masking walls.
     1053    ALLOCATE( wall_flags_0(nzb:nzt,nys:nyn,nxl:nxr) )
     1054    wall_flags_0 = 0
     1055
     1056    IF ( scalar_advec == 'ws-scheme' )  THEN
     1057!
     1058!--    Set flags to steer the degradation of the advection scheme in advec_ws
     1059!--    near topography, inflow- and outflow boundaries as well as bottom and
     1060!--    top of model domain. wall_flags_0 remains zero for all non-prognostic
     1061!--    grid points.
     1062       DO  i = nxl, nxr
     1063          DO  j = nys, nyn
     1064             DO  k = nzb_s_inner(j,i)+1, nzt
     1065!
     1066!--             scalar - x-direction
     1067!--             WS1 (0), WS3 (1), WS5 (2)
     1068                IF ( k <= nzb_s_inner(j,i+1) .OR. ( outflow_l .AND. i == nxl ) &
     1069                     .OR. ( outflow_r .AND. i == nxr ) )  THEN
     1070                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 0 )
     1071                ELSEIF ( k <= nzb_s_inner(j,i+2) .OR. k <= nzb_s_inner(j,i-1)  &
     1072                         .OR. ( outflow_r .AND. i == nxr-1 )                   &
     1073                         .OR. ( outflow_l .AND. i == nxlu  ) )  THEN
     1074                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 1 )
     1075                ELSE
     1076                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 2 )
     1077                ENDIF
     1078!
     1079!--             scalar - y-direction
     1080!--             WS1 (3), WS3 (4), WS5 (5)
     1081                IF ( k <= nzb_s_inner(j+1,i) .OR. ( outflow_s .AND. j == nys ) &
     1082                     .OR. ( outflow_n .AND. j == nyn ) )  THEN
     1083                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 3 )
     1084!--             WS3
     1085                ELSEIF ( k <= nzb_s_inner(j+2,i) .OR. k <= nzb_s_inner(j-1,i)  &
     1086                        .OR. ( outflow_s .AND. j == nysv )                     &
     1087                        .OR. ( outflow_n .AND. j == nyn-1  ) )  THEN
     1088                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 4 )
     1089!--             WS5
     1090                ELSE
     1091                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 5 )
     1092                ENDIF
     1093!
     1094!--             scalar - z-direction
     1095!--             WS1 (6), WS3 (7), WS5 (8)
     1096                flag_set = .FALSE.
     1097                IF ( k == nzb_s_inner(j,i) + 1 .OR. k == nzt )  THEN
     1098                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 6 )
     1099                   flag_set = .TRUE.
     1100                ELSEIF ( k == nzb_s_inner(j,i) + 2 .OR. k == nzt - 1 )  THEN
     1101                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 7 )
     1102                   flag_set = .TRUE.
     1103                ELSEIF ( k > nzb_s_inner(j,i) .AND. .NOT. flag_set )  THEN
     1104                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 8 )
     1105                ENDIF
     1106             ENDDO
     1107          ENDDO
     1108       ENDDO
     1109    ENDIF
     1110
     1111    IF ( momentum_advec == 'ws-scheme' )  THEN
     1112!
     1113!--    Set wall_flags_0 to steer the degradation of the advection scheme in advec_ws
     1114!--    near topography, inflow- and outflow boundaries as well as bottom and
     1115!--    top of model domain. wall_flags_0 remains zero for all non-prognostic
     1116!--    grid points.
     1117       DO  i = nxl, nxr
     1118          DO  j = nys, nyn
     1119             DO  k = nzb_u_inner(j,i)+1, nzt
     1120!
     1121!--             u component - x-direction
     1122!--             WS1 (9), WS3 (10), WS5 (11)
     1123                IF ( k <= nzb_u_inner(j,i+1)                                  &
     1124                     .OR. ( outflow_r .AND. i == nxr ) )  THEN
     1125                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 9 )
     1126                ELSEIF ( k <= nzb_u_inner(j,i+2) .OR. k <= nzb_u_inner(j,i-1) &
     1127                         .OR. ( outflow_r .AND. i == nxr-1 )                  &
     1128                         .OR. ( outflow_l .AND. i == nxlu  ) )  THEN
     1129                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 10 )
     1130                ELSE
     1131                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 11 )
     1132                ENDIF
     1133!
     1134!--             u component - y-direction
     1135!--             WS1 (12), WS3 (13), WS5 (14)
     1136                IF ( k <= nzb_u_inner(j+1,i) .OR. ( outflow_s .AND. j == nys ) &
     1137                    .OR. ( outflow_n .AND. j == nyn ) )  THEN
     1138                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 12 )
     1139                ELSEIF ( k <= nzb_u_inner(j+2,i) .OR. k <= nzb_u_inner(j-1,i)  &
     1140                    .OR. ( outflow_s .AND. j == nysv )                         &
     1141                    .OR. ( outflow_n .AND. j == nyn-1  ) )  THEN
     1142                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 13 )
     1143                ELSE
     1144                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 14 )
     1145                ENDIF
     1146!
     1147!--             u component - z-direction
     1148!--             WS1 (15), WS3 (16), WS5 (17)
     1149                flag_set = .FALSE.
     1150                IF ( k == nzb_u_inner(j,i) + 1 .OR. k == nzt )  THEN
     1151                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 15 )
     1152                   flag_set = .TRUE.
     1153                ELSEIF ( k == nzb_u_inner(j,i) + 2 .OR. k == nzt - 1 )  THEN
     1154                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 16 )
     1155                   flag_set = .TRUE.
     1156                ELSEIF ( k > nzb_u_inner(j,i) .AND. .NOT. flag_set )  THEN
     1157                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 17 )
     1158                ENDIF
     1159
     1160             ENDDO
     1161          ENDDO
     1162       ENDDO
     1163
     1164       DO  i = nxl, nxr
     1165          DO  j = nys, nyn
     1166             DO  k = nzb_v_inner(j,i)+1, nzt
     1167!
     1168!--             v component - x-direction
     1169!--             WS1 (18), WS3 (19), WS5 (20)
     1170                IF ( k <= nzb_v_inner(j,i+1) .OR. ( outflow_l .AND. i == nxl ) &
     1171                    .OR. ( outflow_r .AND. i == nxr ) )  THEN
     1172                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 18 )
     1173!--             WS3
     1174                ELSEIF ( k <= nzb_v_inner(j,i+2) .OR. k <= nzb_v_inner(j,i-1)  &
     1175                    .OR. ( outflow_r .AND. i == nxr-1 )                        &
     1176                    .OR. ( outflow_l .AND. i == nxlu  ) )  THEN
     1177                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 19 )
     1178                ELSE
     1179                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 20 )
     1180                ENDIF
     1181!
     1182!--             v component - y-direction
     1183!--             WS1 (21), WS3 (22), WS5 (23)
     1184                IF ( k <= nzb_v_inner(j+1,i)                                   &
     1185                     .OR. ( outflow_n .AND. j == nyn ) )  THEN
     1186                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 21 )
     1187                ELSEIF ( k <= nzb_v_inner(j+2,i) .OR. k <= nzb_v_inner(j-1,i)  &
     1188                         .OR. ( outflow_s .AND. j == nysv )                    &
     1189                         .OR. ( outflow_n .AND. j == nyn-1  ) )  THEN
     1190                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 22 )
     1191                ELSE
     1192                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 23 )
     1193                ENDIF
     1194!
     1195!--             v component - z-direction
     1196!--             WS1 (24), WS3 (25), WS5 (26)
     1197                flag_set = .FALSE.
     1198                IF ( k == nzb_v_inner(j,i) + 1 .OR. k == nzt )  THEN
     1199                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 24 )
     1200                   flag_set = .TRUE.
     1201                ELSEIF ( k == nzb_v_inner(j,i) + 2 .OR. k == nzt - 1 )  THEN
     1202                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 25 )
     1203                   flag_set = .TRUE.
     1204                ELSEIF ( k > nzb_v_inner(j,i) .AND. .NOT. flag_set )  THEN
     1205                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 26 )
     1206                ENDIF
     1207
     1208             ENDDO
     1209          ENDDO
     1210       ENDDO
     1211       DO  i = nxl, nxr
     1212          DO  j = nys, nyn
     1213             DO  k = nzb_w_inner(j,i)+1, nzt
     1214!
     1215!--             w component - x-direction
     1216!--             WS1 (27), WS3 (28), WS5 (29)
     1217                IF ( k <= nzb_w_inner(j,i+1) .OR. ( outflow_l .AND. i == nxl ) &
     1218                    .OR. ( outflow_r .AND. i == nxr ) )  THEN
     1219                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 27 )
     1220                ELSEIF ( k <= nzb_w_inner(j,i+2) .OR. k <= nzb_w_inner(j,i-1)  &
     1221                         .OR. ( outflow_r .AND. i == nxr-1 )                   &
     1222                         .OR. ( outflow_l .AND. i == nxlu  ) )  THEN
     1223                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 28 )
     1224                ELSE
     1225                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i),29 )
     1226                ENDIF
     1227!
     1228!--             w component - y-direction
     1229!--             WS1 (30), WS3 (31), WS5 (32)
     1230                IF ( k <= nzb_w_inner(j+1,i) .OR. ( outflow_s .AND. j == nys ) &
     1231                     .OR. ( outflow_n .AND. j == nyn ) )  THEN
     1232                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 30 )
     1233                ELSEIF ( k <= nzb_w_inner(j+2,i) .OR. k <= nzb_w_inner(j-1,i)  &
     1234                         .OR. ( outflow_s .AND. j == nysv )                    &
     1235                         .OR. ( outflow_n .AND. j == nyn-1  ) )  THEN
     1236                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 31 )
     1237                ELSE
     1238                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 32 )
     1239                ENDIF
     1240!
     1241!--             w component - z-direction
     1242!--             WS1 (33), WS3 (34), WS5 (35)
     1243                flag_set = .FALSE.
     1244                IF ( k == nzb_w_inner(j,i) .OR. k == nzb_w_inner(j,i) + 1      &
     1245                                           .OR. k == nzt )  THEN
     1246!
     1247!--                Please note, at k == nzb_w_inner(j,i) a flag is explictely
     1248!--                set, although this is not a prognostic level. However,
     1249!--                contrary to the advection of u,v and s this is necessary
     1250!--                because flux_t(nzb_w_inner(j,i)) is used for the tendency
     1251!--                at k == nzb_w_inner(j,i)+1.
     1252                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 33 )
     1253                   flag_set = .TRUE.
     1254                ELSEIF ( k == nzb_w_inner(j,i) + 2 .OR. k == nzt - 1 )  THEN
     1255                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 34 )
     1256                   flag_set = .TRUE.
     1257                ELSEIF ( k > nzb_w_inner(j,i) .AND. .NOT. flag_set )  THEN
     1258                   wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 35 )
     1259                ENDIF
     1260
     1261             ENDDO
     1262          ENDDO
     1263       ENDDO
     1264
     1265    ENDIF
    10381266
    10391267!
Note: See TracChangeset for help on using the changeset viewer.