Changeset 1385 for LMDZ4/trunk/libf/phylmd/geo2atm.F90

Timestamp:

May 7, 2010, 9:15:50 AM (14 years ago)

Author:

Laurent Fairhead

Message:

Correction of bug in wind component transformation at the poles

---

Correction d'un bug dans la transformation du vent aux pôles

OM

File:

• LMDZ4/trunk/libf/phylmd/geo2atm.F90 (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/geo2atm.F90

@@ -39 +39 @@
   ! Value at North Pole
   IF (is_north_pole) THEN
-     pu(:, 1) = pu(1, 1)
-     pv(:, 1) = pv(1, 1) 
-     pr(:, 1) = pr(1, 1)
+     pu(:, 1) = -px (1,1)
+     pv(:, 1) = -py (1,1)
+     pr(:, 1) = 0.0
   ENDIF
   
   ! Value at South Pole     
   IF (is_south_pole) THEN
-     pu(:,jm) = pu(1,jm)
-     pv(:,jm) = pv(1,jm)
-     pr(:,jm) = pr(1,jm)
+     pu(:,jm) = -px (1,jm)
+     pv(:,jm) = -py (1,jm)
+     pr(:,jm) = 0.0
   ENDIF