Changeset 3605 for LMDZ6/branches/Ocean_skin/libf/phylmd/geo2atm.F90

Timestamp:

Nov 21, 2019, 4:43:45 PM (5 years ago)

Author:

lguez

Message:

Merge revisions 3427:3600 of trunk into branch Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/geo2atm.F90 (modified) (3 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/geo2atm.F90

@@ -5 +5 @@
   USE dimphy
   USE mod_phys_lmdz_para
-
+  USE mod_grid_phy_lmdz, only: grid_type, unstructured, regular_lonlat
   IMPLICIT NONE
   INCLUDE 'YOMCST.h'
+  CHARACTER (len = 6)                :: clmodnam
+  CHARACTER (len = 20)               :: modname = 'geo2atm'
+  CHARACTER (len = 80)               :: abort_message
 
 ! Change wind coordinates from cartesian geocentric to local spherical
 ! NB! Fonctionne probablement uniquement en MPI seul (sans OpenMP)
 !
+! Geocentric :
+  ! axe x is eastward : crosses (0N,90E) point.
+  ! axe y  crosses (0N,180E) point.
+  ! axe z is 'up' : crosses north pole.
+!
+! NB! Aux poles, fonctionne probablement uniquement en MPI seul (sans OpenMP)
+
   INTEGER, INTENT (IN)                 :: im, jm
   REAL, DIMENSION (im,jm), INTENT(IN)  :: px, py, pz
@@ -17 +27 @@
   REAL, DIMENSION (im,jm), INTENT(OUT) :: pu, pv, pr
 
-  REAL :: rad
+  REAL :: rad,reps
 
 
   rad = rpi / 180.0E0
-  
+  reps = 1.0e-5
+
   pu(:,:) = &
        - px(:,:) * SIN(rad * plon(:,:)) &
@@ -36 +47 @@
        + pz(:,:) * SIN(rad * plat(:,:))
 
-  ! Value at North Pole
-  IF (is_north_pole_dyn) THEN
-     pu(:, 1) = -px (1,1)
-     pv(:, 1) = -py (1,1)
-     pr(:, 1) = 0.0
-  ENDIF
+  IF (grid_type==regular_lonlat) THEN
+    ! Value at North Pole
+    IF (is_north_pole_dyn) THEN
+       pu(:, 1) = -px (1,1)
+       pv(:, 1) = -py (1,1)
+       pr(:, 1) = 0.0
+    ENDIF
   
-  ! Value at South Pole     
-  IF (is_south_pole_dyn) THEN
-     pu(:,jm) = -px (1,jm)
-     pv(:,jm) = -py (1,jm)
-     pr(:,jm) = 0.0
-  ENDIF
+    ! Value at South Pole     
+    IF (is_south_pole_dyn) THEN
+       pu(:,jm) = -px (1,jm)
+       pv(:,jm) = -py (1,jm)
+       pr(:,jm) = 0.0
+    ENDIF
+
+  ELSE IF (grid_type==unstructured) THEN
+     ! Pole nord pour Dynamico
+     WHERE ( plat(:,:) >= 90.0-reps )
+        pu(:,:) =  py(:,:)
+        pv(:,:) = -px(:,:)
+        pr(:,:) = 0.0e0
+     END WHERE
+
+  ELSE
+     abort_message='Problem: unknown grid type'
+     CALL abort_physic(modname,abort_message,1)
+  END IF
+
+  
+  
   
 END SUBROUTINE geo2atm