Changeset 3092 for trunk/LMDZ.MARS/libf/phymars/writediagfi.F

Timestamp:

Oct 18, 2023, 3:35:24 PM (12 months ago)

Author:

jbclement

Message:

Marc PCM:

• Correction of a bug in "writediagfi.F": the case of using the 1D model with parallelization was not anticipated so that the "diagfi.nc" file was filled with NaNf?;
• Addition of the file "start1D.txt" as an example in the directory deftank/;
• Some "cosmetic" modifications in "improvedclouds_mod.F", "write_output_mod.F90" and "testphys1d.F90".

JBC

File:

• trunk/LMDZ.MARS/libf/phymars/writediagfi.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/writediagfi.F

@@ -312 +312 @@
         if (dim.eq.3) then
 
+          IF (klon_glo>1) THEN ! General case
 #ifdef CPP_PARA
           ! Gather field on a "global" (without redundant longitude) array
@@ -327 +328 @@
 !         Passage variable physique -->  variable dynamique
 !         recast (copy) variable from physics grid to dynamics grid
-          IF (klon_glo>1) THEN ! General case
            DO l=1,nbp_lev
              DO i=1,nbp_lon+1
@@ -341 +341 @@
              ENDDO
            ENDDO
+#endif
           ELSE ! 1D model case
            dx3_1d(1,1:nbp_lev)=px(1,1:nbp_lev)
           ENDIF
-#endif
 !         Ecriture du champs
 
@@ -416 +416 @@
         else if (dim.eq.2) then
 
+          IF (klon_glo>1) THEN ! General case
 #ifdef CPP_PARA
           ! Gather field on a "global" (without redundant longitude) array
@@ -433 +434 @@
 !         Passage variable physique -->  physique dynamique
 !         recast (copy) variable from physics grid to dynamics grid
-          IF (klon_glo>1) THEN ! General case
              DO i=1,nbp_lon+1
                 dx2(i,1)=px(1,1)
@@ -445 +445 @@
                 dx2(nbp_lon+1,j)=dx2(1,j)
              ENDDO
+#endif
           ELSE ! 1D model case
             dx2_1d=px(1,1)
           ENDIF
-#endif
 
           if (is_master) then