Changeset 2350

Timestamp:

Aug 25, 2015, 1:40:19 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Corrections for 1D case where nbp_lon=nbp_lat=1.
EM + MPL

Location:

LMDZ5/trunk/libf/phylmd

Files:

• aaam_bud.F90 (modified) (2 diffs)
• iophy.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/aaam_bud.F90

@@ -6 +6 @@
 
   USE dimphy
-  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
+  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo
   IMPLICIT NONE
   ! ======================================================================
@@ -127 +127 @@
   hadley = 1.E18
   hadday = 1.E18*24.*3600.
-  dlat = xpi/real(nbp_lat-1)
+  IF(klon_glo.EQ.1) THEN
+    dlat = xpi
+  ELSE
+    dlat = xpi/real(nbp_lat-1)
+  ENDIF
   dlon = 2.*xpi/real(nbp_lon)
 

LMDZ5/trunk/libf/phylmd/iophy.F90

@@ -76 +76 @@
     
 !$OMP MASTER  
-    ALLOCATE(io_lat(nbp_lat-1/(nbp_lon*(nbp_lat-1))))
-    io_lat(1)=rlat_glo(1)
-    io_lat(nbp_lat-1/(nbp_lon*(nbp_lat-1)))=rlat_glo(klon_glo)
-    IF ((nbp_lon*nbp_lat) > 1) then
+    ALLOCATE(io_lat(nbp_lat))
+    IF (klon_glo == 1) THEN
+      io_lat(1)=rlat_glo(1)
+    ELSE
+      io_lat(1)=rlat_glo(1)
+      io_lat(nbp_lat)=rlat_glo(klon_glo)
       DO i=2,nbp_lat-1
         io_lat(i)=rlat_glo(2+(i-2)*nbp_lon)