Changeset 154

Timestamp:

11/08/17 16:29:28 (6 years ago)

Author:

aquiquet

Message:

flottab: bug in marais detection

Location:

trunk/SOURCES

Files:

• flottab2-0.7.f90 (modified) (6 diffs)
• spinup_mod.f90 (modified) (1 diff)

trunk/SOURCES/flottab2-0.7.f90

@@ -1335 +1335 @@
 integer,dimension(0:n_ta_max)     :: compt_marais        !< contient les equivalence entre les taches
 integer,dimension(0:n_ta_max)     :: nb_pts_marais !< indique le nombre de points par tache
-logical,dimension(0:n_ta_max)     :: marais      !< T si iceberg, F si calotte posee
+logical,dimension(0:n_ta_max)     :: marais      !< T si flottants entoure de poses, F sinon
  
 ! 1-initialisation
@@ -1381 +1381 @@
            
            !si un des voisins n'est pas glace alors la tache n'est pas un marais
-           do pm=-1,1,2
-              if ( (ice(i+pm,j).eq.0) .or. (ice(i,j+pm).eq.0) ) then
-                 marais(label)=.false.
-              endif
-           enddo
+           if ( (ice(i+1,j).eq.0) .or. (ice(i,j+1).eq.0) .or. (ice(i-1,j).eq.0) .or. (ice(i,j-1).eq.0) ) then
+              marais(label)=.false.
+           endif
 
            ! si 2 taches differentes sont dans le masque, il faut les identifier dans compt
@@ -1404 +1402 @@
            
            !si un des voisins n'est pas glace alors la tache n'est pas un marais
-           do pm=-1,1,2
-              if ( (ice(i+pm,j).eq.0) .or. (ice(i,j+pm).eq.0) ) then
-                 marais(label_max)=.false.
-              endif
-           enddo
+           if ( (ice(i+1,j).eq.0) .or. (ice(i,j+1).eq.0) .or. (ice(i-1,j).eq.0) .or. (ice(i,j-1).eq.0) ) then
+              marais(label_max)=.false.
+           endif
 
            
@@ -1424 +1420 @@
 enddo
 
-
+ marais(0)=.false.
 
 ! On reorganise compt en ecrivant le numero de la tache fondamentale
@@ -1439 +1435 @@
 
 !$OMP PARALLEL
-!$OMP DO REDUCTION(+:nb_pts_tache)
+!$OMP DO REDUCTION(+:nb_pts_marais)
 do j=1,ny
    do i=1,nx
@@ -1457 +1453 @@
 enddo
 
+debug_3D(:,:,64)=real(table_out_marais(:,:))
+
 end subroutine determin_marais
 

trunk/SOURCES/spinup_mod.f90

@@ -606 +606 @@
        end do
     end do
-    debug_3D(:,:,63)=ux_deformation(:,:)
-    debug_3D(:,:,64)=uy_deformation(:,:)
+    !debug_3D(:,:,63)=ux_deformation(:,:)
+    !debug_3D(:,:,64)=uy_deformation(:,:)
 
     if (itracebug.eq.1)  call tracebug(' fin de calc_ubar_def ')