Changeset 6078 for LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

Timestamp:

Feb 16, 2026, 6:01:58 PM (5 days ago)

Author:

evignon

Message:

harmonisation des lectures de z0 dans les .def afin de pouvoir forcer les cas 1D au format standard
en z0m et z0h sur tout type de sous-surface

File:

• LMDZ6/trunk/libf/phylmd/conf_phys_m.f90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -129 +129 @@
     REAL, SAVE          :: f_rugoro_omp   , z0min_omp
     REAL, SAVE          :: z0m_seaice_omp,z0h_seaice_omp
-    REAL, SAVE          :: z0m_landice_omp,ratio_z0hz0m_landice_omp
+    REAL, SAVE          :: z0m_landice_omp,z0h_landice_omp
+    REAL, SAVE          :: z0m_ocean_omp, z0h_ocean_omp
+    REAL, SAVE          :: z0m_land_omp, z0h_land_omp
+    REAL                :: ratio_z0hz0m_land, ratio_z0hz0m_ocean, ratio_z0hz0m_landice, ratio_z0hz0m_seaice
     REAL, SAVE          :: min_wind_speed_omp,f_gust_wk_omp,f_gust_bl_omp,f_qsat_oce_omp, f_z0qh_oce_omp
-    INTEGER, SAVE       :: iflag_gusts_omp,iflag_z0_oce_omp
+    INTEGER, SAVE       :: iflag_gusts_omp,iflag_z0_oce_omp, iflag_z0_ter_omp
 
     REAL :: seuil_inversion
@@ -1004 +1007 @@
     !Config Desc = fraction of shortwave penetrating into the ice
     !Config Def  = 0.3
-    !Config Help = pour le modele de glace de mer
+    !Config Help = pour le modale de glace de mer
     si_pen_frac_omp = 0.3
     CALL getin('si_pen_frac', si_pen_frac_omp)
@@ -1932 +1935 @@
     iflag_z0_oce_omp=0
     CALL getin('iflag_z0_oce',iflag_z0_oce_omp)
+
+    iflag_z0_ter_omp=0
+    CALL getin('iflag_z0_ter',iflag_z0_ter_omp)
+
     !
     iflag_gusts_omp=0
@@ -1939 +1946 @@
     CALL getin('min_wind_speed',min_wind_speed_omp)
 
-    z0m_seaice_omp = 0.002 ; CALL getin('z0m_seaice',z0m_seaice_omp)
-    z0h_seaice_omp = 0.002 ; CALL getin('z0h_seaice',z0h_seaice_omp)
-
-
-    z0m_landice_omp = 0.001 ; CALL getin('z0m_landice',z0m_landice_omp)
-    ratio_z0hz0m_landice_omp = 1. ; CALL getin('ratio_z0hz0m_landice',ratio_z0hz0m_landice_omp)
-
     f_rugoro_omp = 0.
     CALL getin('f_rugoro',f_rugoro_omp)
@@ -1952 +1952 @@
     CALL getin('z0min',z0min_omp)
 
-    ! PARAMETERS FOR SNOW AND ICE MELTING
+
+
+    z0m_seaice_omp = 0.002
+    CALL getin('z0m_seaice',z0m_seaice_omp)
+    ratio_z0hz0m_seaice =  1.
+    CALL getin('ratio_z0hz0m_seaice',ratio_z0hz0m_seaice)
+    z0h_seaice_omp = z0m_seaice_omp * ratio_z0hz0m_seaice
+    CALL getin('z0h_seaice',z0h_seaice_omp)
+
+    z0m_landice_omp = 0.001
+    CALL getin('z0m_landice',z0m_landice_omp)
+    ratio_z0hz0m_landice =  1.
+    CALL getin('ratio_z0hz0m_landice',ratio_z0hz0m_landice)
+    z0h_landice_omp = z0m_landice_omp * ratio_z0hz0m_landice
+    CALL getin('z0h_landice',z0h_landice_omp)
+
+    z0m_ocean_omp = z0min_omp
+    CALL getin('z0m_ocean',z0m_ocean_omp)
+    ratio_z0hz0m_ocean =  1.
+    CALL getin('ratio_z0hz0m_ocean',ratio_z0hz0m_ocean)
+    z0h_ocean_omp = z0m_ocean_omp * ratio_z0hz0m_ocean
+    CALL getin('z0h_ocean',z0h_ocean_omp)
+
+    z0m_land_omp = z0min_omp
+    CALL getin('z0m_land',z0m_land_omp)
+    ratio_z0hz0m_land =  1.
+    CALL getin('ratio_z0hz0m_land',ratio_z0hz0m_land)
+    z0h_land_omp = z0m_land_omp * ratio_z0hz0m_land
+    CALL getin('z0h_land',z0h_land_omp)
+
+
+      ! PARAMETERS FOR SNOW AND ICE MELTING
     chasno_tun_omp=0.15
     CALL getin('chasno_tun',chasno_tun_omp)
@@ -2379 +2410 @@
     iflag_gusts=iflag_gusts_omp
     iflag_z0_oce=iflag_z0_oce_omp
+    iflag_z0_ter=iflag_z0_ter_omp
 
     z0m_seaice=z0m_seaice_omp
     z0h_seaice=z0h_seaice_omp
     z0m_landice=z0m_landice_omp
-    ratio_z0hz0m_landice=ratio_z0hz0m_landice_omp
+    z0h_landice=z0h_landice_omp
+    z0m_ocean=z0m_ocean_omp
+    z0h_ocean=z0h_ocean_omp
+    z0m_land=z0m_land_omp
+    z0h_land=z0h_land_omp
 
     chasno_tun=chasno_tun_omp
@@ -2682 +2718 @@
     WRITE(lunout,*) ' f_rugoro = ',f_rugoro
     WRITE(lunout,*) ' z0min = ',z0min
+    WRITE(lunout,*) ' z0m_landice = ',z0m_landice
+    WRITE(lunout,*) ' z0h_landice = ',z0h_landice
+    WRITE(lunout,*) ' z0m_seaice = ',z0m_seaice
+    WRITE(lunout,*) ' z0h_seaice = ',z0h_seaice
+    WRITE(lunout,*) ' z0m_ocean = ',z0m_ocean
+    WRITE(lunout,*) ' z0h_ocean = ',z0h_ocean
+    WRITE(lunout,*) ' z0m_land = ',z0m_land
+    WRITE(lunout,*) ' z0h_land = ',z0h_land
     WRITE(lunout,*) ' tmax_fonte_cv = ', tmax_fonte_cv
     WRITE(lunout,*) ' iflag_wake = ', iflag_wake