Changeset 58


Ignore:
Timestamp:
04/26/16 14:28:27 (8 years ago)
Author:
aquiquet
Message:

Bug fix in flottab: surface below sea level after isostasy

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SOURCES/flottab2-0.7.f90

    r4 r58  
    594594!!$end if 
    595595 
    596 ! / test aurel 
    597 debug_3d(:,:,43)=table_out(:,:) 
    598  
    599 iceberg(:)=.true. 
    600 do j=1,ny 
    601    do i=1,nx 
    602       if(.not.flot(i,j)) then 
    603          iceberg(table_out(i,j))=.false. 
    604       end if 
    605    end do 
    606 end do 
    607  
    608 do j=1,ny 
    609    do i=1,nx 
    610  
    611       debug_3d(i,j,56)=nb_pts_tache(table_out(i,j)) 
    612       if (iceberg(table_out(i,j))) then 
    613          debug_3d(i,j,57)=1 
    614       else 
    615          debug_3d(i,j,57)=0        
    616       end if 
    617    end do 
    618 end do 
    619 ! fin test aurel / 
    620  
    621  
    622596!ice(:,:)=0   Attention, voir si ca marche toujours pour l'Antarctique et heminord ! 
    623597 
Note: See TracChangeset for help on using the changeset viewer.