Changeset 7811 for branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

Timestamp:

2017-03-20T14:36:18+01:00 (7 years ago)

Author:

jcastill

Message:

Further fix to the calculation of surface roughness for the significant wave height read from a wave model (nn_z0_met=3)

File:

• branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90 (modified) (1 diff)

branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

@@ -200 +200 @@
          zhsro(:,:) = MAX(rsbc_zs2 * ustars2(:,:) * zdep(:,:)**1.5, rn_hsro) ! zhsro = rn_frac_hs * Hsw (eq. 11)
       CASE ( 3 )             ! Roughness given by the wave model (coupled or read in file) 
-         zhsro(:,:) = hsw(:,:)
+         WHERE( hsw == 0._wp ) ! surface roughness length according to Charnock formula when sign. wave height 0
+            zhsro = MAX(rn_charn / grav * ustars2, rn_hsro)
+         ELSEWHERE
+            zhsro = MAX(hsw, rn_hsro)
+         END WHERE
       END SELECT