Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/atm2geo.F

Timestamp:

Aug 17, 2006, 5:41:51 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/atm2geo.F (modified) (3 diffs)

LMDZ4/branches/V3_test/libf/phylmd/atm2geo.F

@@ -4 +4 @@
 C
       SUBROUTINE atm2geo ( im, jm, pte, ptn, plon, plat, pxx, pyy, pzz )
+      USE dimphy
+      IMPLICIT NONE
+      include 'dimensions90.h'
 cc
 cc Change wind local atmospheric coordinates to
 cc geocentric
 cc
-c$$$      INCLUDE 'param.h'
+cxxx      INCLUDE 'param.h'
 c
       INTEGER, INTENT (in)              :: im, jm
@@ -22 +25 @@
 c
       LOGICAL, SAVE :: linit = .FALSE.
+c$OMP THREADPRIVATE(linit)
 c
-c$$$      IF ( .NOT. linit ) THEN 
+cxxx      IF ( .NOT. linit ) THEN 
           zsinlon = SIN (rad * plon)
           zcoslon = COS (rad * plon)
@@ -29 +33 @@
           zcoslat = COS (rad * plat)
           linit = .TRUE.
-c$$$      ENDIF 
+cxxx      ENDIF 
 c
       pxx = - zsinlon * pte - zsinlat * zcoslon * ptn
       pyy =   zcoslon * pte - zsinlat * zsinlon * ptn
       pzz =   zcoslat * ptn
+
 c
 c Value at North Pole
-      pxx ( :,  1) = - ptn ( 1, 1)
-      pyy ( :,  1) = - pte ( 1, 1)
-      pzz ( :,  1) = 0.0
+      IF (jjphy_begin==1) THEN
+        pxx ( :,  1) = - ptn ( 1, 1)
+        pyy ( :,  1) = - pte ( 1, 1)
+        pzz ( :,  1) = 0.0
+      ENDIF
 c Value at South Pole
-      pxx ( :, jm) = + ptn ( 1, jm)
-      pyy ( :, jm) = + pte ( 1, jm)
-      pzz ( :, jm) = 0.0
-c
+      
+      IF (jjphy_end==jjm+1) THEN
+        pxx ( :, jm) = + ptn ( 1, jm)
+        pyy ( :, jm) = + pte ( 1, jm)
+        pzz ( :, jm) = 0.0
+      ENDIF
+ 
       RETURN 
       END SUBROUTINE atm2geo