Changeset 798 for trunk/LMDZ.MARS/libf/dyn3d/gcm.F

Timestamp:

Oct 1, 2012, 9:33:54 AM (12 years ago)

Author:

tnavarro

Message:

modification that allows to have "1+1=2" with dynamical core only.

File:

• trunk/LMDZ.MARS/libf/dyn3d/gcm.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/dyn3d/gcm.F

@@ -189 +189 @@
 c
 
-      CALL pression ( ip1jmp1, ap, bp, ps, p       )
 
       call dump2d(iip1,jjp1,ps,'PRESSION SURFACE')
-      CALL exner_hyb( ip1jmp1, ps, p,beta, pks, pk, pkf )
 
 c
@@ -249 +247 @@
    1  CONTINUE
 c
+c TN 09/2012. To ensure "1+1=2" in dynamical core :
+c update atmospheric pressure IN the main loop
+      CALL pression ( ip1jmp1, ap, bp, ps, p       )
+      CALL exner_hyb( ip1jmp1, ps, p,beta, pks, pk, pkf )
+
       IF( MOD( itau, 10* day_step ).EQ.0 )  THEN
         CALL test_period ( ucov,vcov,teta,q,p,phis )