- Timestamp:
- Jul 15, 2020, 10:14:36 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_driver.f90
r3762 r3763 634 634 else !compress 635 635 nloc = len 636 ktop1 = 0637 precip1 = 0638 cbmf1 = 0639 plfc1 = 0640 wbeff1 = 0641 ptop21 = 0642 sigd1 = 0643 wd1 = 0644 cape1 = Cape_noconv645 cin1 = Cin_noconv646 Plim11 = 0647 plim21 = 0648 supmax01 = 0649 asupmaxmin1 = 0650 epmax_diag1 = 0651 ft1 = 0652 fq1 = 0653 fu1 = 0654 fv1 = 0655 ma1 = 0656 mip1 = 0657 upwd1 = 0658 dnwd1 = 0659 dnwd01 = 0660 qcondc1 = 0661 ftd1 = 0662 fqd1 = 0663 asupmax1 = 0664 da1 = 0665 mp1 = 0666 d1a1 = 0667 dam1 = 0668 qta1 = 0669 evap1 = 0670 ep1 = 0671 eplaMm1 = 0672 wdtrainA1 = 0673 wdtrainS1 = 0674 wdtrainM1 = 0675 qtc1 = 0676 sigt1 = 0677 vprecip1 = 0678 vprecipi1 = 0679 phi1 = 0680 phi21 = 0681 sigij1 = 0682 elij1 = 0683 epmlmMm1 = 0684 636 coef_convective(:) = merge(1, 0, iflag1(:) == 0) 685 637 call cv3a_driver_compressed(nloc, nd, ntra, & … … 810 762 811 763 call driver_log('cv3_undilute2') 764 qta(:,:) = 0 765 ep(:,:) = 0 812 766 CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, & 813 767 tnk, qnk, gznk, hnk, t, q, qs, gz, & 814 768 p, ph, h, tv, lv, lf, pbase, buoybase, plcl, & 815 769 inb, tp, tvp, clw, hp, ep, sigp, buoy, & 816 frac_a, frac_s, qpreca, qta) 817 770 frac_a, frac_s, qpreca, qta) 771 818 772 ! epmax_cape 819 773 ! on recalcule ep et hp 820 774 call driver_log('cv3_epmax_cape') 775 epmax_diag(:) = 0 821 776 call cv3_epmax_fn_cape(nloc, ncum, nd & 822 777 , ep, hp, icb, inb, clw, nk, t, h, hnk, lv, lf, frac_s & … … 828 783 ! ------------------------------------------------------------------- 829 784 call enter_profile("cv3p_mixing") 785 elij(:,:,:) = 0 786 supmax(:,:) = 0 830 787 IF (iflag_mix >= 1) THEN 831 CALL zilch(supmax, nloc*nd)832 788 call driver_log('cv3p_mixing') 833 789 CALL cv3p_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, & … … 836 792 ment, qent, hent, uent, vent, nent, & 837 793 sigij, elij, supmax, ments, qents, traent) 838 ELSE839 CALL zilch(supmax, nloc*nd)840 794 END IF 841 795 call exit_profile("cv3p_mixing") … … 845 799 ! ------------------------------------------------------------------- 846 800 801 cape(:) = -1 802 cin(:) = -100000. 847 803 ptop2(:) = 0 848 804 coef_clos(:) = 1. 805 Plim1(:) = 0 806 plim2(:) = 0 807 supmax0(:) = 0 808 asupmaxmin(:) = 0 809 cbmf(:) = 0 810 plfc(:) = 0 811 wbeff(:) = 0 812 asupmax(:,:) = 0 813 814 ok_inhib = (iflag_mix == 2) 815 849 816 IF (iflag_clos == 0) THEN 850 817 call driver_log('cv3_closure') … … 852 819 pbase, p, ph, tv, buoy, & 853 820 sig, w0, cape, m, iflag) 854 END IF ! iflag_clos==0 855 856 ok_inhib = iflag_mix == 2 821 END IF ! iflag_clos==0 857 822 858 823 IF (iflag_clos == 1) PRINT *, ' pas d appel cv3p_closure' … … 958 923 ENDIF 959 924 925 926 960 927 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 961 928 !--- passive tracers … … 964 931 call driver_log('cv3_tracer') 965 932 call enter_profile("cv3_tracer") 933 934 sigij(:,:,:) = 0 935 vprecip(:,:) = 0 936 ! GLITCHY : vprecip is unused in cv3_tracer and sigij is intent in 966 937 CALL cv3_tracer(nloc, -1, ncum, nd, nd, & 967 938 ment, sigij, da, phi, phi2, d1a, dam, & 968 939 ep, vprecip, elij, clw, epmlmMm, eplaMm, & 969 940 icb, inb) 941 942 vprecipi = 0 970 943 call exit_profile("cv3_tracer") 971 944 END IF ! ncum>0
Note: See TracChangeset
for help on using the changeset viewer.