Changeset 3169

Timestamp:

Jan 3, 2024, 12:41:29 PM (11 months ago)

Author:

emillour

Message:

Mars PCM:
Minor fix for XIOS output of a scalar in write_output_mod.F90.
Also added more helpful error messages in xios_output_mod.F90.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/write_output_mod.F90 (modified) (1 diff)
• libf/phymars/xios_output_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4426 +4426 @@
 == 30/12/2023 == EM
 Fix typo in previous commit
+
+== 03/01/2024 == EM
+Minor fix for XIOS output of a scalar in write_output_mod.F90.
+Also added more helpful error messages in xios_output_mod.F90.

trunk/LMDZ.MARS/libf/phymars/write_output_mod.F90

@@ -38 +38 @@
 #ifdef CPP_XIOS
     ! only send the field to xios if the user asked for it
-    if (xios_is_active_field(field_name)) call send_xios_field(field_name,(/field/))
+    if (xios_is_active_field(field_name)) call send_xios_field(field_name,field)
 #endif
 

trunk/LMDZ.MARS/libf/phymars/xios_output_mod.F90

@@ -243 +243 @@
     REAL,DIMENSION(klon_mpi) :: buffer_omp
     REAL :: Field2d(nbp_lon,jj_nb)
+    
+    CHARACTER(len=128) :: msg
 
     IF (prt_level >= 10) WRITE(lunout,*)'Begin histrwrite2d_xios ',trim(field_name)
-    IF (SIZE(field)/=klon) CALL abort_physic('iophy::histwrite2d_xios','Field first DIMENSION not equal to klon',1)
-    
+    IF (SIZE(field)/=klon) THEN
+      WRITE(msg,*) "Pb with field "//trim(field_name)//&
+                   " : Field first DIMENSION not equal to klon"
+      CALL abort_physic('iophy::histwrite2d_xios',trim(msg),1)
+    ENDIF
     CALL Gather_omp(field,buffer_omp)    
 !$OMP MASTER
@@ -276 +281 @@
     INTEGER :: ip, n, nlev
 
+    CHARACTER(len=128) :: msg
+
   IF (prt_level >= 10) write(lunout,*)'Begin histrwrite3d_xios ',trim(field_name)
 
-    !Et on.... écrit
-    IF (SIZE(field,1)/=klon) CALL abort_physic('iophy::histwrite3d','Field first DIMENSION not equal to klon',1)
+    !And we write...
+    IF (SIZE(field,1)/=klon) THEN
+      WRITE(msg,*) "Pb with field "//trim(field_name)//&
+                   " : Field first DIMENSION not equal to klon"
+      CALL abort_physic('iophy::histwrite3d',trim(msg),1)
+    ENDIF
+    
     nlev=SIZE(field,2)