Changeset 4673

Timestamp:

Sep 4, 2023, 4:49:20 PM (9 months ago)

Author:

evignon

Message:

correction commit precedent

Location:

LMDZ6/trunk/libf

Files:

• phylmd/surf_landice_mod.F90 (modified) (1 diff)
• phylmdiso/surf_landice_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90

@@ -379 +379 @@
            ! snow density increases with snow age and
            ! increases even faster in case of sedimentation of blowing snow or rain
-           tau_dens=max(tau_densmin, & tau_dens0*exp(-abs(precip_bs(j))/pbst_bs-abs(precip_rain(j))/prt_bs) &
-                    *exp(-max(tsurf(j)-272.0,0.)))
+           tau_dens=max(tau_densmin, tau_dens0*exp(-abs(precip_bs(j))/pbst_bs - & 
+                    abs(precip_rain(j))/prt_bs)*exp(-max(tsurf(j)-272.0,0.)))
            rhos(j)=rho0+(rhohard-rho0)*(1.-exp(-agesno(j)*86400.0/tau_dens))
            ! blowing snow flux formula used in MAR

LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90

@@ -468 +468 @@
            ! snow density increases with snow age and
            ! increases even faster in case of sedimentation of blowing snow or rain
-           tau_dens=max(tau_densmin, & tau_dens0*exp(-abs(precip_bs(j))/pbst_bs-abs(precip_rain(j))/prt_bs) &
-                    *exp(-max(tsurf(j)-272.0,0.)))
+           tau_dens=max(tau_densmin, tau_dens0*exp(-abs(precip_bs(j))/pbst_bs &
+                   -abs(precip_rain(j))/prt_bs) *exp(-max(tsurf(j)-272.0,0.)))
            rhos(j)=rho0+(rhohard-rho0)*(1.-exp(-agesno(j)*86400.0/tau_dens))
            ! blowing snow flux formula used in MAR