FULLPOS users' guide : Quick update
from
Aladin cycle 10 & Arpege/Ifs cycle 20
to
Aladin cycle 12 & Arpege/Ifs cycle 22
The name of a subdomain can have up to 10 characters
New subdomains with default setup :
CFPDOM &
NLAT
& NLON & RLATN & RLATS & RLONW
& RLONE
'REUNION05' & 61 &
141 & -5.
& -35. & 30.
& 100.
'OCINDIEN' & 67
& 89 & 33.
& -66. & 0.
& 132.
LFPOLE option : removed
Additive variables are at disposal to control the climatology usage
:
RFPCSAB : Critical sand percentage difference between source
and model climatology for computing relative soil moisture in ISBA ; default
is 0.01
RFPCD2 : Critical soil depth difference between source and model
climatology for computing relative soil moisture in ISBA ; default is 0.001
LFPMOIS : Month allowed for climatology usage :
.FALSE. => month of the model (forecast)
.TRUE. => month of the initial file
Default is .FALSE. ; it is recommended to set LFPMOIS=.TRUE. to enable
in-line post-processing
Ozone, cloud fraction and upper air water related variables :
For ozone, cloud fraction, specific humidity, liquid water or ice,
it is possible to consider the post-processed fields either as gridpoint
or as spectral coefficients,independantly of the state of the model variable
(so it is a way to convert these fields from spectral to gridpoint or vice
versa) :
LFPGPQ : .TRUE. to compute gridpoint humidity rather than
spectral ; default = LECMWF
LFPGPL : .TRUE. to compute gridpoint liquid water rather than
spectral ; default =.TRUE.
LFPGPI : .TRUE. to compute gridpoint ice rather than spectral
; default =.TRUE.
LFPGPA : .TRUE. to compute gridpoint cloud fraction rather
than spectral ; default =.TRUE.
LFPGPO3 : .TRUE. to compute gridpoint ozone rather than spectral
; default = LECMWF
These variables are in the namelist NAMFPC.
Adiabatic post-processing :
To run the post-processing in the adiabatic model, you should (not
ask for physical fields and) carefully remove the physical fields from
the model, by setting the following variables in namelists :
/NAMPHY
LSOLV=.FALSE.,
LFGEL=.FALSE.,
LFGELS=.FALSE.,
LMPHYS=.FALSE.,
LNEBN=.FALSE.,
LREASUR=.TRUE.,
/END
/NAMDPHY
NVSO=0,
NVCLIV=0,
NVRS=0,
NVSF=0,
NVSG=0,
NCSV=0,
NVCLIN=0,
NVCLIP=0,
/END
Notice : for technical reasons, it is important to keep LREASUR=.TRUE.
NFPXFLD namelist parameter :is still valid in the distributed
memory code ;
the specificity is just that the "extraction" process consists in inter-processors
communications.
NFPROMA namelist parameter : has been renamed NFPROMAG
Packing/unpacking :
In the Distributed Memory code, the packing/unpacking of the global
fields can be distributed : in namelist NAMPAR1 :
NSTRIN : number of processors used for unpacking ; the best performances
are obtained when this variable is set like the total number of processors
used in the model (NPROC).
NSTROUT : number of processors used for packing ; the
best performances are obtained for an intermediate value (roughly NPROC/2).
Extension zone for Aladin :
Once you have requested the extension zone, you may set the variable
of the namelist NAMFPEZO :
NSTREFP : number of processors used to compute the extension
zone ; default is 1 (increasing this parameter seems not to bring benefits).
NFPROMEL : Vector length in the computation of the extension
zone ; default is odd and is the maximum possible value.
Notice : the computation of the extension zone is possible only
with Aladin library.
Customized complexion of certain fields :
LFPRH100 : .TRUE. to convert relative humidity in percent ;
default = LECMWF
LFPLOSP : is a variable to handle surface pressure (Ps)
or its logarithm, taking also into account the value of NLOSP :
LFPLOSP =.TRUE. =>
Fill Ps array with Ps or
Ln(Ps), according to NLOSP ;
Fill Ln(Ps) array always
with Ln(Ps)
LFPLOSP =FALSE. =>
Fill Ps array always with
Ps ;
Fill Ln(Ps) array with Ps
or Ln(Ps), according to NLOSP
If (LECMWF) then the default value is .FALSE. ; Else, the default value
is .TRUE. except for the configuration 927 where
it is set internally.
These variables are in the namelist NAMFPC.
There is a new (and post-processable) surface field in historical file : 'SURFRESERV.GLACE',
BOGUSSING :
It is possible to translate a cyclone inside an Aladin file. This procedure
will translate the cyclone to a geographical point given in namelist. The
following variables should be set in the namelist NAMFPG :
NFPTTYP : set NFPTTYP=2 to activate the translation if the output grid
is not global
FPMUCEN : In Aladin : sine of the latitude of the observed cyclone
; default value is sin(ELAT0)
FPLOCEN : In Aladin : longitude of the observed cyclone (in degrees)
; default value is ELON0
Notice : to define the translation vector, the cyclone seen by the
model is located through the minimum of surface pressure in the whole model
area (C+I).
Usually, this procedure is used to re-build the upper-air fields of
an historical file Aladin (using the configuration EE927 in adiabatic mode).
In order not to translate the orography, one should first lower the orography
to zero, then translate, and finally re-set the original orography.
It is possible to "inject" data out of Aladin post-processing inside
a file Arpege : this is called "the bogussing configuration", or
"the configuration 927E" because it works like the configuration
E927 but in a "reverse" way.
To use this configuration, you should use a namelist of configuration
927, but invoking "aladin" instead "arpifs", and set in namelist NAMFPC
the new variable NFPINCR=1 (default value is 0). Furthermore, you have
to provide 3 input files instead of one :
ICMSH${CNMEXP(1:4)}INIT : Aladin bogussed coupling file
ELSCF${CNMEXP(1:4)}ALBC : Aladin coupling file (before bogussing)
BGPX${CNMEXP(1:4)}${CFPDOM} : Arpege background file
Notice : the Arpege background file should contain all the fields
of a normal historical file in gridpoint. Furthermore : in this file surface
pressure should be the true surface pressure and not its logarithm.
Important notice : for the time being, this configuration "927E"
works only in adiabatic mode (i.e. : it is not yet possible to inject physical
surface fields except orography).
New 3D dynamic field available in post-processing :
ozone mixing ratio : name = CNO3MX ; default value ='OZONE' ;
number of bits for packing = N(B/S)O3MX
New surface dynamic field available in post-processing :
Log. of surface pressure : name = CNLNSP ; default value = 'LOG.SURF.PRESS
' ; number of bits for packing = N(B/S)LNSP
New physical surface fields available in post-processing :
Frozen superficial soil wetness : name = CNFSSW ; default value = 'SURFRESERV.GLACE'
;
number of bits for packing = NBFSSW
Frozen deep soil wetness
: name = CNFDSW ; default value = 'PROFRESERV.GLACE' ;
number of bits for packing = NBFDSW
New cumulated fluxes available in post-processing :
Surface downward moon radiation : name = CNCSMR ; default value
= 'SURFRAYT.LUNE.DE' ;
number of bits for packing =NBCSMR
BEWARE of cloudiness CNXC='SURFCUMUL NEBUL ' :
This is the cloudiness out of the subroutine "Araneb", at the level
the closest to 550 hPa !
Redundant fields to be removed from post-processing very soon :
U at bottom level
: CNXUBL='SURFU NIVBAS'
V at bottom level
: CNXVBL='SURFV NIVBAS'
Temperature at bottom level
: CNXTBL ='SURFT NIVBAS'
Specific humidity at bottom level : CNXQBL='SURFQ
NIVBAS'
Geopotential at bottom level
: CNXGBL='SURFGEOP NIVBAS'
Surface temperature
: CNXST='SURFTEMPE SURF'
Deep Soil temperature
: CNXDT='SURFTEMPE PROF'
Surface Water content
: CNXSW='SURFRESER SURF'
Deep Soil Water content
: CNXDW='SURFRESER PROF'
Snow mass
: CNXSNS='SURFNEIGE'
NAMDYN :
This namelist contains variables related to the dynamics of the model,
that may be useful for the post-processing :
NLOSP : To use surface pressure or its logarithm as the primitive
variable
NDLNPR : special discretization of delta(Log(P)) for non-hydrostatism.