PageRenderTime 57ms CodeModel.GetById 20ms RepoModel.GetById 0ms app.codeStats 1ms

/wrfv2_fire/dyn_em/module_bc_em.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 1199 lines | 853 code | 197 blank | 149 comment | 11 complexity | e54e2d5b2cf9585faac56925e27bb28a MD5 | raw file
Possible License(s): AGPL-1.0

Large files files are truncated, but you can click here to view the full file

    

  1. !WRF:MODEL_LAYER:BOUNDARY
    2. 

  2. !
    3. 

  3. MODULE module_bc_em
    4. 

  4. 

  5.    USE module_bc
    6. 

  6.    USE module_configure
    7. 

  7.    USE module_wrf_error
    8. 

  8. 

  9. CONTAINS
    10. 

  10. 

  11. !------------------------------------------------------------------------
    12. 

  12. 

  13.    SUBROUTINE spec_bdyupdate_ph( ph_save, field,      &
    14. 

  14.                                field_tend, mu_tend, muts, dt,     &
    15. 

  15.                                variable_in, config_flags, & 
    16. 

  16.                                spec_zone,                  &
    17. 

  17.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    18. 

  18.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    19. 

  19.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    20. 

  20.                                its,ite, jts,jte, kts,kte )
    21. 

  21. 

  22. !  This subroutine adds the tendencies in the boundary specified region.
    23. 

  23. !  spec_zone is the width of the outer specified b.c.s that are set here.
    24. 

  24. !  (JD August 2000)
    25. 

  25. 

  26.       IMPLICIT NONE
    27. 

  27. 

  28.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    29. 

  29.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    30. 

  30.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    31. 

  31.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    32. 

  32.       INTEGER,      INTENT(IN   )    :: spec_zone
    33. 

  33.       CHARACTER,    INTENT(IN   )    :: variable_in
    34. 

  34.       REAL,         INTENT(IN   )    :: dt
    35. 

  35. 

  36. 

  37.       REAL,  DIMENSION( ims:ime , kms:kme , jms:jme ), INTENT(INOUT) :: field
    38. 

  38.       REAL,  DIMENSION( ims:ime , kms:kme , jms:jme ), INTENT(IN   ) :: field_tend, ph_save
    39. 

  39.       REAL,  DIMENSION( ims:ime , jms:jme ), INTENT(IN   ) :: mu_tend, muts
    40. 

  40.       TYPE( grid_config_rec_type ) config_flags
    41. 

  41. 

  42.       CHARACTER  :: variable
    43. 

  43.       INTEGER    :: i, j, k, ibs, ibe, jbs, jbe, itf, jtf, ktf
    44. 

  44.       INTEGER    :: b_dist, b_limit
    45. 

  45. 

  46. !     Local array
    47. 

  47. 

  48.       REAL,  DIMENSION( its:ite , jts:jte ) :: mu_old
    49. 

  49.       LOGICAL    :: periodic_x
    50. 

  50. 

  51.       periodic_x = config_flags%periodic_x
    52. 

  52. 

  53.       variable = variable_in
    54. 

  54. 

  55.       IF (variable == 'U') variable = 'u'
    56. 

  56.       IF (variable == 'V') variable = 'v'
    57. 

  57.       IF (variable == 'M') variable = 'm'
    58. 

  58.       IF (variable == 'H') variable = 'h'
    59. 

  59. 

  60.       ibs = ids
    61. 

  61.       ibe = ide-1
    62. 

  62.       itf = min(ite,ide-1)
    63. 

  63.       jbs = jds
    64. 

  64.       jbe = jde-1
    65. 

  65.       jtf = min(jte,jde-1)
    66. 

  66.       ktf = kde-1
    67. 

  67.       IF (variable == 'u') ibe = ide
    68. 

  68.       IF (variable == 'u') itf = min(ite,ide)
    69. 

  69.       IF (variable == 'v') jbe = jde
    70. 

  70.       IF (variable == 'v') jtf = min(jte,jde)
    71. 

  71.       IF (variable == 'm') ktf = kte
    72. 

  72.       IF (variable == 'h') ktf = kte
    73. 

  73. 

  74.       IF (jts - jbs .lt. spec_zone) THEN
    75. 

  75. ! Y-start boundary
    76. 

  76.         DO j = jts, min(jtf,jbs+spec_zone-1)
    77. 

  77.           b_dist = j - jbs
    78. 

  78.           b_limit = b_dist
    79. 

  79.           IF(periodic_x)b_limit = 0
    80. 

  80.           DO k = kts, ktf
    81. 

  81.             DO i = max(its,b_limit+ibs), min(itf,ibe-b_limit)
    82. 

  82. 

  83.               mu_old(i,j) = muts(i,j) - dt*mu_tend(i,j)
    84. 

  84. 

  85.               field(i,k,j) = field(i,k,j)*mu_old(i,j)/muts(i,j) + &
    86. 

  86.                    dt*field_tend(i,k,j)/muts(i,j) +               &
    87. 

  87.                    ph_save(i,k,j)*(mu_old(i,j)/muts(i,j) - 1.)
    88. 

  88. 

  89.             ENDDO
    90. 

  90.           ENDDO
    91. 

  91.         ENDDO
    92. 

  92.       ENDIF 
    93. 

  93.       IF (jbe - jtf .lt. spec_zone) THEN 
    94. 

  94. ! Y-end boundary 
    95. 

  95.         DO j = max(jts,jbe-spec_zone+1), jtf 
    96. 

  96.           b_dist = jbe - j 
    97. 

  97.           b_limit = b_dist
    98. 

  98.           IF(periodic_x)b_limit = 0
    99. 

  99.           DO k = kts, ktf 
    100. 

  100.             DO i = max(its,b_limit+ibs), min(itf,ibe-b_limit)
    101. 

  101. 

  102.               mu_old(i,j) = muts(i,j) - dt*mu_tend(i,j)
    103. 

  103. 

  104.               field(i,k,j) = field(i,k,j)*mu_old(i,j)/muts(i,j) + &
    105. 

  105.                    dt*field_tend(i,k,j)/muts(i,j) +               &
    106. 

  106.                    ph_save(i,k,j)*(mu_old(i,j)/muts(i,j) - 1.)
    107. 

  107. 

  108.             ENDDO
    109. 

  109.           ENDDO
    110. 

  110.         ENDDO
    111. 

  111.       ENDIF 
    112. 

  112. 

  113.     IF(.NOT.periodic_x)THEN
    114. 

  114.       IF (its - ibs .lt. spec_zone) THEN
    115. 

  115. ! X-start boundary
    116. 

  116.         DO i = its, min(itf,ibs+spec_zone-1)
    117. 

  117.           b_dist = i - ibs
    118. 

  118.           DO k = kts, ktf
    119. 

  119.             DO j = max(jts,b_dist+jbs+1), min(jtf,jbe-b_dist-1)
    120. 

  120. 

  121.               mu_old(i,j) = muts(i,j) - dt*mu_tend(i,j)
    122. 

  122. 

  123.               field(i,k,j) = field(i,k,j)*mu_old(i,j)/muts(i,j) + &
    124. 

  124.                    dt*field_tend(i,k,j)/muts(i,j) +               &
    125. 

  125.                    ph_save(i,k,j)*(mu_old(i,j)/muts(i,j) - 1.)
    126. 

  126. 

  127.             ENDDO
    128. 

  128.           ENDDO
    129. 

  129.         ENDDO
    130. 

  130.       ENDIF 
    131. 

  131. 

  132.       IF (ibe - itf .lt. spec_zone) THEN
    133. 

  133. ! X-end boundary
    134. 

  134.         DO i = max(its,ibe-spec_zone+1), itf
    135. 

  135.           b_dist = ibe - i
    136. 

  136.           DO k = kts, ktf
    137. 

  137.             DO j = max(jts,b_dist+jbs+1), min(jtf,jbe-b_dist-1)
    138. 

  138. 

  139.               mu_old(i,j) = muts(i,j) - dt*mu_tend(i,j)
    140. 

  140. 

  141.               field(i,k,j) = field(i,k,j)*mu_old(i,j)/muts(i,j) + &
    142. 

  142.                    dt*field_tend(i,k,j)/muts(i,j) +               &
    143. 

  143.                    ph_save(i,k,j)*(mu_old(i,j)/muts(i,j) - 1.)
    144. 

  144. 

  145.             ENDDO
    146. 

  146.           ENDDO
    147. 

  147.         ENDDO
    148. 

  148.       ENDIF 
    149. 

  149.     ENDIF
    150. 

  150. 

  151.    END SUBROUTINE spec_bdyupdate_ph
    152. 

  152. 

  153. !------------------------------------------------------------------------
    154. 

  154. 

  155.    SUBROUTINE relax_bdy_dry ( config_flags,                                    &
    156. 

  156.                               ru_tendf, rv_tendf, ph_tendf, t_tendf,           &
    157. 

  157.                               rw_tendf, mu_tend,                               &
    158. 

  158.                               ru, rv, ph, t,                                   &
    159. 

  159.                               w, mu, mut,                                      &
    160. 

  160.                               u_bxs,u_bxe,u_bys,u_bye,                         &
    161. 

  161.                               v_bxs,v_bxe,v_bys,v_bye,                         &
    162. 

  162.                               ph_bxs,ph_bxe,ph_bys,ph_bye,                     &
    163. 

  163.                               t_bxs,t_bxe,t_bys,t_bye,                         &
    164. 

  164.                               w_bxs,w_bxe,w_bys,w_bye,                         &
    165. 

  165.                               mu_bxs,mu_bxe,mu_bys,mu_bye,                     &
    166. 

  166.                               u_btxs,u_btxe,u_btys,u_btye,                     &
    167. 

  167.                               v_btxs,v_btxe,v_btys,v_btye,                     &
    168. 

  168.                               ph_btxs,ph_btxe,ph_btys,ph_btye,                 &
    169. 

  169.                               t_btxs,t_btxe,t_btys,t_btye,                     &
    170. 

  170.                               w_btxs,w_btxe,w_btys,w_btye,                     &
    171. 

  171.                               mu_btxs,mu_btxe,mu_btys,mu_btye,                 &
    172. 

  172.                               spec_bdy_width, spec_zone, relax_zone,           &
    173. 

  173.                               dtbc, fcx, gcx,             &
    174. 

  174.                               ids,ide, jds,jde, kds,kde,  & ! domain dims
    175. 

  175.                               ims,ime, jms,jme, kms,kme,  & ! memory dims
    176. 

  176.                               ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    177. 

  177.                               its, ite, jts, jte, kts, kte)
    178. 

  178.    IMPLICIT NONE
    179. 

  179. 

  180.    !  Input data.
    181. 

  181.    TYPE( grid_config_rec_type ) config_flags
    182. 

  182. 

  183.    INTEGER ,               INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
    184. 

  184.                                             ims, ime, jms, jme, kms, kme, &
    185. 

  185.                                             ips, ipe, jps, jpe, kps, kpe, & 
    186. 

  186.                                             its, ite, jts, jte, kts, kte
    187. 

  187.    INTEGER ,               INTENT(IN   ) :: spec_bdy_width, spec_zone, relax_zone
    188. 

  188. 

  189.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(IN   ) :: ru,     &
    190. 

  190.                                                                       rv,     &
    191. 

  191.                                                                       ph,     &
    192. 

  192.                                                                       w,      &
    193. 

  193.                                                                       t
    194. 

  194.    REAL , DIMENSION( ims:ime , jms:jme  ) , INTENT(IN   )          :: mu  , &
    195. 

  195.                                                                       mut
    196. 

  196.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(INOUT) :: ru_tendf, &
    197. 

  197.                                                                       rv_tendf, &
    198. 

  198.                                                                       ph_tendf, &
    199. 

  199.                                                                       rw_tendf, &
    200. 

  200.                                                                       t_tendf
    201. 

  201.    REAL , DIMENSION( ims:ime , jms:jme  ) , INTENT(INOUT)          :: mu_tend
    202. 

  202.    REAL , DIMENSION( spec_bdy_width) , INTENT(IN   ) :: fcx, gcx
    203. 

  203. 

  204.    REAL,  DIMENSION( jms:jme , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: u_bxs,u_bxe, &
    205. 

  205.                                                                                v_bxs,v_bxe, &
    206. 

  206.                                                                                ph_bxs,ph_bxe, &
    207. 

  207.                                                                                w_bxs,w_bxe, &
    208. 

  208.                                                                                t_bxs,t_bxe, &
    209. 

  209.                                                                                u_btxs,u_btxe, &
    210. 

  210.                                                                                v_btxs,v_btxe, &
    211. 

  211.                                                                                ph_btxs,ph_btxe, &
    212. 

  212.                                                                                w_btxs,w_btxe, &
    213. 

  213.                                                                                t_btxs,t_btxe
    214. 

  214. 

  215.    REAL,  DIMENSION( ims:ime , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: u_bys,u_bye, &
    216. 

  216.                                                                                v_bys,v_bye, &
    217. 

  217.                                                                                ph_bys,ph_bye, &
    218. 

  218.                                                                                w_bys,w_bye, &
    219. 

  219.                                                                                t_bys,t_bye, &
    220. 

  220.                                                                                u_btys,u_btye, &
    221. 

  221.                                                                                v_btys,v_btye, &
    222. 

  222.                                                                                ph_btys,ph_btye, &
    223. 

  223.                                                                                w_btys,w_btye, &
    224. 

  224.                                                                                t_btys,t_btye
    225. 

  225. 

  226. 

  227.    REAL,  DIMENSION( jms:jme , 1:1     , spec_bdy_width    ), INTENT(IN   ) :: mu_bxs,mu_bxe, &
    228. 

  228.                                                                                mu_btxs,mu_btxe
    229. 

  229. 

  230.    REAL,  DIMENSION( ims:ime , 1:1     , spec_bdy_width    ), INTENT(IN   ) :: mu_bys,mu_bye, &
    231. 

  231.                                                                                mu_btys,mu_btye
    232. 

  232.    REAL, INTENT(IN   ) :: dtbc
    233. 

  233. 

  234. ! changed to tile dimensions, 20090923, JM
    235. 

  235.    REAL , DIMENSION( its-1:ite+1 , kts:kte, jts-1:jte+1  ) :: rfield
    236. 

  236.    INTEGER :: i_start, i_end, j_start, j_end, i, j, k
    237. 

  237. 

  238.            CALL relax_bdytend ( ru, ru_tendf,             &
    239. 

  239.                                u_bxs,u_bxe,u_bys,u_bye,u_btxs,u_btxe,u_btys,u_btye, &
    240. 

  240.                                'u'        , config_flags, &
    241. 

  241.                                spec_bdy_width, spec_zone, relax_zone, &
    242. 

  242.                                dtbc, fcx, gcx,             &
    243. 

  243.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    244. 

  244.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    245. 

  245.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    246. 

  246.                                its,ite, jts,jte, kts,kte )
    247. 

  247.            CALL relax_bdytend ( rv, rv_tendf,             &
    248. 

  248.                                v_bxs,v_bxe,v_bys,v_bye,v_btxs,v_btxe,v_btys,v_btye, &
    249. 

  249.                                'v'        , config_flags, &
    250. 

  250.                                spec_bdy_width, spec_zone, relax_zone, &
    251. 

  251.                                dtbc, fcx, gcx,             &
    252. 

  252.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    253. 

  253.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    254. 

  254.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    255. 

  255.                                its,ite, jts,jte, kts,kte )
    256. 

  256. 

  257. ! rfield will be calculated beyond tile limits because relax_bdytend
    258. 

  258. !   requires a 5-point stencil, and this avoids need for inter-tile/patch 
    259. 

  259. !   communication here
    260. 

  260.            i_start = max(its-1, ids)
    261. 

  261.            i_end = min(ite+1, ide-1)
    262. 

  262.            j_start = max(jts-1, jds)
    263. 

  263.            j_end = min(jte+1, jde-1)
    264. 

  264. 

  265.            DO j=j_start,j_end
    266. 

  266.            DO k=kts,kte
    267. 

  267.            DO i=i_start,i_end
    268. 

  268.               rfield(i,k,j) = ph(i,k,j)*mut(i,j)
    269. 

  269.            ENDDO
    270. 

  270.            ENDDO
    271. 

  271.            ENDDO
    272. 

  272.            CALL relax_bdytend_tile ( rfield, ph_tendf,             &
    273. 

  273.                                ph_bxs,ph_bxe,ph_bys,ph_bye, ph_btxs,ph_btxe,ph_btys,ph_btye,       &
    274. 

  274.                                'h'        , config_flags, &
    275. 

  275.                                spec_bdy_width, spec_zone, relax_zone, &
    276. 

  276.                                dtbc, fcx, gcx,             &
    277. 

  277.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    278. 

  278.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    279. 

  279.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    280. 

  280.                                its,ite, jts,jte, kts,kte,  &
    281. 

  281.                                its-1, ite+1, jts-1,jte+1,kts,kte )  ! dims of first argument
    282. 

  282.            DO j=j_start,j_end
    283. 

  283.            DO k=kts,kte-1
    284. 

  284.            DO i=i_start,i_end
    285. 

  285.               rfield(i,k,j) = t(i,k,j)*mut(i,j)
    286. 

  286.            ENDDO
    287. 

  287.            ENDDO
    288. 

  288.            ENDDO
    289. 

  289.            CALL relax_bdytend_tile ( rfield, t_tendf,              &
    290. 

  290.                                t_bxs,t_bxe,t_bys,t_bye, t_btxs,t_btxe,t_btys,t_btye,       &
    291. 

  291.                                't'        , config_flags, &
    292. 

  292.                                spec_bdy_width, spec_zone, relax_zone, &
    293. 

  293.                                dtbc, fcx, gcx,             &
    294. 

  294.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    295. 

  295.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    296. 

  296.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    297. 

  297.                                its,ite, jts,jte, kts,kte,  &
    298. 

  298.                                its-1, ite+1, jts-1,jte+1,kts,kte )  ! dims of first argument
    299. 

  299.            CALL relax_bdytend ( mu, mu_tend,               &
    300. 

  300.                                mu_bxs,mu_bxe,mu_bys,mu_bye, mu_btxs,mu_btxe,mu_btys,mu_btye,                &
    301. 

  301.                                'm'        , config_flags,  &
    302. 

  302.                                spec_bdy_width, spec_zone, relax_zone, &
    303. 

  303.                                dtbc, fcx, gcx,             &
    304. 

  304.                                ids,ide, jds,jde, 1  ,1  ,  & ! domain dims
    305. 

  305.                                ims,ime, jms,jme, 1  ,1  ,  & ! memory dims
    306. 

  306.                                ips,ipe, jps,jpe, 1  ,1  ,  & ! patch  dims
    307. 

  307.                                its,ite, jts,jte, 1  ,1   )
    308. 

  308. 

  309.          IF( config_flags%nested) THEN
    310. 

  310. 

  311.            i_start = max(its-1, ids)
    312. 

  312.            i_end = min(ite+1, ide-1)
    313. 

  313.            j_start = max(jts-1, jds)
    314. 

  314.            j_end = min(jte+1, jde-1)
    315. 

  315. 

  316.            DO j=j_start,j_end
    317. 

  317.            DO k=kts,kte
    318. 

  318.            DO i=i_start,i_end
    319. 

  319.               rfield(i,k,j) = w(i,k,j)*mut(i,j)
    320. 

  320.            ENDDO
    321. 

  321.            ENDDO
    322. 

  322.            ENDDO
    323. 

  323.            
    324. 

  324.            CALL relax_bdytend_tile ( rfield, rw_tendf,             &
    325. 

  325.                                w_bxs,w_bxe,w_bys,w_bye, w_btxs,w_btxe,w_btys,w_btye,       &
    326. 

  326.                                'h'        , config_flags, &
    327. 

  327.                                spec_bdy_width, spec_zone, relax_zone, &
    328. 

  328.                                dtbc, fcx, gcx,             &
    329. 

  329.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    330. 

  330.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    331. 

  331.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    332. 

  332.                                its,ite, jts,jte, kts,kte,  &
    333. 

  333.                                its-1, ite+1, jts-1,jte+1,kts,kte )  ! dims of first argument
    334. 

  334. 

  335.         END IF
    336. 

  336. 

  337.    END SUBROUTINE relax_bdy_dry 
    338. 

  338. !------------------------------------------------------------------------
    339. 

  339.    SUBROUTINE relax_bdy_scalar ( scalar_tend,                &
    340. 

  340.                                  scalar, mu,                 &
    341. 

  341.                                  scalar_bxs,scalar_bxe,scalar_bys,scalar_bye, &
    342. 

  342.                                  scalar_btxs,scalar_btxe,scalar_btys,scalar_btye, &
    343. 

  343.                                  spec_bdy_width, spec_zone, relax_zone,       &
    344. 

  344.                                  dtbc, fcx, gcx,             &
    345. 

  345.                                  config_flags,               &
    346. 

  346.                                  ids,ide, jds,jde, kds,kde,  & ! domain dims
    347. 

  347.                                  ims,ime, jms,jme, kms,kme,  & ! memory dims
    348. 

  348.                                  ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    349. 

  349.                                  its, ite, jts, jte, kts, kte)
    350. 

  350.    IMPLICIT NONE
    351. 

  351. 

  352.    !  Input data.
    353. 

  353.    TYPE( grid_config_rec_type ) config_flags
    354. 

  354. 

  355.    INTEGER ,               INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
    356. 

  356.                                             ims, ime, jms, jme, kms, kme, &
    357. 

  357.                                             ips, ipe, jps, jpe, kps, kpe, & 
    358. 

  358.                                             its, ite, jts, jte, kts, kte
    359. 

  359.    INTEGER ,               INTENT(IN   ) :: spec_bdy_width, spec_zone, relax_zone
    360. 

  360. 

  361.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(IN   ) :: scalar
    362. 

  362.    REAL , DIMENSION( ims:ime , jms:jme  ) , INTENT(IN   ) :: mu
    363. 

  363.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(INOUT) :: scalar_tend
    364. 

  364.    REAL , DIMENSION( spec_bdy_width) , INTENT(IN   ) :: fcx, gcx
    365. 

  365. 

  366.    REAL,  DIMENSION( jms:jme , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: scalar_bxs,scalar_bxe, &
    367. 

  367.                                                                                scalar_btxs,scalar_btxe
    368. 

  368.    REAL,  DIMENSION( ims:ime , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: scalar_bys,scalar_bye, &
    369. 

  369.                                                                                scalar_btys,scalar_btye
    370. 

  370.    REAL, INTENT(IN   ) :: dtbc
    371. 

  371. !Local
    372. 

  372.    INTEGER :: i,j,k, i_start, i_end, j_start, j_end
    373. 

  373.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) :: rscalar
    374. 

  374. 

  375. ! rscalar will be calculated beyond tile limits because relax_bdytend
    376. 

  376. !   requires a 5-point stencil, and this avoids need for inter-tile/patch 
    377. 

  377. !   communication here
    378. 

  378.            i_start = max(its-1, ids)
    379. 

  379.            i_end = min(ite+1, ide-1)
    380. 

  380.            j_start = max(jts-1, jds)
    381. 

  381.            j_end = min(jte+1, jde-1)
    382. 

  382. 

  383.            DO j=j_start,j_end
    384. 

  384.            DO k=kts,min(kte,kde-1)
    385. 

  385.            DO i=i_start,i_end
    386. 

  386.               rscalar(i,k,j) = scalar(i,k,j)*mu(i,j)
    387. 

  387.            ENDDO
    388. 

  388.            ENDDO
    389. 

  389.            ENDDO
    390. 

  390. 

  391.            CALL relax_bdytend (rscalar, scalar_tend,             &
    392. 

  392.                                scalar_bxs,scalar_bxe,scalar_bys,scalar_bye, scalar_btxs,scalar_btxe,scalar_btys,scalar_btye,       &
    393. 

  393.                                'q'        , config_flags, &
    394. 

  394.                                spec_bdy_width, spec_zone, relax_zone, &
    395. 

  395.                                dtbc, fcx, gcx,             &
    396. 

  396.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    397. 

  397.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    398. 

  398.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    399. 

  399.                                its,ite, jts,jte, kts,kte )
    400. 

  400. 

  401. 

  402.    END SUBROUTINE relax_bdy_scalar 
    403. 

  403. 

  404. !------------------------------------------------------------------------
    405. 

  405.    SUBROUTINE spec_bdy_dry ( config_flags,                        &
    406. 

  406.                              ru_tend, rv_tend, ph_tend, t_tend,   &
    407. 

  407.                              rw_tend, mu_tend,                    &
    408. 

  408.                              u_bxs,u_bxe,u_bys,u_bye,             &
    409. 

  409.                              v_bxs,v_bxe,v_bys,v_bye,             &
    410. 

  410.                              ph_bxs,ph_bxe,ph_bys,ph_bye,         &
    411. 

  411.                              t_bxs,t_bxe,t_bys,t_bye,             &
    412. 

  412.                              w_bxs,w_bxe,w_bys,w_bye,             &
    413. 

  413.                              mu_bxs,mu_bxe,mu_bys,mu_bye,         &
    414. 

  414.                              u_btxs,u_btxe,u_btys,u_btye,         &
    415. 

  415.                              v_btxs,v_btxe,v_btys,v_btye,         &
    416. 

  416.                              ph_btxs,ph_btxe,ph_btys,ph_btye,     &
    417. 

  417.                              t_btxs,t_btxe,t_btys,t_btye,         &
    418. 

  418.                              w_btxs,w_btxe,w_btys,w_btye,         &
    419. 

  419.                              mu_btxs,mu_btxe,mu_btys,mu_btye,     &
    420. 

  420.                              spec_bdy_width, spec_zone,           &
    421. 

  421.                              ids,ide, jds,jde, kds,kde,  & ! domain dims
    422. 

  422.                              ims,ime, jms,jme, kms,kme,  & ! memory dims
    423. 

  423.                              ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    424. 

  424.                              its, ite, jts, jte, kts, kte)
    425. 

  425.    IMPLICIT NONE
    426. 

  426. 

  427.    !  Input data.
    428. 

  428.    TYPE( grid_config_rec_type ) config_flags
    429. 

  429. 

  430. 

  431.    INTEGER ,               INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
    432. 

  432.                                             ims, ime, jms, jme, kms, kme, &
    433. 

  433.                                             ips, ipe, jps, jpe, kps, kpe, & 
    434. 

  434.                                             its, ite, jts, jte, kts, kte
    435. 

  435.    INTEGER ,               INTENT(IN   ) :: spec_bdy_width, spec_zone
    436. 

  436. 

  437.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(OUT  ) :: ru_tend, &
    438. 

  438.                                                                       rv_tend, &
    439. 

  439.                                                                       ph_tend, &
    440. 

  440.                                                                       rw_tend, &
    441. 

  441.                                                                       t_tend
    442. 

  442.    REAL , DIMENSION( ims:ime , jms:jme  ) , INTENT(OUT  )          :: mu_tend
    443. 

  443. 

  444.    REAL,  DIMENSION( jms:jme , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: u_bxs,u_bxe,  &
    445. 

  445.                                                                                v_bxs,v_bxe,  &
    446. 

  446.                                                                                ph_bxs,ph_bxe, &
    447. 

  447.                                                                                w_bxs,w_bxe, &
    448. 

  448.                                                                                t_bxs,t_bxe,  &
    449. 

  449.                                                                                u_btxs,u_btxe, &
    450. 

  450.                                                                                v_btxs,v_btxe, &
    451. 

  451.                                                                                ph_btxs,ph_btxe, &
    452. 

  452.                                                                                w_btxs,w_btxe, &
    453. 

  453.                                                                                t_btxs,t_btxe
    454. 

  454. 

  455.    REAL,  DIMENSION( ims:ime , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: u_bys,u_bye,  &
    456. 

  456.                                                                                v_bys,v_bye,  &
    457. 

  457.                                                                                ph_bys,ph_bye, &
    458. 

  458.                                                                                w_bys,w_bye, &
    459. 

  459.                                                                                t_bys,t_bye,  &
    460. 

  460.                                                                                u_btys,u_btye, &
    461. 

  461.                                                                                v_btys,v_btye, &
    462. 

  462.                                                                                ph_btys,ph_btye, &
    463. 

  463.                                                                                w_btys,w_btye, &
    464. 

  464.                                                                                t_btys,t_btye
    465. 

  465. 

  466.    REAL,  DIMENSION( jms:jme , 1:1 ,     spec_bdy_width    ), INTENT(IN   ) :: mu_bxs,mu_bxe, &
    467. 

  467.                                                                                mu_btxs,mu_btxe
    468. 

  468. 

  469.    REAL,  DIMENSION( ims:ime , 1:1 ,     spec_bdy_width    ), INTENT(IN   ) :: mu_bys,mu_bye, &
    470. 

  470.                                                                                mu_btys,mu_btye
    471. 

  471.          CALL spec_bdytend (   ru_tend,                &
    472. 

  472.                                u_bxs,u_bxe,u_bys,u_bye, u_btxs,u_btxe,u_btys,u_btye,    &
    473. 

  473.                                'u'     , config_flags, &
    474. 

  474.                                spec_bdy_width, spec_zone, &
    475. 

  475.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    476. 

  476.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    477. 

  477.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    478. 

  478.                                its,ite, jts,jte, kts,kte )
    479. 

  479.          CALL spec_bdytend (   rv_tend,                &
    480. 

  480.                                v_bxs,v_bxe,v_bys,v_bye, v_btxs,v_btxe,v_btys,v_btye,    &
    481. 

  481.                                'v'     , config_flags, &
    482. 

  482.                                spec_bdy_width, spec_zone, &
    483. 

  483.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    484. 

  484.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    485. 

  485.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    486. 

  486.                                its,ite, jts,jte, kts,kte )
    487. 

  487.          CALL spec_bdytend (   ph_tend,                &
    488. 

  488.                                ph_bxs,ph_bxe,ph_bys,ph_bye, ph_btxs,ph_btxe,ph_btys,ph_btye,    &
    489. 

  489.                                'h'     , config_flags, &
    490. 

  490.                                spec_bdy_width, spec_zone, &
    491. 

  491.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    492. 

  492.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    493. 

  493.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    494. 

  494.                                its,ite, jts,jte, kts,kte )
    495. 

  495.          CALL spec_bdytend (   t_tend,                &
    496. 

  496.                                t_bxs,t_bxe,t_bys,t_bye, t_btxs,t_btxe,t_btys,t_btye,    &
    497. 

  497.                                't'     , config_flags, &
    498. 

  498.                                spec_bdy_width, spec_zone, &
    499. 

  499.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    500. 

  500.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    501. 

  501.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    502. 

  502.                                its,ite, jts,jte, kts,kte )
    503. 

  503.          CALL spec_bdytend (   mu_tend,                &
    504. 

  504.                                mu_bxs,mu_bxe,mu_bys,mu_bye, mu_btxs,mu_btxe,mu_btys,mu_btye,       &
    505. 

  505.                                'm'     , config_flags, &
    506. 

  506.                                spec_bdy_width, spec_zone, &
    507. 

  507.                                ids,ide, jds,jde, 1  ,1  ,  & ! domain dims
    508. 

  508.                                ims,ime, jms,jme, 1  ,1  ,  & ! memory dims
    509. 

  509.                                ips,ipe, jps,jpe, 1  ,1  ,  & ! patch  dims
    510. 

  510.                                its,ite, jts,jte, 1  ,1   )
    511. 

  511. 

  512.          if(config_flags%nested)                           &
    513. 

  513.          CALL spec_bdytend (   rw_tend,                    &
    514. 

  514.                                w_bxs,w_bxe,w_bys,w_bye, w_btxs,w_btxe,w_btys,w_btye,                  &
    515. 

  515.                                'h'     , config_flags,     &
    516. 

  516.                                spec_bdy_width, spec_zone,  &
    517. 

  517.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    518. 

  518.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    519. 

  519.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    520. 

  520.                                its,ite, jts,jte, kts,kte )
    521. 

  521. 

  522.    END SUBROUTINE spec_bdy_dry 
    523. 

  523. 

  524. !------------------------------------------------------------------------
    525. 

  525.    SUBROUTINE spec_bdy_scalar ( scalar_tend,    &
    526. 

  526.                           scalar_bxs,scalar_bxe,scalar_bys,scalar_bye,  &
    527. 

  527.                           scalar_btxs,scalar_btxe,scalar_btys,scalar_btye, &
    528. 

  528.                           spec_bdy_width, spec_zone,                   &
    529. 

  529.                           config_flags,               &
    530. 

  530.                           ids,ide, jds,jde, kds,kde,  & ! domain dims
    531. 

  531.                           ims,ime, jms,jme, kms,kme,  & ! memory dims
    532. 

  532.                           ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    533. 

  533.                           its, ite, jts, jte, kts, kte)
    534. 

  534.    IMPLICIT NONE
    535. 

  535. 

  536.    !  Input data.
    537. 

  537.    TYPE( grid_config_rec_type ) config_flags
    538. 

  538. 

  539. 

  540.    INTEGER ,               INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
    541. 

  541.                                             ims, ime, jms, jme, kms, kme, &
    542. 

  542.                                             ips, ipe, jps, jpe, kps, kpe, & 
    543. 

  543.                                             its, ite, jts, jte, kts, kte
    544. 

  544.    INTEGER ,               INTENT(IN   ) :: spec_bdy_width, spec_zone
    545. 

  545. 

  546.    REAL , DIMENSION( ims:ime , kms:kme, jms:jme  ) , INTENT(OUT  ) :: scalar_tend
    547. 

  547. 

  548.    REAL,  DIMENSION( jms:jme , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: scalar_bxs,scalar_bxe, &
    549. 

  549.                                                                                scalar_btxs,scalar_btxe
    550. 

  550. 

  551.    REAL,  DIMENSION( ims:ime , kds:kde , spec_bdy_width    ), INTENT(IN   ) :: scalar_bys,scalar_bye, &
    552. 

  552.                                                                                scalar_btys,scalar_btye
    553. 

  553. 

  554. !Local
    555. 

  555.    INTEGER :: i,j,k
    556. 

  556. 

  557. 

  558.          CALL spec_bdytend (   scalar_tend,                &
    559. 

  559.                                scalar_bxs,scalar_bxe,scalar_bys,scalar_bye, scalar_btxs,scalar_btxe,scalar_btys,scalar_btye,    &
    560. 

  560.                                'q'     , config_flags, &
    561. 

  561.                                spec_bdy_width, spec_zone, &
    562. 

  562.                                ids,ide, jds,jde, kds,kde,  & ! domain dims
    563. 

  563.                                ims,ime, jms,jme, kms,kme,  & ! memory dims
    564. 

  564.                                ips,ipe, jps,jpe, kps,kpe,  & ! patch  dims
    565. 

  565.                                its,ite, jts,jte, kts,kte )
    566. 

  566. 

  567. 

  568.    END SUBROUTINE spec_bdy_scalar 
    569. 

  569. 

  570. !------------------------------------------------------------------------
    571. 

  571. 

  572.    SUBROUTINE set_phys_bc_dry_1( config_flags, u_1, u_2, v_1, v_2,   &
    573. 

  573.                                  rw_1, rw_2, w_1, w_2,           &
    574. 

  574.                                  t_1, t_2, tp_1, tp_2, pp, pip,  &
    575. 

  575.                                  ids,ide, jds,jde, kds,kde,      &
    576. 

  576.                                  ims,ime, jms,jme, kms,kme,      &
    577. 

  577.                                  ips,ipe, jps,jpe, kps,kpe,      &
    578. 

  578.                                  its,ite, jts,jte, kts,kte      )
    579. 

  579. 

  580. !
    581. 

  581. !  this is just a wraper to call the boundary condition routines
    582. 

  582. !  for each variable
    583. 

  583. !
    584. 

  584. 

  585.       IMPLICIT NONE
    586. 

  586. 

  587.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    588. 

  588.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    589. 

  589.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    590. 

  590.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    591. 

  591. 

  592.       TYPE( grid_config_rec_type ) config_flags
    593. 

  593. 

  594.       REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: &
    595. 

  595.            u_1,u_2, v_1, v_2, rw_1, rw_2, w_1, w_2,                  &
    596. 

  596.            t_1, t_2, tp_1, tp_2, pp, pip
    597. 

  597. 

  598. 

  599. 

  600.       CALL set_physical_bc3d( u_1  , 'u', config_flags,               &
    601. 

  601.                               ids, ide, jds, jde, kds, kde,       &
    602. 

  602.                               ims, ime, jms, jme, kms, kme,       &
    603. 

  603.                               ips, ipe, jps, jpe, kps, kpe,       &
    604. 

  604.                               its, ite, jts, jte, kts, kte )
    605. 

  605.       CALL set_physical_bc3d( u_2  , 'u', config_flags,               &
    606. 

  606.                               ids, ide, jds, jde, kds, kde,       &
    607. 

  607.                               ims, ime, jms, jme, kms, kme,       &
    608. 

  608.                               ips, ipe, jps, jpe, kps, kpe,       &
    609. 

  609.                               its, ite, jts, jte, kts, kte )
    610. 

  610.       CALL set_physical_bc3d( v_1  , 'v', config_flags,               &
    611. 

  611.                               ids, ide, jds, jde, kds, kde,       &
    612. 

  612.                               ims, ime, jms, jme, kms, kme,       &
    613. 

  613.                               ips, ipe, jps, jpe, kps, kpe,       &
    614. 

  614.                               its, ite, jts, jte, kts, kte )
    615. 

  615.       CALL set_physical_bc3d( v_2  , 'v', config_flags,               &
    616. 

  616.                               ids, ide, jds, jde, kds, kde,       &
    617. 

  617.                               ims, ime, jms, jme, kms, kme,       &
    618. 

  618.                               ips, ipe, jps, jpe, kps, kpe,       &
    619. 

  619.                               its, ite, jts, jte, kts, kte )
    620. 

  620.       CALL set_physical_bc3d( rw_1 , 'w', config_flags,               &
    621. 

  621.                               ids, ide, jds, jde, kds, kde,       &
    622. 

  622.                               ims, ime, jms, jme, kms, kme,       &
    623. 

  623.                               ips, ipe, jps, jpe, kps, kpe,       &
    624. 

  624.                               its, ite, jts, jte, kts, kte )
    625. 

  625.       CALL set_physical_bc3d( rw_2 , 'w', config_flags,               &
    626. 

  626.                               ids, ide, jds, jde, kds, kde,       &
    627. 

  627.                               ims, ime, jms, jme, kms, kme,       &
    628. 

  628.                               ips, ipe, jps, jpe, kps, kpe,       &
    629. 

  629.                               its, ite, jts, jte, kts, kte )
    630. 

  630.       CALL set_physical_bc3d( w_1  , 'w', config_flags,               &
    631. 

  631.                               ids, ide, jds, jde, kds, kde,       &
    632. 

  632.                               ims, ime, jms, jme, kms, kme,       &
    633. 

  633.                               ips, ipe, jps, jpe, kps, kpe,       &
    634. 

  634.                               its, ite, jts, jte, kts, kte )
    635. 

  635.       CALL set_physical_bc3d( w_2  , 'w', config_flags,               &
    636. 

  636.                               ids, ide, jds, jde, kds, kde,       &
    637. 

  637.                               ims, ime, jms, jme, kms, kme,       &
    638. 

  638.                               ips, ipe, jps, jpe, kps, kpe,       &
    639. 

  639.                               its, ite, jts, jte, kts, kte )
    640. 

  640.       CALL set_physical_bc3d( t_1, 'p', config_flags,                 &
    641. 

  641.                               ids, ide, jds, jde, kds, kde,       &
    642. 

  642.                               ims, ime, jms, jme, kms, kme,       &
    643. 

  643.                               ips, ipe, jps, jpe, kps, kpe,       &
    644. 

  644.                               its, ite, jts, jte, kts, kte )
    645. 

  645.       CALL set_physical_bc3d( t_2, 'p', config_flags,                 &
    646. 

  646.                               ids, ide, jds, jde, kds, kde,       &
    647. 

  647.                               ims, ime, jms, jme, kms, kme,       &
    648. 

  648.                               ips, ipe, jps, jpe, kps, kpe,       &
    649. 

  649.                               its, ite, jts, jte, kts, kte )
    650. 

  650.       CALL set_physical_bc3d( tp_1, 'p', config_flags,                &
    651. 

  651.                               ids, ide, jds, jde, kds, kde,       &
    652. 

  652.                               ims, ime, jms, jme, kms, kme,       &
    653. 

  653.                               ips, ipe, jps, jpe, kps, kpe,       &
    654. 

  654.                               its, ite, jts, jte, kts, kte )
    655. 

  655.       CALL set_physical_bc3d( tp_2, 'p', config_flags,                &
    656. 

  656.                               ids, ide, jds, jde, kds, kde,       &
    657. 

  657.                               ims, ime, jms, jme, kms, kme,       &
    658. 

  658.                               ips, ipe, jps, jpe, kps, kpe,       &
    659. 

  659.                               its, ite, jts, jte, kts, kte )
    660. 

  660.       CALL set_physical_bc3d( pp , 'p', config_flags,                 &
    661. 

  661.                               ids, ide, jds, jde, kds, kde,       &
    662. 

  662.                               ims, ime, jms, jme, kms, kme,       &
    663. 

  663.                               ips, ipe, jps, jpe, kps, kpe,       &
    664. 

  664.                               its, ite, jts, jte, kts, kte )
    665. 

  665.       CALL set_physical_bc3d( pip , 'p', config_flags,                &
    666. 

  666.                               ids, ide, jds, jde, kds, kde,       &
    667. 

  667.                               ims, ime, jms, jme, kms, kme,       &
    668. 

  668.                               ips, ipe, jps, jpe, kps, kpe,       &
    669. 

  669.                               its, ite, jts, jte, kts, kte )
    670. 

  670. 

  671.   END SUBROUTINE set_phys_bc_dry_1
    672. 

  672. 

  673. !--------------------------------------------------------------
    674. 

  674. 

  675.    SUBROUTINE set_phys_bc_dry_2( config_flags,                     &
    676. 

  676.                                  u_1, u_2, v_1, v_2, w_1, w_2,     &
    677. 

  677.                                  t_1, t_2, ph_1, ph_2, mu_1, mu_2, &
    678. 

  678.                                  ids,ide, jds,jde, kds,kde,        &
    679. 

  679.                                  ims,ime, jms,jme, kms,kme,        &
    680. 

  680.                                  ips,ipe, jps,jpe, kps,kpe,        &
    681. 

  681.                                  its,ite, jts,jte, kts,kte        )
    682. 

  682. 

  683. !
    684. 

  684. !  this is just a wraper to call the boundary condition routines
    685. 

  685. !  for each variable
    686. 

  686. !
    687. 

  687. 

  688.       IMPLICIT NONE
    689. 

  689. 

  690.       TYPE( grid_config_rec_type ) config_flags
    691. 

  691. 

  692.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    693. 

  693.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    694. 

  694.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    695. 

  695.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    696. 

  696. 

  697.       REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: &
    698. 

  698.          u_1, u_2, v_1, v_2, w_1, w_2,                       &
    699. 

  699.          t_1, t_2, ph_1, ph_2
    700. 

  700. 

  701.       REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: &
    702. 

  702.                              mu_1, mu_2
    703. 

  703. 

  704. 

  705.       CALL set_physical_bc3d( u_1, 'U', config_flags,           &
    706. 

  706.                               ids, ide, jds, jde, kds, kde, &
    707. 

  707.                               ims, ime, jms, jme, kms, kme, &
    708. 

  708.                               ips, ipe, jps, jpe, kps, kpe, &
    709. 

  709.                               its, ite, jts, jte, kts, kte )
    710. 

  710. 

  711.       CALL set_physical_bc3d( u_2, 'U', config_flags,           &
    712. 

  712.                               ids, ide, jds, jde, kds, kde, &
    713. 

  713.                               ims, ime, jms, jme, kms, kme, &
    714. 

  714.                               ips, ipe, jps, jpe, kps, kpe, &
    715. 

  715.                               its, ite, jts, jte, kts, kte )
    716. 

  716. 

  717.       CALL set_physical_bc3d( v_1 , 'V', config_flags,          &
    718. 

  718.                               ids, ide, jds, jde, kds, kde, &
    719. 

  719.                               ims, ime, jms, jme, kms, kme, &
    720. 

  720.                               ips, ipe, jps, jpe, kps, kpe, &
    721. 

  721.                               its, ite, jts, jte, kts, kte )
    722. 

  722.       CALL set_physical_bc3d( v_2 , 'V', config_flags,          &
    723. 

  723.                               ids, ide, jds, jde, kds, kde, &
    724. 

  724.                               ims, ime, jms, jme, kms, kme, &
    725. 

  725.                               ips, ipe, jps, jpe, kps, kpe, &
    726. 

  726.                               its, ite, jts, jte, kts, kte )
    727. 

  727. 

  728.       CALL set_physical_bc3d( w_1, 'w', config_flags,           &
    729. 

  729.                               ids, ide, jds, jde, kds, kde, &
    730. 

  730.                               ims, ime, jms, jme, kms, kme, &
    731. 

  731.                               ips, ipe, jps, jpe, kps, kpe, &
    732. 

  732.                               its, ite, jts, jte, kts, kte )
    733. 

  733.       CALL set_physical_bc3d( w_2, 'w', config_flags,           &
    734. 

  734.                               ids, ide, jds, jde, kds, kde, &
    735. 

  735.                               ims, ime, jms, jme, kms, kme, &
    736. 

  736.                               ips, ipe, jps, jpe, kps, kpe, &
    737. 

  737.                               its, ite, jts, jte, kts, kte )
    738. 

  738. 

  739.       CALL set_physical_bc3d( t_1, 'p', config_flags,           &
    740. 

  740.                               ids, ide, jds, jde, kds, kde, &
    741. 

  741.                               ims, ime, jms, jme, kms, kme, &
    742. 

  742.                               ips, ipe, jps, jpe, kps, kpe, &
    743. 

  743.                               its, ite, jts, jte, kts, kte )
    744. 

  744. 

  745.       CALL set_physical_bc3d( t_2, 'p', config_flags,           &
    746. 

  746.                               ids, ide, jds, jde, kds, kde, &
    747. 

  747.                               ims, ime, jms, jme, kms, kme, &
    748. 

  748.                               ips, ipe, jps, jpe, kps, kpe, &
    749. 

  749.                               its, ite, jts, jte, kts, kte )
    750. 

  750. 

  751.       CALL set_physical_bc3d( ph_1 , 'w', config_flags,         &
    752. 

  752.                               ids, ide, jds, jde, kds, kde, &
    753. 

  753.                               ims, ime, jms, jme, kms, kme, &
    754. 

  754.                               ips, ipe, jps, jpe, kps, kpe, &
    755. 

  755.                               its, ite, jts, jte, kts, kte )
    756. 

  756. 

  757.       CALL set_physical_bc3d( ph_2 , 'w', config_flags,         &
    758. 

  758.                               ids, ide, jds, jde, kds, kde, &
    759. 

  759.                               ims, ime, jms, jme, kms, kme, &
    760. 

  760.                               ips, ipe, jps, jpe, kps, kpe, &
    761. 

  761.                               its, ite, jts, jte, kts, kte )
    762. 

  762. 

  763.       CALL set_physical_bc2d( mu_1, 't', config_flags, &
    764. 

  764.                               ids, ide, jds, jde,  &
    765. 

  765.                               ims, ime, jms, jme,  &
    766. 

  766.                               ips, ipe, jps, jpe,  &
    767. 

  767.                               its, ite, jts, jte  )
    768. 

  768. 

  769.       CALL set_physical_bc2d( mu_2, 't', config_flags, &
    770. 

  770.                               ids, ide, jds, jde,  &
    771. 

  771.                               ims, ime, jms, jme,  &
    772. 

  772.                               ips, ipe, jps, jpe,  &
    773. 

  773.                               its, ite, jts, jte  )
    774. 

  774. 

  775.    END SUBROUTINE set_phys_bc_dry_2
    776. 

  776. 

  777. !------------------------------------------------------------------------
    778. 

  778. 

  779.    SUBROUTINE set_phys_bc_smallstep_1( config_flags, ru_1, du, rv_1, dv,   &
    780. 

  780.                                        ids,ide, jds,jde, kds,kde,      &
    781. 

  781.                                        ims,ime, jms,jme, kms,kme,      &
    782. 

  782.                                        ips,ipe, jps,jpe, kps,kpe,      &
    783. 

  783.                                        its,ite, jts,jte, kts,kte      )
    784. 

  784. 

  785. !
    786. 

  786. !  this is just a wraper to call the boundary condition routines
    787. 

  787. !  for each variable
    788. 

  788. !
    789. 

  789. 

  790.       IMPLICIT NONE
    791. 

  791. 

  792.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    793. 

  793.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    794. 

  794.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    795. 

  795.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    796. 

  796. 

  797.       TYPE( grid_config_rec_type ) config_flags
    798. 

  798. 

  799.       REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: &
    800. 

  800.            ru_1,du, rv_1, dv
    801. 

  801. 

  802.       CALL set_physical_bc3d( ru_1  , 'u', config_flags,              &
    803. 

  803.                               ids, ide, jds, jde, kds, kde,       &
    804. 

  804.                               ims, ime, jms, jme, kms, kme,       &
    805. 

  805.                               ips, ipe, jps, jpe, kps, kpe,       &
    806. 

  806.                               its, ite, jts, jte, kts, kde )
    807. 

  807.       CALL set_physical_bc3d( du , 'u', config_flags,                 &
    808. 

  808.                               ids, ide, jds, jde, kds, kde,       &
    809. 

  809.                               ims, ime, jms, jme, kms, kme,       &
    810. 

  810.                               ips, ipe, jps, jpe, kps, kpe,       &
    811. 

  811.                               its, ite, jts, jte, kts, kde )
    812. 

  812.       CALL set_physical_bc3d( rv_1  , 'v', config_flags,              &
    813. 

  813.                               ids, ide, jds, jde, kds, kde,       &
    814. 

  814.                               ims, ime, jms, jme, kms, kme,       &
    815. 

  815.                               ips, ipe, jps, jpe, kps, kpe,       &
    816. 

  816.                               its, ite, jts, jte, kts, kde )
    817. 

  817.       CALL set_physical_bc3d( dv  , 'v', config_flags,                &
    818. 

  818.                               ids, ide, jds, jde, kds, kde,       &
    819. 

  819.                               ims, ime, jms, jme, kms, kme,       &
    820. 

  820.                               ips, ipe, jps, jpe, kps, kpe,       &
    821. 

  821.                               its, ite, jts, jte, kts, kde )
    822. 

  822. 

  823.   END SUBROUTINE set_phys_bc_smallstep_1
    824. 

  824. 

  825. !-------------------------------------------------------------------
    826. 

  826. 

  827.    SUBROUTINE rk_phys_bc_dry_1( config_flags, u, v, rw, w,  &
    828. 

  828.                                 muu, muv, mut, php, alt, p, &
    829. 

  829.                                 ids,ide, jds,jde, kds,kde,  &
    830. 

  830.                                 ims,ime, jms,jme, kms,kme,  &
    831. 

  831.                                 ips,ipe, jps,jpe, kps,kpe,  &
    832. 

  832.                                 its,ite, jts,jte, kts,kte  )
    833. 

  833. 

  834. !
    835. 

  835. !  this is just a wraper to call the boundary condition routines
    836. 

  836. !  for each variable
    837. 

  837. !
    838. 

  838. 

  839.       IMPLICIT NONE
    840. 

  840. 

  841.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    842. 

  842.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    843. 

  843.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    844. 

  844.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    845. 

  845. 

  846.       TYPE( grid_config_rec_type ) config_flags
    847. 

  847. 

  848.       REAL, DIMENSION( ims:ime, kms:kme, jms:jme ),                    &
    849. 

  849.                                 INTENT(INOUT) ::  u, v, rw, w, php, alt, p
    850. 

  850.       REAL, DIMENSION( ims:ime, jms:jme ),                             &
    851. 

  851.                                 INTENT(INOUT) ::    muu, muv, mut
    852. 

  852. 

  853.       CALL set_physical_bc3d( u  , 'u', config_flags,             &
    854. 

  854.                               ids, ide, jds, jde, kds, kde,       &
    855. 

  855.                               ims, ime, jms, jme, kms, kme,       &
    856. 

  856.                               ips, ipe, jps, jpe, kps, kpe,       &
    857. 

  857.                               its, ite, jts, jte, kts, kte )
    858. 

  858.       CALL set_physical_bc3d( v  , 'v', config_flags,             &
    859. 

  859.                               ids, ide, jds, jde, kds, kde,       &
    860. 

  860.                               ims, ime, jms, jme, kms, kme,       &
    861. 

  861.                               ips, ipe, jps, jpe, kps, kpe,       &
    862. 

  862.                               its, ite, jts, jte, kts, kte )
    863. 

  863.       CALL set_physical_bc3d(rw , 'w', config_flags,              &
    864. 

  864.                               ids, ide, jds, jde, kds, kde,       &
    865. 

  865.                               ims, ime, jms, jme, kms, kme,       &
    866. 

  866.                               ips, ipe, jps, jpe, kps, kpe,       &
    867. 

  867.                               its, ite, jts, jte, kts, kte )
    868. 

  868.       CALL set_physical_bc3d( w , 'w', config_flags,              &
    869. 

  869.                               ids, ide, jds, jde, kds, kde,       &
    870. 

  870.                               ims, ime, jms, jme, kms, kme,       &
    871. 

  871.                               ips, ipe, jps, jpe, kps, kpe,       &
    872. 

  872.                               its, ite, jts, jte, kts, kte )
    873. 

  873.       CALL set_physical_bc3d( php , 'w', config_flags,            &
    874. 

  874.                               ids, ide, jds, jde, kds, kde,       &
    875. 

  875.                               ims, ime, jms, jme, kms, kme,       &
    876. 

  876.                               ips, ipe, jps, jpe, kps, kpe,       &
    877. 

  877.                               its, ite, jts, jte, kts, kte )
    878. 

  878.       CALL set_physical_bc3d( alt, 't', config_flags,             &
    879. 

  879.                               ids, ide, jds, jde, kds, kde,       &
    880. 

  880.                               ims, ime, jms, jme, kms, kme,       &
    881. 

  881.                               ips, ipe, jps, jpe, kps, kpe,       &
    882. 

  882.                               its, ite, jts, jte, kts, kte )
    883. 

  883. 

  884.       CALL set_physical_bc3d( p, 'p', config_flags,               &
    885. 

  885.                               ids, ide, jds, jde, kds, kde,       &
    886. 

  886.                               ims, ime, jms, jme, kms, kme,       &
    887. 

  887.                               ips, ipe, jps, jpe, kps, kpe,       &
    888. 

  888.                               its, ite, jts, jte, kts, kte )
    889. 

  889. 

  890.       CALL set_physical_bc2d( muu, 'u', config_flags,  &
    891. 

  891.                               ids, ide, jds, jde,      &
    892. 

  892.                               ims, ime, jms, jme,      &
    893. 

  893.                               ips, ipe, jps, jpe,      &
    894. 

  894.                               its, ite, jts, jte  )
    895. 

  895. 

  896.       CALL set_physical_bc2d( muv, 'v', config_flags,  &
    897. 

  897.                               ids, ide, jds, jde,      &
    898. 

  898.                               ims, ime, jms, jme,      &
    899. 

  899.                               ips, ipe, jps, jpe,      &
    900. 

  900.                               its, ite, jts, jte  )
    901. 

  901. 

  902.       CALL set_physical_bc2d( mut, 't', config_flags,  &
    903. 

  903.                               ids, ide, jds, jde,      &
    904. 

  904.                               ims, ime, jms, jme,      &
    905. 

  905.                               ips, ipe, jps, jpe,      &
    906. 

  906.                               its, ite, jts, jte  )
    907. 

  907. 

  908.   END SUBROUTINE rk_phys_bc_dry_1
    909. 

  909. 

  910. !------------------------------------------------------------------------
    911. 

  911. 

  912.   SUBROUTINE rk_phys_bc_dry_2( config_flags, u, v, w,      &
    913. 

  913.                                t, ph, mu,                  &
    914. 

  914.                                ids,ide, jds,jde, kds,kde,  &
    915. 

  915.                                ims,ime, jms,jme, kms,kme,  &
    916. 

  916.                                ips,ipe, jps,jpe, kps,kpe,  &
    917. 

  917.                                its,ite, jts,jte, kts,kte  )
    918. 

  918. 

  919. !
    920. 

  920. !  this is just a wraper to call the boundary condition routines
    921. 

  921. !  for each variable
    922. 

  922. !
    923. 

  923. 

  924.       IMPLICIT NONE
    925. 

  925. 

  926.       INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
    927. 

  927.       INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
    928. 

  928.       INTEGER,      INTENT(IN   )    :: ips,ipe, jps,jpe, kps,kpe
    929. 

  929.       INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
    930. 

  930. 

  931.       TYPE( grid_config_rec_type ) config_flags
    932. 

  932. 

  933.       REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: &
    934. 

  934.                              u, v, w, t, ph
    935. 

  935. 

  936.       REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: &
    937. 

  937.                              mu
    938. 

  938. 

  939.       CALL set_physical_bc3d( u   , 'U', config_flags,            &
    940. 

  940.                               ids, ide, jds, jde, kds, kde,       &
    941. 

  941.                               ims, ime, jms, jme, kms, kme,       &
    942. 

  942.                               ips, ipe, jps, jpe, kps, kpe,       &
    943. 

  943.                               its, ite, jts, jte, kts, kte )
    944. 

  944.       CALL set_physical_bc3d( v   , 'V', config_flags,            &
    945. 

  945.                               ids, ide, jds, jde, kds, kde,       &
    946. 

  946.                               ims, ime, jms, jme, kms, kme,       &
    947. 

  947.                               ips, ipe, jps, jpe, kps, kpe,       &
    948. 

  948.                               its, ite, jts, jte, kts, kte )
    949. 

  949.       CALL set_physical_bc3d( w  , 'w', config_flags,             &
    950. 

  950.                               ids, ide, jds, jde, kds, kde,       &
    951. 

  951.                               ims, ime, jms, jme, kms, kme,       &
    952. 

  952.                               ips, ipe, jps, jpe, kps, kpe,       &
    953. 

  953.                               its, ite, jts, jte, kts, kte )
    954. 

  954.       CALL set_physical_bc3d( t, 'p', config_flags,               &
    955. 

  955.                               ids, ide, jds, jde, kds, kde,       &
    956. 

  956.                               ims, ime, jms, jme, kms, kme,       &
    957. 

  957.                               ips, ipe, jps, jpe, kps, kpe,       &
    958. 

  958.                               its, ite, jts, jte, kts, kte )
    959. 

  959.       CALL set_physical_bc3d( ph  , 'w', config_flags,            &
    960. 

  960.                               ids, ide, jds, jde, kds, kde,       &
    961. 

  961.                               ims, ime, jms, jme, kms, kme…
    962.

Large files files are truncated, but you can click here to view the full file