Changeset 84

Timestamp:

07/13/16 09:59:57 (8 years ago)

Author:

aquiquet

Message:

Small bug correction in dragging modules for gfortran compatibility

Location:

trunk/SOURCES

Files:

• Greeneem_files/dragging_neem_mod.f90 (modified) (1 diff)
• dragging_neff_slope_mod.f90 (modified) (2 diffs)

trunk/SOURCES/Greeneem_files/dragging_neem_mod.f90

@@ -31 +31 @@
 
 implicit none
-logical,dimension(nx,ny) :: fleuvemx         !< fleuves sont les tableaux courants (dep. time)   
-logical,dimension(nx,ny) :: fleuvemy
+!logical,dimension(nx,ny) :: fleuvemx         !< fleuves sont les tableaux courants (dep. time)   
+!logical,dimension(nx,ny) :: fleuvemy
 logical,dimension(nx,ny) :: fleuve      
 logical,dimension(nx,ny) :: cote    

trunk/SOURCES/dragging_neff_slope_mod.f90

@@ -228 +228 @@
 
 
-!!$
-!!$fleuve_maj: do j=2,ny-1  
-!!$ifleuve:   do i=2,nx-1                     
-!!$
-!!$cote_detect :  if (cote(i,j)) then
-!!$         idep=i
-!!$         jdep=j
-!!$
-!!$         if (socle_cry(i,j).lt.0.) then    ! dans une vallee
-!!$            fleuve(i,j)=.true. 
-!!$         else
-!!$            cote(i,j)=.false.
-!!$            cycle ifleuve
-!!$         endif
-!!$
-!!$suit : do l=1,lmax         ! debut de la boucle de suivi, lmax longueur maxi des fleuves
-!!$           i_moins1=max(idep-1,2)
-!!$           j_moins1=max(jdep-1,1)
-!!$           i_plus1=min(idep+1,nx)
-!!$           j_plus1=min(jdep+1,ny)
-!!$           
-!!$! recherche du max en suivant le socle le plus profond
-!!$! * en excluant les points flottants
-!!$! * et ceux qui sont deja tagges fleuves
-!!$
-!!$           valmax=1000.
-!!$
-!!$           do jloc=j_moins1,j_plus1
-!!$              do iloc=i_moins1,i_plus1
-!!$ 
-!!$                 if ((B(iloc,jloc).lt.valmax)      & 
-!!$                      .and.(.not.flot(iloc,jloc))  &
-!!$                      .and.(.not.fleuve(iloc,jloc)).and.(socle_cry(iloc,jloc).lt.cry_lim)) then
-!!$                    imax=iloc
-!!$                    jmax=jloc
-!!$                    valmax=B(iloc,jloc)
-!!$                 endif
-!!$              end do
-!!$           end do
-!!$
-!!$         if ((hwater(imax,jmax).gt.hwatstream).and.(socle_cry(i,j).lt.cry_lim)) then
-!!$            fleuve(imax,jmax)=.true.
-!!$            idep=imax
-!!$            jdep=jmax
-!!$         else
-!!$            fleuve(imax,jmax)=.false.
-!!$            exit suit
-!!$         end if
-!!$ 
-!!$      end do suit
-!!$
-!!$   end if cote_detect
-!!$
-!!$end do ifleuve
-!!$end do fleuve_maj
-
 ! aurel, we add the neff threshold:
 where ((neff(:,:).le.seuil_neff).and.(.not.flot(:,:)).and.(H(:,:).gt.1.)) fleuve(:,:)=.true.
@@ -311 +255 @@
    do i=1,nx
       !if ((not(flot(i,j))).and.(hwater(i,j).gt.0.1)) slowssa(i,j)=.true.
-      if ((not(flot(i,j))).and.(ibase(i,j).ne.1).and.(H(i,j).gt.1.)) slowssa(i,j)=.true.
+      if ((.not.(flot(i,j))).and.(ibase(i,j).ne.1).and.(H(i,j).gt.1.)) slowssa(i,j)=.true.
    end do
 end do