Changeset 3757
- Timestamp:
- Jul 15, 2020, 10:14:33 PM (4 years ago)
- Location:
- LMDZ6/branches/Optimisation_LMDZ/libf/phylmd
- Files:
-
- 1 edited
- 1 copied
-
cv3a_driver.f90 (copied) (copied from LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cva_driver.F90) (2 diffs)
-
cva_driver.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_driver.f90
r3755 r3757 1 2 ! $Id$ 3 4 SUBROUTINE cv a_driver(len, nd, ndp1, ntra, nloc, k_upper, &1 module cv3a_driver_mod 2 contains 3 4 SUBROUTINE cv3a_driver(len, nd, ndp1, ntra, nloc, k_upper, & 5 5 iflag_con, iflag_mix, iflag_ice_thermo, iflag_clos, ok_conserv_q, & 6 6 !! delt, t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, & ! jyg … … 1381 1381 1382 1382 RETURN 1383 END SUBROUTINE cva_driver 1383 END SUBROUTINE cv3a_driver 1384 1385 END MODULE -
LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cva_driver.F90
r3709 r3757 1 2 ! $Id$3 4 1 SUBROUTINE cva_driver(len, nd, ndp1, ntra, nloc, k_upper, & 2 iflag_con, iflag_mix, iflag_ice_thermo, iflag_clos, ok_conserv_q, & 3 delt, comp_threshold, & 4 t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, & 5 u1, v1, tra1, & 6 p1, ph1, & 7 Ale1, Alp1, omega1, & 8 sig1feed1, sig2feed1, wght1, & 9 iflag1, ft1, fq1, fu1, fv1, ftra1, & 10 precip1, kbas1, ktop1, & 11 cbmf1, plcl1, plfc1, wbeff1, & 12 sig1, w01, & 13 ptop21, sigd1, & 14 ma1, mip1, Vprecip1, Vprecipi1, upwd1, dnwd1, dnwd01, & 15 qcondc1, wd1, & 16 cape1, cin1, tvp1, & 17 ftd1, fqd1, & 18 Plim11, Plim21, asupmax1, supmax01, asupmaxmin1, & 19 lalim_conv1, & 20 da1, phi1, mp1, phi21, d1a1, dam1, sigij1, wghti1, & 21 qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & 22 wdtrainA1, wdtrainS1, wdtrainM1, qtc1, sigt1, tau_cld_cv, & 23 coefw_cld_cv, & 24 epmax_diag1) 25 use cv3a_driver_mod 26 IMPLICIT NONE 27 INTEGER, INTENT (IN) :: len 28 INTEGER, INTENT (IN) :: nd 29 INTEGER, INTENT (IN) :: ndp1 30 INTEGER, INTENT (IN) :: ntra 31 INTEGER, INTENT(IN) :: nloc ! (nloc=len) pour l'instant 32 INTEGER, INTENT (IN) :: k_upper 33 INTEGER, INTENT (IN) :: iflag_con 34 INTEGER, INTENT (IN) :: iflag_mix 35 INTEGER, INTENT (IN) :: iflag_ice_thermo 36 INTEGER, INTENT (IN) :: iflag_clos 37 LOGICAL, INTENT (IN) :: ok_conserv_q 38 REAL, INTENT (IN) :: tau_cld_cv 39 REAL, INTENT (IN) :: coefw_cld_cv 40 REAL, INTENT (IN) :: delt 41 REAL, INTENT (IN) :: comp_threshold 42 REAL, DIMENSION (len, nd), INTENT (IN) :: t1 43 REAL, DIMENSION (len, nd), INTENT (IN) :: q1 44 REAL, DIMENSION (len, nd), INTENT (IN) :: qs1 45 REAL, DIMENSION (len, nd), INTENT (IN) :: t1_wake 46 REAL, DIMENSION (len, nd), INTENT (IN) :: q1_wake 47 REAL, DIMENSION (len, nd), INTENT (IN) :: qs1_wake 48 REAL, DIMENSION (len), INTENT (IN) :: s1_wake 49 REAL, DIMENSION (len, nd), INTENT (IN) :: u1 50 REAL, DIMENSION (len, nd), INTENT (IN) :: v1 51 REAL, DIMENSION (len, nd, ntra), INTENT (IN) :: tra1 52 REAL, DIMENSION (len, nd), INTENT (IN) :: p1 53 REAL, DIMENSION (len, ndp1), INTENT (IN) :: ph1 54 REAL, DIMENSION (len), INTENT (IN) :: Ale1 55 REAL, DIMENSION (len), INTENT (IN) :: Alp1 56 REAL, DIMENSION (len, nd), INTENT (IN) :: omega1 57 REAL, INTENT (IN) :: sig1feed1 ! pressure at lower bound of feeding layer 58 REAL, INTENT (IN) :: sig2feed1 ! pressure at upper bound of feeding layer 59 REAL, DIMENSION (nd), INTENT (IN) :: wght1 ! weight density determining the feeding mixture 60 INTEGER, DIMENSION (len), INTENT (IN) :: lalim_conv1 61 62 ! Input/Output 63 REAL, DIMENSION (len, nd), INTENT (INOUT) :: sig1 64 REAL, DIMENSION (len, nd), INTENT (INOUT) :: w01 65 66 ! Output 67 INTEGER, DIMENSION (len), INTENT (OUT) :: iflag1 68 REAL, DIMENSION (len, nd), INTENT (OUT) :: ft1 69 REAL, DIMENSION (len, nd), INTENT (OUT) :: fq1 70 REAL, DIMENSION (len, nd), INTENT (OUT) :: fu1 71 REAL, DIMENSION (len, nd), INTENT (OUT) :: fv1 72 REAL, DIMENSION (len, nd, ntra), INTENT (OUT) :: ftra1 73 REAL, DIMENSION (len), INTENT (OUT) :: precip1 74 INTEGER, DIMENSION (len), INTENT (OUT) :: kbas1 75 INTEGER, DIMENSION (len), INTENT (OUT) :: ktop1 76 REAL, DIMENSION (len), INTENT (OUT) :: cbmf1 77 REAL, DIMENSION (len), INTENT (OUT) :: plcl1 78 REAL, DIMENSION (len), INTENT (OUT) :: plfc1 79 REAL, DIMENSION (len), INTENT (OUT) :: wbeff1 80 REAL, DIMENSION (len), INTENT (OUT) :: ptop21 81 REAL, DIMENSION (len), INTENT (OUT) :: sigd1 82 REAL, DIMENSION (len, nd), INTENT (OUT) :: ma1 ! adiab. asc. mass flux (staggered grid) 83 REAL, DIMENSION (len, nd), INTENT (OUT) :: mip1 ! mass flux shed from adiab. ascent (extensive) 84 REAL, DIMENSION (len, ndp1), INTENT (OUT) :: vprecip1 ! tot precipitation flux (staggered grid) 85 REAL, DIMENSION (len, ndp1), INTENT (OUT) :: vprecipi1 ! ice precipitation flux (staggered grid) 86 REAL, DIMENSION (len, nd), INTENT (OUT) :: upwd1 ! upwd sat. mass flux (staggered grid) 87 REAL, DIMENSION (len, nd), INTENT (OUT) :: dnwd1 ! dnwd sat. mass flux (staggered grid) 88 REAL, DIMENSION (len, nd), INTENT (OUT) :: dnwd01 ! unsat. mass flux (staggered grid) 89 REAL, DIMENSION (len, nd), INTENT (OUT) :: qcondc1 ! max cloud condensate (intensive) ! cld 90 REAL, DIMENSION (len), INTENT (OUT) :: wd1 ! gust 91 REAL, DIMENSION (len), INTENT (OUT) :: cape1 92 REAL, DIMENSION (len), INTENT (OUT) :: cin1 93 REAL, DIMENSION (len, nd), INTENT (OUT) :: tvp1 ! Virt. temp. in the adiab. ascent 94 REAL, DIMENSION (len, nd), INTENT (OUT) :: ftd1 ! Temp. tendency due to the sole unsat. drafts 95 REAL, DIMENSION (len, nd), INTENT (OUT) :: fqd1 ! Moist. tendency due to the sole unsat. drafts 96 REAL, DIMENSION (len), INTENT (OUT) :: Plim11 97 REAL, DIMENSION (len), INTENT (OUT) :: Plim21 98 REAL, DIMENSION (len, nd), INTENT (OUT) :: asupmax1 ! Highest mixing fraction of mixed updraughts 99 REAL, DIMENSION (len), INTENT (OUT) :: supmax01 100 REAL, DIMENSION (len), INTENT (OUT) :: asupmaxmin1 101 REAL, DIMENSION (len, nd), INTENT (OUT) :: qtc1 ! in cloud water content (intensive) ! cld 102 REAL, DIMENSION (len, nd), INTENT (OUT) :: sigt1 ! fract. cloud area (intensive) ! cld 103 104 ! RomP >>> 105 REAL, DIMENSION (len, nd), INTENT (OUT) :: wdtrainA1, wdtrainS1, wdtrainM1 ! precipitation sources (extensive) 106 REAL, DIMENSION (len, nd), INTENT (OUT) :: mp1 ! unsat. mass flux (staggered grid) 107 REAL, DIMENSION (len, nd), INTENT (OUT) :: da1 ! detrained mass flux of adiab. asc. air (extensive) 108 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi1 ! mass flux of envt. air in mixed draughts (extensive) 109 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: epmlmMm1 ! (extensive) 110 REAL, DIMENSION (len, nd), INTENT (OUT) :: eplaMm1 ! (extensive) 111 REAL, DIMENSION (len, nd), INTENT (OUT) :: evap1 ! evaporation rate in precip. downdraft. (intensive) 112 REAL, DIMENSION (len, nd), INTENT (OUT) :: ep1 113 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: sigij1 ! mass fraction of env. air in mixed draughts (intensive) 114 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: elij1! cond. water per unit mass of mixed draughts (intensive) 115 REAL, DIMENSION (len, nd), INTENT (OUT) :: qta1 ! total water per unit mass of the adiab. asc. (intensive) 116 REAL, DIMENSION (len, nd), INTENT (OUT) :: clw1 ! cond. water per unit mass of the adiab. asc. (intensive) 117 !JYG,RL 118 REAL, DIMENSION (len, nd), INTENT (OUT) :: wghti1 ! final weight of the feeding layers (extensive) 119 !JYG,RL 120 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi21 ! (extensive) 121 REAL, DIMENSION (len, nd), INTENT (OUT) :: d1a1 ! (extensive) 122 REAL, DIMENSION (len, nd), INTENT (OUT) :: dam1 ! (extensive) 123 ! RomP <<< 124 REAL, DIMENSION (len ), INTENT (OUT) :: epmax_diag1 125 126 if(iflag_con == 3) then 127 call cv3a_driver(len, nd, ndp1, ntra, nloc, k_upper, & 128 iflag_con, iflag_mix, iflag_ice_thermo, iflag_clos, ok_conserv_q, & 129 delt, comp_threshold, & 130 t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, & 131 u1, v1, tra1, & 132 p1, ph1, & 133 Ale1, Alp1, omega1, & 134 sig1feed1, sig2feed1, wght1, & 135 iflag1, ft1, fq1, fu1, fv1, ftra1, & 136 precip1, kbas1, ktop1, & 137 cbmf1, plcl1, plfc1, wbeff1, & 138 sig1, w01, & 139 ptop21, sigd1, & 140 ma1, mip1, Vprecip1, Vprecipi1, upwd1, dnwd1, dnwd01, & 141 qcondc1, wd1, & 142 cape1, cin1, tvp1, & 143 ftd1, fqd1, & 144 Plim11, Plim21, asupmax1, supmax01, asupmaxmin1, & 145 lalim_conv1, & 146 da1, phi1, mp1, phi21, d1a1, dam1, sigij1, wghti1, & 147 qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & 148 wdtrainA1, wdtrainS1, wdtrainM1, qtc1, sigt1, tau_cld_cv, & 149 coefw_cld_cv, & 150 epmax_diag1) 151 elseif(iflag_con == 4) then 152 call cv4a_driver(len, nd, ndp1, ntra, nloc, k_upper, & 153 iflag_con, iflag_mix, iflag_ice_thermo, iflag_clos, ok_conserv_q, & 154 delt, comp_threshold, & 155 t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, & 156 u1, v1, tra1, & 157 p1, ph1, & 158 Ale1, Alp1, omega1, & 159 sig1feed1, sig2feed1, wght1, & 160 iflag1, ft1, fq1, fu1, fv1, ftra1, & 161 precip1, kbas1, ktop1, & 162 cbmf1, plcl1, plfc1, wbeff1, & 163 sig1, w01, & 164 ptop21, sigd1, & 165 ma1, mip1, Vprecip1, Vprecipi1, upwd1, dnwd1, dnwd01, & 166 qcondc1, wd1, & 167 cape1, cin1, tvp1, & 168 ftd1, fqd1, & 169 Plim11, Plim21, asupmax1, supmax01, asupmaxmin1, & 170 lalim_conv1, & 171 da1, phi1, mp1, phi21, d1a1, dam1, sigij1, wghti1, & 172 qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & 173 wdtrainA1, wdtrainS1, wdtrainM1, qtc1, sigt1, tau_cld_cv, & 174 coefw_cld_cv, & 175 epmax_diag1) 176 else 177 call abort_physic("cva_driver", "iflag_con is not compatible", 1) 178 endif 179 180 END SUBROUTINE 181 182 SUBROUTINE cv4a_driver(len, nd, ndp1, ntra, nloc, k_upper, & 5 183 iflag_con, iflag_mix, iflag_ice_thermo, iflag_clos, ok_conserv_q, & 6 184 !! delt, t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, & ! jyg … … 1381 1559 1382 1560 RETURN 1383 END SUBROUTINE cv a_driver1561 END SUBROUTINE cv4a_driver
Note: See TracChangeset
for help on using the changeset viewer.
