Changeset 1562 for palm/trunk
- Timestamp:
- Mar 9, 2015 1:10:37 PM (10 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/advec_ws.f90 (modified) (31 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_ws.f90
r1561 r1562 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! Bugfix: REAL constants provided with KIND-attribute in call of 23 ! intrinsic function like MAX and MIN 23 24 ! 24 25 ! Former revisions: … … 757 758 !-- Calculate ratio of upwind gradients. Note, Min/Max is just to 758 759 !-- avoid if statements. 759 rl = ( MAX( 0.0 , u(k,j,i) - u_gtrans ) *&760 rl = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) * & 760 761 ABS( ( sk(k,j,i-1) - sk(k,j,i-2) ) / & 761 762 ( sk(k,j,i) - sk(k,j,i-1) + 1E-20_wp ) & 762 763 ) + & 763 MIN( 0.0 , u(k,j,i) - u_gtrans ) *&764 MIN( 0.0_wp, u(k,j,i) - u_gtrans ) * & 764 765 ABS( ( sk(k,j,i) - sk(k,j,i+1) ) / & 765 766 ( sk(k,j,i-1) - sk(k,j,i) + 1E-20_wp ) & … … 767 768 ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp ) 768 769 769 rr = ( MAX( 0.0 , u(k,j,i+1) - u_gtrans ) *&770 rr = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 770 771 ABS( ( sk(k,j,i) - sk(k,j,i-1) ) / & 771 772 ( sk(k,j,i+1) - sk(k,j,i) + 1E-20_wp ) & 772 773 ) + & 773 MIN( 0.0 , u(k,j,i+1) - u_gtrans ) *&774 MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 774 775 ABS( ( sk(k,j,i+1) - sk(k,j,i+2) ) / & 775 776 ( sk(k,j,i) - sk(k,j,i+1) + 1E-20_wp ) & … … 777 778 ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp ) 778 779 779 rs = ( MAX( 0.0 , v(k,j,i) - v_gtrans ) *&780 rs = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) * & 780 781 ABS( ( sk(k,j-1,i) - sk(k,j-2,i) ) / & 781 782 ( sk(k,j,i) - sk(k,j-1,i) + 1E-20_wp ) & 782 783 ) + & 783 MIN( 0.0 , v(k,j,i) - v_gtrans ) *&784 MIN( 0.0_wp, v(k,j,i) - v_gtrans ) * & 784 785 ABS( ( sk(k,j,i) - sk(k,j+1,i) ) / & 785 786 ( sk(k,j-1,i) - sk(k,j,i) + 1E-20_wp ) & … … 787 788 ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp ) 788 789 789 rn = ( MAX( 0.0 , v(k,j+1,i) - v_gtrans ) *&790 rn = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 790 791 ABS( ( sk(k,j,i) - sk(k,j-1,i) ) / & 791 792 ( sk(k,j+1,i) - sk(k,j,i) + 1E-20_wp ) & 792 793 ) + & 793 MIN( 0.0 , v(k,j+1,i) - v_gtrans ) *&794 MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 794 795 ABS( ( sk(k,j+1,i) - sk(k,j+2,i) ) / & 795 796 ( sk(k,j,i) - sk(k,j+1,i) + 1E-20_wp ) & … … 803 804 k_mmm = k - 3 * ibit8 804 805 805 rd = ( MAX( 0.0 , w(k-1,j,i) ) *&806 rd = ( MAX( 0.0_wp, w(k-1,j,i) ) * & 806 807 ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i) ) / & 807 808 ( sk(k-1,j,i) - sk(k_mm,j,i) + 1E-20_wp ) & 808 809 ) + & 809 MIN( 0.0 , w(k-1,j,i) ) *&810 MIN( 0.0_wp, w(k-1,j,i) ) * & 810 811 ABS( ( sk(k-1,j,i) - sk(k,j,i) ) / & 811 812 ( sk(k_mm,j,i) - sk(k-1,j,i) + 1E-20_wp ) & … … 813 814 ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp ) 814 815 815 rt = ( MAX( 0.0 , w(k,j,i) ) *&816 rt = ( MAX( 0.0_wp, w(k,j,i) ) * & 816 817 ABS( ( sk(k,j,i) - sk(k-1,j,i) ) / & 817 818 ( sk(k+1,j,i) - sk(k,j,i) + 1E-20_wp ) & 818 819 ) + & 819 MIN( 0.0 , w(k,j,i) ) *&820 MIN( 0.0_wp, w(k,j,i) ) * & 820 821 ABS( ( sk(k+1,j,i) - sk(k_pp,j,i) ) / & 821 822 ( sk(k,j,i) - sk(k+1,j,i) + 1E-20_wp ) & … … 986 987 !-- Calculate ratio of upwind gradients. Note, Min/Max is just to 987 988 !-- avoid if statements. 988 rl = ( MAX( 0.0 , u(k,j,i) - u_gtrans ) *&989 rl = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) * & 989 990 ABS( ( sk(k,j,i-1) - sk(k,j,i-2) ) / & 990 991 ( sk(k,j,i) - sk(k,j,i-1) + 1E-20_wp ) & 991 992 ) + & 992 MIN( 0.0 , u(k,j,i) - u_gtrans ) *&993 MIN( 0.0_wp, u(k,j,i) - u_gtrans ) * & 993 994 ABS( ( sk(k,j,i) - sk(k,j,i+1) ) / & 994 995 ( sk(k,j,i-1) - sk(k,j,i) + 1E-20_wp ) & … … 996 997 ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp ) 997 998 998 rr = ( MAX( 0.0 , u(k,j,i+1) - u_gtrans ) *&999 rr = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 999 1000 ABS( ( sk(k,j,i) - sk(k,j,i-1) ) / & 1000 1001 ( sk(k,j,i+1) - sk(k,j,i) + 1E-20_wp ) & 1001 1002 ) + & 1002 MIN( 0.0 , u(k,j,i+1) - u_gtrans ) *&1003 MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 1003 1004 ABS( ( sk(k,j,i+1) - sk(k,j,i+2) ) / & 1004 1005 ( sk(k,j,i) - sk(k,j,i+1) + 1E-20_wp ) & … … 1006 1007 ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp ) 1007 1008 1008 rs = ( MAX( 0.0 , v(k,j,i) - v_gtrans ) *&1009 rs = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) * & 1009 1010 ABS( ( sk(k,j-1,i) - sk(k,j-2,i) ) / & 1010 1011 ( sk(k,j,i) - sk(k,j-1,i) + 1E-20_wp ) & 1011 1012 ) + & 1012 MIN( 0.0 , v(k,j,i) - v_gtrans ) *&1013 MIN( 0.0_wp, v(k,j,i) - v_gtrans ) * & 1013 1014 ABS( ( sk(k,j,i) - sk(k,j+1,i) ) / & 1014 1015 ( sk(k,j-1,i) - sk(k,j,i) + 1E-20_wp ) & … … 1016 1017 ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp ) 1017 1018 1018 rn = ( MAX( 0.0 , v(k,j+1,i) - v_gtrans ) *&1019 rn = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 1019 1020 ABS( ( sk(k,j,i) - sk(k,j-1,i) ) / & 1020 1021 ( sk(k,j+1,i) - sk(k,j,i) + 1E-20_wp ) & 1021 1022 ) + & 1022 MIN( 0.0 , v(k,j+1,i) - v_gtrans ) *&1023 MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 1023 1024 ABS( ( sk(k,j+1,i) - sk(k,j+2,i) ) / & 1024 1025 ( sk(k,j,i) - sk(k,j+1,i) + 1E-20_wp ) & … … 1032 1033 k_mmm = k - 3 * ibit8 1033 1034 1034 rd = ( MAX( 0.0 , w(k-1,j,i) ) *&1035 rd = ( MAX( 0.0_wp, w(k-1,j,i) ) * & 1035 1036 ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i) ) / & 1036 1037 ( sk(k-1,j,i) - sk(k_mm,j,i) + 1E-20_wp ) & 1037 1038 ) + & 1038 MIN( 0.0 , w(k-1,j,i) ) *&1039 MIN( 0.0_wp, w(k-1,j,i) ) * & 1039 1040 ABS( ( sk(k-1,j,i) - sk(k,j,i) ) / & 1040 1041 ( sk(k_mm,j,i) - sk(k-1,j,i) + 1E-20_wp ) & … … 1042 1043 ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp ) 1043 1044 1044 rt = ( MAX( 0.0 , w(k,j,i) ) *&1045 rt = ( MAX( 0.0_wp, w(k,j,i) ) * & 1045 1046 ABS( ( sk(k,j,i) - sk(k-1,j,i) ) / & 1046 1047 ( sk(k+1,j,i) - sk(k,j,i) + 1E-20_wp ) & 1047 1048 ) + & 1048 MIN( 0.0 , w(k,j,i) ) *&1049 MIN( 0.0_wp, w(k,j,i) ) * & 1049 1050 ABS( ( sk(k+1,j,i) - sk(k_pp,j,i) ) / & 1050 1051 ( sk(k,j,i) - sk(k+1,j,i) + 1E-20_wp ) & … … 2893 2894 !-- Calculate ratio of upwind gradients. Note, Min/Max is just 2894 2895 !-- to avoid if statements. 2895 rl = ( MAX( 0.0 , u(k,j,i) - u_gtrans ) *&2896 rl = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) * & 2896 2897 ABS( ( sk(k,j,i-1) - sk(k,j,i-2) ) /& 2897 2898 ( sk(k,j,i) - sk(k,j,i-1) + 1E-20_wp ) & 2898 2899 ) + & 2899 MIN( 0.0 , u(k,j,i) - u_gtrans ) *&2900 MIN( 0.0_wp, u(k,j,i) - u_gtrans ) * & 2900 2901 ABS( ( sk(k,j,i) - sk(k,j,i+1) ) /& 2901 2902 ( sk(k,j,i-1) - sk(k,j,i) + 1E-20_wp ) & … … 2903 2904 ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp ) 2904 2905 2905 rr = ( MAX( 0.0 , u(k,j,i+1) - u_gtrans ) *&2906 rr = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 2906 2907 ABS( ( sk(k,j,i) - sk(k,j,i-1) ) /& 2907 2908 ( sk(k,j,i+1) - sk(k,j,i) + 1E-20_wp ) & 2908 2909 ) + & 2909 MIN( 0.0 , u(k,j,i+1) - u_gtrans ) *&2910 MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 2910 2911 ABS( ( sk(k,j,i+1) - sk(k,j,i+2) ) /& 2911 2912 ( sk(k,j,i) - sk(k,j,i+1) + 1E-20_wp ) & … … 2913 2914 ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp ) 2914 2915 2915 rs = ( MAX( 0.0 , v(k,j,i) - v_gtrans ) *&2916 rs = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) * & 2916 2917 ABS( ( sk(k,j-1,i) - sk(k,j-2,i) ) /& 2917 2918 ( sk(k,j,i) - sk(k,j-1,i) + 1E-20_wp ) & 2918 2919 ) + & 2919 MIN( 0.0 , v(k,j,i) - v_gtrans ) *&2920 MIN( 0.0_wp, v(k,j,i) - v_gtrans ) * & 2920 2921 ABS( ( sk(k,j,i) - sk(k,j+1,i) ) /& 2921 2922 ( sk(k,j-1,i) - sk(k,j,i) + 1E-20_wp ) & … … 2923 2924 ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp ) 2924 2925 2925 rn = ( MAX( 0.0 , v(k,j+1,i) - v_gtrans ) *&2926 rn = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 2926 2927 ABS( ( sk(k,j,i) - sk(k,j-1,i) ) /& 2927 2928 ( sk(k,j+1,i) - sk(k,j,i) + 1E-20_wp ) & 2928 2929 ) + & 2929 MIN( 0.0 , v(k,j+1,i) - v_gtrans ) *&2930 MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 2930 2931 ABS( ( sk(k,j+1,i) - sk(k,j+2,i) ) /& 2931 2932 ( sk(k,j,i) - sk(k,j+1,i) + 1E-20_wp ) & … … 2939 2940 k_mmm = k - 3 * ibit8 2940 2941 2941 rd = ( MAX( 0.0 , w(k-1,j,i) ) *&2942 rd = ( MAX( 0.0_wp, w(k-1,j,i) ) * & 2942 2943 ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i) ) / & 2943 2944 ( sk(k-1,j,i) - sk(k_mm,j,i) + 1E-20_wp ) & 2944 2945 ) + & 2945 MIN( 0.0 , w(k-1,j,i) ) *&2946 MIN( 0.0_wp, w(k-1,j,i) ) * & 2946 2947 ABS( ( sk(k-1,j,i) - sk(k,j,i) ) / & 2947 2948 ( sk(k_mm,j,i) - sk(k-1,j,i) + 1E-20_wp ) & … … 2949 2950 ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp ) 2950 2951 2951 rt = ( MAX( 0.0 , w(k,j,i) ) *&2952 rt = ( MAX( 0.0_wp, w(k,j,i) ) * & 2952 2953 ABS( ( sk(k,j,i) - sk(k-1,j,i) ) / & 2953 2954 ( sk(k+1,j,i) - sk(k,j,i) + 1E-20_wp ) & 2954 2955 ) + & 2955 MIN( 0.0 , w(k,j,i) ) *&2956 MIN( 0.0_wp, w(k,j,i) ) * & 2956 2957 ABS( ( sk(k+1,j,i) - sk(k_pp,j,i) ) / & 2957 2958 ( sk(k,j,i) - sk(k+1,j,i) + 1E-20_wp ) & … … 3118 3119 !-- Calculate ratio of upwind gradients. Note, Min/Max is just 3119 3120 !-- to avoid if statements. 3120 rl = ( MAX( 0.0 , u(k,j,i) - u_gtrans ) *&3121 rl = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) * & 3121 3122 ABS( ( sk(k,j,i-1) - sk(k,j,i-2) ) /& 3122 3123 ( sk(k,j,i) - sk(k,j,i-1) + 1E-20_wp ) & 3123 3124 ) + & 3124 MIN( 0.0 , u(k,j,i) - u_gtrans ) *&3125 MIN( 0.0_wp, u(k,j,i) - u_gtrans ) * & 3125 3126 ABS( ( sk(k,j,i) - sk(k,j,i+1) ) /& 3126 3127 ( sk(k,j,i-1) - sk(k,j,i) + 1E-20_wp ) & … … 3128 3129 ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp ) 3129 3130 3130 rr = ( MAX( 0.0 , u(k,j,i+1) - u_gtrans ) *&3131 rr = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 3131 3132 ABS( ( sk(k,j,i) - sk(k,j,i-1) ) /& 3132 3133 ( sk(k,j,i+1) - sk(k,j,i) + 1E-20_wp ) & 3133 3134 ) + & 3134 MIN( 0.0 , u(k,j,i+1) - u_gtrans ) *&3135 MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 3135 3136 ABS( ( sk(k,j,i+1) - sk(k,j,i+2) ) /& 3136 3137 ( sk(k,j,i) - sk(k,j,i+1) + 1E-20_wp ) & … … 3138 3139 ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp ) 3139 3140 3140 rs = ( MAX( 0.0 , v(k,j,i) - v_gtrans ) *&3141 rs = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) * & 3141 3142 ABS( ( sk(k,j-1,i) - sk(k,j-2,i) ) /& 3142 3143 ( sk(k,j,i) - sk(k,j-1,i) + 1E-20_wp ) & 3143 3144 ) + & 3144 MIN( 0.0 , v(k,j,i) - v_gtrans ) *&3145 MIN( 0.0_wp, v(k,j,i) - v_gtrans ) * & 3145 3146 ABS( ( sk(k,j,i) - sk(k,j+1,i) ) /& 3146 3147 ( sk(k,j-1,i) - sk(k,j,i) + 1E-20_wp ) & … … 3148 3149 ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp ) 3149 3150 3150 rn = ( MAX( 0.0 , v(k,j+1,i) - v_gtrans ) *&3151 rn = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 3151 3152 ABS( ( sk(k,j,i) - sk(k,j-1,i) ) /& 3152 3153 ( sk(k,j+1,i) - sk(k,j,i) + 1E-20_wp ) & 3153 3154 ) + & 3154 MIN( 0.0 , v(k,j+1,i) - v_gtrans ) *&3155 MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 3155 3156 ABS( ( sk(k,j+1,i) - sk(k,j+2,i) ) /& 3156 3157 ( sk(k,j,i) - sk(k,j+1,i) + 1E-20_wp ) & … … 3164 3165 k_mmm = k - 3 * ibit8 3165 3166 3166 rd = ( MAX( 0.0 , w(k-1,j,i) ) *&3167 rd = ( MAX( 0.0_wp, w(k-1,j,i) ) * & 3167 3168 ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i) ) / & 3168 3169 ( sk(k-1,j,i) - sk(k_mm,j,i) + 1E-20_wp ) & 3169 3170 ) + & 3170 MIN( 0.0 , w(k-1,j,i) ) *&3171 MIN( 0.0_wp, w(k-1,j,i) ) * & 3171 3172 ABS( ( sk(k-1,j,i) - sk(k,j,i) ) / & 3172 3173 ( sk(k_mm,j,i) - sk(k-1,j,i) + 1E-20_wp ) & … … 3174 3175 ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp ) 3175 3176 3176 rt = ( MAX( 0.0 , w(k,j,i) ) *&3177 rt = ( MAX( 0.0_wp, w(k,j,i) ) * & 3177 3178 ABS( ( sk(k,j,i) - sk(k-1,j,i) ) / & 3178 3179 ( sk(k+1,j,i) - sk(k,j,i) + 1E-20_wp ) & 3179 3180 ) + & 3180 MIN( 0.0 , w(k,j,i) ) *&3181 MIN( 0.0_wp, w(k,j,i) ) * & 3181 3182 ABS( ( sk(k+1,j,i) - sk(k_pp,j,i) ) / & 3182 3183 ( sk(k,j,i) - sk(k+1,j,i) + 1E-20_wp ) & … … 3631 3632 !-- Calculate ratio of upwind gradients. Note, Min/Max is just 3632 3633 !-- to avoid if statements. 3633 rl = ( MAX( 0.0 , u(k,j,i) - u_gtrans ) *&3634 rl = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) * & 3634 3635 ABS( ( sk(k,j,i-1) - sk(k,j,i-2) ) /& 3635 3636 ( sk(k,j,i) - sk(k,j,i-1) + 1E-20_wp ) & 3636 3637 ) + & 3637 MIN( 0.0 , u(k,j,i) - u_gtrans ) *&3638 MIN( 0.0_wp, u(k,j,i) - u_gtrans ) * & 3638 3639 ABS( ( sk(k,j,i) - sk(k,j,i+1) ) /& 3639 3640 ( sk(k,j,i-1) - sk(k,j,i) + 1E-20_wp ) & … … 3641 3642 ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp ) 3642 3643 3643 rr = ( MAX( 0.0 , u(k,j,i+1) - u_gtrans ) *&3644 rr = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 3644 3645 ABS( ( sk(k,j,i) - sk(k,j,i-1) ) /& 3645 3646 ( sk(k,j,i+1) - sk(k,j,i) + 1E-20_wp ) & 3646 3647 ) + & 3647 MIN( 0.0 , u(k,j,i+1) - u_gtrans ) *&3648 MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) * & 3648 3649 ABS( ( sk(k,j,i+1) - sk(k,j,i+2) ) /& 3649 3650 ( sk(k,j,i) - sk(k,j,i+1) + 1E-20_wp ) & … … 3651 3652 ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp ) 3652 3653 3653 rs = ( MAX( 0.0 , v(k,j,i) - v_gtrans ) *&3654 rs = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) * & 3654 3655 ABS( ( sk(k,j-1,i) - sk(k,j-2,i) ) /& 3655 3656 ( sk(k,j,i) - sk(k,j-1,i) + 1E-20_wp ) & 3656 3657 ) + & 3657 MIN( 0.0 , v(k,j,i) - v_gtrans ) *&3658 MIN( 0.0_wp, v(k,j,i) - v_gtrans ) * & 3658 3659 ABS( ( sk(k,j,i) - sk(k,j+1,i) ) /& 3659 3660 ( sk(k,j-1,i) - sk(k,j,i) + 1E-20_wp ) & … … 3661 3662 ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp ) 3662 3663 3663 rn = ( MAX( 0.0 , v(k,j+1,i) - v_gtrans ) *&3664 rn = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 3664 3665 ABS( ( sk(k,j,i) - sk(k,j-1,i) ) /& 3665 3666 ( sk(k,j+1,i) - sk(k,j,i) + 1E-20_wp ) & 3666 3667 ) + & 3667 MIN( 0.0 , v(k,j+1,i) - v_gtrans ) *&3668 MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) * & 3668 3669 ABS( ( sk(k,j+1,i) - sk(k,j+2,i) ) /& 3669 3670 ( sk(k,j,i) - sk(k,j+1,i) + 1E-20_wp ) & … … 3677 3678 k_mmm = k - 3 * ibit8 3678 3679 3679 rd = ( MAX( 0.0 , w(k-1,j,i) ) *&3680 rd = ( MAX( 0.0_wp, w(k-1,j,i) ) * & 3680 3681 ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i) ) / & 3681 3682 ( sk(k-1,j,i) - sk(k_mm,j,i) + 1E-20_wp ) & 3682 3683 ) + & 3683 MIN( 0.0 , w(k-1,j,i) ) *&3684 MIN( 0.0_wp, w(k-1,j,i) ) * & 3684 3685 ABS( ( sk(k-1,j,i) - sk(k,j,i) ) / & 3685 3686 ( sk(k_mm,j,i) - sk(k-1,j,i) + 1E-20_wp ) & … … 3687 3688 ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp ) 3688 3689 3689 rt = ( MAX( 0.0 , w(k,j,i) ) *&3690 rt = ( MAX( 0.0_wp, w(k,j,i) ) * & 3690 3691 ABS( ( sk(k,j,i) - sk(k-1,j,i) ) / & 3691 3692 ( sk(k+1,j,i) - sk(k,j,i) + 1E-20_wp ) & 3692 3693 ) + & 3693 MIN( 0.0 , w(k,j,i) ) *&3694 MIN( 0.0_wp, w(k,j,i) ) * & 3694 3695 ABS( ( sk(k+1,j,i) - sk(k_pp,j,i) ) / & 3695 3696 ( sk(k,j,i) - sk(k+1,j,i) + 1E-20_wp ) &
Note: See TracChangeset
for help on using the changeset viewer.
