Changeset 3213 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 12, 2024, 3:05:15 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:
Correction of a bug in "writediagsoil.F90": the case of using the 1D model with parallelization was not anticipated so that the "diagsoil.nc" file was filled with meaningless data.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/writediagsoil.F90 (modified) (6 diffs)

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

@@ -206 +206 @@
 if (dimpx.eq.3) then ! Case of a 3D variable
   ! A. Recast data along 'dynamics' grid
+  if (klon_glo>1) then ! General case
 #ifdef CPP_PARA
   ! gather field on a "global" (without redundant longitude) array
@@ -219 +220 @@
 !$OMP BARRIER
 #else
-  if (klon_glo>1) then ! General case
    do l=1,nsoilmx
     ! handle the poles
@@ -235 +235 @@
     enddo
    enddo
+#endif
   else ! 1D model case
    data3_1d(1,1:nsoilmx)=px(1,1:nsoilmx)
   endif
-#endif
   
   ! B. Write (append) the variable to the NetCDF file
@@ -292 +292 @@
 
   ! A. Recast data along 'dynamics' grid
+  if (klon_glo>1) then ! General case
 #ifdef CPP_PARA
   ! gather field on a "global" (without redundant longitude) array
@@ -306 +307 @@
 !$OMP BARRIER
 #else
-  if (klon_glo>1) then ! general case
     ! handle the poles
     do i=1,nbp_lon+1
@@ -320 +320 @@
       data2(nbp_lon+1,j)=data2(1,j) ! extra (modulo) longitude
     enddo
+#endif
   else ! 1D model case
     data2_1d=px(1,1)
   endif
-#endif
 
   ! B. Write (append) the variable to the NetCDF file