Changeset 1673 for LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F
- Timestamp:
- Oct 27, 2012, 4:23:07 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F
r1632 r1673 549 549 C PRINT*,dyq(1) 550 550 C PRINT*,dyqv(iip1+1) 551 C ap n=abs(dyq(1)/dyqv(iip1+1))551 C appn=abs(dyq(1)/dyqv(iip1+1)) 552 552 C PRINT*,dyq(ip1jm+1) 553 553 C PRINT*,dyqv(ip1jm-iip1+1) 554 C ap s=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1))554 C apps=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1)) 555 555 C DO ij=2,iim 556 C ap n=amax1(abs(dyq(ij)/dyqv(ij)),apn)557 C ap s=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),aps)556 C appn=amax1(abs(dyq(ij)/dyqv(ij)),appn) 557 C apps=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),apps) 558 558 C ENDDO 559 C ap n=min(pente_max/apn,1.)560 C ap s=min(pente_max/aps,1.)559 C appn=min(pente_max/appn,1.) 560 C apps=min(pente_max/apps,1.) 561 561 C 562 562 C … … 564 564 C 565 565 C IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.) 566 C & ap n=0.566 C & appn=0. 567 567 C IF(dyqv(ismax(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1)* 568 568 C & dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.) 569 C & ap s=0.569 C & apps=0. 570 570 C 571 571 C limitation des pentes aux poles 572 572 C DO ij=1,iip1 573 C dyq(ij)=ap n*dyq(ij)574 C dyq(ip1jm+ij)=ap s*dyq(ip1jm+ij)573 C dyq(ij)=appn*dyq(ij) 574 C dyq(ip1jm+ij)=apps*dyq(ip1jm+ij) 575 575 C ENDDO 576 576 C
Note: See TracChangeset
for help on using the changeset viewer.