Changeset 3488 for LMDZ6/trunk

Timestamp:

Apr 26, 2019, 1:41:36 PM (6 years ago)

Author:

Laurent Fairhead

Message:

Some modification needed for the outputs for the unstructured grid case
AC

File:

• LMDZ6/trunk/libf/phylmd/iophy.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/iophy.F90

@@ -1067 +1067 @@
             ENDIF 
           ELSE IF (grid_type==unstructured) THEN
-            CALL xios_send_field(var%name, buffer_omp)
+            IF (SIZE(field) == klon) then
+              CALL xios_send_field(var%name, buffer_omp)
+            ELSE
+               CALL xios_send_field(var%name, field)
+            ENDIF
+
           ENDIF
           IF (prt_level >= 10) THEN
@@ -1094 +1099 @@
                     ENDIF
                   ELSE IF (grid_type==unstructured) THEN
-                    CALL xios_send_field(var%name, buffer_omp)
+                    IF (SIZE(field) == klon) then
+                      CALL xios_send_field(var%name, buffer_omp)
+                    ELSE
+                      CALL xios_send_field(var%name, field)
+                    ENDIF
                   ENDIF
 
@@ -1248 +1257 @@
             ENDIF
           ELSE IF (grid_type==unstructured) THEN
-            CALL xios_send_field(var%name, buffer_omp(:,1:nlevx))
+            IF (SIZE(field,1) == klon) then
+               CALL xios_send_field(var%name, buffer_omp(:,1:nlevx))
+            ELSE
+               CALL xios_send_field(var%name, field)
+            ENDIF
         ENDIF
 
@@ -1274 +1287 @@
                   ENDIF
                 ELSE IF (grid_type==unstructured) THEN
-                  CALL xios_send_field(var%name, buffer_omp(:,1:nlevx))
+                  IF (SIZE(field,1) == klon) then
+                     CALL xios_send_field(var%name, buffer_omp(:,1:nlevx))
+                  ELSE
+                     CALL xios_send_field(var%name, field)
+                  ENDIF
                 ENDIF
 
@@ -1362 +1379 @@
     !Et sinon on.... écrit
     IF (SIZE(field)/=klon .AND. SIZE(field)/=klev .AND. SIZE(field)/=klev+1) CALL abort_physic('iophy::histwrite2d_xios','Field first DIMENSION not equal to klon/klev',1)
-    
     IF (SIZE(field) == klev .OR. SIZE(field) == klev+1) then
 !$OMP MASTER
+
         CALL xios_send_field(field_name,field)
 !$OMP END MASTER