Changeset 3838

Timestamp:

Jul 7, 2025, 3:11:29 PM (2 weeks ago)

Author:

jbclement

Message:

Mars PCM 1D:
Fixing formula for prescribing atmospheric water profile in the 1D model (already under discretization) + Addition of more explicit comments.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/testphys1d.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4894 +4894 @@
 
 == 04/07/2025 == JBC
-Imposing the prescribed atmospheric water profile though 'dq' instead of 'q' directly (more correct).
+Imposing the prescribed atmospheric water profile though 'dq' instead of 'q' directly for the 1D model (more correct).
+
+== 04/07/2025 == JBC
+More explicit code and comments about the prescribtion of the atmospheric water profile in the 1D model.
+
+== 07/07/2025 == JBC
+Fixing formula for prescribing atmospheric water profile in the 1D model (already under discretization) + Addition of more explicit comments.

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

@@ -195 +195 @@
     ! 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 +206 @@
     !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 +245 @@
     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
     ! ----------------------------------