Changeset 469 for branches/GRISLIv3/SOURCES/tab-litho-0.3.f90
- Timestamp:
- 01/26/24 17:34:15 (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/GRISLIv3/SOURCES/tab-litho-0.3.f90
r446 r469 88 88 WE(:,:) = WE(:,:)/(som*ROM*g) 89 89 90 91 92 ! ATTENTION DESSOUS VALABLE SEULEMENT EN 40 km93 94 ! la partie ci-dessous doit etre mise dans une routine specifique Philippe95 ! pour faire exactement la meme isostasie que Philippe96 !if ((icouple.eq.2).or.(icouple.eq.4)) then97 98 ! LBLOC=1099 100 ! do I=-LBLOC,LBLOC101 ! do J=-LBLOC,LBLOC102 ! We(i,j)=0. ! mise a 0 partout103 ! end do104 ! end do105 ! open(num_phil,file='../INPUT/iso-philippe')106 107 ! read(num_phil,*)108 ! read(num_phil,*)109 ! read(num_phil,*)110 ! read(num_phil,*)111 ! do J=-10,10112 ! read(num_phil,*)(We(i,j),i=-10,10)113 ! end do114 ! close(num_phil)115 ! do i=-LBLOC,LBLOC116 ! do j=-LBLOC,LBLOC117 ! we(i,j)=we(i,j)/11030/g/rom118 ! end do119 ! end do120 !endif121 122 90 return 123 91 end subroutine tab_litho
Note: See TracChangeset
for help on using the changeset viewer.