Read standard parameter ancillary data in NetCDF (#324)

# CABLE
## Description

Fixes #281 

Aside from `climate` parameters, BIOS ancillary data is being read from
the gridinfo. They are supposed to be read via a separate BLAZE NetCDF
file #392

The original 0.5 degree gridinfo file
([`gridinfo_CSIRO_CRU05x05_4tiles.nc`](https://github.com/CABLE-LSM/bios_testing/blob/b45e0736cfd7fd1de1014319eb551b924542a96c/cable.nml#L4C63-L4C96))
went through 2 stages to make it compatible with reading to ancillary
inputs

TODO: Once these PRs are merged, point to these specific files
1. Subsampling at 0.05 degrees using script
CABLE-LSM/gists#3
2. Add BIOS coastal remapping script for gridinfo to have the same mask
for all input parameters. CABLE-LSM/gists#4

## Type of change

New feature

## Testing

Run `act9test` and `1000pts-configuration` (TODO: provide links) to
`BIOS_ancillary_netcdf` and `CABLE_POP-TRENDY` with `derived_parameters`
and `check_parameter_values` being called after `cable_load_bios_params`
instead of within `load_parameters`.

Author: abhaasgoyal

offline/cable_LUC_EXPT.F90

@@ -3,6 +3,7 @@ MODULE CABLE_LUC_EXPT
   use cable_common_module,  only: is_leapyear, leap_day, handle_err, handle_iostat,&
                                   get_dimid, LatNames, LonNames, TimeNames
   use cable_io_vars_module, only: logn, land_x, land_y, landpt, latitude, longitude
+  USE cable_common_module, ONLY: cable_user
   use cable_def_types_mod,  only: mland
 
   implicit none
@@ -67,13 +68,14 @@ MODULE CABLE_LUC_EXPT
 
   ! ------------------------------------------------------------------
 
-  SUBROUTINE LUC_EXPT_INIT(LUC_EXPT)
+  SUBROUTINE LUC_EXPT_INIT(inMVG, LUC_EXPT)
 
     use netcdf,                    only: nf90_open, nf90_nowrite, nf90_inq_varid, nf90_inq_dimid, &
          nf90_inquire_dimension, nf90_inq_varid, nf90_get_att, nf90_get_var, nf90_close
 
     implicit none
 
+    INTEGER, INTENT(IN) :: inMVG(:, :)
     type(luc_expt_type), intent(inout) :: luc_expt
 
     REAL :: tmp
@@ -379,6 +381,12 @@ SUBROUTINE LUC_EXPT_INIT(LUC_EXPT)
        PrimOnly_fID = -1
     ENDIF
 
+    IF (TRIM(cable_user%MetType) .EQ. "bios") THEN
+       DO k = 1, mland
+          LUC_EXPT%biome(k) = inMVG(landpt(k)%ilon,landpt(k)%ilat) 
+       ENDDO
+    END IF
+
     ! Determine woody fraction (forest and shrub cover).
     call get_woody_fraction(LUC_EXPT)
 
@@ -457,14 +465,12 @@ end subroutine luc_expt_zero
   subroutine get_woody_fraction(LUC_EXPT) 
     ! Determine woody fraction (forest and shrub cover) from ancillary data.
 
-    use cable_bios_met_obs_params, only: cable_bios_load_biome
     use cable_common_module,       only: cable_user
 
     implicit none
 
     type(LUC_EXPT_type), intent(inout) :: LUC_EXPT
     real    :: CPC(mland)
-    integer :: MVG(mland)
     real    :: projection_factor
 
     if (TRIM(cable_user%MetType) .EQ. "bios") then
@@ -473,10 +479,8 @@ subroutine get_woody_fraction(LUC_EXPT)
       ! Woody fraction (woodfrac) is then calculated from CPC.
 
       ! read bios parameter file to NVIS Major Vegetation Group "biomes"
-      call cable_bios_load_biome(MVG)
 
       ! adjust fraction woody cover based on Major Vegetation Group
-      LUC_EXPT%biome = MVG
       LUC_EXPT%ivegp = 2
       projection_factor = 0.65
       WHERE (LUC_EXPT%biome .eq. 1)
@@ -665,21 +669,17 @@ END SUBROUTINE LUC_EXPT_SET_TILES
   ! ------------------------------------------------------------------
 
 
-  SUBROUTINE LUC_EXPT_SET_TILES_BIOS(inVeg, inPfrac, LUC_EXPT )
-
-    USE cable_bios_met_obs_params, ONLY: cable_bios_load_fracC4
+  SUBROUTINE LUC_EXPT_SET_TILES_BIOS(inVeg, inPfrac, infracC4, LUC_EXPT )
 
     IMPLICIT NONE
 
     INTEGER,              INTENT(INOUT) :: inVeg(:,:,:)
     REAL,                 INTENT(INOUT) :: inPFrac(:,:,:)
-    TYPE (LUC_EXPT_TYPE), INTENT(INOUT) :: LUC_EXPT
+    REAL,                 INTENT(IN)    :: infracC4(:,:)
+    TYPE (LUC_EXPT_TYPE), INTENT(IN) :: LUC_EXPT
 
-    REAL :: fracC4(mland)
     INTEGER :: k, m, n
 
-    CALL cable_bios_load_fracC4(fracC4)
-
     DO k=1, mland
        m = landpt(k)%ilon
        n = landpt(k)%ilat
@@ -691,7 +691,7 @@ SUBROUTINE LUC_EXPT_SET_TILES_BIOS(inVeg, inPfrac, LUC_EXPT )
           inPFrac(m,n,2:3) = 0.0
           inPFrac(m,n,1)   = 1.0
           if ( LUC_EXPT%grass(k) .gt. 0.01 ) then
-             if (fracC4(k).gt.0.5) then
+             if (infracC4(m, n).gt.0.5) then
                 inVeg(m,n,2) = 7 ! C4 grass
              else
                 inVeg(m,n,2) = 6 ! C3 grass
@@ -704,7 +704,7 @@ SUBROUTINE LUC_EXPT_SET_TILES_BIOS(inVeg, inPfrac, LUC_EXPT )
 
           inVeg(m,n,1) = LUC_EXPT%ivegp(k)
           inVeg(m,n,2) = LUC_EXPT%ivegp(k)
-          if (fracC4(k).gt.0.5) then
+          if (infracC4(m, n).gt.0.5) then
              inVeg(m,n,3) = 7 ! C4 grass
           else
              inVeg(m,n,3) = 6 ! C3 grass