Changeset 3200

Timestamp:

Feb 2, 2024, 6:06:10 PM (10 months ago)

Author:

jbclement

Message:

Mars PCM:
Small update following r3188 to keep the CO2 mixing ration constant if "CO2cond_ps = 0" + Some cleanings.
JBC

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/changelog.txt

@@ -4456 +4456 @@
 
 == 29/01/2024 == JBC
-Addition of the parameter 'CO2cond_ps' (= 1 by default) for 1D. This coefficient controls the surface pressure change. If 'relaxcoef_ps = 1', then surface pressure varies normally. If 'relaxcoef_ps = 0', then surface pressure is kept constant. The ratio of polar cap surface over planetary surface is a typical value (8.3e-4) for tests. To be defined in "callphys.def" so that both PCM and PEM can read it.
+Addition of the parameter 'CO2cond_ps' (= 1 by default) for 1D. This coefficient controls the surface pressure change. If 'CO2cond_ps = 1', then surface pressure varies normally. If 'CO2cond_ps = 0', then surface pressure is kept constant. The ratio of polar cap surface over planetary surface is a typical value (8.3e-4) for tests. To be defined in "callphys.def" so that both PCM and PEM can read it.
 
 == 30/01/2024 == CS
 Update of the reference to default values of coeff_detrainment, coeff_injection, ti_injection and tf_injection (following old revision r2639) in callphys.def.MCD6
+
+== 02/02/2024 == JBC
+Small update following r3188 to keep the CO2 mixing ration constant if "CO2cond_ps = 0" + Some cleanings.

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

@@ -78 +78 @@
 real, dimension(1),                  intent(out) :: latitude, longitude, cell_area
 real,                                intent(out) :: atm_wat_profile, atm_wat_tau ! Force atmospheric water profiles
-real,                                intent(out) :: CO2cond_ps                 ! Relaxation coefficient for psurf
+real,                                intent(out) :: CO2cond_ps                   ! Coefficient to control the surface pressure change
 
 !=======================================================================
@@ -530 +530 @@
 endif !(.not. therestartfi)
 
-! emissivity
+! Emissivity
 ! ----------
 if (.not. therestartfi) then

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

@@ -270 +270 @@
 
     ! Increment tracers
-    q(1,:,:) = q(1,:,:) + dttestphys*dq(1,:,:)
+    if (abs(CO2cond_ps) < 1.e-10) then
+        do iq = 1,nq
+            if (iq == igcm_co2) cycle
+            q(1,:,iq) = q(1,:,iq) + dttestphys*dq(1,:,iq)
+        enddo
+    else
+        q = q + dttestphys*dq
+    endif
 enddo ! End of time stepping loop (idt=1,ndt)
 

trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90

@@ -6 +6 @@
 !   author: LL, 06/2023
 !   ------
-!   
+!
 !=======================================================================
 
-    IMPLICIT NONE
+implicit none
 
-    LOGICAL, SAVE :: paleoclimate    ! False by default, is activate  for paleoclimates specific processes (e.g., lag layer)
-                                     ! is initialized in conf_phys
+logical, save :: paleoclimate ! False by default, is activate  for paleoclimates specific processes (e.g., lag layer)
+                              ! is initialized in conf_phys
+
 !$OMP THREADPRIVATE(paleoclimate)
-
-    real, save, allocatable :: h2o_ice_depth(:,:)  ! Thickness of the lag before H2O ice [m]
-    real, save, allocatable :: lag_co2_ice(:,:)  ! Thickness of the lag before CO2 ice [m]
-    real, save :: albedo_perennialco2             ! Albedo for perennial co2 ice [1]
-    real, save, allocatable :: d_coef(:,:)  ! Diffusion coeficent
-    LOGICAL,SAVE :: lag_layer ! does lag layer is present?
+    real,    save, allocatable, dimension(:,:) :: h2o_ice_depth       ! Thickness of the lag before H2O ice [m]
+    real,    save, allocatable, dimension(:,:) :: lag_co2_ice         ! Thickness of the lag before CO2 ice [m]
+    real,    save, allocatable, dimension(:,:) :: d_coef              ! Diffusion coeficent
+    real,    save                              :: albedo_perennialco2 ! Albedo for perennial co2 ice [1]
+    logical, save                              :: lag_layer           ! Does lag layer is present?
 !$OMP THREADPRIVATE(h2o_ice_depth,d_coef,lag_co2_ice,albedo_perennialco2)
 
-    CONTAINS
+!=======================================================================
+contains
+!=======================================================================
 
+SUBROUTINE ini_paleoclimate_h(ngrid,nslope)
 
-  subroutine ini_paleoclimate_h(ngrid,nslope)
+implicit none
 
-  implicit none
-  integer,intent(in) :: ngrid  ! number of atmospheric columns
-  integer,intent(in) :: nslope ! number of slope within a mesh 
+integer, intent(in) :: ngrid  ! number of atmospheric columns
+integer, intent(in) :: nslope ! number of slope within a mesh
 
-    allocate(h2o_ice_depth(ngrid,nslope))
-    allocate(lag_co2_ice(ngrid,nslope))
-    allocate(d_coef(ngrid,nslope))
-  end subroutine ini_paleoclimate_h
+allocate(h2o_ice_depth(ngrid,nslope))
+allocate(lag_co2_ice(ngrid,nslope))
+allocate(d_coef(ngrid,nslope))
 
-  subroutine end_paleoclimate_h
+END SUBROUTINE ini_paleoclimate_h
 
-  implicit none
-    if (allocated(d_coef)) deallocate(d_coef)
-    if (allocated(h2o_ice_depth)) deallocate(h2o_ice_depth)
-    if (allocated(lag_co2_ice)) deallocate(lag_co2_ice)
-  end subroutine end_paleoclimate_h
+!=======================================================================
+SUBROUTINE end_paleoclimate_h
 
+implicit none
 
-    END MODULE paleoclimate_mod
+if (allocated(d_coef)) deallocate(d_coef)
+if (allocated(h2o_ice_depth)) deallocate(h2o_ice_depth)
+if (allocated(lag_co2_ice)) deallocate(lag_co2_ice)
+
+END SUBROUTINE end_paleoclimate_h
+
+END MODULE paleoclimate_mod