Changeset 5188 for LMDZ6/trunk/libf

Timestamp:

Sep 14, 2024, 2:23:57 PM (3 months ago)

Author:

evignon

Message:

changes in the roughness length parameterization over landice.
Now account for a 2-day average temperature instead of instantaneous temperature
Valentin Wiener

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• surf_landice_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -400 +400 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: tpot, tpote, ue, uq, uwat, ve, vq, vwat, zxffonte
 !$OMP THREADPRIVATE(tpot, tpote, ue, uq, uwat, ve, vq, vwat, zxffonte)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxustartlic, zxrhoslic, zxqsaltlic
-!$OMP THREADPRIVATE(zxustartlic, zxrhoslic, zxqsaltlic)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxustartlic, zxrhoslic, zxqsaltlic, tempsmoothlic
+!$OMP THREADPRIVATE(zxustartlic, zxrhoslic, zxqsaltlic, tempsmoothlic)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfqcalving
 !$OMP THREADPRIVATE(zxfqcalving)
@@ -1061 +1061 @@
       ALLOCATE(zxrunofflic(klon), runoff_diag(klon))
       runoff_diag(:)=0.
-      ALLOCATE(zxustartlic(klon), zxrhoslic(klon), zxqsaltlic(klon))
-      zxustartlic(:)=0. ; zxrhoslic(:)=0. ; zxqsaltlic(:)=0.
+      ALLOCATE(zxustartlic(klon), zxrhoslic(klon), zxqsaltlic(klon), tempsmoothlic(klon))
+      zxustartlic(:)=0. ; zxrhoslic(:)=0. ; zxqsaltlic(:)=0. ; tempsmoothlic(:)=0.
       ALLOCATE(rain_lsc(klon))
       ALLOCATE(rain_num(klon))
@@ -1483 +1483 @@
 ! SN runoff_diag
       DEALLOCATE(zxrunofflic, runoff_diag)
-      DEALLOCATE(zxustartlic, zxrhoslic, zxqsaltlic)
+      DEALLOCATE(zxustartlic, zxrhoslic, zxqsaltlic, tempsmoothlic)
       DEALLOCATE(zxtsol, snow_lsc, zxfqfonte, zxqsurf)
       DEALLOCATE(rain_lsc)

LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90

@@ -36 +36 @@
     USE cpl_mod,          ONLY : cpl_send_landice_fields
     USE calcul_fluxs_mod
-    USE phys_local_var_mod, ONLY : zxrhoslic, zxustartlic, zxqsaltlic
+    USE phys_local_var_mod, ONLY : zxrhoslic, zxustartlic, zxqsaltlic, tempsmoothlic
     USE phys_output_var_mod, ONLY : snow_o,zfra_o
 #ifdef ISO   
@@ -185 +185 @@
 
     REAL,DIMENSION(klon) :: alb1,alb2
+    REAL                 :: time_tempsmooth,coef_tempsmooth
     REAL,DIMENSION(klon) :: precip_totsnow, evap_totsnow
     REAL, DIMENSION (klon,6) :: alb6
@@ -230 +231 @@
            PRINT*, 'alb_nir_sno_lic',alb_nir_sno_lic
   
+  DO j=1,knon
+       i = knindex(j)
+       tempsmoothlic(i) = temp_air(j)
+  ENDDO
   firstcall=.false.
   ENDIF
@@ -445 +450 @@
     z0h(1:knon) = z0h_landice
 else
-    ! parameterization of z0=f(T) following measurements in Adelie Land by Amory et al 2018
+    ! parameterization of z0=f(T) following measurements in Adelie Land by Amory et al 2017
     coefa = 0.1658 !0.1862 !Ant
     coefb = -50.3869 !-55.7718 !Ant
@@ -456 +461 @@
     coefc = log(z03/z02)/(ta3-ta2)
     coefd = log(z03)-coefc*ta3
+    time_tempsmooth=2.*86400.
+    coef_tempsmooth=min(1.,dtime/time_tempsmooth)
+    !coef_tempsmooth=0.
     do j=1,knon 
-      if (temp_air(j) .lt. ta1) then
+      i=knindex(j)
+      print*, "tempsmoothlic", tempsmoothlic(i)
+      tempsmoothlic(i)=temp_air(j)*coef_tempsmooth+tempsmoothlic(i)*(1.-coef_tempsmooth)
+      if (tempsmoothlic(i) .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
+      else if (tempsmoothlic(i).ge.ta1 .and. tempsmoothlic(i).lt.ta2) then
+        z0m(j) = exp(coefa*tempsmoothlic(i) + coefb)
+      else if (tempsmoothlic(i).ge.ta2 .and. tempsmoothlic(i).lt.ta3) then
         ! if st > 0, melting induce smooth surface
-        z0m(j) = exp(coefc*temp_air(j) + coefd)
+        z0m(j) = exp(coefc*tempsmoothlic(i) + coefd)
       else
         z0m(j) = z03