Changeset 5079 for LMDZ6/trunk/libf/dyn3d_common/advzp.F
- Timestamp:
- Jul 19, 2024, 11:28:59 AM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d_common/advzp.F
r2600 r5079 135 135 C Conversion des flux de masses en kg 136 136 137 DO 500l = 1,llm138 DO 500j = 1,jjp1139 DO 500i = 1,iip1137 DO l = 1,llm 138 DO j = 1,jjp1 139 DO i = 1,iip1 140 140 wgri (i,j,llm+1-l) = w (i,j,l) 141 500 CONTINUE 141 END DO 142 END DO 143 END DO 142 144 do j=1,jjp1 143 145 do i=1,iip1 … … 154 156 C boucle sur les latitudes 155 157 C 156 DO 1K=1,LAT158 DO K=1,LAT 157 159 C 158 160 C place limits on appropriate moments before transport … … 161 163 IF(.NOT.LIMIT) GO TO 101 162 164 C 163 DO 10JV=1,NTRA164 DO 10L=1,NIV165 DO 100I=1,LON166 IF(S0(I,K,L,JV) .GT.0.) THEN165 DO JV=1,NTRA 166 DO L=1,NIV 167 DO I=1,LON 168 IF(S0(I,K,L,JV)>0.) THEN 167 169 SLPMAX=S0(I,K,L,JV) 168 170 S1MAX =1.5*SLPMAX … … 180 182 SYZ(I,K,L,JV)=0. 181 183 ENDIF 182 100 CONTINUE 183 10 CONTINUE 184 END DO 185 END DO 186 END DO 184 187 C 185 188 101 CONTINUE … … 193 196 C 2- reajusts moments remaining in the box 194 197 C 195 DO 11L=1,NIV-1198 DO L=1,NIV-1 196 199 LP=L+1 197 200 C 198 DO 110I=1,LON199 C 200 IF(WGRI(I,K,L) .LT.0.) THEN201 DO I=1,LON 202 C 203 IF(WGRI(I,K,L)<0.) THEN 201 204 FM(I,L)=-WGRI(I,K,L)*DTZ 202 205 ALF(I)=FM(I,L)/SM(I,K,LP) … … 215 218 ALF4 (I)=ALF1(I)*ALF1Q(I) 216 219 C 217 110 CONTINUE218 C 219 DO 111JV=1,NTRA220 DO 1110I=1,LON221 C 222 IF(WGRI(I,K,L) .LT.0.) THEN220 END DO 221 C 222 DO JV=1,NTRA 223 DO I=1,LON 224 C 225 IF(WGRI(I,K,L)<0.) THEN 223 226 C 224 227 F0 (I,L,JV)=ALF (I)* ( S0(I,K,LP,JV)-ALF1(I)* … … 273 276 ENDIF 274 277 C 275 1110 CONTINUE276 111 CONTINUE277 C 278 11 CONTINUE278 END DO 279 END DO 280 C 281 END DO 279 282 C 280 283 C puts the temporary moments Fi into appropriate neighboring boxes 281 284 C 282 DO 12L=1,NIV-1285 DO L=1,NIV-1 283 286 LP=L+1 284 287 C 285 DO 120I=1,LON286 C 287 IF(WGRI(I,K,L) .LT.0.) THEN288 DO I=1,LON 289 C 290 IF(WGRI(I,K,L)<0.) THEN 288 291 SM(I,K,L)=SM(I,K,L)+FM(I,L) 289 292 ALF(I)=FM(I,L)/SM(I,K,L) … … 299 302 ALF3(I)=ALF1(I)-ALF(I) 300 303 C 301 120 CONTINUE302 C 303 DO 121JV=1,NTRA304 DO 1210I=1,LON305 C 306 IF(WGRI(I,K,L) .LT.0.) THEN304 END DO 305 C 306 DO JV=1,NTRA 307 DO I=1,LON 308 C 309 IF(WGRI(I,K,L)<0.) THEN 307 310 C 308 311 TEMPTM=-ALF(I)*S0(I,K,L,JV)+ALF1(I)*F0(I,L,JV) … … 342 345 ENDIF 343 346 C 344 1210 CONTINUE345 121 CONTINUE346 C 347 12 CONTINUE347 END DO 348 END DO 349 C 350 END DO 348 351 C 349 352 C fin de la boucle principale sur les latitudes 350 353 C 351 1 CONTINUE354 END DO 352 355 C 353 356 DO l = 1,llm
Note: See TracChangeset
for help on using the changeset viewer.