- Timestamp:
- Jul 24, 2024, 4:23:34 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90
r5116 r5117 539 539 do j=2,jmp1-1 540 540 PRINT*,'avant if ',cos(rlatu(j)),coslat0 541 if(cos(rlatu(j))<coslat0) THEN541 IF (cos(rlatu(j))<coslat0) THEN 542 542 ! nb de pts affectes par le filtrage de part et d'autre du pt 543 543 ifiltre=(coslat0/cos(rlatu(j))-1.)/2. … … 548 548 wwf(ifiltre+1)=(coslat0/cos(rlatu(j))-1.)/2.-ifiltre 549 549 do i=1,imp1-1 550 if(masque(i,j)>0.9) THEN550 IF (masque(i,j)>0.9) THEN 551 551 ssz=phis(i,j) 552 552 do ifi=1,ifiltre+1 553 553 ii=i+ifi 554 if(ii>imp1-1) ii=ii-imp1+1554 IF (ii>imp1-1) ii=ii-imp1+1 555 555 ssz=ssz+wwf(ifi)*phis(ii,j) 556 556 ii=i-ifi 557 if(ii<1) ii=ii+imp1-1557 IF (ii<1) ii=ii+imp1-1 558 558 ssz=ssz+wwf(ifi)*phis(ii,j) 559 559 enddo
Note: See TracChangeset
for help on using the changeset viewer.