Ignore:
Timestamp:
Feb 18, 2009 5:50:38 PM (13 years ago)
Author:
letzel
Message:
  • External pressure gradient (check_parameters, init_3d_model, header, modules, parin, prognostic_equations)
  • New topography case 'single_street_canyon'
File:
1 edited

Legend:

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

    r198 r240  
    44! Actual revisions:
    55! -----------------
    6 !
     6! external pressure gradient
    77!
    88! Former revisions:
     
    173173!--       Drag by plant canopy
    174174          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
     175
     176!
     177!--       External pressure gradient
     178          IF ( dp_external )  THEN
     179             DO  k = dp_level_ind_b+1, nzt
     180                tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
     181             ENDDO
     182          ENDIF
     183
    175184          CALL user_actions( i, j, 'u-tendency' )
    176185
     
    244253!--       Drag by plant canopy
    245254          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )     
     255
     256!
     257!--       External pressure gradient
     258          IF ( dp_external )  THEN
     259             DO  k = dp_level_ind_b+1, nzt
     260                tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
     261             ENDDO
     262          ENDIF
    246263
    247264          CALL user_actions( i, j, 'v-tendency' )
     
    878895             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
    879896
     897!
     898!--          External pressure gradient
     899             IF ( dp_external )  THEN
     900                DO  k = dp_level_ind_b+1, nzt
     901                   tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
     902                ENDDO
     903             ENDIF
     904
    880905             CALL user_actions( i, j, 'u-tendency' )
    881906
     
    931956!--          Drag by plant canopy
    932957             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )       
     958
     959!
     960!--          External pressure gradient
     961             IF ( dp_external )  THEN
     962                DO  k = dp_level_ind_b+1, nzt
     963                   tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
     964                ENDDO
     965             ENDIF
    933966
    934967             CALL user_actions( i, j, 'v-tendency' )
     
    13471380    IF ( plant_canopy )  CALL plant_canopy_model( 1 )
    13481381
     1382!
     1383!-- External pressure gradient
     1384    IF ( dp_external )  THEN
     1385       DO  i = nxlu, nxr
     1386          DO  j = nys, nyn
     1387             DO  k = dp_level_ind_b+1, nzt
     1388                tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
     1389             ENDDO
     1390          ENDDO
     1391       ENDDO
     1392    ENDIF
     1393
    13491394    CALL user_actions( 'u-tendency' )
    13501395
     
    14221467!-- Drag by plant canopy
    14231468    IF ( plant_canopy )  CALL plant_canopy_model( 2 )
     1469
     1470!
     1471!-- External pressure gradient
     1472    IF ( dp_external )  THEN
     1473       DO  i = nxl, nxr
     1474          DO  j = nysv, nyn
     1475             DO  k = dp_level_ind_b+1, nzt
     1476                tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
     1477             ENDDO
     1478          ENDDO
     1479       ENDDO
     1480    ENDIF
     1481
    14241482    CALL user_actions( 'v-tendency' )
    14251483
Note: See TracChangeset for help on using the changeset viewer.