Commit a4ffcd62 authored by Gesa Meyer's avatar Gesa Meyer Committed by Joe Melton
Browse files

Correct monthly outputs and modify job-options-template comment

parent ed90d081
Loading
Loading
Loading
Loading
+3 −2
Original line number Diff line number Diff line
@@ -85,8 +85,9 @@
                            ! useTracer -> '13C' [Not implemented yet] means the tracer is 13C and will then call a 13C fractionation scheme. 
        tracerCO2file = ''  ! Tracer CO2 file, this file needs to be 
                            ! correctly chosen for the useTracer option. It uses the transientCO2 and fixedYearCO2 switches to determine how it is read in.
        tracertimeAdjust = 0 ,! This adjusts the read in year for the tracer. E.g. if you start your run in year 1900 for readMetStartYear
                            ! and tracertimeAdjust -> -1000 then it will read in year 900 from the tracerCO2file. Default is 0.
        tracertimeAdjust = 0 ,  ! This adjusts the read in year for the tracer. Default is 0.
                                ! E.g. if you start your run in year 1900 for readMetStartYear
                                ! and tracertimeAdjust -> -1000 then it will read in year 900 from the tracerCO2file. 

    !CO2 switches:
        transientCO2 = .false. , ! Read in time varying CO2 concentration from CO2File or if set to false then use only the year of fixedYearCO2 value
+5 −3
Original line number Diff line number Diff line
@@ -595,7 +595,7 @@ contains
    implicit none

    character( * ), intent(in)  :: timeFreq
    integer :: totsteps, totyrs, i, j, m, n
    integer :: totsteps, totyrs, i, j, m, n, l
    integer :: lastDOY, numsteps, totdays, dct, lastday, k
    real :: runningDays
    logical :: leapnow
@@ -635,15 +635,17 @@ contains

      allocate(timeVect(totyrs * 12))

      l = 1
      do i = 1, runyrs
        if (c_switch%runStartYear + i-1 >= c_switch%jmosty) then 
          if (c_switch%leap) call findLeapYears(metYearOrder(i),leapnow,lastDOY)
          do m = 1, 12
            j = ((i - 1) * 12) + m
            j = ((l - 1) * 12) + m
            timeVect(j) = runningDays + real(monthend(m + 1))
          end do
          runningDays = runningDays + real(lastDOY)    
          l = l + 1  
        end if 
        runningDays = runningDays + real(lastDOY) 
      end do
      
    case ("daily")