Changeset 37
- Timestamp:
- 01/21/16 17:24:56 (8 years ago)
- Location:
- trunk/SOURCES
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/3D-physique-gen_mod.f90
r11 r37 192 192 real :: toblim !< 193 193 real :: moteurmax !< 194 real :: betamax !< 194 real :: betamax !< (Pa) frottement maxi sous les streams 195 195 real :: deltaacc !< 196 196 real :: pvimin !< valeur de pvi pour les noeuds fictifs -
trunk/SOURCES/dragging_hwat_sedim_mod.f90
r14 r37 26 26 27 27 real :: tostick ! pour la glace posee 28 real :: tobmax ! tobmax : (Pa) frottement maxi sous les streams29 28 real :: tob_ile ! pour les iles 30 29 real :: cry_lim=50. ! courbure limite pour le suivi des fleuves … … 63 62 implicit none 64 63 65 namelist/drag_hwat_sedim/file_sedim,hwatstream,cf, tobmax,toblim,seuil_hwater,seuil_sedim,seuil_neff,coef_sedim,coef_gz,coef_ile64 namelist/drag_hwat_sedim/file_sedim,hwatstream,cf,betamax,toblim,seuil_hwater,seuil_sedim,seuil_neff,coef_sedim,coef_gz,coef_ile 66 65 67 66 if (itracebug.eq.1) write(num_tracebug,*)' dragging avec hwatermax' 68 67 68 iter_beta=0 69 69 70 70 ! formats pour les ecritures dans 42 … … 83 83 write(num_rep_42,*) 'hwatstream = ', hwatstream 84 84 write(num_rep_42,*) 'cf = ', cf 85 write(num_rep_42,*) ' tobmax = ', tobmax85 write(num_rep_42,*) 'betamax = ', betamax 86 86 write(num_rep_42,*) 'toblim = ', toblim 87 87 write(num_rep_42,*) 'seuil_hwater = ', seuil_hwater … … 97 97 write(num_rep_42,428) '! cf coefficient de la loi de frottement fonction Neff' 98 98 write(num_rep_42,428) '! seulement pour les points cotiers' 99 write(num_rep_42,428) '! tobmax : (Pa) frottement maxi sous les streams '99 write(num_rep_42,428) '! betamax : (Pa) frottement maxi sous les streams ' 100 100 write(num_rep_42,428) '! toblim : (Pa) pour les iles ' 101 101 write(num_rep_42,428) '! seuil_hwater (m) : seuil hwater pour avoir glissement sur zone sediment' … … 108 108 109 109 tostick=1.e5 ! valeurs pour les points non flgzmx 110 tob_ile= tobmax/2.110 tob_ile=betamax/2. 111 111 moteurmax=toblim 112 112 … … 129 129 130 130 ! lecture du fichier sediment : 131 open(20,file=TRIM(DIRNAMEINP)//trim(file_sedim)) ! 'sediment_ij_hemin40.dat') 132 do j=1,ny 133 do i=1,nx 134 read(20,*),toto,toto,h_sedim(i,j) 135 enddo 136 enddo 137 h_sedim(:,:)=max(h_sedim(:,:),0.) ! pas d'epaisseurs de sediments negatives 138 131 !open(20,file=TRIM(DIRNAMEINP)//trim(file_sedim)) ! 'sediment_ij_hemin40.dat') 132 !do j=1,ny 133 ! do i=1,nx 134 ! read(20,*),toto,toto,h_sedim(i,j) 135 ! enddo 136 !enddo 137 !h_sedim(:,:)=max(h_sedim(:,:),0.) ! pas d'epaisseurs de sediments negatives 138 h_sedim(:,:)=10. 139 139 140 return 140 141 end subroutine init_dragging
Note: See TracChangeset
for help on using the changeset viewer.