Changeset 4161 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Apr 1, 2026, 3:17:47 PM (6 weeks ago)

Author:

jmauxion

Message:

Mars PCM:

• Fix improvedclouds_mod by adding a missing "else" condition on zdq.
• Add the possibility to output watercaptag in physiq_mod.
• Fix write_output to handle logical arrays with an unstructured grid.
• Update some of the utils (localtime/zrecast) to handle phisfi.

JM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• improvedclouds_mod.F90 (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)
• write_output_mod.F90 (modified) (1 diff)

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

@@ -443 +443 @@
         IF (spenttime.ne.0) then
           zdq=(dMicetot/spenttime)!*(ptimestep-spenttime)
+        ELSE
+          ! Initialization for spenttime=0
+          zdq=zpotcond(ig,l)*((ptimestep-spenttime)/ptimestep)
         ENDIF
         zdq=abs(zdq)

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

@@ -3337 +3337 @@
          call write_output("phisfi","Surface geopotential",
      &                    "m2s-2",phisfi(:))
+         call write_output("watercaptag","Watercap tag"
+     &         ,"Boolean",watercaptag(:))
          if (grid_type == regular_lonlat) then
            call write_output("area","Mesh area","m2",

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

@@ -285 +285 @@
 logical, dimension(:), intent(in) :: field
 ! Local argument used to convert logical to real
-real, dimension(ngrid) :: field_real
-logical :: is_active ! For XIOS, should this field be sent or not
-
+real, allocatable, dimension(:) :: field_real
+logical :: is_active ! For XIOS, should this field be sent or not
+
+allocate(field_real(size(field,1)))
 field_real = 0.
 where (field) field_real = 1.