Changeset 4947 for LMDZ6

Timestamp:

May 21, 2024, 9:36:17 PM (6 months ago)

Author:

evignon

Message:

ajout de la parameterisation de rugosite de la neige sur les calottes de Amory et al. 2017

File:

• LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90

@@ -146 +146 @@
     REAL, DIMENSION(klon)  :: fluxbs_1, fluxbs_2, bsweight_fresh
     LOGICAL, DIMENSION(klon) :: ok_remaining_freshsnow
+    REAL  :: ta1, ta2, ta3, z01, z02, z03, coefa, coefb, coefc, coefd
+
 
 ! End definition
@@ -350 +352 @@
 !
 !****************************************************************************************
-    z0m = z0m_landice
-    z0h = z0h_landice
-    !z0m = SQRT(z0m**2+rugoro**2)
-
+
+if (z0m_landice .GT. 0.) then
+    z0m(1:knon) = z0m_landice
+    z0h(1:knon) = z0h_landice
+else
+    ! parameterization of z0=f(T) following measurements in Adelie Land by Amory et al 2018
+    coefa = 0.1658 !0.1862 !Ant
+    coefb = -50.3869 !-55.7718 !Ant
+    ta1 = 253.15 !255. Ant
+    ta2 = 273.15
+    ta3 = 273.15+3
+    z01 = exp(coefa*ta1 + coefb) !~0.2 ! ~0.25 mm
+    z02 = exp(coefa*ta2 + coefb) !~6  !~7 mm
+    z03 = z01
+    coefc = log(z03/z02)/(ta3-ta2)
+    coefd = log(z03)-coefc*ta3
+    do j=1,knon 
+      if (temp_air(j) .lt. ta1) then
+        z0m(j) = z01
+      else if (temp_air(j).ge.ta1 .and. temp_air(j).lt.ta2) then
+        z0m(j) = exp(coefa*temp_air(j) + coefb)
+      else if (temp_air(j).ge.ta2 .and. temp_air(j).lt.ta3) then
+        ! if st > 0, melting induce smooth surface
+        z0m(j) = exp(coefc*temp_air(j) + coefd)
+      else
+        z0m(j) = z03
+      endif
+      z0h(j)=z0m(j)
+    enddo
+
+endif    
+ 
 
 !****************************************************************************************
@@ -366 +396 @@
    if (ok_bs) then
        fluxbs(:)=0.
-       do j=1,klon
+       do j=1,knon
           ws1(j)=(u1(j)**2+v1(j)**2)**0.5
           ustar(j)=(cdragm(j)*(u1(j)**2+v1(j)**2))**0.5