Changeset 3333 for trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90

Timestamp:

May 17, 2024, 11:30:17 AM (7 months ago)

Author:

llange

Message:

Mars PCM
Fixing a bug in vdif_cd: a "residual", used as criterion to enter an iterative loop, was wrongly initialized. Hence, for some points, the algorithm does not go into the loop, and a wrong value of Cd, Ch was computed.
Also some cleaning/small fixing with save variables with OMP.

LL

File:

• trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90 (modified) (1 diff)

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

@@ -20 +20 @@
     real,    save                              :: albedo_perennialco2     ! Albedo for perennial co2 ice [1]
     logical, save                              :: lag_layer               ! Does lag layer is present?
-    logical, save                              :: include_waterbuoyancy    ! Include the effect of water buoyancy when computing the sublimation of water ice ?
-!$OMP THREADPRIVATE(h2o_ice_depth,d_coef,lag_co2_ice,albedo_perennialco2,include_waterbuoyancy)
+    logical, save                              :: include_waterbuoyancy   ! Include the effect of water buoyancy when computing the sublimation of water ice ?
+!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,albedo_perennialco2,lag_layer,include_waterbuoyancy)
 
 !=======================================================================