Changeset 38 for branches

Timestamp:

01/21/16 17:26:45 (8 years ago)

Author:

dumas

Message:

Merge avec la r37 du trunk

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/3D-physique-gen_mod.f90 (modified) (1 diff)
• SOURCES/dragging_hwat_sedim_mod.f90 (modified) (6 diffs)

branches/iLoveclim/SOURCES/3D-physique-gen_mod.f90

@@ -192 +192 @@
   real ::  toblim                       !< 
   real ::  moteurmax                    !< 
-  real ::  betamax                      !< 
+  real ::  betamax                      !< (Pa) frottement maxi sous les streams
   real ::  deltaacc                     !< 
   real ::  pvimin                       !< valeur de pvi pour les noeuds fictifs

branches/iLoveclim/SOURCES/dragging_hwat_sedim_mod.f90

@@ -26 +26 @@
 
 real :: tostick   ! pour la glace posee
-real :: tobmax    ! tobmax : (Pa) frottement maxi sous les streams
 real :: tob_ile    ! pour les iles
 real :: cry_lim=50.  ! courbure limite pour le suivi des fleuves
@@ -63 +62 @@
 implicit none
  
-namelist/drag_hwat_sedim/file_sedim,hwatstream,cf,tobmax,toblim,seuil_hwater,seuil_sedim,seuil_neff,coef_sedim,coef_gz,coef_ile    
+namelist/drag_hwat_sedim/file_sedim,hwatstream,cf,betamax,toblim,seuil_hwater,seuil_sedim,seuil_neff,coef_sedim,coef_gz,coef_ile    
 
 if (itracebug.eq.1) write(num_tracebug,*)' dragging avec hwatermax'
 
+iter_beta=0
 
 ! formats pour les ecritures dans 42
@@ -83 +83 @@
 write(num_rep_42,*) 'hwatstream     = ', hwatstream 
 write(num_rep_42,*) 'cf             = ', cf
-write(num_rep_42,*) 'tobmax         = ', tobmax
+write(num_rep_42,*) 'betamax        = ', betamax
 write(num_rep_42,*) 'toblim         = ', toblim
 write(num_rep_42,*) 'seuil_hwater   = ', seuil_hwater
@@ -97 +97 @@
 write(num_rep_42,428) '! cf coefficient de la loi de frottement fonction Neff'
 write(num_rep_42,428) '! seulement pour les points cotiers'
-write(num_rep_42,428) '! tobmax : (Pa) frottement maxi sous les streams '
+write(num_rep_42,428) '! betamax : (Pa) frottement maxi sous les streams '
 write(num_rep_42,428) '! toblim : (Pa)  pour les iles '
 write(num_rep_42,428) '! seuil_hwater (m) : seuil hwater pour avoir glissement sur zone sediment'
@@ -108 +108 @@
 
 tostick=1.e5   ! valeurs pour les points non flgzmx
-tob_ile=tobmax/2.
+tob_ile=betamax/2.
 moteurmax=toblim
 
@@ -129 +129 @@
 
 ! lecture du fichier sediment :
-open(20,file=TRIM(DIRNAMEINP)//trim(file_sedim)) !   'sediment_ij_hemin40.dat')
-do j=1,ny
-   do i=1,nx
-      read(20,*),toto,toto,h_sedim(i,j)
-   enddo
-enddo
-h_sedim(:,:)=max(h_sedim(:,:),0.)  ! pas d'epaisseurs de sediments negatives
-   
+!open(20,file=TRIM(DIRNAMEINP)//trim(file_sedim)) !   'sediment_ij_hemin40.dat')
+!do j=1,ny
+!   do i=1,nx
+!      read(20,*),toto,toto,h_sedim(i,j)
+!   enddo
+!enddo
+!h_sedim(:,:)=max(h_sedim(:,:),0.)  ! pas d'epaisseurs de sediments negatives
+h_sedim(:,:)=10. 
+  
       return
       end subroutine init_dragging