- Timestamp:
- Oct 21, 2020 3:09:15 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/bulk_cloud_model_mod.f90
r4754 r4755 1881 1881 !-- Here the the simple kessler scheme is used. 1882 1882 IF( microphysics_kessler ) THEN 1883 CALL autoconversion_kessler_ij( i, j )1884 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j )1883 CALL autoconversion_kessler_ij( i, j ) 1884 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j ) 1885 1885 ! 1886 1886 !-- Here the seifert beheng scheme is used. Cloud concentration is assumed to a constant value 1887 1887 !-- an qc a diagnostic value. 1888 1888 ELSEIF ( microphysics_seifert .AND. .NOT. microphysics_morrison ) THEN 1889 CALL adjust_cloud_ij( i, j )1889 CALL adjust_cloud_ij( i, j ) 1890 1890 IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics ) THEN 1891 CALL adjust_ice_ij( i, j )1892 CALL nucleation_ice_ij( i, j )1893 CALL homogeneous_freezing_cloud_ij( i, j )1894 CALL deposition_ice_ij( i, j )1891 CALL adjust_ice_ij( i, j ) 1892 CALL nucleation_ice_ij( i, j ) 1893 CALL homogeneous_freezing_cloud_ij( i, j ) 1894 CALL deposition_ice_ij( i, j ) 1895 1895 IF ( snow .AND. graupel ) THEN 1896 CALL deposition_snow_ij( i, j )1897 CALL deposition_graupel_ij( i, j )1898 CALL selfcollection_ice_ij( i, j )1899 CALL selfcollection_snow_ij( i, j )1900 CALL selfcollection_graupel_ij( i, j )1901 CALL collection_graupel_ice_ij( i, j )1902 CALL collection_snow_ice_ij( i, j )1903 CALL collection_graupel_snow_ij( i, j )1904 CALL riming_graupel_cloud_ij( i, j )1905 CALL riming_graupel_rain_ij( i, j )1906 CALL riming_ice_cloud_ij( i, j )1907 CALL riming_ice_rain_ij( i, j )1908 CALL riming_snow_cloud_ij( i, j )1909 CALL riming_snow_rain_ij( i, j )1910 CALL heterogeneous_freezing_rain_ij( i, j )1911 CALL melting_ice_ij( i, j )1912 CALL melting_snow_ij( i, j )1913 CALL melting_graupel_ij( i, j )1914 CALL evaporation_graupel_ij( i, j )1915 CALL evaporation_snow_ij( i, j )1896 CALL deposition_snow_ij( i, j ) 1897 CALL deposition_graupel_ij( i, j ) 1898 CALL selfcollection_ice_ij( i, j ) 1899 CALL selfcollection_snow_ij( i, j ) 1900 CALL selfcollection_graupel_ij( i, j ) 1901 CALL collection_graupel_ice_ij( i, j ) 1902 CALL collection_snow_ice_ij( i, j ) 1903 CALL collection_graupel_snow_ij( i, j ) 1904 CALL riming_graupel_cloud_ij( i, j ) 1905 CALL riming_graupel_rain_ij( i, j ) 1906 CALL riming_ice_cloud_ij( i, j ) 1907 CALL riming_ice_rain_ij( i, j ) 1908 CALL riming_snow_cloud_ij( i, j ) 1909 CALL riming_snow_rain_ij( i, j ) 1910 CALL heterogeneous_freezing_rain_ij( i, j ) 1911 CALL melting_ice_ij( i, j ) 1912 CALL melting_snow_ij( i, j ) 1913 CALL melting_graupel_ij( i, j ) 1914 CALL evaporation_graupel_ij( i, j ) 1915 CALL evaporation_snow_ij( i, j ) 1916 1916 ENDIF 1917 1917 ENDIF 1918 CALL autoconversion_ij( i, j )1919 CALL accretion_ij( i, j )1920 CALL selfcollection_breakup_rain_ij( i, j )1921 CALL evaporation_rain_ij( i, j )1922 CALL sedimentation_rain_ij( i, j )1923 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i ,j )1918 CALL autoconversion_ij( i, j ) 1919 CALL accretion_ij( i, j ) 1920 CALL selfcollection_breakup_rain_ij( i, j ) 1921 CALL evaporation_rain_ij( i, j ) 1922 CALL sedimentation_rain_ij( i, j ) 1923 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i ,j ) 1924 1924 IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics ) THEN 1925 IF ( ice_crystal_sedimentation ) CALL sedimentation_ice_ij( i, j )1926 IF ( snow_sedimentation .AND. snow ) CALL sedimentation_snow_ij( i, j )1927 IF ( graupel_sedimentation .AND. graupel ) CALL sedimentation_graupel_ij( i, j )1925 IF ( ice_crystal_sedimentation ) CALL sedimentation_ice_ij( i, j ) 1926 IF ( snow_sedimentation .AND. snow ) CALL sedimentation_snow_ij( i, j ) 1927 IF ( graupel_sedimentation .AND. graupel ) CALL sedimentation_graupel_ij( i, j ) 1928 1928 ENDIF 1929 1929 ! … … 1931 1931 !-- allocated. 1932 1932 ELSEIF ( microphysics_morrison_no_rain .AND. .NOT. microphysics_seifert ) THEN 1933 CALL activation_cloud_ij( i, j )1934 CALL condensation_cloud_ij( i, j )1935 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j )1933 CALL activation_cloud_ij( i, j ) 1934 CALL condensation_cloud_ij( i, j ) 1935 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j ) 1936 1936 ! 1937 1937 !-- Here the full morrison scheme is used and all processes of Seifert and Beheng are included 1938 1938 ELSEIF ( microphysics_morrison .AND. microphysics_seifert ) THEN 1939 CALL adjust_cloud_ij( i, j )1940 CALL activation_cloud_ij( i, j )1941 CALL condensation_cloud_ij( i, j )1939 CALL adjust_cloud_ij( i, j ) 1940 CALL activation_cloud_ij( i, j ) 1941 CALL condensation_cloud_ij( i, j ) 1942 1942 IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics ) THEN 1943 CALL adjust_ice_ij( i, j )1944 CALL nucleation_ice_ij( i, j )1945 CALL homogeneous_freezing_cloud_ij( i, j )1946 CALL deposition_ice_ij( i, j )1943 CALL adjust_ice_ij( i, j ) 1944 CALL nucleation_ice_ij( i, j ) 1945 CALL homogeneous_freezing_cloud_ij( i, j ) 1946 CALL deposition_ice_ij( i, j ) 1947 1947 IF ( snow .AND. graupel ) THEN 1948 CALL deposition_snow_ij ( i,j )1949 CALL deposition_graupel_ij( i, j )1950 CALL selfcollection_ice_ij( i, j )1951 CALL selfcollection_snow_ij( i, j )1952 CALL selfcollection_graupel_ij( i, j )1953 CALL collection_graupel_ice_ij( i, j )1954 CALL collection_snow_ice_ij( i, j )1955 CALL collection_graupel_snow_ij( i, j )1956 CALL riming_graupel_cloud_ij( i, j )1957 CALL riming_graupel_rain_ij( i, j )1958 CALL riming_ice_cloud_ij( i, j )1959 CALL riming_ice_rain_ij( i, j )1960 CALL riming_snow_cloud_ij( i, j )1961 CALL riming_snow_rain_ij( i, j )1962 CALL heterogeneous_freezing_rain_ij( i, j )1963 CALL melting_ice_ij( i, j )1964 CALL melting_snow_ij( i, j )1965 CALL melting_graupel_ij( i, j )1966 CALL evaporation_graupel_ij( i, j )1967 CALL evaporation_snow_ij( i, j )1948 CALL deposition_snow_ij( i, j ) 1949 CALL deposition_graupel_ij( i, j ) 1950 CALL selfcollection_ice_ij( i, j ) 1951 CALL selfcollection_snow_ij( i, j ) 1952 CALL selfcollection_graupel_ij( i, j ) 1953 CALL collection_graupel_ice_ij( i, j ) 1954 CALL collection_snow_ice_ij( i, j ) 1955 CALL collection_graupel_snow_ij( i, j ) 1956 CALL riming_graupel_cloud_ij( i, j ) 1957 CALL riming_graupel_rain_ij( i, j ) 1958 CALL riming_ice_cloud_ij( i, j ) 1959 CALL riming_ice_rain_ij( i, j ) 1960 CALL riming_snow_cloud_ij( i, j ) 1961 CALL riming_snow_rain_ij( i, j ) 1962 CALL heterogeneous_freezing_rain_ij( i, j ) 1963 CALL melting_ice_ij( i, j ) 1964 CALL melting_snow_ij( i, j ) 1965 CALL melting_graupel_ij( i, j ) 1966 CALL evaporation_graupel_ij( i, j ) 1967 CALL evaporation_snow_ij( i, j ) 1968 1968 ENDIF 1969 1969 ENDIF 1970 CALL autoconversion_ij( i, j )1971 CALL accretion_ij( i, j )1972 CALL selfcollection_breakup_rain_ij( i, j )1973 CALL evaporation_rain_ij( i, j )1974 CALL sedimentation_rain_ij( i, j )1975 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j )1970 CALL autoconversion_ij( i, j ) 1971 CALL accretion_ij( i, j ) 1972 CALL selfcollection_breakup_rain_ij( i, j ) 1973 CALL evaporation_rain_ij( i, j ) 1974 CALL sedimentation_rain_ij( i, j ) 1975 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i, j ) 1976 1976 IF ( microphysics_ice_phase .AND. simulated_time > start_ice_microphysics ) THEN 1977 IF ( ice_crystal_sedimentation ) CALL sedimentation_ice_ij( i, j )1978 IF ( snow_sedimentation .AND. snow ) CALL sedimentation_snow_ij( i, j )1979 IF ( graupel_sedimentation .AND. graupel ) CALL sedimentation_graupel_ij( i, j )1980 ENDIF 1981 1982 ENDIF 1983 1984 CALL calc_precipitation_amount_ij( i, j )1977 IF ( ice_crystal_sedimentation ) CALL sedimentation_ice_ij( i, j ) 1978 IF ( snow_sedimentation .AND. snow ) CALL sedimentation_snow_ij( i, j ) 1979 IF ( graupel_sedimentation .AND. graupel ) CALL sedimentation_graupel_ij( i, j ) 1980 ENDIF 1981 1982 ENDIF 1983 1984 CALL calc_precipitation_amount_ij( i, j ) 1985 1985 1986 1986 ENDIF … … 7452 7452 ! Calculates melting of snow after SB2006 7453 7453 !--------------------------------------------------------------------------------------------------! 7454 SUBROUTINE melting_snow_ij ( i,j )7454 SUBROUTINE melting_snow_ij( i, j ) 7455 7455 7456 7456 INTEGER(iwp) :: i !< loop index … … 7612 7612 ! Calculates evaporation of melting graupel particles after SB2006 7613 7613 !--------------------------------------------------------------------------------------------------! 7614 SUBROUTINE evaporation_graupel_ij ( i,j )7614 SUBROUTINE evaporation_graupel_ij( i, j ) 7615 7615 7616 7616 INTEGER(iwp) :: i !< loop index … … 7739 7739 ! Calculates evaporation of melting snow particles after SB2006 7740 7740 !--------------------------------------------------------------------------------------------------! 7741 SUBROUTINE evaporation_snow_ij ( i,j )7741 SUBROUTINE evaporation_snow_ij( i, j ) 7742 7742 7743 7743 INTEGER(iwp) :: i !< loop index … … 7863 7863 !> Heteorogenous freezing of rain droplets (Seifert and Beheng, 2006). 7864 7864 !--------------------------------------------------------------------------------------------------! 7865 SUBROUTINE heterogeneous_freezing_rain_ij( i, j )7865 SUBROUTINE heterogeneous_freezing_rain_ij( i, j ) 7866 7866 7867 7867 INTEGER(iwp) :: i !< loop index … … 8009 8009 !> Homogeneous freezing of cloud droplets (Seifert and Beheng, 2006). 8010 8010 !--------------------------------------------------------------------------------------------------! 8011 SUBROUTINE homogeneous_freezing_cloud_ij( i, j )8011 SUBROUTINE homogeneous_freezing_cloud_ij( i, j ) 8012 8012 8013 8013 INTEGER(iwp) :: i !< loop index
Note: See TracChangeset
for help on using the changeset viewer.