After discussion allocatable was more appropriate than optional. Also…

… added check for variables being allocated

src/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F

@@ -606,8 +606,8 @@ subroutine readland_hrldas(wrfinput_flnm,                            &
     real,    dimension(xstart:xend,ystart:yend), intent(out) :: seaice
     real,    dimension(xstart:xend,ystart:yend), intent(out) :: msftx
     real,    dimension(xstart:xend,ystart:yend), intent(out) :: msfty
-    real,    dimension(xstart:xend,ystart:yend), intent(inout), optional :: glact
-    real,    dimension(xstart:xend,ystart:yend), intent(inout), optional :: vis_icealb
+    real,    dimension(:,:), intent(inout), allocatable :: glact
+    real,    dimension(:,:), intent(inout), allocatable :: vis_icealb
 
     character(len=256) :: units
     integer :: ierr
@@ -699,6 +699,8 @@ subroutine readland_hrldas(wrfinput_flnm,                            &
           stop "FATAL ERROR: In module_hrldas_netcdf_io.F readland_hrldas()"// &
                " - crocus_opt on but netcdf read of glacinfo returned non-zero"
        end if
+
+       if (allocated(glact) .eqv. .false.) error stop "FATAL ERROR: glact not initialized but crocus is on"
        ! Get Glacier thickness (glact)
        call get_2d_netcdf("glacier_thickness", ncid, glact,   units, xstart, xend, ystart, yend, NOT_FATAL, ierr)
        if (ierr /= 0) then
@@ -707,6 +709,7 @@ subroutine readland_hrldas(wrfinput_flnm,                            &
                ". Setting glacier thickness to 0.0"
        end if
 
+       if (allocated(vis_icealb) .eqv. .false.) error stop "FATAL ERROR: vis_icealb not initialized but crocus is on"
        ! Get visible ice albedo map (vis_icealb)
        call get_2d_netcdf("min_vis", ncid, vis_icealb,   units, xstart, xend, ystart, yend, NOT_FATAL, ierr)
        if (ierr /= 0) then