Ignore:
Timestamp:
Jun 20, 2017 9:51:42 AM (4 years ago)
Author:
schwenkel
Message:

implementation of new bulk microphysics scheme

File:
1 edited

Legend:

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

    r2101 r2292  
    2525! -----------------
    2626! $Id$
     27! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
     28! includes two more prognostic equations for cloud drop concentration (nc) 
     29! and cloud water content (qc).
     30!
     31! 2101 2017-01-05 16:42:31Z suehring
    2732!
    2833! 2000 2016-08-20 18:09:15Z knoop
     
    10121017          ENDIF
    10131018
    1014        ELSEIF ( sk_char == 'qr' )  THEN
    1015 
    1016 !
    1017 !--       Rain water content boundary condition at the bottom boundary:
     1019       ELSEIF ( sk_char == 'qc' )  THEN
     1020
     1021!
     1022!--       Cloud water content boundary condition at the bottom boundary:
    10181023!--       Dirichlet (fixed surface rain water content).
    10191024          DO  i = nxl, nxr
     
    10261031
    10271032!
     1033!--       Cloud water content boundary condition at the top boundary: Dirichlet
     1034          DO  i = nxl, nxr
     1035             DO  j = nys, nyn
     1036                sk_p(nzt+2,j,i)   = sk_p(nzt+1,j,i)
     1037                sk_p(nzt+3,j,i)   = sk_p(nzt+1,j,i)
     1038             ENDDO
     1039          ENDDO
     1040
     1041       ELSEIF ( sk_char == 'qr' )  THEN
     1042
     1043!
     1044!--       Rain water content boundary condition at the bottom boundary:
     1045!--       Dirichlet (fixed surface rain water content).
     1046          DO  i = nxl, nxr
     1047             DO  j = nys, nyn
     1048                sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
     1049                sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
     1050                sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
     1051             ENDDO
     1052          ENDDO
     1053
     1054!
    10281055!--       Rain water content boundary condition at the top boundary: Dirichlet
     1056          DO  i = nxl, nxr
     1057             DO  j = nys, nyn
     1058                sk_p(nzt+2,j,i)   = sk_p(nzt+1,j,i)
     1059                sk_p(nzt+3,j,i)   = sk_p(nzt+1,j,i)
     1060             ENDDO
     1061          ENDDO
     1062
     1063       ELSEIF ( sk_char == 'nc' )  THEN
     1064
     1065!
     1066!--       Cloud drop concentration boundary condition at the bottom boundary:
     1067!--       Dirichlet (fixed surface cloud drop concentration).
     1068          DO  i = nxl, nxr
     1069             DO  j = nys, nyn
     1070                sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
     1071                sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
     1072                sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
     1073             ENDDO
     1074          ENDDO
     1075
     1076!
     1077!--       Cloud drop concentration boundary condition at the top boundary: Dirichlet
    10291078          DO  i = nxl, nxr
    10301079             DO  j = nys, nyn
Note: See TracChangeset for help on using the changeset viewer.