Skip to content
GitLab
Select Git revision
  • 191-speed-up-runmodel_lag-when-it-is-used-alone
  • dev default
  • master
  • 171-check-error-direct-call-of-as-data-frame-posixct-is-deprecated-use-as-data-frame-vector-or-as
  • 161-runmodel_lag-proposal-for-a-better-transformation-and-screening-of-the-parameters
  • 165-add-missing-isintstore-argument-in-createcaliboptions
  • 33-merge-airgr-and-airgrplus-functions
  • 115-add-pcict-date-format-managment
  • rmOptions
  • 136-decreased-performance-of-calibration-execution-time
  • 118-plot-outputsmodel-remove-case-sensitivity-to-the-which-argument
  • 96-parallelize-the-grid-screening-step-during-calibration
  • 60-migrate-transfoparam-functions-and-createcaliboptions-to-s3-methods
  • airGRplus
  • CRAN_v1.7.6
  • CRAN_v1.7.5
  • CRAN_v1.7.4
  • CRAN_v1.7.3
  • CRAN_v1.7.2
  • CRAN_v1.7.1
  • CRAN_v1.7.0
  • CRAN_v1.6.12
  • CRAN_v1.6.11
  • CRAN_v1.6.10.4
  • CRAN_v1.6.9.27
  • CRAN_v1.6.9.24
  • CRAN_v1.6.9.23
  • CRAN_v1.6.9.21
  • CRAN_v1.4.3.65
  • CRAN_v1.4.3.60
  • CRAN_v1.4.3.52
  • CRAN_v1.3.2.42
  • CRAN_v1.3.2.41
  • CRAN_v1.3.2.40
Find file BlameHistoryPermalink
frun_GR1A.f 4.48 KiB
1
2
3 
      SUBROUTINE frun_GR1A(
4 
                                 !inputs
5 
     &                             LInputs      , ! [integer] length of input and output series
6 
     &                             InputsPrecip , ! [double] input series of total precipitation [mm/year]
7 
     &                             InputsPE     , ! [double] input series of potential evapotranspiration (PE) [mm/year]
8 
     &                             NParam       , ! [integer] number of model parameters
9 
     &                             Param        , ! [double] parameter set
10 
     &                             NStates      , ! [integer] number of state variables
11 
     &                             StateStart   , ! [double] state variables used when the model run starts (none here)
12 
     &                             NOutputs     , ! [integer] number of output series
13 
     &                             IndOutputs   , ! [integer] indices of output series
14 
                                 !outputs
15 
     &                             Outputs      , ! [double] output series
16 
     &                             StateEnd     ) ! [double] state variables at the end of the model run (none here)
17
18
19 
      !DEC$ ATTRIBUTES DLLEXPORT :: frun_GR1A
20
21
22 
      Implicit None
23 
      !### input and output variables
24 
      integer, intent(in) :: LInputs,NParam,NStates,NOutputs
25 
      doubleprecision, dimension(LInputs)  :: InputsPrecip
26 
      doubleprecision, dimension(LInputs)  :: InputsPE
27 
      doubleprecision, dimension(NParam)   :: Param
28 
      doubleprecision, dimension(NStates)  :: StateStart
29 
      doubleprecision, dimension(NStates)  :: StateEnd
30 
      integer, dimension(NOutputs) :: IndOutputs
31 
      doubleprecision, dimension(LInputs,NOutputs) :: Outputs
32
33 
      !parameters, internal states and variables
34 
      integer NMISC
35 
      parameter (NMISC=3)
36 
      doubleprecision MISC(NMISC)
37 
      doubleprecision P0,P1,E1,Q
38 
      integer I,K
39
40 
      !--------------------------------------------------------------
41 
      !Initializations
42 
      !--------------------------------------------------------------
43
44 
      !parameter values
45 
      !Param(1) : PE adjustment factor [-]
46
47 
      !initialization of model outputs
48 
      Q = -999.999
49 
      MISC = -999.999
50 
c      Outputs = -999.999  !initialization made in R
51
52 
      StateStart = -999.999  ! CRAN-compatibility updates
53 
      StateEnd = -999.999  ! CRAN-compatibility updates
54
55
56 
      !--------------------------------------------------------------
57 
      !Time loop
58 
      !--------------------------------------------------------------
59 
      DO k=2,LInputs
60 
        P0=InputsPrecip(k-1)
61 
        P1=InputsPrecip(k)
62 
        E1=InputsPE(k)
63 
c        Q = -999.999
64 
c        MISC = -999.999
65 
        !model run on one time step
66 
        CALL MOD_GR1A(Param,P0,P1,E1,Q,MISC)
67 
        !storage of outputs
68 
        DO I=1,NOutputs
69 
        Outputs(k,I)=MISC(IndOutputs(I))
70 
        ENDDO
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124
ENDDO RETURN ENDSUBROUTINE c################################################################################################################################ C********************************************************************** SUBROUTINE MOD_GR1A(Param,P0,P1,E1,Q,MISC) C Calculation of streamflow on a single time step with the GR1A model C Inputs: C Param Vector of model parameters (Param(1) [-]) C P0 Value of rainfall during the previous time step [mm/year] C P1 Value of rainfall during the current time step [mm/year] C E1 Value of potential evapotranspiration during the current time step [mm/year] C Outputs: C Q Value of simulated flow at the catchment outlet for the current time step [mm/year] C MISC Vector of model outputs for the time step [mm/year] C********************************************************************** Implicit None INTEGER NMISC,NParam PARAMETER (NMISC=3) PARAMETER (NParam=1) DOUBLEPRECISION Param(NParam) DOUBLEPRECISION MISC(NMISC) DOUBLEPRECISION P0,P1,E1,Q DOUBLEPRECISION tt ! speed-up C Runoff ! speed-up tt = (0.7*P1+0.3*P0)/Param(1)/E1 Q=P1*(1.-1./SQRT(1.+tt*tt)) ! Q=P1*(1.-1./(1.+((0.7*P1+0.3*P0)/Param(1)/E1)**2.)**0.5) ! fin speed-up C Variables storage MISC( 1)=E1 ! PE ! [numeric] observed potential evapotranspiration [mm/year] MISC( 2)=P1 ! Precip ! [numeric] observed total precipitation [mm/year] MISC( 3)=Q ! Qsim ! [numeric] simulated outflow at catchment outlet [mm/year] ENDSUBROUTINE