Changeset 13995

Timestamp:

2020-12-02T14:46:49+01:00 (4 years ago)

Author:

laurent

Message:

Fix to access z0 without "Cdn_oce" in "ABL/ablmod.F90"!

File:

• NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL/ablmod.F90 (modified) (2 diffs)

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL/ablmod.F90

@@ -18 +18 @@
    USE dom_oce, ONLY  : tmask
    USE sbc_oce, ONLY  : ght_abl, ghw_abl, e3t_abl, e3w_abl, jpka, jpkam1, rhoa
-   USE sbcblk         ! use rn_efac, cdn_oce
+   USE sbcblk         ! use rn_efac
    USE sbc_phy        ! Catalog of functions for physical/meteorological parameters in the marine boundary layer
    !
@@ -135 +135 @@
          ustar2(ji,jj) = zzoce
 #endif
-         zrough(ji,jj) = ght_abl(2) * EXP( - vkarmn / SQRT( MAX( Cdn_oce(ji,jj), 1.e-4 ) ) ) !<-- recover the value of z0 from Cdn_oce
+         !#LB: sorry Cdn_oce is gone:
+         !zrough(ji,jj) = ght_abl(2) * EXP( - vkarmn / SQRT( MAX( Cdn_oce(ji,jj), 1.e-4 ) ) ) !<-- recover the value of z0 from Cdn_oce
       END_2D
+
+      zrough(:,:) = z0_from_Cd( ght_abl(2), pCd_du(:,:) / MAX( pwndm(:,:), 0.5_wp ) ! #LB: z0_from_Cd is define in sbc_phy.F90...
+         
       !
       !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<