| Directory: | ./ |
|---|---|
| File: | phys/cv3a_uncompress.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 77 | 138 | 55.8% |
| Branches: | 15 | 222 | 6.8% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | 127337 | SUBROUTINE cv3a_uncompress(nloc, len, ncum, nd, ntra, idcum, compress, & | |
| 2 | iflag, kbas, ktop, & | ||
| 3 | 240 | precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, & | |
| 4 | ft, fq, fu, fv, ftra, & | ||
| 5 | sigd, ma, mip, vprecip, vprecipi, upwd, dnwd, dnwd0, & | ||
| 6 | qcondc, wd, cape, cin, & | ||
| 7 | tvp, & | ||
| 8 | ftd, fqd, & | ||
| 9 | plim1, plim2, asupmax, supmax0, & | ||
| 10 | asupmaxmin, & | ||
| 11 | 240 | da, phi, mp, phi2, d1a, dam, sigij, & ! RomP+AC+jyg | |
| 12 | qta, clw, elij, evap, ep, epmlmMm, eplaMm, & ! RomP+jyg | ||
| 13 | wdtrainA, wdtrainS, wdtrainM, & ! RomP | ||
| 14 | qtc, sigt, & | ||
| 15 | epmax_diag, & ! epmax_cape | ||
| 16 | iflag1, kbas1, ktop1, & | ||
| 17 | 240 | precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, ptop21, & | |
| 18 | ft1, fq1, fu1, fv1, ftra1, & | ||
| 19 | sigd1, ma1, mip1, vprecip1, vprecipi1, upwd1, dnwd1, dnwd01, & | ||
| 20 | qcondc1, wd1, cape1, cin1, & | ||
| 21 | tvp1, & | ||
| 22 | ftd1, fqd1, & | ||
| 23 | plim11, plim21, asupmax1, supmax01, & | ||
| 24 | asupmaxmin1, & | ||
| 25 | 240 | da1, phi1, mp1, phi21, d1a1, dam1, sigij1, & ! RomP+AC+jyg | |
| 26 | qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP+jyg | ||
| 27 | wdtrainA1, wdtrainS1, wdtrainM1, & ! RomP | ||
| 28 | qtc1, sigt1, & | ||
| 29 | epmax_diag1) ! epmax_cape | ||
| 30 | |||
| 31 | ! ************************************************************** | ||
| 32 | ! * | ||
| 33 | ! CV3A_UNCOMPRESS * | ||
| 34 | ! * | ||
| 35 | ! * | ||
| 36 | ! written by : Sandrine Bony-Lena , 17/05/2003, 11.22.15 * | ||
| 37 | ! modified by : Jean-Yves Grandpeix, 23/06/2003, 10.36.17 * | ||
| 38 | ! ************************************************************** | ||
| 39 | |||
| 40 | IMPLICIT NONE | ||
| 41 | |||
| 42 | include "cv3param.h" | ||
| 43 | |||
| 44 | ! inputs: | ||
| 45 | INTEGER, INTENT (IN) :: nloc, len, ncum, nd, ntra | ||
| 46 | INTEGER, DIMENSION (nloc), INTENT (IN) :: idcum(nloc) | ||
| 47 | !jyg< | ||
| 48 | LOGICAL, INTENT (IN) :: compress | ||
| 49 | !>jyg | ||
| 50 | INTEGER, DIMENSION (nloc), INTENT (IN) ::iflag, kbas, ktop | ||
| 51 | REAL, DIMENSION (nloc), INTENT (IN) :: precip, cbmf, plcl, plfc | ||
| 52 | REAL, DIMENSION (nloc), INTENT (IN) :: wbeff | ||
| 53 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: sig, w0 | ||
| 54 | REAL, DIMENSION (nloc), INTENT (IN) :: ptop2 | ||
| 55 | REAL, DIMENSION (nloc), INTENT (IN) :: epmax_diag | ||
| 56 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: ft, fq, fu, fv | ||
| 57 | REAL, DIMENSION (nloc, nd, ntra), INTENT (IN) :: ftra | ||
| 58 | REAL, DIMENSION (nloc), INTENT (IN) :: sigd | ||
| 59 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: ma, mip | ||
| 60 | REAL, DIMENSION (nloc, nd+1), INTENT (IN) :: vprecip | ||
| 61 | REAL, DIMENSION (nloc, nd+1), INTENT (IN) :: vprecipi | ||
| 62 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: upwd, dnwd, dnwd0 | ||
| 63 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: qcondc | ||
| 64 | REAL, DIMENSION (nloc), INTENT (IN) :: wd, cape, cin | ||
| 65 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: tvp | ||
| 66 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: ftd, fqd | ||
| 67 | REAL, DIMENSION (nloc), INTENT (IN) :: plim1, plim2 | ||
| 68 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: asupmax | ||
| 69 | REAL, DIMENSION (nloc), INTENT (IN) :: supmax0, asupmaxmin | ||
| 70 | |||
| 71 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: da | ||
| 72 | REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: phi !AC! | ||
| 73 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: mp !RomP | ||
| 74 | REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: phi2 !RomP | ||
| 75 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: d1a, dam !RomP | ||
| 76 | REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: sigij !RomP | ||
| 77 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: qta !jyg | ||
| 78 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: clw !RomP | ||
| 79 | REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: elij !RomP | ||
| 80 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: evap, ep !RomP | ||
| 81 | REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: epmlmMm !RomP+jyg | ||
| 82 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: eplamM !RomP+jyg | ||
| 83 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: qtc, sigt !RomP | ||
| 84 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: wdtrainA, wdtrainS, wdtrainM !RomP | ||
| 85 | |||
| 86 | ! outputs: | ||
| 87 | INTEGER, DIMENSION (len), INTENT (OUT) :: iflag1, kbas1, ktop1 | ||
| 88 | REAL, DIMENSION (len), INTENT (OUT) :: precip1, cbmf1, plcl1, plfc1 | ||
| 89 | REAL, DIMENSION (len), INTENT (OUT) :: wbeff1 | ||
| 90 | REAL, DIMENSION (len, nd), INTENT (OUT) :: sig1, w01 | ||
| 91 | REAL, DIMENSION (len), INTENT (OUT) :: epmax_diag1 ! epmax_cape | ||
| 92 | REAL, DIMENSION (len), INTENT (OUT) :: ptop21 | ||
| 93 | REAL, DIMENSION (len, nd), INTENT (OUT) :: ft1, fq1, fu1, fv1 | ||
| 94 | REAL, DIMENSION (len, nd, ntra), INTENT (OUT) :: ftra1 | ||
| 95 | REAL, DIMENSION (len), INTENT (OUT) :: sigd1 | ||
| 96 | REAL, DIMENSION (len, nd), INTENT (OUT) :: ma1, mip1 | ||
| 97 | REAL, DIMENSION (len, nd+1), INTENT (OUT) :: vprecip1 | ||
| 98 | REAL, DIMENSION (len, nd+1), INTENT (OUT) :: vprecipi1 | ||
| 99 | REAL, DIMENSION (len, nd), INTENT (OUT) :: upwd1, dnwd1, dnwd01 | ||
| 100 | REAL, DIMENSION (len, nd), INTENT (OUT) :: qcondc1 | ||
| 101 | REAL, DIMENSION (len), INTENT (OUT) :: wd1, cape1, cin1 | ||
| 102 | REAL, DIMENSION (len, nd), INTENT (OUT) :: tvp1 | ||
| 103 | REAL, DIMENSION (len, nd), INTENT (OUT) :: ftd1, fqd1 | ||
| 104 | REAL, DIMENSION (len), INTENT (OUT) :: plim11, plim21 | ||
| 105 | REAL, DIMENSION (len, nd), INTENT (OUT) :: asupmax1 | ||
| 106 | REAL, DIMENSION (len), INTENT (OUT) :: supmax01, asupmaxmin1 | ||
| 107 | |||
| 108 | REAL, DIMENSION (len, nd), INTENT (OUT) :: da1 | ||
| 109 | REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi1 !AC! | ||
| 110 | REAL, DIMENSION (len, nd), INTENT (OUT) :: mp1 !RomP | ||
| 111 | REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi21 !RomP | ||
| 112 | REAL, DIMENSION (len, nd), INTENT (OUT) :: d1a1, dam1 !RomP !RomP | ||
| 113 | REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: sigij1 !RomP | ||
| 114 | REAL, DIMENSION (len, nd), INTENT (OUT) :: qta1 !jyg | ||
| 115 | REAL, DIMENSION (len, nd), INTENT (OUT) :: clw1 !RomP | ||
| 116 | REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: elij1 !RomP | ||
| 117 | REAL, DIMENSION (len, nd), INTENT (OUT) :: evap1, ep1 !RomP | ||
| 118 | REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: epmlmMm1 !RomP+jyg | ||
| 119 | REAL, DIMENSION (len, nd), INTENT (OUT) :: eplamM1 !RomP+jyg | ||
| 120 | REAL, DIMENSION (len, nd), INTENT (OUT) :: qtc1, sigt1 !RomP | ||
| 121 | REAL, DIMENSION (len, nd), INTENT (OUT) :: wdtrainA1, wdtrainS1, wdtrainM1 !RomP | ||
| 122 | |||
| 123 | |||
| 124 | ! local variables: | ||
| 125 | INTEGER i, k, j | ||
| 126 | INTEGER jdcum | ||
| 127 | ! c integer k1,k2 | ||
| 128 | |||
| 129 | !jyg< | ||
| 130 |
1/2✓ Branch 0 taken 240 times.
✗ Branch 1 not taken.
|
240 | IF (compress) THEN |
| 131 | !>jyg | ||
| 132 |
2/2✓ Branch 0 taken 127097 times.
✓ Branch 1 taken 240 times.
|
127337 | DO i = 1, ncum |
| 133 | 127097 | sig1(idcum(i), nd) = sig(i, nd) | |
| 134 | 127097 | ptop21(idcum(i)) = ptop2(i) | |
| 135 | 127097 | sigd1(idcum(i)) = sigd(i) | |
| 136 | 127097 | precip1(idcum(i)) = precip(i) | |
| 137 | 127097 | cbmf1(idcum(i)) = cbmf(i) | |
| 138 | 127097 | plcl1(idcum(i)) = plcl(i) | |
| 139 | 127097 | plfc1(idcum(i)) = plfc(i) | |
| 140 | 127097 | wbeff1(idcum(i)) = wbeff(i) | |
| 141 | 127097 | iflag1(idcum(i)) = iflag(i) | |
| 142 | 127097 | kbas1(idcum(i)) = kbas(i) | |
| 143 | 127097 | ktop1(idcum(i)) = ktop(i) | |
| 144 | 127097 | wd1(idcum(i)) = wd(i) | |
| 145 | 127097 | cape1(idcum(i)) = cape(i) | |
| 146 | 127097 | cin1(idcum(i)) = cin(i) | |
| 147 | 127097 | plim11(idcum(i)) = plim1(i) | |
| 148 | 127097 | plim21(idcum(i)) = plim2(i) | |
| 149 | 127097 | supmax01(idcum(i)) = supmax0(i) | |
| 150 | 127097 | asupmaxmin1(idcum(i)) = asupmaxmin(i) | |
| 151 | 127337 | epmax_diag1(idcum(i)) = epmax_diag(i) | |
| 152 | END DO | ||
| 153 | |||
| 154 |
2/2✓ Branch 0 taken 6480 times.
✓ Branch 1 taken 240 times.
|
6720 | DO k = 1, nl |
| 155 |
2/2✓ Branch 0 taken 3431619 times.
✓ Branch 1 taken 6480 times.
|
3438339 | DO i = 1, ncum |
| 156 | 3431619 | sig1(idcum(i), k) = sig(i, k) | |
| 157 | 3431619 | w01(idcum(i), k) = w0(i, k) | |
| 158 | 3431619 | ft1(idcum(i), k) = ft(i, k) | |
| 159 | 3431619 | fq1(idcum(i), k) = fq(i, k) | |
| 160 | 3431619 | fu1(idcum(i), k) = fu(i, k) | |
| 161 | 3431619 | fv1(idcum(i), k) = fv(i, k) | |
| 162 | 3431619 | ma1(idcum(i), k) = ma(i, k) | |
| 163 | 3431619 | mip1(idcum(i), k) = mip(i, k) | |
| 164 | 3431619 | vprecip1(idcum(i), k) = vprecip(i, k) | |
| 165 | 3431619 | vprecipi1(idcum(i), k) = vprecipi(i, k) | |
| 166 | 3431619 | upwd1(idcum(i), k) = upwd(i, k) | |
| 167 | 3431619 | dnwd1(idcum(i), k) = dnwd(i, k) | |
| 168 | 3431619 | dnwd01(idcum(i), k) = dnwd0(i, k) | |
| 169 | 3431619 | qcondc1(idcum(i), k) = qcondc(i, k) | |
| 170 | 3431619 | tvp1(idcum(i), k) = tvp(i, k) | |
| 171 | 3431619 | ftd1(idcum(i), k) = ftd(i, k) | |
| 172 | 3431619 | fqd1(idcum(i), k) = fqd(i, k) | |
| 173 | 3431619 | asupmax1(idcum(i), k) = asupmax(i, k) | |
| 174 | |||
| 175 | 3431619 | da1(idcum(i), k) = da(i, k) !AC! | |
| 176 | 3431619 | mp1(idcum(i), k) = mp(i, k) !RomP | |
| 177 | 3431619 | d1a1(idcum(i), k) = d1a(i, k) !RomP | |
| 178 | 3431619 | dam1(idcum(i), k) = dam(i, k) !RomP | |
| 179 | 3431619 | qta1(idcum(i), k) = qta(i, k) !jyg | |
| 180 | 3431619 | clw1(idcum(i), k) = clw(i, k) !RomP | |
| 181 | 3431619 | evap1(idcum(i), k) = evap(i, k) !RomP | |
| 182 | 3431619 | ep1(idcum(i), k) = ep(i, k) !RomP | |
| 183 | 3431619 | eplamM1(idcum(i), k) = eplamM(i, k) !RomP+jyg | |
| 184 | 3431619 | wdtrainA1(idcum(i), k) = wdtrainA(i, k) !RomP | |
| 185 | 3431619 | wdtrainS1(idcum(i), k) = wdtrainS(i, k) !RomP | |
| 186 | 3431619 | wdtrainM1(idcum(i), k) = wdtrainM(i, k) !RomP | |
| 187 | 3431619 | qtc1(idcum(i), k) = qtc(i, k) | |
| 188 | 3438099 | sigt1(idcum(i), k) = sigt(i, k) | |
| 189 | |||
| 190 | END DO | ||
| 191 | END DO | ||
| 192 | |||
| 193 | ! Fluxes are defined on a staggered grid and extend up to nl+1 | ||
| 194 |
2/2✓ Branch 0 taken 240 times.
✓ Branch 1 taken 127097 times.
|
127337 | DO i = 1, ncum |
| 195 | 127097 | ma1(idcum(i), nlp) = 0. | |
| 196 | 127097 | vprecip1(idcum(i), nlp) = 0. | |
| 197 | 127097 | vprecipi1(idcum(i), nlp) = 0. | |
| 198 | 127097 | upwd1(idcum(i), nlp) = 0. | |
| 199 | 127097 | dnwd1(idcum(i), nlp) = 0. | |
| 200 | 127337 | dnwd01(idcum(i), nlp) = 0. | |
| 201 | END DO | ||
| 202 | |||
| 203 | ! AC! do 2100 j=1,ntra | ||
| 204 | ! AC!c oct3 do 2110 k=1,nl | ||
| 205 | ! AC! do 2110 k=1,nd ! oct3 | ||
| 206 | ! AC! do 2120 i=1,ncum | ||
| 207 | ! AC! ftra1(idcum(i),k,j)=ftra(i,k,j) | ||
| 208 | ! AC! 2120 continue | ||
| 209 | ! AC! 2110 continue | ||
| 210 | ! AC! 2100 continue | ||
| 211 | |||
| 212 | ! AC! | ||
| 213 | !jyg< | ||
| 214 | ! Essais pour gagner du temps en diminuant l'adressage indirect | ||
| 215 | !! DO j = 1, nd | ||
| 216 | !! DO k = 1, nd | ||
| 217 | !! DO i = 1, ncum | ||
| 218 | !! phi1(idcum(i), k, j) = phi(i, k, j) !AC! | ||
| 219 | !! phi21(idcum(i), k, j) = phi2(i, k, j) !RomP | ||
| 220 | !! sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP | ||
| 221 | !! elij1(idcum(i), k, j) = elij(i, k, j) !RomP | ||
| 222 | !! epmlmMm(idcum(i), k, j) = epmlmMm(i, k, j) !RomP+jyg | ||
| 223 | !! END DO | ||
| 224 | !! END DO | ||
| 225 | !! END DO | ||
| 226 | |||
| 227 | !! DO i = 1, ncum | ||
| 228 | !! jdcum=idcum(i) | ||
| 229 | !! phi1 (jdcum, 1:nl+1, 1:nl+1) = phi (i, 1:nl+1, 1:nl+1) !AC! | ||
| 230 | !! phi21 (jdcum, 1:nl+1, 1:nl+1) = phi2 (i, 1:nl+1, 1:nl+1) !RomP | ||
| 231 | !! sigij1 (jdcum, 1:nl+1, 1:nl+1) = sigij (i, 1:nl+1, 1:nl+1) !RomP | ||
| 232 | !! elij1 (jdcum, 1:nl+1, 1:nl+1) = elij (i, 1:nl+1, 1:nl+1) !RomP | ||
| 233 | !! epmlmMm1(jdcum, 1:nl+1, 1:nl+1) = epmlmMm(i, 1:nl+1, 1:nl+1) !RomP+jyg | ||
| 234 | !! END DO | ||
| 235 | ! These tracer associated arrays are defined up to nl, not nl+1 | ||
| 236 |
2/2✓ Branch 0 taken 127097 times.
✓ Branch 1 taken 240 times.
|
127337 | DO i = 1, ncum |
| 237 | 127097 | jdcum=idcum(i) | |
| 238 |
2/2✓ Branch 0 taken 3431619 times.
✓ Branch 1 taken 127097 times.
|
3558956 | DO k = 1,nl |
| 239 |
2/2✓ Branch 0 taken 92653713 times.
✓ Branch 1 taken 3431619 times.
|
96212429 | DO j = 1,nl |
| 240 | 92653713 | phi1 (jdcum, j, k) = phi (i, j, k) !AC! | |
| 241 | 92653713 | phi21 (jdcum, j, k) = phi2 (i, j, k) !RomP | |
| 242 | 92653713 | sigij1 (jdcum, j, k) = sigij (i, j, k) !RomP | |
| 243 | 92653713 | elij1 (jdcum, j, k) = elij (i, j, k) !RomP | |
| 244 | 96085332 | epmlmMm1(jdcum, j, k) = epmlmMm(i, j, k) !RomP+jyg | |
| 245 | END DO | ||
| 246 | ENDDO | ||
| 247 | ENDDO | ||
| 248 | !>jyg | ||
| 249 | ! AC! | ||
| 250 | |||
| 251 | |||
| 252 | ! do 2220 k2=1,nd | ||
| 253 | ! do 2210 k1=1,nd | ||
| 254 | ! do 2200 i=1,ncum | ||
| 255 | ! ment1(idcum(i),k1,k2) = ment(i,k1,k2) | ||
| 256 | ! sigij1(idcum(i),k1,k2) = sigij(i,k1,k2) | ||
| 257 | ! 2200 enddo | ||
| 258 | ! 2210 enddo | ||
| 259 | ! 2220 enddo | ||
| 260 | ! | ||
| 261 | !jyg< | ||
| 262 | ELSE !(compress) | ||
| 263 | ! | ||
| 264 | ✗ | sig1(:,nd) = sig(:,nd) | |
| 265 | ✗ | ptop21(:) = ptop2(:) | |
| 266 | ✗ | sigd1(:) = sigd(:) | |
| 267 | ✗ | precip1(:) = precip(:) | |
| 268 | ✗ | cbmf1(:) = cbmf(:) | |
| 269 | ✗ | plcl1(:) = plcl(:) | |
| 270 | ✗ | plfc1(:) = plfc(:) | |
| 271 | ✗ | wbeff1(:) = wbeff(:) | |
| 272 | ✗ | iflag1(:) = iflag(:) | |
| 273 | ✗ | kbas1(:) = kbas(:) | |
| 274 | ✗ | ktop1(:) = ktop(:) | |
| 275 | ✗ | wd1(:) = wd(:) | |
| 276 | ✗ | cape1(:) = cape(:) | |
| 277 | ✗ | cin1(:) = cin(:) | |
| 278 | ✗ | plim11(:) = plim1(:) | |
| 279 | ✗ | plim21(:) = plim2(:) | |
| 280 | ✗ | supmax01(:) = supmax0(:) | |
| 281 | ✗ | asupmaxmin1(:) = asupmaxmin(:) | |
| 282 | ! | ||
| 283 | ✗ | sig1(:, 1:nl) = sig(:, 1:nl) | |
| 284 | ✗ | w01(:, 1:nl) = w0(:, 1:nl) | |
| 285 | ✗ | ft1(:, 1:nl) = ft(:, 1:nl) | |
| 286 | ✗ | fq1(:, 1:nl) = fq(:, 1:nl) | |
| 287 | ✗ | fu1(:, 1:nl) = fu(:, 1:nl) | |
| 288 | ✗ | fv1(:, 1:nl) = fv(:, 1:nl) | |
| 289 | ✗ | ma1(:, 1:nl) = ma(:, 1:nl) | |
| 290 | ✗ | mip1(:, 1:nl) = mip(:, 1:nl) | |
| 291 | ✗ | vprecip1(:, 1:nl) = vprecip(:, 1:nl) | |
| 292 | ✗ | vprecipi1(:, 1:nl) = vprecipi(:, 1:nl) | |
| 293 | ✗ | upwd1(:, 1:nl) = upwd(:, 1:nl) | |
| 294 | ✗ | dnwd1(:, 1:nl) = dnwd(:, 1:nl) | |
| 295 | ✗ | dnwd01(:, 1:nl) = dnwd0(:, 1:nl) | |
| 296 | ✗ | qcondc1(:, 1:nl) = qcondc(:, 1:nl) | |
| 297 | ✗ | tvp1(:, 1:nl) = tvp(:, 1:nl) | |
| 298 | ✗ | ftd1(:, 1:nl) = ftd(:, 1:nl) | |
| 299 | ✗ | fqd1(:, 1:nl) = fqd(:, 1:nl) | |
| 300 | ✗ | asupmax1(:, 1:nl) = asupmax(:, 1:nl) | |
| 301 | |||
| 302 | ✗ | da1(:, 1:nl) = da(:, 1:nl) !AC! | |
| 303 | ✗ | mp1(:, 1:nl) = mp(:, 1:nl) !RomP | |
| 304 | ✗ | d1a1(:, 1:nl) = d1a(:, 1:nl) !RomP | |
| 305 | ✗ | dam1(:, 1:nl) = dam(:, 1:nl) !RomP | |
| 306 | ✗ | qta1(:, 1:nl) = qta(:, 1:nl) !jyg | |
| 307 | ✗ | clw1(:, 1:nl) = clw(:, 1:nl) !RomP | |
| 308 | ✗ | evap1(:, 1:nl) = evap(:, 1:nl) !RomP | |
| 309 | ✗ | ep1(:, 1:nl) = ep(:, 1:nl) !RomP | |
| 310 | ✗ | eplamM1(:, 1:nl) = eplamM(:, 1:nl) !RomP+jyg | |
| 311 | ✗ | wdtrainA1(:, 1:nl) = wdtrainA(:, 1:nl) !RomP | |
| 312 | ✗ | wdtrainS1(:, 1:nl) = wdtrainS(:, 1:nl) !RomP | |
| 313 | ✗ | wdtrainM1(:, 1:nl) = wdtrainM(:, 1:nl) !RomP | |
| 314 | ✗ | qtc1(:, 1:nl) = qtc(:, 1:nl) | |
| 315 | ✗ | sigt1(:, 1:nl) = sigt(:, 1:nl) | |
| 316 | ! | ||
| 317 | ✗ | ma1(:, nlp) = 0. | |
| 318 | ✗ | vprecip1(:, nlp) = 0. | |
| 319 | ✗ | vprecipi1(:, nlp) = 0. | |
| 320 | ✗ | upwd1(:, nlp) = 0. | |
| 321 | ✗ | dnwd1(:, nlp) = 0. | |
| 322 | ✗ | dnwd01(:, nlp) = 0. | |
| 323 | |||
| 324 | ! | ||
| 325 | ✗ | phi1 (:, 1:nl, 1:nl) = phi (:, 1:nl, 1:nl) !AC! | |
| 326 | ✗ | phi21 (:, 1:nl, 1:nl) = phi2 (:, 1:nl, 1:nl) !RomP | |
| 327 | ✗ | sigij1 (:, 1:nl, 1:nl) = sigij (:, 1:nl, 1:nl) !RomP | |
| 328 | ✗ | elij1 (:, 1:nl, 1:nl) = elij (:, 1:nl, 1:nl) !RomP | |
| 329 | ✗ | epmlmMm1(:, 1:nl, 1:nl) = epmlmMm(:, 1:nl, 1:nl) !RomP+jyg | |
| 330 | ENDIF !(compress) | ||
| 331 | !>jyg | ||
| 332 | |||
| 333 | 240 | RETURN | |
| 334 | END SUBROUTINE cv3a_uncompress | ||
| 335 | |||
| 336 |