Changeset 3912 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

Timestamp:

Sep 5, 2025, 2:46:19 PM (4 months ago)

Author:

jbclement

Message:

Mars PCM 1D:
Prescription of the atmospheric water ice profile in 1D. The boolean 'ctrl_h2oice' activates the option, the file "profile_ref_h2ovap" defines the prescribed profile and the value 'relaxtime_h2oice', if positive, activates the relaxation and gives the constant of time.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90 (modified) (4 diffs)

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

@@ -6 +6 @@
 use phyredem,            only: physdem0, physdem1
 use watersat_mod,        only: watersat
-use tracer_mod,          only: igcm_h2o_vap, noms
+use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, noms
 use comcstfi_h,          only: pi, g, rcp, cpp
 use time_phylmdz_mod,    only: daysec
@@ -83 +83 @@
 logical :: startfiles_1D, therestart1D, therestartfi, there
 
-! Prescription of atmospheric water profile
-logical                         :: prescribed_h2ovap
-real                            :: h2ovap_relax_time
-real, dimension(nlayer)         :: q_prescribed_h2o_vap
+! Prescription of atmospheric water/ice profiles
+logical                         :: ctrl_h2ovap, ctrl_h2oice
+real                            :: relaxtime_h2ovap, relaxtime_h2oice
+real, dimension(nlayer)         :: qref_h2ovap, qref_h2oice
 real, dimension(:), allocatable :: zqsat
 !=======================================================================
@@ -142 +142 @@
 call init_testphys1d('start1D.txt','startfi.nc',therestart1D,therestartfi,ngrid,nlayer,odpref,nq,q, &
                      time,psurf,u,v,temp,ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w, &
-                     play,plev,latitude,longitude,cell_area,prescribed_h2ovap,h2ovap_relax_time,q_prescribed_h2o_vap)
+                     play,plev,latitude,longitude,cell_area,                                        &
+                     ctrl_h2ovap,relaxtime_h2ovap,qref_h2ovap,ctrl_h2oice,relaxtime_h2oice,qref_h2oice)
 
 ! Write a "startfi" file
@@ -208 +209 @@
     ! 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
+    ! The physics computes the tendency on q for igcm_h2o_vap/igcm_h2o_ice, here we can mimic an effect
+    ! of the "3D dynamics" either by forcing a profile or by relaxing towards a prescribed profile
+    if (water) then
+        if (ctrl_h2ovap) then ! If atmospheric water vapour profile is prescribed
+            if (relaxtime_h2ovap < 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(:),qref_h2ovap(:)) - q(1,:,igcm_h2o_vap))/dttestphys
+                dq(1,:,igcm_h2o_vap) = (qref_h2ovap(:) - q(1,:,igcm_h2o_vap))/dttestphys
+            else ! Relaxation
+                dq(1,:,igcm_h2o_vap) = dq(1,:,igcm_h2o_vap) - (q(1,:,igcm_h2o_vap) - qref_h2ovap(:))/relaxtime_h2ovap
+            endif
+        endif
+        if (ctrl_h2oice) then ! If atmospheric water ice profile is prescribed
+            if (relaxtime_h2oice < 0.) then ! Forcing
+                dq(1,:,igcm_h2o_ice) = (qref_h2oice(:) - q(1,:,igcm_h2o_ice))/dttestphys
+            else ! Relaxation
+                dq(1,:,igcm_h2o_ice) = dq(1,:,igcm_h2o_ice) - (q(1,:,igcm_h2o_ice) - qref_h2oice(:))/relaxtime_h2oice
+            endif
         endif
      endif