Index: /trunk/LMDZ.MARS/changelog.txt
===================================================================
--- /trunk/LMDZ.MARS/changelog.txt	(revision 3106)
+++ /trunk/LMDZ.MARS/changelog.txt	(revision 3107)
@@ -4273,4 +4273,11 @@
 Prettyfied solarlong.F and made it a module. Likewise for conf_phys.F
 
+Some code tidying:
+Made pi in module comcstfi_h a parameter (and not a variable recomputed at
+various points by various routines) and added module routine init_comcstfi_h
+to cleanly initialize module variables.
+Moved iniorbit.F to be a module routine of planete_h since it initializes
+(some of ) the module variables it contains.
+
 == 23/10/2023 == EV
 We added the to module vdifc the possibily of subsurface intercation, mostly to have the option of buried glaicer that can lose ice and create polar layers.
@@ -4288,9 +4295,4 @@
 qeq, SSI water vapor quantity
 
-
-Some code tidying:
-Made pi in module comcstfi_h a parameter (and not a variable recomputed at
-various points by various routines) and added module routine init_comcstfi_h
-to cleanly initialize module variables.
-Moved iniorbit.F to be a module routine of planete_h since it initializes
-(some of ) the module variables it contains.
+== 26/10/2023 == JBC
+Few small fixes following r3098.
Index: /trunk/LMDZ.MARS/deftank/field_def_physics_mars.xml
===================================================================
--- /trunk/LMDZ.MARS/deftank/field_def_physics_mars.xml	(revision 3106)
+++ /trunk/LMDZ.MARS/deftank/field_def_physics_mars.xml	(revision 3107)
@@ -293,11 +293,11 @@
 
             <field id="fluxtop_lw"
-                   long_name="Longwave radiation to space at the top of the atmosphere" 
+                   long_name="Longwave radiation at the top of the atmosphere" 
                    unit="W.m-2" />
             <field id="fluxtop_dn_sw"
-                   long_name="Incoming Solar radiative flux from space" 
+                   long_name="Incoming shortwave radiation at the top of the atmosphere" 
                    unit="W.m-2" />
             <field id="fluxtop_up_sw"
-                   long_name="Outgoing Solar radiative flux to space" 
+                   long_name="Upward shortwave radiation at the top of the atmosphere" 
                    unit="W.m-2" />
 
@@ -325,7 +325,4 @@
             <field id="fluxsurf_dn_sw_slope07"
                    long_name="Incoming shortwave radiation at the surface on slope 07"  
-                   unit="W.m-2" />
-            <field id="fluxtop_up_sw"
-                   long_name="Upward shortwave radiation at the top of the atmosphere" 
                    unit="W.m-2" />
             
@@ -1209,4 +1206,16 @@
                    long_name="Waterdensity_soil of slope 07"
                    unit="kg.m-3" />
+            <field id="zdqsdif_ssi_frost"
+                   long_name="Flux between frost and subsurface"
+                   unit="kg.m-2.s-1" />
+            <field id="subtimestep"
+                   long_name="vdifc substimestep length"
+                   unit="s" />
+            <field id="zdq_subtimestep"
+                   long_name="Actual flux zdqsdif*subtimestep"
+                   unit="kg.m-2" />
+            <field id="zdq_end"
+                   long_name="Flux after all contributions"
+                   unit="kg.m-2.s-1" />
 
         </field_group>
Index: /trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3106)
+++ /trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3107)
@@ -1110,7 +1110,8 @@
               ENDIF
               call write_output('zdqsdif_ssi_frost',
-     &                '','',zdqsdif_ssi_frost(ig,1))
+     &                'Flux between frost and subsurface','kg.m-2.s-1',
+     &                zdqsdif_ssi_frost(ig,1))
               call write_output('subtimestep',
-     &                '','',subtimestep)
+     &                'vdifc substimestep length','s',subtimestep)
              ! ENDDO !subsurface ice subslope
 
@@ -1284,12 +1285,12 @@
 
          call write_output("surf_h2o_lh",
-     &                          "Ground ice latent heat flux",
-     &                               "W.m-2",surf_h2o_lh(:,iflat))
-                       call write_output('zdq_subtimestep',
-     &                '','',zdqsdif(:)*subtimestep)
-                       call write_output('zdq_end',
-     &                '','',zdqsdif(:))
-                       call write_output('vdifc_subtimestep',
-     &                '','',subtimestep)
+     &                     "Ground ice latent heat flux",
+     &                     "W.m-2",surf_h2o_lh(:,iflat))
+         call write_output('zdq_subtimestep',
+     &                     'Actual flux zdqsdif*subtimestep',
+     &                     'kg.m-2',zdqsdif(:)*subtimestep)
+         call write_output('zdq_end',
+     &                     'Flux after all contributions',
+     &                     'kg.m-2.s-1',zdqsdif(:))
 C       Diagnostic output for HDO
 !        if (hdo) then