ReportNETCDF.f90

Go to the documentation of this file.

!==============================================================================
!===============================================================
SUBROUTINE reportnetcdf(year, nspc)
    use, intrinsic :: iso_fortran_env
    use, intrinsic :: ieee_arithmetic
    use all_par, only: maxlat, maxlon, writeoutput, reportintervals, simustartyear, maxlc, maxhc
    use filelistmodule, only: biodiv_nc, biomass_nc, number_nc,seed_nc, ingrowth_nc,npp_nc,basalarea_nc,&
                              antagonist_nc,pollen_nc,lai_nc,light_nc,&
                              heightstruct_nc, nc_file_id
 
    use netcdf, only: nf90_put_var
 
    IMPLICIT NONE
    
    INTEGER, INTENT(in) :: year, nspc
    INTEGER :: ReportYear = 0
    INTEGER :: i, rpti, yr
    REAL :: fillval
        
    fillval = ieee_value(fillval, ieee_quiet_nan)   ! some kind of NA values
    yr = year + simustartyear                       ! actual year
    rpti = 1
    do i=1,10  ! determine current report interval
        if (yr >= reportintervals(i,2) .and.  yr <= reportintervals(i,3) )  then
            rpti=reportintervals(i,1)
            exit
        end if 
    end do
    if (mod(year, rpti) == 0 .and. writeoutput%netcdf) then ! only if output for this year and netcdf output in general is whished 
        reportyear = reportyear+1 ! increase report year every call.
        if(writeoutput%light) then
            call m_check(nf90_put_var(nc_file_id, light_nc%varId, light_nc%values4D, &
                start = (/ reportyear, 1, 1 /),     count = (/ 1, maxhc+1, maxlc, maxlon, maxlat /))&
                ,__line__, __file__)
            light_nc%values4D = fillval
        end if
        if(writeoutput%biodiv) then
            call m_check(nf90_put_var(nc_file_id, biodiv_nc%varId, biodiv_nc%values2D, &
                start = (/ reportyear, 1, 1 /),     count = (/ 1, maxlon, maxlat /))&
                ,__line__, __file__)
            biodiv_nc%values2D = fillval
        end if
  ! as example described: 
  ! The values of biomass have been put to Biomass_NC%values3D in Report.
  ! count of them are   now put to the variable Biomass_NC%varId in the file NC_FILE_ID at the position ReportYear 
  ! it is checked whether this assignment is ok (m_check)  
  ! then the variable is set to NA        
        if(writeoutput%biomass) then
            call m_check(nf90_put_var(nc_file_id, biomass_nc%varId, biomass_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /),   count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            biomass_nc%values3D = fillval
        end if
        if(writeoutput%number) then
            call m_check(nf90_put_var(nc_file_id, number_nc%varId, number_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            number_nc%values3D = fillval
        end if
        if(writeoutput%seeds) then
            call m_check(nf90_put_var(nc_file_id, seed_nc%varId, seed_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            seed_nc%values3D = fillval
        end if
        if(writeoutput%ingrowth) then
            call m_check(nf90_put_var(nc_file_id, ingrowth_nc%varId, ingrowth_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            ingrowth_nc%values3D = fillval
        end if
        if(writeoutput%NPP) then
            call m_check(nf90_put_var(nc_file_id, npp_nc%varId, npp_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            npp_nc%values3D = fillval
        end if
        if(writeoutput%basalArea) then
            call m_check(nf90_put_var(nc_file_id, basalarea_nc%varId, basalarea_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            basalarea_nc%values3D = fillval
        end if
        if(writeoutput%antagonists) then        
            call m_check(nf90_put_var(nc_file_id, antagonist_nc%varId, antagonist_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            antagonist_nc%values3D = fillval
        end if
        if(writeoutput%pollen) then     
            call m_check(nf90_put_var(nc_file_id, pollen_nc%varId, pollen_nc%values3D, &
                start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                ,__line__, __file__)
            pollen_nc%values3D = fillval
        end if
        if(writeoutput%lai) then
                call m_check(nf90_put_var(nc_file_id, lai_nc%varId, lai_nc%values3D, &
                    start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                    ,__line__, __file__)
                lai_nc%values3D = fillval
        end if
 
        if(writeoutput%hstruct) then
            do i=1,size(heightstruct_nc)
                call m_check(nf90_put_var(nc_file_id, heightstruct_nc(i)%varId, heightstruct_nc(i)%values3D, &
                    start = (/ reportyear, 1, 1, 1 /), count = (/ 1, nspc+1, maxlon, maxlat /))&
                    ,__line__, __file__)
                heightstruct_nc(i)%values3D = fillval
            end do
        end if
    end if
end SUBROUTINE