General description of the file outcdf.out
General information
First one needs to define at least one type of grid, one set of levels, and one set of steps. These will then get used in the sortie (for dynamic fields) resp. sortie_p (for physics fields) declarations.
Example outcfg.out
grid=1,model;
levels=1,eta,-1;
levels=2,pres,[1000.,850.];
steps=1,hour,0.;
steps=2,step,<0.,6.,1.>;
sortie ([ME,PN,P0,PT] ,grid, 1, levels, 1, steps, 2)
sortie_p([MF,MY,MG,ALVIS,DLAT,DLON,WFLUX] ,grid, 1, levels, 1, steps, 1)
Output GRID definitions
Examples:
grid=1,model; grid=2,core; grid=3,free;
General syntax:
grid=gridid,model/core/free;
gridid - number to identify gridset to relate to sortie statement model - total grid of the model core - in LAM: model grid minus the pilot region free - in LAM: only the free area, excluding blending (Davies sponge) and pilot area
IMPORTANT NOTE:
Limit the number of definitions for "grid" to improve the efficiency in the output routines. The maximum number of definitions is 4.
Output LEVEL definitions
Examples:
evels=1,eta,-1; levels=2,eta,0; levels=6,eta,<1,12,2>; levels=5,pres,[1000.,850.,700.,500.];
General syntax
levels=levelid ,pres/eta,[list]/<list>
levelid - number to identify levelset to relate to sortie statement eta - model levels (eta) The integer corresponds to the "nth" defined model level thus make sure that the largest digit <= # of elements defined in hyb '-1' will give all model levels '0' will give the diagnostic level (2 m for temperature & humidity; 10 m for winds) pres - pressure levels
lists formats: [a,b,c] means level a,b and c are requested <a,b,c> means levels a to b, every c are requested
IMPORTANT NOTE:
Limit the number of definitions for "level" to improve the efficiency the output routines. The maximum number of definitions is 32. The maximum number of levels defined in each "level" command is 201
Output FREQUENCY definitions
Examples:
steps=1,step,-1; steps=2,hour,0.,3.,6.,12.,24.,48.; steps=3,hour,<0.,48.,3.>; steps=4,init,hour,<0.,6.,1.>; steps=5,hour,<0.,24.,3.>,<24.,48.,6.>; steps=6,step,<0,5,1>,[6,12,18];
General syntax
steps=stepid, [init],step/hour,[list]/<list>;
stepid - number to identify stepset to relate to sortie statement 'init' - means this command only applies to output during the initialization period 'hour' - output in hours 'step' - output in timesteps in the model '-1' will give every timestep of the model
list formats: [a,b,c] means a,b and c are requested <a,b,c> means a to b, incrementing every c are requested
IMPORTANT NOTE:
Limit the number of definitions for "steps" to improve the efficiency the output routines. It can really make a huge difference!!! The maximum number of definitions is 32. The maximum number of step/hour defined in each "steps" command is 50000
Accumulators & averages/min/max fields
When running GEM there are two ways to trigger the resetting of accumulators and min/max fields. They cannot be mixed! You need to use either one way or the other for all variables.
1) Use moyhr & acchr
One can set 'moyhr' and 'acchr' in the physics_cfgs namelist of the gem_settings.nml to the time in hours at which one would like to reset all(!) average/min/max fields resp. all(!) accumulators.
2) Set flags in outcfg.out
Averages, accumulators, averages of accumulators, minimum and maximum fields can now be defined as such in the file ‘outcfg.out’.
When running the UQAM version of GEM these fields have to be reset this way.
To define one of the above fields as such one has to add the specific “type” to ‘sortie_p’. There are five different options:
average : Average (divides field by number of time steps since last output and sets it to zero after each output) accum : Accumulator (sets field to zero after each output) avgacc : Average of an accumulator (divides field by the output interval in seconds and sets it to zero after each output) min : Minimum (sets field to +1.E+36 after each output) max : Maximum (sets field to -1.E+36 after each output)
The difference between ‘accum’ and ‘avgacc’ is that with ‘avgacc’ the accumulator gets divided by the output time step Δt. The unit of the accumulators will therefore an an additional [1/s]. For example, then using ‘accum’ PR will be in [m], but when using ‘avgacc’ PR will be in '[m/s].
Examples:
sortie_p([PR,SN,TRAF] , grid, 1, levels, 1, steps, 2, avgacc) sortie_p([T5] , grid, 1, levels, 2, steps, 3, min) sortie_p([T9] , grid, 1, levels, 2, steps, 3, max) sortie_p([TSAV] , grid, 1, levels, 2, steps, 3, average) sortie_p([TJ,TT2M,TSKN] , grid, 1, levels, 2, steps, 2)
Ajouter un commentaire