PRO check_plasma, inp, grid, plasma
    ;+#check_plasma
    ;+Checks if plasma paramters structure is valid
    ;+***
    ;+##Input Arguments
    ;+     **inputs**: Input structure
    ;+
    ;+     **grid**: Interpolation grid structure
    ;+
    ;+     **plasma**: Plasma parameters structure
    ;+
    ;+##Example Usage
    ;+```idl
    ;+IDL> check_plasma, inputs, grid, plasma
    ;+```
    err_status=0
    info,'Checking plasma parameters...'

    nr = grid.nr
    nz = grid.nz
    w = where("nphi" eq strlowcase(TAG_NAMES(grid)),nw)
    if nw eq 0 then begin
        info,'"nphi" is missing from the plasma, assuming axisymmetry'
        nphi = 1
    endif else begin
        nphi = grid.nphi
    endelse

    nthermal = plasma.nthermal
    zero_int = {dims:0, type:'INT'}
    zero_string = {dims:0, type:'STRING'}
    zero_double = {dims:0, type:'DOUBLE'}
    nthermal_double = {dims:[nthermal], type:'DOUBLE'}
    nrnznphi_double = {dims:[nr,nz,nphi], type:'DOUBLE'}
    ntnrnznphi_double = {dims:[nthermal,nr,nz,nphi], type:'DOUBLE'}
    nrnznphi_int    = {dims:[nr,nz,nphi], type:'INT'}
    schema = {time:zero_double, $
              nthermal:zero_int,$
              impurity_charge:zero_int,$
              species_mass:nthermal_double,$
              vr:nrnznphi_double, $
              vt:nrnznphi_double, $
              vz:nrnznphi_double, $
              dene:nrnznphi_double, $
              denimp:nrnznphi_double, $
              deni:ntnrnznphi_double, $
              denn:nrnznphi_double, $
	      ti:nrnznphi_double, $
              te:nrnznphi_double, $
              zeff:nrnznphi_double, $
              mask:nrnznphi_int, $
              data_source:zero_string}

    check_struct_schema,schema,plasma,err_status,desc = "plasma parameters"
    if err_status eq 1 then begin
        goto, GET_OUT
    endif

    if plasma.data_source eq '' then begin
        error, 'Invalid data source. An empty string is not a data source.'
        err_status = 1
    endif

    ;;Electron density
    plasma.dene = plasma.dene > 0. ;[1/cm^3]

    ;;Impurity density
    plasma.denimp = plasma.denimp > 0. ;[1/cm^3]

    ;;Ion Density
    plasma.deni = plasma.deni > 0. ;[1/cm^3]

    ;;Neutral density
    plasma.denn = plasma.denn > 0. ;[1/cm^3]

    ;;Zeff
    plasma.zeff = plasma.impurity_charge < plasma.zeff > 1.0

    ;;Electron temperature
    plasma.te = plasma.te > 0. ;[keV]

    ;;Ion temperature
    plasma.ti = plasma.ti > 0. ;[keV]

    if abs(plasma.time - inp.time) gt 0.02 then begin
        warn,'Plasma time and input time do not match'
        print,'Input time: ',inp.time
        print,'Plasma time: ',plasma.time
    endif

    plasma = create_struct(plasma, grid)

    GET_OUT:
    if err_status ne 0 then begin
        error,'Invalid plasma parameters. Exiting...',/halt
    endif else begin
        success,'Plasma parameters are valid'
    endelse
END