Changeset 1385 for LMDZ4/trunk


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/atm2geo.F90

    r1146 r1385  
    3434! Value at North Pole 
    3535  IF (is_north_pole) THEN
    36      pxx(:, 1) = pxx(1,1)
    37      pyy(:, 1) = pyy(1,1)
     36     pxx(:, 1) = - pte (1, 1)
     37     pyy(:, 1) = - ptn (1, 1)
    3838     pzz(:, 1) = pzz(1,1)
    3939  ENDIF
  • LMDZ4/trunk/libf/phylmd/geo2atm.F90

    r1146 r1385  
    3939  ! Value at North Pole
    4040  IF (is_north_pole) THEN
    41      pu(:, 1) = pu(1, 1)
    42      pv(:, 1) = pv(1, 1)
    43      pr(:, 1) = pr(1, 1)
     41     pu(:, 1) = -px (1,1)
     42     pv(:, 1) = -py (1,1)
     43     pr(:, 1) = 0.0
    4444  ENDIF
    4545 
    4646  ! Value at South Pole     
    4747  IF (is_south_pole) THEN
    48      pu(:,jm) = pu(1,jm)
    49      pv(:,jm) = pv(1,jm)
    50      pr(:,jm) = pr(1,jm)
     48     pu(:,jm) = -px (1,jm)
     49     pv(:,jm) = -py (1,jm)
     50     pr(:,jm) = 0.0
    5151  ENDIF
    5252 
Note: See TracChangeset for help on using the changeset viewer.