- Timestamp:
- May 7, 2018 2:38:42 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r3004 r3011 25 25 ! ----------------- 26 26 ! $Id$ 27 ! removed redundant if statement 28 ! 29 ! 3004 2018-04-27 12:33:25Z Giersch 27 30 ! precipitation_rate removed 28 31 ! … … 691 694 ALLOCATE( q(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 692 695 q_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 693 tq_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 696 tq_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 697 vpt(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 694 698 #else 695 699 ALLOCATE( q_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 696 700 q_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 697 q_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 701 q_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 702 vpt_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 698 703 #endif 699 704 700 ! 701 !-- 3D-arrays needed for humidity 702 IF ( humidity ) THEN 703 #if defined( __nopointer ) 704 ALLOCATE( vpt(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 705 #else 706 ALLOCATE( vpt_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 707 #endif 708 709 IF ( cloud_physics ) THEN 705 IF ( cloud_physics ) THEN 710 706 ! 711 707 !-- Liquid water content 712 708 #if defined( __nopointer ) 713 709 ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 714 710 #else 715 711 ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 716 712 #endif 717 713 718 714 ! 719 !-- 720 715 !-- 3D-cloud water content 716 IF ( .NOT. microphysics_morrison ) THEN 721 717 #if defined( __nopointer ) 722 718 ALLOCATE( qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 723 719 #else 724 720 ALLOCATE( qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 725 721 #endif 726 727 ! 728 !-- 729 730 731 ! 732 !-- 733 734 735 736 ! 737 !-- 722 ENDIF 723 ! 724 !-- Precipitation amount and rate (only needed if output is switched) 725 ALLOCATE( precipitation_amount(nysg:nyng,nxlg:nxrg) ) 726 727 ! 728 !-- 3d-precipitation rate 729 ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 730 731 IF ( microphysics_morrison ) THEN 732 ! 733 !-- 3D-cloud drop water content, cloud drop concentration arrays 738 734 #if defined( __nopointer ) 739 ALLOCATE( nc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&740 nc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&741 qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&742 qc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&743 tnc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&744 735 ALLOCATE( nc(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 736 nc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 737 qc(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 738 qc_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 739 tnc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 740 tqc_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 745 741 #else 746 ALLOCATE( nc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&747 nc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&748 nc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&749 qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&750 qc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&751 742 ALLOCATE( nc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 743 nc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 744 nc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 745 qc_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 746 qc_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 747 qc_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 752 748 #endif 753 754 755 756 ! 757 !-- 749 ENDIF 750 751 IF ( microphysics_seifert ) THEN 752 ! 753 !-- 3D-rain water content, rain drop concentration arrays 758 754 #if defined( __nopointer ) 759 ALLOCATE( nr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&760 nr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&761 qr(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&762 qr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&763 tnr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&764 755 ALLOCATE( nr(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 756 nr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 757 qr(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 758 qr_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 759 tnr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 760 tqr_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 765 761 #else 766 ALLOCATE( nr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&767 nr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&768 nr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&769 qr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&770 qr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&771 762 ALLOCATE( nr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 763 nr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 764 nr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 765 qr_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 766 qr_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 767 qr_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 772 768 #endif 773 774 775 776 777 778 ! 779 !-- 769 ENDIF 770 771 ENDIF 772 773 IF ( cloud_droplets ) THEN 774 ! 775 !-- Liquid water content, change in liquid water content 780 776 #if defined( __nopointer ) 781 ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&782 777 ALLOCATE ( ql(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 778 ql_c(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 783 779 #else 784 ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),&785 780 ALLOCATE ( ql_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 781 ql_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 786 782 #endif 787 783 ! 788 !-- Real volume of particles (with weighting), volume of particles 789 ALLOCATE ( ql_v(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 790 ql_vp(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 791 ENDIF 792 793 ENDIF 794 795 ENDIF 796 784 !-- Real volume of particles (with weighting), volume of particles 785 ALLOCATE ( ql_v(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 786 ql_vp(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 787 ENDIF 788 789 ENDIF 797 790 798 791 IF ( passive_scalar ) THEN
Note: See TracChangeset
for help on using the changeset viewer.