Changeset 2429 for LMDZ5/trunk/libf/dynphy_lonlat/mod_interface_dyn_phys.F90

Timestamp:

Jan 27, 2016, 1:43:09 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Correction on the calculation of the surface of the grid at the poles (problem was introduced
in r2222).
Due to the different distribution of OMP tasks in the dynamics and the physics, had to
introduce 2 new logical variables, is_pole_north_phy and is_pole_south_phy, and so decided
to rename the old is_north_pole/is_south_pole to is_north_pole_dyn/is_south_pole_dyn to
stay coherent and, hopefully, clear things up a bit.
LF

File:

• LMDZ5/trunk/libf/dynphy_lonlat/mod_interface_dyn_phys.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/dynphy_lonlat/mod_interface_dyn_phys.F90

@@ -22 +22 @@
     
     k=1
-    IF (is_north_pole) THEN
+    IF (is_north_pole_dyn) THEN
       index_i(k)=1
       index_j(k)=1
@@ -28 +28 @@
     ELSE
       DO i=ii_begin,iim
-        index_i(k)=i
-        index_j(k)=jj_begin
-        k=k+1
-       ENDDO
+        index_i(k)=i
+        index_j(k)=jj_begin
+        k=k+1
+      ENDDO
     ENDIF
     
     DO j=jj_begin+1,jj_end-1
       DO i=1,iim
-        index_i(k)=i
-        index_j(k)=j
-        k=k+1
+        index_i(k)=i
+        index_j(k)=j
+        k=k+1
       ENDDO
     ENDDO
     
-    IF (is_south_pole) THEN
+    IF (is_south_pole_dyn) THEN
       index_i(k)=1
       index_j(k)=jj_end
     ELSE
       DO i=1,ii_end
-        index_i(k)=i
-        index_j(k)=jj_end
-        k=k+1
+        index_i(k)=i
+        index_j(k)=jj_end
+        k=k+1
        ENDDO
     ENDIF