Index: trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90	(revision 3836)
+++ trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90	(revision 3838)
@@ -195,6 +195,5 @@
     ! Wind increment: specific for 1D
     ! -------------------------------
-    ! The physics compute the tendencies on u and v,
-    ! here we just add Coriolos effect
+    ! The physics computes the tendencies on u and v, here we just add Coriolos effect
     !do ilayer = 1,nlayer
     !    du(ilayer) = du(ilayer) + ptif*(v(ilayer) - grv)
@@ -207,4 +206,19 @@
     !endif
 
+    ! Water tracer increment: specific for 1D
+    ! ---------------------------------------
+    ! The physics computes the tendency on q(1,:,igcm_h2o_vap), here we can mimic an effect of the "3D dynamics"
+    ! either by forcing a profile or by relaxing towards a prescribed profile
+    if (water .and. prescribed_h2ovap) then ! If atmospheric water profile is prescribed
+        if (h2ovap_relax_time < 0.) then ! Forcing
+            ! For some tests: unless it reaches saturation (maximal value)
+            !call watersat(nlayer,temp,play,zqsat)
+            !dq(1,:,igcm_h2o_vap) = (min(zqsat(:),q_prescribed_h2o_vap(:)) - q(1,:,igcm_h2o_vap))/dttestphys
+            dq(1,:,igcm_h2o_vap) = (q_prescribed_h2o_vap(:) - q(1,:,igcm_h2o_vap))/dttestphys
+        else ! Relaxation
+            dq(1,:,igcm_h2o_vap) = dq(1,:,igcm_h2o_vap) - (q(1,:,igcm_h2o_vap) - q_prescribed_h2o_vap(:))/h2ovap_relax_time
+        endif
+     endif
+
     ! Compute time for next time step
     ! -------------------------------
@@ -231,15 +245,4 @@
     play = aps + psurf*bps
 
-    ! Prescription of atmospheric water if asked
-    ! ------------------------------------------
-    if (water .and. prescribed_h2ovap) then ! If atmospheric water profile is prescribed
-        if (h2ovap_relax_time < 0.) then ! Forcing
-            !call watersat(nlayer,temp,play,zqsat)
-            dq(1,:,igcm_h2o_vap) = (q_prescribed_h2o_vap(:) - q(1,:,igcm_h2o_vap))/dttestphys
-        else ! Relaxation
-            dq(1,:,igcm_h2o_vap) = (q_prescribed_h2o_vap(:) + (q(1,:,igcm_h2o_vap) - q_prescribed_h2o_vap(:))*dexp(-dttestphys/h2ovap_relax_time) - q(1,:,igcm_h2o_vap))/dttestphys
-        endif
-     endif
-
     ! Compute tracers for next time step
     ! ----------------------------------