Changeset 3381 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Jun 13, 2024, 3:38:38 PM (6 months ago)

Author:

llange

Message:

Mars PEM
Fixing bug to update the CO2 ice tendency: the case when a glacier diseapper was not handled
LL

File:

• trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope_mod.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope_mod.F90

@@ -52 +52 @@
                       - (beta_clap_co2/(alpha_clap_co2-log(vmr_co2_pem(i,t)*ps_PCM_2(i,t)*(global_avg_press_new/global_avg_press_PCM)/100.)))**4
             enddo
+            if (ave < 1e-4) ave = 0.
+            tendencies_co2_ice_phys(i,islope) = tendencies_co2_ice_phys_ini(i,islope) - coef*ave/timelen
         endif
-        if (ave < 1e-4) ave = 0.
-        tendencies_co2_ice_phys(i,islope) = tendencies_co2_ice_phys_ini(i,islope) - coef*ave/timelen
     enddo
 enddo