Changeset 5363

Timestamp:

Dec 3, 2024, 8:26:45 AM (43 hours ago)

Author:

evignon

Message:

Valentin Wiener, Etienne Vignon: ustar dans les sorties n'est plus un diag ad-hoc mais la variable
directement reliée à la tension de vent en surface

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -849 +849 @@
        CALL histwrite_phy(o_sicf, zx_tmp_fi2d)
        CALL histwrite_phy(o_q2m, zq2m)
-       IF (vars_defined) zx_tmp_fi2d = zustar
+       !IF (vars_defined) zx_tmp_fi2d = zustar
+       !set ustar output variable as directly related to actual surface stress
+       IF (vars_defined) THEN
+          zx_tmp_fi2d=0.
+          DO nsrf=1,nbsrf
+             zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*((t_seri(:,1)*RD/pplay(:,1))**0.5)*(fluxu(:,1,nsrf)**2+fluxv(:,1,nsrf)**2)**0.25
+          ENDDO           
+       ENDIF
        CALL histwrite_phy(o_ustar, zx_tmp_fi2d)
        CALL histwrite_phy(o_u10m, zu10m)