Changeset 3100 for trunk/LMDZ.GENERIC/libf/aeronostd

Timestamp:

Oct 25, 2023, 10:34:54 AM (20 months ago)

Author:

bhatnags

Message:

Generic-PCM

This commit updates the slab ocean module to a parallelisable dynamic slab ocean module. This is particularly relevant if you want to be able to use oceanic heat transport in parallel mode.

It has the following features:

(a) Computes sea ice creation and evolution.
(b) Snow has thermodynamic properties.
(c) Computes oceanic horizontal transport (diffusion & surface-wind driven Ekman transport).
(d) Can be used in parallel mode.

Required callphys.def flags:
The slab ocean and its dependencies can be activated with the following flags (already added to deftank):
## Ocean options
## ~
# Model slab-ocean (Main flag for slab ocean)
ok_slab_ocean = .true.
# The following flags can only be set to true if ok_slab_ocean is true
# Ekman transport
slab_ekman = .true.
# Ekman zonal advection
slab_ekman_zonadv = .true.
# Horizontal diffusion (default coef_hdiff=25000., can be changed)
slab_hdiff = .true.
# Slab-ocean timestep (in physics timesteps)
cpl_pas = 1
# Gent-McWilliams? Scheme (can only be true if slab_ekman is true)
slab_gm = .true.

Notes:
In the current state, the model crashes if moistadjustment = .true. Unsure whether this is due to the slab or is an inherent issue with moistadj (under investigation).

SB and EM

File:

• trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F

@@ -1172 +1172 @@
 
       use chimiedata_h,  only: albedo_snow_chim, albedo_co2_ice_chim
-      use slab_ice_h,    only: h_alb_ice, 
-     &                         alb_ice_min, alb_ice_max
-      use tracer_h,      only: igcm_h2o_ice, igcm_co2_ice
+!      use slab_ice_h,    only: h_alb_ice, alb_ice_min, alb_ice_max
+      use ocean_slab_mod, only: h_alb_ice
+      use ocean_slab_mod, only: alb_ice_min
+      use ocean_slab_mod, only: alb_ice_max
       use callkeys_mod,  only: ok_slab_ocean, co2cond, alb_ocean
       use phys_state_var_mod, only: albedo_bareground,