Data Files¶
head.CONF¶
- A default file head.CONF is automatically copied onto inputs directory
- DO NOT MODIFY ITS FORMAT (no more space or less space !)
- The first 4 letters must be equal to the variable nrac in routine parameters.F90
- The following number is the number of the rank
- The 4 following reals define the extension of the domain
- Northern latitude (at jmax)
- Southern latitude (at jmin+1)
- Eastern longitude (at imax)
- Western longitude (at imin+1)
- The 2 following reals set the step of latitude and longitude respectively (in degrees)
- The 2 following reals set the averaged step of latitude and longitude respectively (in meter) – indicative values, not used by the model
- 1st integer = imax
- 2nd integer = jmax
- The last 4 integers are the location of the open boundaries of the following rank (child grid) in the grid of the current rank : jobcmax,jobcmin,iobcmax,iobcmin.
- The last 4 letters is the CONF name of the child grid (following rank or rank of another configuration as a unique rank is able to save 2d open bounary files for several configurations at a same time)
PS : to save several 2d obc files, duplicate the line and modify the last 4 integers and the last 4 letters
Set up output files : output.dat¶
The file output.dat for setting up of output files is located under $RDIR/CONF/CONF-CASE(/rank*)
- THREE REQUIREMENTS :
- THE MAIN OUPUT MUST BE DEFINED AT FIRST
- The user can choose l_demer=.true. or l_out_precdate=.true. (both at the same time is not allowed)
- The choice l_demer=.true. or l_out_precdate=.true. is possible just once inside the file output.dat.
Definition of parameters are given at the beginning of the output.dat file See below the list of parameters :
Champs | the name of the output file will be champs_suffix.nc, suffix being defining in paracom.txt |
06/03/2005 00:00:00 | date_startout (date of the beginning output) |
09/03/2005 00:00:00 | date_endout (date to the end output) |
1.0d0 | pasor (time step of saving) |
.true. | l_posit (.true. If horizontal extent given with indexes i & j ;.false. if given with longitude/latitudes) |
0 421 0 501 | iextmin,iextmax,jextmin,jextmax (horizontal extent) |
1 30 | kextmin,kextmax (between 1 and kmax) (vertical extent) |
1 1 1 | i-step,j-step,k-step (space step for i,j, k) |
.false. | l_filebydate (a file for each date (if.true) or all records in the same (if false) |
.false. | l_out_off (saving of 3D velocities in conservative flux form (for offline purpose) |
.false. | l_out_precdate (saving at precise dates) |
.false. | l_demer (filtering fields using Demerliac filter) |
.true. | l_out_ssh (saving of sea surface height) |
.true. | l_out_ubt (saving of barotropic (2D) velocities) |
.false. | l_out_u3d (saving of total (3D) velocities) |
.true. | l_out_sal (saving of salinities) |
.false. | l_out_temp (saving of temperature) |
.false. | l_out_wind (saving of wind stress) |
.false. | l_out_visc (saving of viscosity coefficients) |
.false. | l_out_turb (saving of turbulence coefficients) |
.false. | l_out_buoy (saving of buoyancy) |
.false. | l_out_atmsflx (saving of cumulative atmospheric fluxes) |
.false. | l_out_wave (saving of waves parameters) |
River Discharges : river.dat¶
The name of the file river.dat is set in :ref:`doc.paraspec` (file_river).
+-----------------------------------------------+-----------------------------------------------------------------+
| *********************************************| ***************************************************** |
| Seine | river name |
| .true. | location in i,j,k (.true.) or in lon,lat,k (.false.) |
| 142 45 12 | location |
| 100 | thickness (m) on which the outflow is distributed |
| .false. | .true. for a stationary river discharge |
| 450.0 | value of the stationary river discharge |
| ../../inputs/debits/deb-seine-nc.dat | river discharges data |
| .true. | set to .true. for reading the river temperature |
| ../../inputs/concentrations/temp-seine.dat | data of temp |
| 10.0 6.0 | mean temp of river (degrees) & annual amplitude (degrees) |
| *********************************************| ***************************************************** |
+-----------------------------------------------+-----------------------------------------------------------------+
- River thickness : estimated downwards from the interface of the layer k
- If the river thickness is set to 0., the outflow is completed at the layer k
- If the river thickness is higher than depth or set to -1., the outflow is distributed over the total depth
- To prescribe the river discharge at the bottom, set k=1 and river_thickness=0.
MPI job : domain decomposition¶
PREPARATION : Create a mpi.txt file if you want to divide your domain along X and Y
- Use decoupe program under $HOMEMARS/../TOOLS/MPI2D_DOMAIN (see Readme.txt)
- cp bathy_CONF.txt $RDIR/CONF/CONF-CASE[/rankX]/mpi.txt
MPI2D management :
- If none mpi.txt, imax > jmax : domain decomposition through vertical stripes
- If none mpi.txt, jmax > imax : domain decomposition through vertical stripes
- If mpi.txt exist (issued from decoupe program) : optimized 2D domain decomposition