Changeset 4245 for LMDZ6/trunk/libf/phylmd

Timestamp:

Sep 14, 2022, 2:53:51 PM (2 years ago)

Author:

evignon

Message:

controle des z0 sur landice dans les .def et nettoyage surf_landice_mod
pour Valentin

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys_m.F90 (modified) (3 diffs)
• surf_landice_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/clesphys.h

@@ -54 +54 @@
        REAL min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce
        REAL z0m_seaice,z0h_seaice
+       REAL z0m_landice, z0h_landice
        INTEGER iflag_gusts,iflag_z0_oce
 
@@ -113 +114 @@
      &     , f_cdrag_ter,f_cdrag_oce,f_rugoro,z0min,tau_gl              &
      &     , min_wind_speed,f_gust_wk,f_gust_bl,f_qsat_oce,f_z0qh_oce   &
-     &     , z0m_seaice,z0h_seaice                                      &
+     &     , z0m_seaice,z0h_seaice,z0m_landice,z0h_landice              &
      &     , freq_outNMC, freq_calNMC                                   &
      &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -133 +133 @@
     REAL, SAVE          :: f_rugoro_omp   , z0min_omp
     REAL, SAVE          :: z0m_seaice_omp,z0h_seaice_omp
+    REAL, SAVE          :: z0m_landice_omp,z0h_landice_omp
     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
@@ -2123 +2124 @@
     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)
+    z0h_landice_omp = 0.001 ; CALL getin('z0h_landice',z0h_landice_omp)
 
     f_rugoro_omp = 0.
@@ -2666 +2671 @@
     z0m_seaice=z0m_seaice_omp
     z0h_seaice=z0h_seaice_omp
+    z0m_landice=z0m_landice_omp
+    z0h_landice=z0h_landice_omp
 
     f_rugoro=f_rugoro_omp

LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90

@@ -156 +156 @@
            PRINT*, 'alb_nir_sno_lic',alb_nir_sno_lic
   
-!  z0m=1.e-3
-!  z0h = z0m
   firstcall=.false.
   ENDIF
@@ -351 +349 @@
 !
 !****************************************************************************************
-    z0m=1.e-3
-    z0h = z0m
-    z0m = SQRT(z0m**2+rugoro**2)
+    z0m = z0m_landice
+    z0h = z0h_landice
+    !z0m = SQRT(z0m**2+rugoro**2)