Changeset 1707 for LMDZ5/branches/testing/libf/dyn3dmem/vlsplt_loc.F
- Timestamp:
- Jan 11, 2013, 10:19:19 AM (11 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 1670-1692,1694-1703,1705-1706
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/dyn3dmem/vlsplt_loc.F
r1669 r1707 1 ! 2 ! $Id$ 3 ! 1 4 SUBROUTINE vlx_loc(q,pente_max,masse,u_m,ijb_x,ije_x) 2 5 … … 372 375 REAL qbyv(ijb_v:ije_v,llm) 373 376 374 REAL qpns,qpsn,ap n,aps,dyn1,dys1,dyn2,dys2,newmasse,fn,fs377 REAL qpns,qpsn,appn,apps,dyn1,dys1,dyn2,dys2,newmasse,fn,fs 375 378 c REAL newq,oldmasse 376 379 Logical extremum,first,testcpu … … 543 546 C PRINT*,dyq(1) 544 547 C PRINT*,dyqv(iip1+1) 545 C ap n=abs(dyq(1)/dyqv(iip1+1))548 C appn=abs(dyq(1)/dyqv(iip1+1)) 546 549 C PRINT*,dyq(ip1jm+1) 547 550 C PRINT*,dyqv(ip1jm-iip1+1) 548 C ap s=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1))551 C apps=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1)) 549 552 C DO ij=2,iim 550 C ap n=amax1(abs(dyq(ij)/dyqv(ij)),apn)551 C ap s=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),aps)553 C appn=amax1(abs(dyq(ij)/dyqv(ij)),appn) 554 C apps=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),apps) 552 555 C ENDDO 553 C ap n=min(pente_max/apn,1.)554 C ap s=min(pente_max/aps,1.)556 C appn=min(pente_max/appn,1.) 557 C apps=min(pente_max/apps,1.) 555 558 C 556 559 C … … 558 561 C 559 562 C IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.) 560 C & ap n=0.563 C & appn=0. 561 564 C IF(dyqv(ismax(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1)* 562 565 C & dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.) 563 C & ap s=0.566 C & apps=0. 564 567 C 565 568 C limitation des pentes aux poles 566 569 C DO ij=1,iip1 567 C dyq(ij)=ap n*dyq(ij)568 C dyq(ip1jm+ij)=ap s*dyq(ip1jm+ij)570 C dyq(ij)=appn*dyq(ij) 571 C dyq(ip1jm+ij)=apps*dyq(ip1jm+ij) 569 572 C ENDDO 570 573 C
Note: See TracChangeset
for help on using the changeset viewer.