Changeset 8 for trunk/libf/dyn3dpar/caldyn_p.F

Timestamp:

Nov 2, 2010, 12:38:23 PM (14 years ago)

Author:

emillour

Message:

Debut de mise a jour de la dynamique parallele par rapport aux modifs dans la partie sequentielle.

Mais NON TESTE , car pas (encore) possibilite de compiler et faire tourner cas simple (type newtonien sans physique).

Voir commit_v8.log pour les details.

Ehouarn

File:

• trunk/libf/dyn3dpar/caldyn_p.F (modified) (5 diffs)

trunk/libf/dyn3dpar/caldyn_p.F

@@ -8 +8 @@
 
       SUBROUTINE caldyn_p
-     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
+     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,tsurpk,phis ,
      $  phi,conser,du,dv,dteta,dp,w,pbaru,pbarv,time )
       USE parallel
@@ -45 +45 @@
       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
       REAL ps(ip1jmp1),phis(ip1jmp1)
-      REAL pk(iip1,jjp1,llm),pkf(ip1jmp1,llm)
+      REAL pk(ip1jmp1,llm),pkf(ip1jmp1,llm)
+      REAL tsurpk(ip1jmp1,llm)
       REAL,SAVE :: vcont(ip1jm,llm),ucont(ip1jmp1,llm)
       REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
@@ -66 +67 @@
       REAL,SAVE :: massebxy(ip1jm,llm)
       REAL,SAVE :: convm(ip1jmp1,llm)
+!      REAL,SAVE :: temp(ip1jmp1,llm)
       INTEGER   ij,l,ijb,ije,ierr
 
@@ -129 +131 @@
       CALL enercin_p ( vcov   , ucov  , vcont     , ucont  , ecin  )
       CALL bernoui_p ( ip1jmp1, llm   , phi       , ecin   , bern  )
-      CALL dudv2_p   ( teta   , pkf   , bern      , du     , dv    )
+      CALL dudv2_p   ( tsurpk , pkf   , bern      , du     , dv    )
 
 #ifdef DEBUG_IO
@@ -184 +186 @@
 c ym ---> exige communication collective ( aussi dans advect)
         CALL sortvarc
-     $ ( itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov )
+     $ (itau,ucov,tsurpk,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov)
 
       ENDIF