From eb3580a89c2322b29b76a40715bac53b5a4edad3 Mon Sep 17 00:00:00 2001 From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr> Date: Wed, 13 Feb 2019 12:36:55 +0100 Subject: [PATCH] Philippe 13/02/2019: initialize XALBUV even if no radiation (needed in CH_INTERP_JVALUES) --- src/MNH/ini_modeln.f90 | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/MNH/ini_modeln.f90 b/src/MNH/ini_modeln.f90 index 1bbb5cf60..443ca3adf 100644 --- a/src/MNH/ini_modeln.f90 +++ b/src/MNH/ini_modeln.f90 @@ -278,7 +278,8 @@ END MODULE MODI_INI_MODEL_n !! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O !! V. Vionnet : 18/07/2017 : add blowing snow scheme !! 01/18 J.Colin Add DRAG -!! P.Wautelet 29/01/2019: bug: add missing zero-size allocations +! P. Wautelet 29/01/2019: bug: add missing zero-size allocations +! P. Wautelet 13/02/2019: initialize XALBUV even if no radiation (needed in CH_INTERP_JVALUES) !--------------------------------------------------------------------------------- ! !* 0. DECLARATIONS @@ -460,6 +461,8 @@ TYPE(TFILEDATA), INTENT(IN) :: TPINIFILE ! Initial file ! !* 0.2 declarations of local variables ! +REAL, PARAMETER :: NALBUV_DEFAULT = 0.01 ! Arbitrary low value for XALBUV +! INTEGER :: JSV ! Loop index INTEGER :: IRESP ! Return code of FM routines INTEGER :: ILUOUT ! Logical unit number of output-listing @@ -1325,6 +1328,7 @@ ALLOCATE(XLW_BANDS (NLWB_MNH)) ALLOCATE(XZENITH (IIU,IJU)) ALLOCATE(XAZIM (IIU,IJU)) ALLOCATE(XALBUV (IIU,IJU)) +XALBUV(:,:) = NALBUV_DEFAULT !Set to an arbitrary low value (XALBUV is needed in CH_INTERP_JVALUES even if no radiation) ALLOCATE(XDIRSRFSWD(IIU,IJU,NSWB_MNH)) ALLOCATE(XSCAFLASWD(IIU,IJU,NSWB_MNH)) ALLOCATE(XFLALWD (IIU,IJU)) -- GitLab