Changeset 1385

Timestamp:

May 7, 2010, 9:15:50 AM (15 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

Location:

LMDZ4/trunk/libf/phylmd

Files:

• atm2geo.F90 (modified) (1 diff)
• geo2atm.F90 (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/atm2geo.F90

@@ -34 +34 @@
 ! Value at North Pole  
   IF (is_north_pole) THEN
-     pxx(:, 1) = pxx(1,1)
-     pyy(:, 1) = pyy(1,1)
+     pxx(:, 1) = - pte (1, 1)
+     pyy(:, 1) = - ptn (1, 1) 
      pzz(:, 1) = pzz(1,1)
   ENDIF

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