Changeset 108 for trunk/libf/dyn3dpar/top_bound_p.F

Timestamp:

Apr 12, 2011, 11:16:02 AM (14 years ago)

Author:

slebonnois

Message:

Sebastien Lebonnois: sponge layer et dissip horizontale.

File:

• trunk/libf/dyn3dpar/top_bound_p.F (modified) (9 diffs)

trunk/libf/dyn3dpar/top_bound_p.F

@@ -64 +64 @@
              rdamp(llm-3)=tau_top_bound/8.
          else if (iflag_top_bound == 2) then
-! couce eponge dans toutes les couches de pression plus faible que
+! couche eponge dans toutes les couches de pression plus faible que
 ! 100 fois la pression de la derniere couche
              rdamp(:)=tau_top_bound
@@ -70 +70 @@
          endif
          first=.false.
-         print*,'TOP_BOUND rdamp=',rdamp
+         print*,'TOP_BOUND mode',mode_top_bound
+         print*,'Coeffs pour la couche eponge a l equateur'
+         print*,'p (Pa)  z(km)  tau (s)   dt*rdamp'
+         do l=1,llm
+           if (rdamp(l).ne.0.) then
+            zkm        = phi(iip1/2,jjp1/2,l)/(1000*g)
+          print*,presnivs(l),zkm,
+     .          1./rdamp(l),
+     .          dt*rdamp(l)
+           endif
+         enddo
 c$OMP END MASTER
 c$OMP BARRIER
@@ -77 +87 @@
 
       CALL massbar_p(masse,massebx,masseby)
+
+c   mode = 0 : pas de sponge
+c   mode = 1 : u et v -> 0
+c   mode = 2 : u et v -> moyenne zonale
+c   mode = 3 : u, v et h -> moyenne zonale
+
+C POUR V
+
 C  CALCUL DES CHAMPS EN MOYENNE ZONALE:
 
@@ -84 +102 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      do l=1,llm
+      if (mode_top_bound.ge.2) then
+       do l=1,llm
         do j=jjb,jje
           zm=0.
@@ -97 +116 @@
           vzon(j,l)=vzon(j,l)/zm
         enddo
-      enddo
+       enddo
+      else
+       do l=1,llm
+        do j=jjb,jje
+          vzon(j,l)=0.
+        enddo
+       enddo
+      endif
 c$OMP END DO NOWAIT   
 
-c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      do l=1,llm
+C   AMORTISSEMENTS LINEAIRES:
+
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
+      if (mode_top_bound.ge.1) then
+       do l=1,llm
         do j=jjb,jje
           do i=1,iip1
-            dv(i,j,l)=dv(i,j,l)-rdamp(l)*(vcov(i,j,l)-vzon(j,l))
-          enddo
-        enddo
-      enddo
-c$OMP END DO NOWAIT
+            dv(i,j,l)= -rdamp(l)*(vcov(i,j,l)-vzon(j,l))
+          enddo
+        enddo
+       enddo
+      endif
+c$OMP END DO NOWAIT
+
+C POUR U ET H
+
+C  CALCUL DES CHAMPS EN MOYENNE ZONALE:
 
       jjb=jj_begin
@@ -116 +150 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      do l=1,llm
+      if (mode_top_bound.ge.2) then
+       do l=1,llm
         do j=jjb,jje
           uzon(j,l)=0.
@@ -126 +161 @@
           uzon(j,l)=uzon(j,l)/zm
         enddo
-      enddo
+       enddo
+      else
+       do l=1,llm
+        do j=jjb,jje
+          uzon(j,l)=0.
+        enddo
+       enddo
+      endif
 c$OMP END DO NOWAIT
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
-      do l=1,llm
+      if (mode_top_bound.ge.3) then
+       do l=1,llm
         do j=jjb,jje
           zm=0.
@@ -140 +183 @@
           tzon(j,l)=tzon(j,l)/zm
         enddo
-      enddo
+       enddo
+      endif
 c$OMP END DO NOWAIT
 
@@ -146 +190 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      do l=1,llm
+      if (mode_top_bound.ge.1) then
+       do l=1,llm
         do j=jjb,jje
           do i=1,iip1
-            du(i,j,l)=du(i,j,l)
-     s               -rdamp(l)*(ucov(i,j,l)-cu(i,j)*uzon(j,l))
-            dh(i,j,l)=dh(i,j,l)-rdamp(l)*(teta(i,j,l)-tzon(j,l))
-          enddo
-       enddo
-      enddo
+            du(i,j,l)= -rdamp(l)*(ucov(i,j,l)-cu(i,j)*uzon(j,l))
+          enddo
+        enddo
+       enddo
+      endif
+c$OMP END DO NOWAIT
+      
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
+      if (mode_top_bound.ge.3) then
+       do l=1,llm
+        do j=jjb,jje
+          do i=1,iip1
+            dh(i,j,l)= -rdamp(l)*(teta(i,j,l)-tzon(j,l))
+          enddo
+        enddo
+       enddo
+      endif
 c$OMP END DO NOWAIT