Changeset 2222 for LMDZ5/trunk/libf/dyn3d/gcm.F

Timestamp:

Mar 10, 2015, 10:58:04 AM (10 years ago)

Author:

Ehouarn Millour

Message:

Bug fix: Poles are single points on physics grid. True mesh area there is thus the sum of corresponding dynamics "polar meshes" areas.
Note that this also implies that some extra work is needed in physics to generate correct areas at "polar points" for outputs (especially for zoomed grids).
EM

File:

• LMDZ5/trunk/libf/dyn3d/gcm.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/gcm.F

@@ -435 +435 @@
          zcvfi(ngridmx) = cv(ip1jm-iim)
          CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
-         WRITE(lunout,*)
-     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
+         ! Poles are single points on physics grid
+         airefi(1)=sum(aire(1:iim))
+         airefi(ngridmx)=sum(aire(ip1jmp1-(iim+1)+1:ip1jmp1-1))
+!         WRITE(lunout,*)
+!     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
 ! Physics:
 #ifdef CPP_PHYS