SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_assim.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 ! ##################
6  MODULE modd_assim
7 ! ##################
8 !
9 !!**** *MODD_ASSIM - declaration of keys for assimilation schemes
10 !!
11 !! PURPOSE
12 !! -------
13 !
14 !!
15 !!** IMPLICIT ARGUMENTS
16 !! ------------------
17 !! None
18 !!
19 !! REFERENCE
20 !! ---------
21 !!
22 !! AUTHOR
23 !! ------
24 !! L. Jarlan *Meteo France*
25 !!
26 !! MODIFICATIONS
27 !! -------------
28 !! Original 23/02/05
29 !
30 ! Inclusion of OI constants 21/05/09 (J.-F. Mahfouf)
31 !! Add all assim keys 04/2012 T.Aspelien
32 !
33 !* 0. DECLARATIONS
34 ! ------------
35 !
36 IMPLICIT NONE
37 
38 !-------------------------------------------------------------------------------
39 !
40 ! Assimilation Scheme Options:
41 !
42  LOGICAL :: LASSIM ! Assimilation or not
43  !'.TRUE.'
44  !'.FALSE.'
45  LOGICAL :: LREAD_ALL = .TRUE.
46  LOGICAL :: LAROME ! If reading AROME fields
47  LOGICAL :: LECSST ! Use ECMWF SST
48  LOGICAL :: LAESST ! SST analysis performed
49  LOGICAL :: LAESNM ! Update snow analysis
50  LOGICAL :: LALADSURF
51  LOGICAL :: LREAD_SST_FROM_FILE ! Read SST from file
52  CHARACTER(LEN=6) :: CFILE_FORMAT_SST ! Format of the SST file ASCII/FA
53  CHARACTER(LEN=6) :: CFILE_FORMAT_FG ! Format of the first guess file ASCII/FA (OI)
54  CHARACTER(LEN=6) :: CFILE_FORMAT_LSM ! Format of the LSM file ASCII/FA (OI/extrapolations)
55  CHARACTER(LEN=6) :: CFILE_FORMAT_OBS ! Format of the observations file ASCII/FA
56  CHARACTER(LEN=6) :: CFILE_FORMAT_CLIM ! Format of the climate file ASCII/FA (OI)
57  LOGICAL :: LWATERTG2 ! Use deep soil temperature as lake temp.
58  LOGICAL :: LEXTRAP_SEA ! Extrapolation of sea points
59  LOGICAL :: LEXTRAP_WATER ! Extrapolation of inland water points
60  LOGICAL :: LEXTRAP_NATURE ! Extrapolation of nature points
61  LOGICAL :: LPRT ! Running VARASSIM in a perturbation mode
62  LOGICAL :: LSIM ! Running VARASSIM in a perturbation mode
63  LOGICAL :: LBEV ! Running VARASSIM to evolve B
64  LOGICAL :: LBFIXED
65  LOGICAL :: LOBSHEADER
66  LOGICAL :: LOBSNAT
67 
68  INTEGER, PARAMETER :: NOBSMAX = 5
69  INTEGER, PARAMETER :: NVARMAX = 5
70  INTEGER,DIMENSION(NOBSMAX) :: NNCO ! Select the type of observations to be assimilated
71  INTEGER,DIMENSION(NVARMAX) :: NNCV ! Select the control variables to be used
72  INTEGER :: NOBSTYPE
73  INTEGER :: NOBS
74  INTEGER :: NIPERT
75  INTEGER :: NIFIC
76  INTEGER :: NIVAR ! counter for ctnrl vars
77  INTEGER :: NVAR ! number of cntrl vars
78  INTEGER :: NBOUTPUT
79  INTEGER :: NPRINTLEV ! Verbosity
80 
81  CHARACTER(LEN=3),DIMENSION(NVARMAX) :: CVAR_M ! X is ctrl
82  ! 'PLUS ' (default)
83  ! 'AVERA'
84  ! '2DVAR'
85  CHARACTER(LEN=100),DIMENSION(NVARMAX) :: CPREFIX_M ! The prefix of the control variables (in PREP.txt file) (max dim)
86  CHARACTER(LEN=10),DIMENSION(:), ALLOCATABLE :: COBS ! Identifier for simulated observations
87  CHARACTER(LEN=3),DIMENSION(:), ALLOCATABLE :: CVAR ! Identifier for control variable
88  CHARACTER(LEN=12) :: CBIO ! Name of Biomass variable
89  CHARACTER(LEN=100) :: CPREFIX_BIO ! The prefix of the Biomass variable
90  CHARACTER(LEN=5) :: CASSIM_ISBA ! OI/EKF
91  CHARACTER(LEN=5) :: CASSIM ! type of correction
92 
93  REAL,DIMENSION(12) :: XALPH
94  REAL,DIMENSION(NVARMAX) :: XTPRT_M ! The perturbation amplitude (max dim)
95  REAL,DIMENSION(NVARMAX) :: XSIGMA_M ! covariance of background errors if B is fixed (max dim)
96 ! ! covariance of model errors if B evolving (max dim)
97  REAL,DIMENSION(NOBSMAX) :: XERROBS_M ! Observational standard deviation
98  REAL,DIMENSION(NOBSMAX) :: XQCOBS_M
99  REAL,DIMENSION(:,:,:,:),ALLOCATABLE :: XF_PATCH ! vector of model observations (for each pacth)
100  REAL,DIMENSION(:,:,:,:),ALLOCATABLE :: XF ! Vector of forecast control variables
101  REAL,DIMENSION(:,:,:),ALLOCATABLE :: XI
102  REAL,DIMENSION(:,:), ALLOCATABLE :: XYO ! vector of observations
103  REAL,DIMENSION(:,:),ALLOCATABLE :: XLAI_PASS
104  REAL,DIMENSION(:,:),ALLOCATABLE :: XBIO_PASS
105  REAL,DIMENSION(:,:),ALLOCATABLE :: XAT2M_ISBA
106  REAL,DIMENSION(:,:),ALLOCATABLE :: XAHU2M_ISBA
107  REAL,DIMENSION(:,:),ALLOCATABLE :: XAZON10M_ISBA
108  REAL,DIMENSION(:,:),ALLOCATABLE :: XAMER10M_ISBA
109  REAL,DIMENSION(:),ALLOCATABLE :: XAT2M_TEB
110  REAL,DIMENSION(:),ALLOCATABLE :: XTPRT ! The perturbation amplitude
111  REAL,DIMENSION(:),ALLOCATABLE :: XSIGMA ! covariance of background errors if B is fixed
112  ! covariance of model errors if B evolving
113  REAL,DIMENSION(:),ALLOCATABLE :: XERROBS
114  REAL,DIMENSION(:),ALLOCATABLE :: XQCOBS
115  REAL :: XSCALE_Q ! scaling factor of Q matrix w.r.t. the initial B
116  REAL :: XSCALE_QLAI
117 !
118 INTEGER :: NENS
119 INTEGER :: NIE
120 REAL :: XASSIM_WINH
121 REAL, DIMENSION(NVARMAX) :: XINFL_M
122 REAL, DIMENSION(NVARMAX) :: XADDINFL_M
123 REAL, DIMENSION(NVARMAX) :: XADDTIMECORR_M
124 REAL, DIMENSION(:), ALLOCATABLE :: XINFL
125 REAL, DIMENSION(:), ALLOCATABLE :: XADDINFL
126 REAL, DIMENSION(:), ALLOCATABLE :: XADDTIMECORR
127 LOGICAL :: LENKF
128 LOGICAL :: LDENKF
129 LOGICAL :: LENS_GEN
130 LOGICAL :: LPB_CORRELATIONS
131 LOGICAL :: LPERTURBATION_RUN
132 LOGICAL :: LBIAS_CORRECTION
133 !
134 ! Constants and options of the soil OI analysis
135 !
136  LOGICAL :: LHUMID, LIMVEG, LISSEW, L_SM_WP, LFGEL, LCLIM, LLDHMT, &
137  LOBSWG, LOBS2M
138  INTEGER :: NMINDJ, NNEBUL, NNEIGT, NNEIGW, NR_SM_WP, NECHGU, NTVGLA, &
139  NSEAICE, NLISSEW, NIDJ, NITRAD
140  REAL :: XANEBUL, XRCLIMN, XRCLIMTP, XRCLIMTS, XRCLIMV, XRCLIMWP, XRCLIMWS, &
141  XSCOEFH, XSCOEFT, XSEVAP, XSIGH2MO, XSIGT2MO, XSNEIGT, XSNEIGW, &
142  XSPRECIP, XSWFC, XV10MX, XRD1, XRTINER, XWCRIN, XWPMX, &
143  XWSMX, XTMERGL, XRZHZ0G, XRCLIMCA, XRCLISST, XRWPIA, XRWPIB, &
144  XRSNSA, XRSNSB, XSALBM, XSALBB, XSEMIB, XSZZ0B, XSMU0, &
145  XSICE, XSEMIM, XRA_SM_WP, XRSCALDW, XSPRECIP2, &
146  XREPSM, XRCDTR, XSIGHP1, XSIGT2MR, XSIGH2MR, XRSABR, &
147  XRARGR, XGWFC, XEWFC, XGWWILT, XEWWILT, XG1WSAT, XG2WSAT, &
148  XREPS1, XREPS2, XREPS3, XADWR, XSODELX(0:9), &
149  XSIGWGO, XSIGWGB, XSIGW2B, XRTHR_QC, XSIGWGO_MAX, XRSCAL_JAC
150 !
151 END MODULE modd_assim