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

Timestamp:

Apr 25, 2025, 6:26:55 PM (3 months ago)

Author:

jbclement

Message:

1D Mars PCM:

• The value for atmospheric water profile is now imposed as a uniform mixing ratio and not.
• When atmospheric water profile is monitored, ice value is not put to 0 anymore.
• The user can now prescribe a non-uniform atmospheric water profile with a file "profile_def_h2o_vap".

JBC

File:

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

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

@@ -81 +81 @@
 real, dimension(1)      :: latitude, longitude, cell_area
 
-! RV & JBC: Use of starting files for 1D
+! Use of starting files for 1D
 logical :: startfiles_1D, therestart1D, therestartfi, there
 
-! JN & JBC: Force atmospheric water profiles
+! Force atmospheric water profiles
 real                            :: atm_wat_profile, atm_wat_tau
-real, dimension(:), allocatable :: zqsat ! useful for (atm_wat_profile=2)
+real, dimension(nlayer)         :: q_def_h2o_vap
+real, dimension(:), allocatable :: zqsat
 !=======================================================================
 
@@ -149 +150 @@
 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,atm_wat_profile,atm_wat_tau)
+                     play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,q_def_h2o_vap)
 
 ! Write a "startfi" file
@@ -196 +197 @@
 
     ! Modify atmospheric water if asked
-    ! Added "atm_wat_profile" flag (JN + JBC)
-    ! ---------------------------------------
+    ! Added "atm_wat_profile" flag
+    ! ----------------------------
     if (water) then
         call watersat(nlayer,temp,play,zqsat)
-        if (0. <= atm_wat_profile) then ! If atmospheric water is monitored
+        if (0. <= atm_wat_profile .and. atm_wat_profile <= 1.) then ! If atmospheric water is monitored with a uniform value
             if (atm_wat_tau < 0.) then ! Case for prescribed atm_wat_profile: wet if >0, dry if =0
-                q(1,:,igcm_h2o_vap) = min(zqsat,atm_wat_profile*g/psurf)
-                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
-            else !  Case for relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
-                q(1,:,igcm_h2o_vap) = atm_wat_profile*g/psurf + (q(1,:,igcm_h2o_vap) - atm_wat_profile*g/psurf)*dexp(-dttestphys/atm_wat_tau)
-                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
+                q(1,:,igcm_h2o_vap) = min(zqsat,atm_wat_profile)
+            else ! Case for relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
+                q(1,:,igcm_h2o_vap) = atm_wat_profile + (q(1,:,igcm_h2o_vap) - atm_wat_profile)*dexp(-dttestphys/atm_wat_tau)
+            endif
+        else if (abs(atm_wat_profile + 2.) < 1.e-15) then ! If atmospheric water is monitored with a user-defined water profile (with file "profile_def_h2o_vap")
+            if (atm_wat_tau < 0.) then ! Case for prescribed atm_wat_profile: wet if >0, dry if =0
+                q(1,:,igcm_h2o_vap) = min(zqsat(:),q_def_h2o_vap(:))
+            else ! Case for relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
+                q(1,:,igcm_h2o_vap) = q_def_h2o_vap(:) + (q(1,:,igcm_h2o_vap) - q_def_h2o_vap(:))*dexp(-dttestphys/atm_wat_tau)
             endif
         endif