Changeset 129 for trunk

Timestamp:

09/29/17 15:34:54 (7 years ago)

Author:

dumas

Message:

calving_frange : edge test with H | flottab2 : bug in openmp | steps_time_loop : iter_visco not imposed at 10.

Location:

trunk/SOURCES

Files:

• calving_frange.f90 (modified) (2 diffs)
• flottab2-0.7.f90 (modified) (3 diffs)
• steps_time_loop.f90 (modified) (1 diff)

trunk/SOURCES/calving_frange.f90

@@ -97 +97 @@
      *(hcoup_abysses-hcoup_plateau)+hcoup_plateau               &
      , hcoup_plateau), hcoup_abysses )
-
+    
 if (meth_hcoup.eq.1) then
    Hcoup(:,:)=coefbmshelf*Hcoup(:,:)
@@ -149 +149 @@
 
 !ifint:    if((front(i,j).gt.0).and.(front(i,j).lt.4)) then
-        ifint:    if(front(i,j).lt.4) then
+!cdc pb avec front    ifint:    if((front(i,j).lt.4).or.((front(i-1,j)+front(i+1,j)+front(i,j-1)+front(i,j+1)).lt.16)) then
+        ifint:    if((H(i-1,j).lt.2.).or.(H(i+1,j).lt.2.).or.(H(i,j-1).lt.2.).or.(H(i,j+1).lt.2)) then ! si on est au bord avec test sur H
 ! ifint: le point doit avoir au - 1 voisin mais ne pas etre entouré de glace
 ! ce qui evite la formation des polynies dans les shelfs

trunk/SOURCES/flottab2-0.7.f90

@@ -1270 +1270 @@
 end do
 !$OMP END DO 
-
+!$OMP BARRIER
 ! traitement des coins
 
@@ -1279 +1279 @@
 
 ! Les points à plus d'un point de grille du bord sont front=-1
+!$OMP WORKSHARE
 nofront(:,:)=0
+!$OMP END WORKSHARE
+!$OMP BARRIER
 
 !$OMP DO
@@ -1290 +1293 @@
 enddo
 !$OMP END DO
+!$OMP BARRIER
+
+!$OMP WORKSHARE
 where (nofront(:,:).eq.-1)
         front(:,:)=-1
 endwhere
-
+!$OMP END WORKSHARE
 !$OMP END PARALLEL              
 

trunk/SOURCES/steps_time_loop.f90

@@ -444 +444 @@
         end if
  
-        iter_visco= 10                                 ! warning, test sur l'impact iteration dragging et visco 
+!cdc        iter_visco= 10                                 ! warning, test sur l'impact iteration dragging et visco 
                                                        ! Cat 18 janv 2013