Changeset 2298 for LMDZ5/branches/testing/libf/phylmd/cv3a_uncompress.F90
- Timestamp:
- Jun 14, 2015, 9:13:32 PM (9 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2238-2257,2259-2271,2273,2277-2282,2284-2288,2290-2291
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/cv3a_uncompress.F90
r2220 r2298 1 SUBROUTINE cv3a_uncompress(nloc, len, ncum, nd, ntra, idcum, iflag, kbas, & 2 ktop, precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, ft, fq, fu, fv, & 3 ftra, sigd, ma, mip, vprecip, upwd, dnwd, dnwd0, qcondc, wd, cape, cin, & 4 tvp, ftd, fqd, plim1, plim2, asupmax, supmax0, asupmaxmin & 5 , da, phi, mp, phi2, d1a, dam, sigij & ! RomP+AC+jyg 6 , clw, elij, evap, ep, epmlmmm, eplamm & ! RomP 7 , wdtraina, wdtrainm & ! RomP 8 , qtc, sigt & 9 10 , iflag1, kbas1, ktop1, precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, & 11 ptop21, ft1, fq1, fu1, fv1, ftra1, sigd1, ma1, mip1, vprecip1, upwd1, & 12 dnwd1, dnwd01, qcondc1, wd1, cape1, cin1, tvp1, ftd1, fqd1, plim11, & 13 plim21, asupmax1, supmax01, asupmaxmin1 & 14 , da1, phi1, mp1, phi21, d1a1, dam1, sigij1 & ! RomP+AC+jyg 15 , clw1, elij1, evap1, ep1, epmlmmm1, eplamm1 & ! RomP 16 , wdtraina1, wdtrainm1 & ! RomP 17 , qtc1, sigt1) 1 SUBROUTINE cv3a_uncompress(nloc, len, ncum, nd, ntra, idcum, compress, & 2 iflag, kbas, ktop, & 3 precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, & 4 ft, fq, fu, fv, ftra, & 5 sigd, ma, mip, vprecip, upwd, dnwd, dnwd0, & 6 qcondc, wd, cape, cin, & 7 tvp, & 8 ftd, fqd, & 9 plim1, plim2, asupmax, supmax0, & 10 asupmaxmin, & 11 da, phi, mp, phi2, d1a, dam, sigij, & ! RomP+AC+jyg 12 clw, elij, evap, ep, epmlmMm, eplaMm, & ! RomP 13 wdtrainA, wdtrainM, & ! RomP 14 qtc, sigt, & 15 16 iflag1, kbas1, ktop1, & 17 precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, ptop21, & 18 ft1, fq1, fu1, fv1, ftra1, & 19 sigd1, ma1, mip1, vprecip1, upwd1, dnwd1, dnwd01, & 20 qcondc1, wd1, cape1, cin1, & 21 tvp1, & 22 ftd1, fqd1, & 23 plim11, plim21, asupmax1, supmax01, & 24 asupmaxmin1, & 25 da1, phi1, mp1, phi21, d1a1, dam1, sigij1, & ! RomP+AC+jyg 26 clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP 27 wdtrainA1, wdtrainM1, & ! RomP 28 qtc1, sigt1) 18 29 19 30 ! ************************************************************** … … 31 42 32 43 ! inputs: 33 INTEGER nloc, len, ncum, nd, ntra 34 INTEGER idcum(nloc) 35 INTEGER iflag(nloc), kbas(nloc), ktop(nloc) 36 REAL precip(nloc), cbmf(nloc), plcl(nloc), plfc(nloc) 37 REAL wbeff(len) 38 REAL sig(nloc, nd), w0(nloc, nd), ptop2(nloc) 39 REAL ft(nloc, nd), fq(nloc, nd), fu(nloc, nd), fv(nloc, nd) 40 REAL ftra(nloc, nd, ntra) 41 REAL sigd(nloc) 42 REAL ma(nloc, nd), mip(nloc, nd), vprecip(nloc, nd+1) 43 REAL upwd(nloc, nd), dnwd(nloc, nd), dnwd0(nloc, nd) 44 REAL qcondc(nloc, nd) 45 REAL wd(nloc), cape(nloc), cin(nloc) 46 REAL tvp(nloc, nd) 47 REAL ftd(nloc, nd), fqd(nloc, nd) 48 REAL plim1(nloc), plim2(nloc) 49 REAL asupmax(nloc, nd), supmax0(nloc) 50 REAL asupmaxmin(nloc) 51 52 REAL da(nloc, nd), phi(nloc, nd, nd) !AC! 53 REAL mp(nloc, nd) !RomP 54 REAL phi2(nloc, nd, nd) !RomP 55 REAL d1a(nloc, nd), dam(nloc, nd) !RomP 56 REAL sigij(nloc, nd, nd) !RomP 57 REAL clw(nloc, nd), elij(nloc, nd, nd) !RomP 58 REAL evap(nloc, nd), ep(nloc, nd) !RomP 59 REAL epmlmmm(nloc, nd, nd), eplamm(nloc, nd) !RomP+jyg 60 REAL qtc(nloc, nd), sigt(nloc, nd) !RomP 61 REAL wdtraina(nloc, nd), wdtrainm(nloc, nd) !RomP 44 INTEGER, INTENT (IN) :: nloc, len, ncum, nd, ntra 45 INTEGER, DIMENSION (nloc), INTENT (IN) :: idcum(nloc) 46 !jyg< 47 LOGICAL, INTENT (IN) :: compress 48 !>jyg 49 INTEGER, DIMENSION (nloc), INTENT (IN) ::iflag, kbas, ktop 50 REAL, DIMENSION (nloc), INTENT (IN) :: precip, cbmf, plcl, plfc 51 REAL, DIMENSION (nloc), INTENT (IN) :: wbeff 52 REAL, DIMENSION (nloc, nd), INTENT (IN) :: sig, w0 53 REAL, DIMENSION (nloc), INTENT (IN) :: ptop2 54 REAL, DIMENSION (nloc, nd), INTENT (IN) :: ft, fq, fu, fv 55 REAL, DIMENSION (nloc, nd, ntra), INTENT (IN) :: ftra 56 REAL, DIMENSION (nloc), INTENT (IN) :: sigd 57 REAL, DIMENSION (nloc, nd), INTENT (IN) :: ma, mip 58 REAL, DIMENSION (nloc, nd+1), INTENT (IN) :: vprecip 59 REAL, DIMENSION (nloc, nd), INTENT (IN) :: upwd, dnwd, dnwd0 60 REAL, DIMENSION (nloc, nd), INTENT (IN) :: qcondc 61 REAL, DIMENSION (nloc), INTENT (IN) :: wd, cape, cin 62 REAL, DIMENSION (nloc, nd), INTENT (IN) :: tvp 63 REAL, DIMENSION (nloc, nd), INTENT (IN) :: ftd, fqd 64 REAL, DIMENSION (nloc), INTENT (IN) :: plim1, plim2 65 REAL, DIMENSION (nloc, nd), INTENT (IN) :: asupmax 66 REAL, DIMENSION (nloc), INTENT (IN) :: supmax0, asupmaxmin 67 68 REAL, DIMENSION (nloc, nd), INTENT (IN) :: da 69 REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: phi !AC! 70 REAL, DIMENSION (nloc, nd), INTENT (IN) :: mp !RomP 71 REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: phi2 !RomP 72 REAL, DIMENSION (nloc, nd), INTENT (IN) :: d1a, dam !RomP 73 REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: sigij !RomP 74 REAL, DIMENSION (nloc, nd), INTENT (IN) :: clw !RomP 75 REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: elij !RomP 76 REAL, DIMENSION (nloc, nd), INTENT (IN) :: evap, ep !RomP 77 REAL, DIMENSION (nloc, nd, nd), INTENT (IN) :: epmlmMm !RomP+jyg 78 REAL, DIMENSION (nloc, nd), INTENT (IN) :: eplamM !RomP+jyg 79 REAL, DIMENSION (nloc, nd), INTENT (IN) :: qtc, sigt !RomP 80 REAL, DIMENSION (nloc, nd), INTENT (IN) :: wdtrainA, wdtrainM !RomP 62 81 63 82 ! outputs: 64 INTEGER iflag1(len), kbas1(len), ktop1(len) 65 REAL precip1(len), cbmf1(len), plcl1(nloc), plfc1(nloc) 66 REAL wbeff1(len) 67 REAL sig1(len, nd), w01(len, nd), ptop21(len) 68 REAL ft1(len, nd), fq1(len, nd), fu1(len, nd), fv1(len, nd) 69 REAL ftra1(len, nd, ntra) 70 REAL sigd1(len) 71 REAL ma1(len, nd), mip1(len, nd), vprecip1(len, nd+1) 72 REAL upwd1(len, nd), dnwd1(len, nd), dnwd01(len, nd) 73 REAL qcondc1(len, nd) 74 REAL wd1(len), cape1(len), cin1(len) 75 REAL tvp1(len, nd) 76 REAL ftd1(len, nd), fqd1(len, nd) 77 REAL plim11(len), plim21(len) 78 REAL asupmax1(len, nd), supmax01(len) 79 REAL asupmaxmin1(len) 80 81 REAL da1(nloc, nd), phi1(nloc, nd, nd) !AC! 82 REAL mp1(nloc, nd) !RomP 83 REAL phi21(nloc, nd, nd) !RomP 84 REAL d1a1(nloc, nd), dam1(nloc, nd) !RomP 85 REAL sigij1(len, nd, nd) !RomP 86 REAL clw1(len, nd), elij1(len, nd, nd) !RomP 87 REAL evap1(len, nd), ep1(len, nd) !RomP 88 REAL epmlmmm1(len, nd, nd), eplamm1(len, nd) !RomP+jyg 89 REAL qtc1(len, nd), sigt1(len, nd) !RomP 90 REAL wdtraina1(len, nd), wdtrainm1(len, nd) !RomP 83 INTEGER, DIMENSION (len), INTENT (OUT) :: iflag1, kbas1, ktop1 84 REAL, DIMENSION (len), INTENT (OUT) :: precip1, cbmf1, plcl1, plfc1 85 REAL, DIMENSION (len), INTENT (OUT) :: wbeff1 86 REAL, DIMENSION (len, nd), INTENT (OUT) :: sig1, w01 87 REAL, DIMENSION (len), INTENT (OUT) :: ptop21 88 REAL, DIMENSION (len, nd), INTENT (OUT) :: ft1, fq1, fu1, fv1 89 REAL, DIMENSION (len, nd, ntra), INTENT (OUT) :: ftra1 90 REAL, DIMENSION (len), INTENT (OUT) :: sigd1 91 REAL, DIMENSION (len, nd), INTENT (OUT) :: ma1, mip1 92 REAL, DIMENSION (len, nd+1), INTENT (OUT) :: vprecip1 93 REAL, DIMENSION (len, nd), INTENT (OUT) :: upwd1, dnwd1, dnwd01 94 REAL, DIMENSION (len, nd), INTENT (OUT) :: qcondc1 95 REAL, DIMENSION (len), INTENT (OUT) :: wd1, cape1, cin1 96 REAL, DIMENSION (len, nd), INTENT (OUT) :: tvp1 97 REAL, DIMENSION (len, nd), INTENT (OUT) :: ftd1, fqd1 98 REAL, DIMENSION (len), INTENT (OUT) :: plim11, plim21 99 REAL, DIMENSION (len, nd), INTENT (OUT) :: asupmax1 100 REAL, DIMENSION (len), INTENT (OUT) :: supmax01, asupmaxmin1 101 102 REAL, DIMENSION (len, nd), INTENT (OUT) :: da1 103 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi1 !AC! 104 REAL, DIMENSION (len, nd), INTENT (OUT) :: mp1 !RomP 105 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: phi21 !RomP 106 REAL, DIMENSION (len, nd), INTENT (OUT) :: d1a1, dam1 !RomP !RomP 107 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: sigij1 !RomP 108 REAL, DIMENSION (len, nd), INTENT (OUT) :: clw1 !RomP 109 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: elij1 !RomP 110 REAL, DIMENSION (len, nd), INTENT (OUT) :: evap1, ep1 !RomP 111 REAL, DIMENSION (len, nd, nd), INTENT (OUT) :: epmlmMm1 !RomP+jyg 112 REAL, DIMENSION (len, nd), INTENT (OUT) :: eplamM1 !RomP+jyg 113 REAL, DIMENSION (len, nd), INTENT (OUT) :: qtc1, sigt1 !RomP 114 REAL, DIMENSION (len, nd), INTENT (OUT) :: wdtrainA1, wdtrainM1 !RomP 91 115 92 116 93 117 ! local variables: 94 118 INTEGER i, k, j 119 INTEGER jdcum 95 120 ! c integer k1,k2 96 121 97 DO i = 1, ncum 98 ptop21(idcum(i)) = ptop2(i) 99 sigd1(idcum(i)) = sigd(i) 100 precip1(idcum(i)) = precip(i) 101 cbmf1(idcum(i)) = cbmf(i) 102 plcl1(idcum(i)) = plcl(i) 103 plfc1(idcum(i)) = plfc(i) 104 wbeff1(idcum(i)) = wbeff(i) 105 iflag1(idcum(i)) = iflag(i) 106 kbas1(idcum(i)) = kbas(i) 107 ktop1(idcum(i)) = ktop(i) 108 wd1(idcum(i)) = wd(i) 109 cape1(idcum(i)) = cape(i) 110 cin1(idcum(i)) = cin(i) 111 plim11(idcum(i)) = plim1(i) 112 plim21(idcum(i)) = plim2(i) 113 supmax01(idcum(i)) = supmax0(i) 114 asupmaxmin1(idcum(i)) = asupmaxmin(i) 115 END DO 116 117 DO k = 1, nd 122 !jyg< 123 IF (compress) THEN 124 !>jyg 118 125 DO i = 1, ncum 119 sig1(idcum(i), k) = sig(i, k) 120 w01(idcum(i), k) = w0(i, k) 121 ft1(idcum(i), k) = ft(i, k) 122 fq1(idcum(i), k) = fq(i, k) 123 fu1(idcum(i), k) = fu(i, k) 124 fv1(idcum(i), k) = fv(i, k) 125 ma1(idcum(i), k) = ma(i, k) 126 mip1(idcum(i), k) = mip(i, k) 127 vprecip1(idcum(i), k) = vprecip(i, k) 128 upwd1(idcum(i), k) = upwd(i, k) 129 dnwd1(idcum(i), k) = dnwd(i, k) 130 dnwd01(idcum(i), k) = dnwd0(i, k) 131 qcondc1(idcum(i), k) = qcondc(i, k) 132 tvp1(idcum(i), k) = tvp(i, k) 133 ftd1(idcum(i), k) = ftd(i, k) 134 fqd1(idcum(i), k) = fqd(i, k) 135 asupmax1(idcum(i), k) = asupmax(i, k) 136 137 da1(idcum(i), k) = da(i, k) !AC! 138 mp1(idcum(i), k) = mp(i, k) !RomP 139 d1a1(idcum(i), k) = d1a(i, k) !RomP 140 dam1(idcum(i), k) = dam(i, k) !RomP 141 clw1(idcum(i), k) = clw(i, k) !RomP 142 evap1(idcum(i), k) = evap(i, k) !RomP 143 ep1(idcum(i), k) = ep(i, k) !RomP 144 eplamm(idcum(i), k) = eplamm(i, k) !RomP+jyg 145 wdtraina1(idcum(i), k) = wdtraina(i, k) !RomP 146 wdtrainm1(idcum(i), k) = wdtrainm(i, k) !RomP 147 qtc1(idcum(i), k) = qtc(i, k) 148 sigt1(idcum(i), k) = sigt(i, k) 149 126 sig1(idcum(i), nd) = sig(i, nd) 127 ptop21(idcum(i)) = ptop2(i) 128 sigd1(idcum(i)) = sigd(i) 129 precip1(idcum(i)) = precip(i) 130 cbmf1(idcum(i)) = cbmf(i) 131 plcl1(idcum(i)) = plcl(i) 132 plfc1(idcum(i)) = plfc(i) 133 wbeff1(idcum(i)) = wbeff(i) 134 iflag1(idcum(i)) = iflag(i) 135 kbas1(idcum(i)) = kbas(i) 136 ktop1(idcum(i)) = ktop(i) 137 wd1(idcum(i)) = wd(i) 138 cape1(idcum(i)) = cape(i) 139 cin1(idcum(i)) = cin(i) 140 plim11(idcum(i)) = plim1(i) 141 plim21(idcum(i)) = plim2(i) 142 supmax01(idcum(i)) = supmax0(i) 143 asupmaxmin1(idcum(i)) = asupmaxmin(i) 150 144 END DO 151 END DO 152 153 DO i = 1, ncum 154 sig1(idcum(i), nd) = sig(i, nd) 155 END DO 156 157 158 ! AC! do 2100 j=1,ntra 159 ! AC!c oct3 do 2110 k=1,nl 160 ! AC! do 2110 k=1,nd ! oct3 161 ! AC! do 2120 i=1,ncum 162 ! AC! ftra1(idcum(i),k,j)=ftra(i,k,j) 163 ! AC! 2120 continue 164 ! AC! 2110 continue 165 ! AC! 2100 continue 166 167 ! AC! 168 DO j = 1, nd 169 DO k = 1, nd 145 146 DO k = 1, nl+1 170 147 DO i = 1, ncum 171 phi1(idcum(i), k, j) = phi(i, k, j) !AC! 172 phi21(idcum(i), k, j) = phi2(i, k, j) !RomP 173 sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP 174 elij1(idcum(i), k, j) = elij(i, k, j) !RomP 175 epmlmmm(idcum(i), k, j) = epmlmmm(i, k, j) !RomP+jyg 148 sig1(idcum(i), k) = sig(i, k) 149 w01(idcum(i), k) = w0(i, k) 150 ft1(idcum(i), k) = ft(i, k) 151 fq1(idcum(i), k) = fq(i, k) 152 fu1(idcum(i), k) = fu(i, k) 153 fv1(idcum(i), k) = fv(i, k) 154 ma1(idcum(i), k) = ma(i, k) 155 mip1(idcum(i), k) = mip(i, k) 156 vprecip1(idcum(i), k) = vprecip(i, k) 157 upwd1(idcum(i), k) = upwd(i, k) 158 dnwd1(idcum(i), k) = dnwd(i, k) 159 dnwd01(idcum(i), k) = dnwd0(i, k) 160 qcondc1(idcum(i), k) = qcondc(i, k) 161 tvp1(idcum(i), k) = tvp(i, k) 162 ftd1(idcum(i), k) = ftd(i, k) 163 fqd1(idcum(i), k) = fqd(i, k) 164 asupmax1(idcum(i), k) = asupmax(i, k) 165 166 da1(idcum(i), k) = da(i, k) !AC! 167 mp1(idcum(i), k) = mp(i, k) !RomP 168 d1a1(idcum(i), k) = d1a(i, k) !RomP 169 dam1(idcum(i), k) = dam(i, k) !RomP 170 clw1(idcum(i), k) = clw(i, k) !RomP 171 evap1(idcum(i), k) = evap(i, k) !RomP 172 ep1(idcum(i), k) = ep(i, k) !RomP 173 eplamM1(idcum(i), k) = eplamM(i, k) !RomP+jyg 174 wdtrainA1(idcum(i), k) = wdtrainA(i, k) !RomP 175 wdtrainM1(idcum(i), k) = wdtrainM(i, k) !RomP 176 qtc1(idcum(i), k) = qtc(i, k) 177 sigt1(idcum(i), k) = sigt(i, k) 178 176 179 END DO 177 180 END DO 178 END DO 179 ! AC! 180 181 182 ! do 2220 k2=1,nd 183 ! do 2210 k1=1,nd 184 ! do 2200 i=1,ncum 185 ! ment1(idcum(i),k1,k2) = ment(i,k1,k2) 186 ! sigij1(idcum(i),k1,k2) = sigij(i,k1,k2) 187 ! 2200 enddo 188 ! 2210 enddo 189 ! 2220 enddo 181 182 ! AC! do 2100 j=1,ntra 183 ! AC!c oct3 do 2110 k=1,nl 184 ! AC! do 2110 k=1,nd ! oct3 185 ! AC! do 2120 i=1,ncum 186 ! AC! ftra1(idcum(i),k,j)=ftra(i,k,j) 187 ! AC! 2120 continue 188 ! AC! 2110 continue 189 ! AC! 2100 continue 190 191 ! AC! 192 !jyg< 193 ! Essais pour gagner du temps en diminuant l'adressage indirect 194 !! DO j = 1, nd 195 !! DO k = 1, nd 196 !! DO i = 1, ncum 197 !! phi1(idcum(i), k, j) = phi(i, k, j) !AC! 198 !! phi21(idcum(i), k, j) = phi2(i, k, j) !RomP 199 !! sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP 200 !! elij1(idcum(i), k, j) = elij(i, k, j) !RomP 201 !! epmlmMm(idcum(i), k, j) = epmlmMm(i, k, j) !RomP+jyg 202 !! END DO 203 !! END DO 204 !! END DO 205 DO i = 1, ncum 206 jdcum=idcum(i) 207 phi1 (jdcum, 1:nl+1, 1:nl+1) = phi (i, 1:nl+1, 1:nl+1) !AC! 208 phi21 (jdcum, 1:nl+1, 1:nl+1) = phi2 (i, 1:nl+1, 1:nl+1) !RomP 209 sigij1 (jdcum, 1:nl+1, 1:nl+1) = sigij (i, 1:nl+1, 1:nl+1) !RomP 210 elij1 (jdcum, 1:nl+1, 1:nl+1) = elij (i, 1:nl+1, 1:nl+1) !RomP 211 epmlmMm1(jdcum, 1:nl+1, 1:nl+1) = epmlmMm(i, 1:nl+1, 1:nl+1) !RomP+jyg 212 END DO 213 !>jyg 214 ! AC! 215 216 217 ! do 2220 k2=1,nd 218 ! do 2210 k1=1,nd 219 ! do 2200 i=1,ncum 220 ! ment1(idcum(i),k1,k2) = ment(i,k1,k2) 221 ! sigij1(idcum(i),k1,k2) = sigij(i,k1,k2) 222 ! 2200 enddo 223 ! 2210 enddo 224 ! 2220 enddo 225 ! 226 !jyg< 227 ELSE !(compress) 228 ! 229 sig1(:,nd) = sig(:,nd) 230 ptop21(:) = ptop2(:) 231 sigd1(:) = sigd(:) 232 precip1(:) = precip(:) 233 cbmf1(:) = cbmf(:) 234 plcl1(:) = plcl(:) 235 plfc1(:) = plfc(:) 236 wbeff1(:) = wbeff(:) 237 iflag1(:) = iflag(:) 238 kbas1(:) = kbas(:) 239 ktop1(:) = ktop(:) 240 wd1(:) = wd(:) 241 cape1(:) = cape(:) 242 cin1(:) = cin(:) 243 plim11(:) = plim1(:) 244 plim21(:) = plim2(:) 245 supmax01(:) = supmax0(:) 246 asupmaxmin1(:) = asupmaxmin(:) 247 ! 248 sig1(:, 1:nl+1) = sig(:, 1:nl+1) 249 w01(:, 1:nl+1) = w0(:, 1:nl+1) 250 ft1(:, 1:nl+1) = ft(:, 1:nl+1) 251 fq1(:, 1:nl+1) = fq(:, 1:nl+1) 252 fu1(:, 1:nl+1) = fu(:, 1:nl+1) 253 fv1(:, 1:nl+1) = fv(:, 1:nl+1) 254 ma1(:, 1:nl+1) = ma(:, 1:nl+1) 255 mip1(:, 1:nl+1) = mip(:, 1:nl+1) 256 vprecip1(:, 1:nl+1) = vprecip(:, 1:nl+1) 257 upwd1(:, 1:nl+1) = upwd(:, 1:nl+1) 258 dnwd1(:, 1:nl+1) = dnwd(:, 1:nl+1) 259 dnwd01(:, 1:nl+1) = dnwd0(:, 1:nl+1) 260 qcondc1(:, 1:nl+1) = qcondc(:, 1:nl+1) 261 tvp1(:, 1:nl+1) = tvp(:, 1:nl+1) 262 ftd1(:, 1:nl+1) = ftd(:, 1:nl+1) 263 fqd1(:, 1:nl+1) = fqd(:, 1:nl+1) 264 asupmax1(:, 1:nl+1) = asupmax(:, 1:nl+1) 265 266 da1(:, 1:nl+1) = da(:, 1:nl+1) !AC! 267 mp1(:, 1:nl+1) = mp(:, 1:nl+1) !RomP 268 d1a1(:, 1:nl+1) = d1a(:, 1:nl+1) !RomP 269 dam1(:, 1:nl+1) = dam(:, 1:nl+1) !RomP 270 clw1(:, 1:nl+1) = clw(:, 1:nl+1) !RomP 271 evap1(:, 1:nl+1) = evap(:, 1:nl+1) !RomP 272 ep1(:, 1:nl+1) = ep(:, 1:nl+1) !RomP 273 eplamM1(:, 1:nl+1) = eplamM(:, 1:nl+1) !RomP+jyg 274 wdtrainA1(:, 1:nl+1) = wdtrainA(:, 1:nl+1) !RomP 275 wdtrainM1(:, 1:nl+1) = wdtrainM(:, 1:nl+1) !RomP 276 qtc1(:, 1:nl+1) = qtc(:, 1:nl+1) 277 sigt1(:, 1:nl+1) = sigt(:, 1:nl+1) 278 ! 279 phi1 (:, 1:nl+1, 1:nl+1) = phi (:, 1:nl+1, 1:nl+1) !AC! 280 phi21 (:, 1:nl+1, 1:nl+1) = phi2 (:, 1:nl+1, 1:nl+1) !RomP 281 sigij1 (:, 1:nl+1, 1:nl+1) = sigij (:, 1:nl+1, 1:nl+1) !RomP 282 elij1 (:, 1:nl+1, 1:nl+1) = elij (:, 1:nl+1, 1:nl+1) !RomP 283 epmlmMm1(:, 1:nl+1, 1:nl+1) = epmlmMm(:, 1:nl+1, 1:nl+1) !RomP+jyg 284 ENDIF !(compress) 285 !>jyg 190 286 191 287 RETURN
Note: See TracChangeset
for help on using the changeset viewer.