:============================================================================= : ENREGISTREMENT DES PARAMETRES SOUS FORME DE NAMELIST : INPUT PARAMETERS NAMELIST :============================================================================= : : namdate3d : date du debut du 3D et date du debut du zoom Agrif : namdate3d : 3D computation starting date and AGRIF nest starting date : : nambathy : parametres de bathymetrie : nambathy : bathymetry related parameters : : namrestart : gestion de la reprise : namrestart : restart parameters : : namsaverestart : gestion de la sauvegarde : namsaverestart : saving parameters : : namoutput : gestion des sorties de resultats : namoutput : output parameters : : namvisc : parametre de viscosite numerique horizontale : namvisc : numerical horizontal viscosity parameters : : namturb : parametre de turbulence sur la verticale (3D) : namturb : vertical turbulence parameters (3D) : : namgridsig : repartition des sigmas (3D) : specification de grille speciale : namgridsig : sigma coordinate distribution (for 3D computations) : special grid specification (rotated) : : namriver : parametres des fleuves : namriver : river related parameters : : namparanum : parametres numeriques : namparanum : numerical parameters : : namthermo : parametres de thermodynamique de l eau de mer : namthermo : sea water thermodynamics parameters : : namtraj : parametres relatifs au calcul de trajectoires : namtraj : trajectory computation-related parameters : : namdiag : controle des diagnostics : namdiag : diagnostics control : : namobc : type de conditions aux limites en niveau : namobc : level boundary conditions type : :----------------------------------------------------------------------------- : namdate3d : dates 3d et zoom : namdate3d : 3D and zoom dates :----------------------------------------------------------------------------- : : date_start3d : date du debut de la simulation 3D (chaine de 19 caracteres : JJ/MM/AAAA hh:mm:ss). Date posterieure ou simultanee a : date_start. Parametre 3D. Type fortran character(len=19). : date_start3d : 3d run starting date (19 character string : DD/MM/AAAA hh:mm:ss). Equal to or later than date_start. : 3D Parameter. Type fortran character(len=19). : : date_startagrif : date du debut du zoom AGRIF (chaine de 19 caracteres : JJ/MM/AAAA hh:mm:ss). Type fortran character(len=19). : date_startagrif : AGRIF zoom starting date (19 character string : DD/MM/AAAA hh:mm:ss). Type fortran character(len=19). : : date_startshift : date_startshiftage en jour du depart du modele par rapport a date_start : pour eliminer la phase de spin-up de la fenetre : date_start-date_end. Type fortran real. : date_startshift : time lag (in days) between the run starting time and date_start : to avoid any spin-up time within the date_start-date_end window : (date_startshift > 0 : the run actually starts date_startshift days after date_start). : (date_startshift < 0 : the run actually starts date_startshift days before date_start). : Type fortran real. : :----------------------------------------------------------------------------- : &namdate3d date_start3d = '04/01/2015 23:08:00' date_startagrif = '25/01/1998 11:12:00' date_startshift = 0.0 / : :----------------------------------------------------------------------------- : nambathy : parametres de bathymetrie. : nambathy : bathymetry-related parameters. :------------------------------------------------------------------------------ : : file_bathy : nom du fichier de topographie et niveau moyen en format netCDF. : Type fortran character(len=lchain) : file_bathy : topography and mean level file name in netCDF format. : Type fortran character(len=lchain) : : l_bathy_meanlev : logique positionne a .true. si la bathymetrie est referencee par : rapport au zero hydrographique et pour activer la lecture de correction : afin de ramener les sondes au niveau moyen. Ces corrections ne concernent : a priori que les zones macrotidales. Type fortran logical. : l_bathy_meanlev : boolean set to .true. if the reference for bathymetric levels is : the level of lowest astronomical tides (hydrographic zero). In : that case, a correction is applied in order to change the vertical : reference to mean sea level. This option usually only : applies to macrotidal areas. Type fortran logical. : : mslshift : cote du niveau moyen (i.e. profondeur autour de laquelle oscille la : mer). A n utiliser que si le reseau de nivellement utilise n est pas : le zero hydrographique (ex: reseau IGN 69, NGF ...) : Type fortran real(kind=rsh). : mslshift : mean sea level shift (i.e. level about which the sea level oscillates) : Useful only when the model vertical reference derives from : a preexisting levelling network (IGN 69, NGF ...) : Type fortran real(kind=rsh). : : name_hx : nom de la variable hx dans le fichier bathymetrique : name_hx : name of hx variable in bathymetry file : : name_hy : nom de la variable hy dans le fichier bathymetrique : name_hy : name of hy variable in bathymetry file : : name_h0 : nom de la variable h0 dans le fichier bathymetrique : name_h0 : name of h0 variable in bathymetry file : : name_bathy_meanlev : nom de la variable "niveau moyen" dans le fichier bathymetrique : name_bathy_meanlev : name of "mean water level" variable in bathymetry file : : bathy_nbsmooth : nombre de boucles de lissage de la : bathymetrie aux endroits de forts gradients de pente. : Type fortran integer. : bathy_nbsmooth : number of bathymetry smoothing iterations where bathymetric : slope is steep. : Type fortran integer. : : l_bathy_comb : elimination des mailles fermees par trois cotes : (contacter V. Garnier si volonte d utiliser cette fonctionalite) : l_bathy_comb : suppression of cells closed on 3 sides : (contact V. Garnier if you want to use this option) : : l_bathy_pit : elimination des cuvettes (h0 < hx,hy) : (contacter V. Garnier si volonte d utiliser cette fonctionalite) : l_bathy_pit : suppression of local pits (h0 < hx,hy) : (contact V. Garnier if you want to use this option) : : hminim : hauteur en metre en deca de laquelle on considere que l on : est sur une maille de terre donc jamais recouverte. La : profondeur en ces mailles est ensuite abaissee dans la : subroutine bathy a -valmanq. Type fortran real(kind=rsh). : hminim : height in meters above which a cell is considered as always dry. : Cells with a depth inferior to hminim (i.e. above this level) : are prescribed to the valmanq default value. Reminder: water depths : are positive downwards. : Type fortran real(kind=rsh). : : hrdef : profondeur en deca de laquelle on reestime hx et hx comme : la moyenne spatiale de h0. : hrdef : depth below which hx and hx are re-computed as the spatial : average of surrounding h0s. : : h0fond : fond de maille residuel (metre). Ce fond de maille est ajoute : pour eviter l assechement complet d une maille et eviter : par exemple dans une maille de rejet la perte de matiere. : h0fond : residual water thickness (in m). This thickness is : added to the computed water level (xe) to avoid possible mass : loss during drying (due to the fact that the scheme for water : continuity is not positive) : : dhj : hauteur minimale d ecoulement sur les bancs decouvrants : exprimee m. Ce parametre est destine a controler les trop : fortes accelerations liees a la reduction (pratiquement a zero) : de la section d ecoulement lors du decouvrement. : dhj : minimum water level on tidal flats (in m). This parameter : prevents irrealistic accelerations induced by the decrease (down : to almost zero) of the flowing section during drying. : : hminkxky : hauteur minimale de la colonne d eau utilisee pour : l evaluation de la diffusion horizontale. : hminkxky : minimum water level used to compute horizontal diffusivity : : l_tssubadv_over_smalldepth : si vrai alors on advecte en 2D les traceurs pour les fonds : inferieurs a hm : si faux, aucune advection pour les fonds inferieurs a hm : l_tssubadv_over_smalldepth : if .TRUE., tracers are advected in 2D if water level below hm : if .FALSE., no advection if water level below hm : : hm : hauteur d eau en deca de laquelle on ne fait plus du : 3D mais du 2D. : hm : water level below which computations are 2D and not 3D. : : fwetmin : Faire tendre fwetmin vers 0 (0.01) pour ameliorer le realisme : du decouvrement en prenant en compte la fraction de cellule qui est : reellement mouillee. : Cette fonctionalite est non utilisee si fwetmin=1.0 : : fwetmin : As fwetmin tends to 0 (0.01), the real fraction of the wet part : of the cell is estimated and used in the continuity, momentum : and tracer equations to improve the realism of wet-drying. : This functionality is unplugged if fwetmin=1.0 : : l_bathy_save : sauvegarde de la bathymetrie (H0,HX,HY) reellement utilisee : lors de la simulation dans le fichier bathy_verif.nc : : l_bathy_save : saving of the bathymetry (H0,HX,HY) used during the : simulation into the bathy_verif.nc file : : file_bathy_child : nom du fichier de topographie de la grille fille (rank avec grille rotationnee uniquement) : (with -Dkey_tide_saveobct) Type fortran character(len=lchain) : file_bathy_child : topography and mean level file name of the child grid in netCDF format. : Usefull for rotated child grid only (with -Dkey_tide_saveobct). Type fortran character(len=lchain) : :----------------------------------------------------------------------------- : &nambathy file_bathy='../inputs/bathy_new_16nov18.nc' l_bathy_meanlev=.false. mslshift=-0.1 !0.0 name_bathy_hx='HX' ! Names of variables in bathymetric files name_bathy_hy='HY' name_bathy_h0='H0' name_bathy_meanlev='niv_moy' bathy_nbsmooth=0 l_bathy_comb=.false. l_bathy_pit=.false. hminim=-7.0 hrdef=30.0 h0fond=0.05 dhj=0.1 hminkxky=0.5 l_tssubadv_over_smalldepth=.true. hm=1.0 fwetmin=1.0 l_bathy_save=.true. file_bathy_child='./' / : :---------------------------------------------------------------------------- : namrestart : gestion de la reprise d une simulation interrompue : namrestart : restart of an interrupted simulation :----------------------------------------------------------------------------- : : l_initfromfile : logique positionne a .true. en cas de reprise. La simulation : sera initialisee a l aide de l etat du modele qui a ete : sauvegarde au cours d une simulation precedente dans le fichier : file_init. Type fortran logical. : l_initfromfile : boolean set to .true. in case of a restart. The new run will : be initialized with parameters/fields saved during a previous : run in the file_init file. Type fortran logical. : : file_init : nom du fichier dans lequel a ete sauve l etat du modele a partir : duquel on souhaite totalement ou partiellement (voir les parametres : suivants) demarrer la simulation. Ce fichier est au format : netcdf. Type fortran character(len=lchain). : file_init : file in which parameters/fields from a previous run are saved : to create the initial conditions of a restart run (see below : which parameters/fields are saved). : netcdf format. Type fortran character(len=lchain). : : l_init_restart : logique positionne a .true. permettre la reproductibilite : des resultats si interruption du run et reprise depuis un : restart file. : l_init_restart boolean set to .true. to get exact results when restarting from : a restart file after a break in the simulation. : : l_init_rtime : logique positionne a .true. pour la reprise de la date (en secondes) : du debut de simulation 3D : l_init_rtime : boolean set to .true. if the starting time of the 3D : simulation is read from the restart file. : : l_init_rdt : logique positionne a .true. pour la reprise du pas de temps : l_init_rdt : boolean set to .true. if the time step is read from the restart file. : : l_init_rssh : logique positionne a .true. pour la reprise des champs : de niveaux de la surface libre : l_init_rssh : boolean set to .true. if the water level is read from the restart file. : : l_init_rbtvel : logique positionne a .true. pour la reprise des champs : de courants barotropes : l_init_rbtvel : boolean set to .true. if barotropic currents are read from : the restart file. : : l_init_r3dvel : logique positionne a .true. pour la reprise des champs : de courants 3D (barotropes+baroclines) : l_init_r3dvel : boolean set to .true. if 3D (barotropic+baroclinic) currents are : read from the restart file. : : l_init_rwz : logique positionne a .true. pour la reprise des champs : de vitesses verticales : l_init_rwz : boolean set to .true. if vertical velocities are read from restarts. : : l_init_rturb : logique positionne a .true. pour la reprise des champs : de viscosite dans le cas de l utilisation d une procedure : de type Smagorinski, et des champs de diffusivite et : de viscosite sur la verticale. : l_init_rturb : boolean set to .true. if viscosity fields (in case of a : Smagorinsky procedure) and vertical diffusivity and viscosity : fields are saved for restarts. : : l_init_rsal : logique positionne a .true. pour la reprise des champs : de salinites. : l_init_rsal : boolean set to .true. if salinity fields are saved for restarts. : : l_init_rtemp : logique positionne a .true. pour la reprise des champs : de temperatures. : l_init_rtemp : boolean set to .true. if temperature fields are saved for restarts. : :----------------------------------------------------------------------------- : &namrestart l_initfromfile = .true. file_init = 'save_2671230_VLtot_Ero_Seine_wssal_20151227000023.nc' l_init_restart=.false. l_init_rtime=.false. l_init_rdt=.false. l_init_rssh=.true. l_init_rbtvel=.true. l_init_r3dvel=.true. l_init_rwz=.true. l_init_rturb=.true. l_init_rsal=.true. l_init_rtemp=.true. / : :------------------------------------------------------------------------------ : namobc : type de conditions aux limites en vitesse. : namobc : velocity boundary conditions :------------------------------------------------------------------------------ : : l_obc_mars : logique positionne a .true. pour appliquer un gradient nul de : vitesse aux frontieres ouvertes. : l_obc_mars : boolean set to .true. to apply zero gradient of velocity at the : open boundaries. : : l_obc_diri : logique positionne a .true. pour utiliser la methode de Dirichlet : aux frontieres ouvertes. : l_obc_diri : boolean set to .true. to use Dirichlet s method at the : open boundaries. : : l_obc_char : logique positionne a .true. pour utiliser la methode des : characteristiques aux frontieres ouvertes. : l_obc_char : boolean set to .true. to use the characteristics method at the : open boundaries. : : l_obc_cycl : logique positionne a .true. pour appliquer des conditions aux : limites cycliques. : l_obc_cycl : boolean set to .true. to apply cyclic open boundaries. : : l_obc_cycl_x : si .true. conditions aux limites periodiques selon l axe des x : l_obc_cycl_x : boolean set to .true. for periodic boundary conditions along x : l_obc_cycl must be set to .true. : : l_obc_cycl_y : si .true. conditions aux limites periodiques selon l axe des y : l_obc_cycl_y : boolean set to .true. for periodic boundary conditions along y : : l_obc_tide : logique positionne a .true. pour simuler le phenomene de maree : et l appliquer aux frontieres ouvertes. Le signal de maree provient : - soit d une recomposition harmonique, : - soit du fichier de conditions aux limites sauve par le rang precedent. : : l_obc_tide : boolean set to .true. to simulate tides and apply then at open boundaries. : The tidal signal comes from : - an harmonic tidal composition, : - or the open boundary file saved by the previous rank. : : l_obc_ogcm : logique positionne a .true. pour appliquer des conditions aux : limites issues d un modele de plus grande echelle ne simulant pas : le phenomene de maree. : l_obc_ogcm : boolean set to .true. to apply open boundaries provived by an : ocean global model that does not simulate the tides. : : l_tide_harmcompo : positionne a .true. pour imposer aux limites une composition : harmonique de maree pour le forcage en niveau. Pratique pour : faire des modeles de detail qui sont habituellement emboites : et que l on peut "desemboiter" ainsi simplement. : Type fortran logical. : l_tide_harmcompo : boolean set to .true. to use harmonic tidal composition for : water levels. Convenient for detailed models which are usually : nested and can be easily un-nested. Type fortran logical. : l_obc_tide must be set to .true. : : l_obc2drank_save : logique positionne a .true. pour ecrire des conditions aux limites : pour un rang suivant. Type fortran logical. : l_obc2drank_save : boolean set to .true. in order to write boundary conditions for : nested computations. Type fortran logical. : : l_obc_south : logique positionne a .true. si presence de frontiere ouverte au sud. : Type fortran logical. : l_obc_south : boolean set to .true. if open boundaries to the south. : Type fortran logical. : : file_obc_ogcm_s : name of file for southern open boundary conditions : : l_obc_north : logique positionne a .true. si presence de frontiere ouverte au nord. : Type fortran logical. : l_obc_north : boolean set to .true. if open boundaries to the north. : Type fortran logical. : : file_obc_ogcm_n : name of file for northern open boundary conditions : : l_obc_east : logique positionne a .true. si presence de frontiere ouverte a l est. : Type fortran logical. : l_obc_east : boolean set to .true. if open boundaries to the east. : Type fortran logical. : : file_obc_ogcm_e : name of file for eastern open boundary conditions : file_obc_ogcm_e : name of file for eastern open boundary conditions : : l_obc_west : logique positionne a .true. si presence de frontiere ouverte a l ouest. : Type fortran logical. : l_obc_west : boolean set to .true. if open boundaries to the west. : Type fortran logical. : : file_obc_ogcm_w : name of file for western open boundary conditions : : l_obc_ogcm_rssh : logique positionne a .true. si la condition aux limites pour la : hauteur de surface de la mer doit etre lue depuis un fichier de : conditions aux limites (defini ci-dessus). : l_obc_ogcm_rssh : boolean set to .true. if the sea surface height boundary condition : has to be read from an open boundary file (defined above) : : l_obc_ogcm_rt : logique positionne a .true. si la condition aux limites pour la : temperature doit etre lue depuis un fichier : l_obc_ogcm_rt : boolean set to .true. if the temperature boundary condition : has to be read from an open boundary file (defined above) : : l_obc_ogcm_rs : logique positionne a .true. si la condition aux limites pour la : salinite doit etre lue depuis un fichier de : conditions aux limites (defini ci-dessus). : l_obc_ogcm_rs : boolean set to .true. if the salinity boundary condition : has to be read from an open boundary file (defined above) : : l_obc_ogcm_ruv : logique positionne a .true. si la condition aux limites pour les : vitesses 2D doit etre lue depuis un fichier de : conditions aux limites (defini ci-dessus). : l_obc_ogcm_ruv : boolean set to .true. if the 2D velocity boundary condition : has to be read from an open boundary file (defined above) : : l_obc_ogcm_ruvz : logique positionne a .true. si la condition aux limites pour les : vitesses 3D doit etre lue depuis un fichier de : conditions aux limites (defini ci-dessus). : l_obc_ogcm_ruv : boolean set to .true. if the 3D velocity boundary condition : has to be read from an open boundary file (defined above) : : l_obc_ogcm_relax : relaxation vers T et S "parent" le long de la frontiere ouverte : (a l interieur du domaine). Le coefficient obc_coefrel doit etre different de 0 (ex : 0.7e-06). : l_obc_ogcm_relax : boolean set to .true. if there is a T and S relaxation towards : their values along the open boundary (inside the calculated domain) : from the previous nest. The obc_coefrel coefficient must be different from 0 (ex : 0.7e-06). : : obc_coefrel : coefficient de relaxation en s^-1. : Si egal a 0, les valeurs de T, S du run "parent" sont prescrites a la frontiere ouverte. : Autrement, relaxation vers T et S "parent" a la frontiere ouverte. : obc_coefrel : relaxation coefficient in s^-1. : If set to 0, the T,S values of the previous nest are prescribed. : Otherwise, there is a relaxation towards the T,S values of the previous nest at the open boundary. : :------------------------------------------------------------------------------ : &namobc l_obc_mars=.true. l_obc_diri=.false. l_obc_char=.false. l_obc_cycl=.false. l_obc_cycl_x=.false. l_obc_cycl_y=.false. l_obc_tide=.true. l_obc_ogcm=.true. l_tide_harmcompo=.true. l_obc2drank_save=.false. l_obc_south=.true. file_obc_ogcm_s='/home/datawork-dhysed-s/PROJETS/AMORAD/FORCAGES/NEW_BATHY/gironde_obc_south_2014_2017.nc' l_obc_north=.false. file_obc_ogcm_n='./' l_obc_west=.true. file_obc_ogcm_w='/home/datawork-dhysed-s/PROJETS/AMORAD/FORCAGES/NEW_BATHY/gironde_obc_west_2014_2017.nc' l_obc_east=.false. file_obc_ogcm_e='./' l_obc_ogcm_rssh=.true. l_obc_ogcm_rt=.false. l_obc_ogcm_rs=.false. l_obc_ogcm_ruv=.false. l_obc_ogcm_ruvz=.false. l_obc_ogcm_relax=.false. obc_coefrel=0.7e-06 / : :----------------------------------------------------------------------------- : nammeteo : parametrisation de la meteo. : meteorological-related parameters :------------------------------------------------------------------------------ : : imeteo_dragtype : type of surface drag coefficient : = 0 -> constant value = cds if imeteo_exchtype=0 : = 1 -> Large and Pond (1981) (wind dependence) : = 2 -> Smith and Banke (1975) (wind dependence) : = 3 -> Geernaert et al. (1986) (wind dependence) : = 4 -> Charnock s relation (1955) (wind dependence) : - not used if l_meteo_hom=T or l_meteo_stat=T : - not used if imeteo_exchtype=1 (uses bilinear interpolation) : : imeteo_exchtype : dependence on air-sea temperature difference for the surface : drag coefficient and the thermal exchange coefficient : = 0 -> no dependence on tdif : = 1 -> as a function of tdif (uses bilinear interpolation) : - not used if l_meteo_hom=T or l_meteo_stat=T : : cds : surface drag coefficient (usi) : - used if imeteo_dragtype=0 and imeteo_exchtype=0 : - used also if l_meteo_hom=T or l_meteo_stat=T : : paref : reference atmospherical pressure : : l_sflx_rpa : .true. to read the atmospherical pressure from a meteorological file : : name_sflx_pa : name of the variable relative to the atmospherical pressure : in the meteorological file : : l_sflx_rwind : .true. to read the wind at 10 m (or wind stress) from a meteorological file : : name_sflx_wind : name of the variable relative to wind in the meteorological file : name_sflx_wtau : name of the variable relative to wind stress in the meteorological file : : l_sflx_rsat : .true. to read the air temperature at 2m (or 10m) from a meteorological file : : name_sflx_sat : name of the variable relative to the air temperature : in the meteorological file : : l_sflx_rrh : .true. to read the relative humidity from a meteorological file : : name_sflx_rh : name of the variable relative to the relative humidity : in the meteorological file : : l_sflx_rsolar : .true. to read the solar flux from a meteorological file : : name_sflx_solar : name of the variable relative to the solar flux : in the meteorological file : : l_sflx_rir : .true. to read the infra-red flux from a meteorological file : : name_sflx_ir : name of the variable relative to the infra-red flux : in the meteorological file : : l_sflx_rflux : .true. to read the atmospherical fluxes in flux form : (turbulent heat and momentum fluxes are not calculated anymore through bulk formulae) : and to read non solar radiation flux (LW + lat + sens) from a meteorological file : : name_sflx_heat : name of the variable relative to the downward non solar radiation flux : (LW + lat + sens) in the meteorological file : : l_sflx_rrain : .true. to read the rain from a meteorological file : : name_sflx_rain : name of the variable relative to the rain in the meteorological file : name_sflx_emp : name of the variable relative to the evaporation minus precipication : in the meteorological file : : l_sflx_rfcc : .true. to read the fractional cloud cover from a meteorological file : : name_sflx_fcc : name of the variable relative to the fractional cloud cover : in the meteorological file : : l_sflx_rpa : .true. to read the atmospherical pressure from a meteorological file : : name_sflx_pa : name of the variable relative to the atmospherical pressure : in the meteorological file : : l_sflx_rath2m : specify if the air temperature or relative humidity is : available at 2m or 10 m. .True. if 2 m. : :----------------------------------------------------------------------------- &nammeteo imeteo_dragtype=0 imeteo_exchtype=0 cds = 0.0016 paref=101500.0 l_sflx_rpa=.true. l_sflx_rwind=.true. l_sflx_rsat=.true. l_sflx_rrh=.true. l_sflx_rsolar=.false. l_sflx_rir=.false. l_sflx_rrain=.false. l_sflx_rfcc=.true. l_sflx_rath2m=.true. l_sflx_rflux=.false. name_sflx_sat='t2m' ! 2-meter air temperature name_sflx_pa='pmer' ! atmospherical pressure name_sflx_rh='rh' ! relative humidity name_sflx_windx='u10m' ! 10-meter U Component name_sflx_windy='v10m' ! 10-meter V Component name_sflx_solar='rad' ! downward solar radiation flux (short wave heat flux) name_sflx_ir ='hflw' ! infra-red flux (long wave heat flux) name_sflx_rain ='rain' ! rain name_sflx_fcc ='fcc' ! cloud cover name_sflx_wtaux='sozotaux' ! zonal wind stress name_sflx_wtauy='sometauy' ! meridional wind stress name_sflx_heat ='flxsurf' ! downward non solar radiation flux (LW + lat + sens) name_sflx_emp ='sowaflup' / ! evaporation minus precipication : :---------------------------------------------------------------------------- : namsaverestart : gestion des Sauvegardes. : namsaverestart : saving parameters :----------------------------------------------------------------------------- : : l_saverestart_1file : logique positionne a .true. pour sauver a intervalles reguliers (longs : de saverestart_step) l etat du modele dans un fichier de sauvegarde, destine ensuite : a la reprise. Type fortran logical. : l_saverestart_1file : boolean set to .true. in order to save the model state at regular : timesteps (of saverestart_step) in a saving file to be used for restarts. : Type fortran logical. : : file_saverestart : nom du fichier dans lequel est sauvee la situation du modele : lorsque l on ne demande pas la sauvegarde sequentielle : (i.e. sauvseq=.false.). Type fortran character(len=lchain). : file_saverestart : file where the model state is saved when sequential saving is : not requested (i.e. sauvseq=.false.). : Type fortran character(len=lchain). : : l_saverestart_bydate : logique positionne a .true. pour effectuer une procedure de sauvegarde : sequentielle : procedure au cours de laquelle on cree differents fichiers : de sauvegarde dont le nom contient la date (jjmmaaaahhmm). A contrario : sauvseq=.false. inhibe cette sauvegarde sequentielle et les sauvegardes : sont effectuees toujours dans le meme fichier file_saverestart. : Type fortran logical. : l_saverestart_bydate : boolean set to .true. in order to plan sequential saving : : several saving files are written, with names which include : the saving time (ddmmaaaahhmm). : sauvseq=.false. inhibits sequential saving and saving is always : written in one and only file called file_saverestart. : Type fortran logical. : : saverestart_step : intervalle de temps entre deux sauvegardes. : Les sauvergardes sont effectuees dans le fichier au format netcdf : intitule file_saverestart. Cet intervalle de temps est compte : en jour. Type fortran real(kind=rsh) : saverestart_step : time step (in days) between consecutive savings. : Saving files are written in a netcdf filed called file_saverestart. : Type fortran real(kind=rsh) : :----------------------------------------------------------------------------- : &namsaverestart l_saverestart_1file = .false. l_saverestart_bydate = .true. file_saverestart= 'save.nc' saverestart_step = 60.0d0 / : :----------------------------------------------------------------------------- : namoutput : gestion des sorties de resultats : namoutput : result outputs :----------------------------------------------------------------------------- : : iscreen : 6 : affichage ecran; 12 : affichage dans le fichier listing : iscreen : 6 : dispay on the screen; 12 : display in the listing file : : iscreenlog : 6 : affichage ecran; 18 : affichage dans le fichier simu.log : iscreenlog : 6 : dispay on the screen; 18 : display in simu.log file : : ierrorlog : 6 : affichage ecran; 17 : affichage dans le fichier error.log : ierrorlog : 6 : dispay on the screen; 17 : display in error.log file : : iwarnlog : 6 : affichage ecran; 16 : affichage dans le fichier warning.log : iwarnlog : 6 : dispay on the screen; 16 : display in warning.log file : : out_dir : chemin du repertoire ou seront sauves les fichiers de sortie : out_dir : path of directory where the output files are written : : file_output : nom du fichier definissant les types de fichiers resultats. : file_output : file name in which output files are defined. : : l_out_nc4par : logique pour activation des entrees/sorties en parallele avec NetCDF-4 : l_out_nc4par : logical to turn the parallel NetCDF-4 writting on. : if .True. : each cpu writes its part in a unique global file : if .False. : each cpu writes its part in its specific file. : : l_out_pack : logique pour activation de la compression en short avec : scale_factor et add_offset : l_out_pack : logical to turn the packing of data values on and to reduce the size of the output file. : if .True., the values are saved in short type according to the following transformation : packed_data_value = nint((unpacked_data_value - add_offset) / scale_factor) :----------------------------------------------------------------------------- : &namoutput iscreen=12 iscreenlog=18 ierrorlog=17 iwarnlog=16 out_dir='/home1/scratch/mdiaz/MARS/RUN_MARS/GIRONDE_CURVI/results' file_output= './output.dat' l_out_nc4par = .true. l_out_pack = .false. / : :----------------------------------------------------------------------------- : namvisc :----------------------------------------------------------------------------- : : fvisc : coefficient dans la formulation de la viscosite : turbulente horizontale nu = fvisc(i)*0.01*(dx**1.15). La : plage de variation de fvisc est : 1 < fvisc < 17 : fvisc : coefficient in the horizontal turbulent viscosity formulation : nu = fvisc(i)*0.01*(dx**1.15). : interval : 1 < fvisc < 17 : : l_smagor : logique positionne a .true. pour evaluer le coefficient de: : viscosite laterale par une procedure de type Smagorinsky. : Positionne a .false. pour avoir une viscosite constante et homogene : (exceptee dans les couches eponges). Type fortran logical.: : l_smagor : boolean set to .true. to compute the lateral viscosity coefficient : according to a Smagorinsky-type procedure. : set to .false. for a constant and homogenous viscosity : (except in sponge layers). Type fortran logical. : : vismin : viscosite minimale m2.s-1. En general ce coefficient est choisi : superieur ou egal a 0.01. Type fortran real(kind=rsh). : vismin : minimum value of viscosity in m2.s-1. Generally greater than or : equal to 0.01. Type fortran real(kind=rsh). : : vismax : viscosite maximale m2.s-1. Ce coefficient depend de la taille : de maille. Un ordre de grandeur est donne par la formule :: : 0.015*dx**1.15 ou dx est la taille de maille du modele. : vismax : maximum value of viscosity in m2.s-1. This coefficient depends : on the grid cell size. The magnitude can be approached by the following : formulation :: : 0.015*dx**1.15 where dx is the model grid cell size. : : cosmag : coefficient de controle de Smagorinsky (sans unite). : la plage de variation est : 0.05 < cosmag < 0.3 : cosmag : Smagorinsky control coefficient (dimensionless). : values in the interval : 0.05 < cosmag < 0.3 : : sponge_nbcell : longueur en nombre de mailles de la couche eponge. : sponge_nbcell : length of the sponge layer (number of cells). : : vismul : coefficient multiplicateur de la viscosite a l interieur de la couche : eponge. Type fortran real(kind=rsh). : vismul : viscosity multiplicative coefficient within the sponge layer. : Type fortran real(kind=rsh). : : l_spongs : logique positionne a .true. pour mettre une couche eponge au sud. : Type fortran logical. : l_spongs : boolean set to .true. to add a sponge layer to the South. : Type fortran logical. : : l_spongn : logique positionne a .true. pour mettre une couche eponge au nord. : Type fortran logical. : l_spongn : boolean set to .true. to add a sponge layer to the North. : Type fortran logical. : : l_sponge : logique positionne a .true. pour mettre une couche eponge a l est. : Type fortran logical. : l_sponge : boolean set to .true. to add a sponge layer to the East. : Type fortran logical. : : l_spongw : logique positionne a .true. pour mettre une couche eponge a l ouest. : Type fortran logical. : l_spongw : boolean set to .true. to add a sponge layer to the West. : Type fortran logical. : :----------------------------------------------------------------------------- : &namvisc fvisc=5.0 l_smagor=.false. vismin=0.2 !0.1 vismax=70.0 !10.0 cosmag=0.27 sponge_nbcell=5 vismul=2.0 l_spongs=.false. l_spongn=.false. l_sponge=.false. l_spongw=.false. / : :----------------------------------------------------------------------------- : namturb :----------------------------------------------------------------------------- : : z0bot : Longueur de rugosite pour l estimation du coefficient de trainee : sur le fond (en m). Plage de variation autour de quelques millimetres. : Parametre du 3D. : z0bot : Roughness length to estimate the bottom drag coefficient (in m). : Of the order of a few mm. 3D model parameter. : : z0surf: Longueur de rugosite pour l estimation du coefficient de trainee : sous la surface (m). Plage de variation autour de quelques metres car : z0surf est egal a 1.6 * Hs (Hauteur significative des vagues) : Parametre du 3D. : z0surf: roughness length to estimate the surface drag coefficient (in m). : Of the order of a few meters since z0surf is equal to 1.6 * Hs (the significant wave height). : (Terray et al., 1996, 2000) : 3D model parameter. : : l_z0b_hom : Logique positionner a false pour lire un fichier netcdf contenant : une longueur de rugosité spatialise, sinon, la valeur de z0bot est utilisee : l_z0b_hom : Boolean set to false to read a netcdf file containing a spatialized roughness length : otherwise the z0bot value is used : : file_z0b : nom du fichier netcdf contenant la longueur de rugosité : file_z0b : filename of netcdf file containing roughness length : : hminfrot : hauteur minimale de la colonne d eau utilisee pour : l evaluation de la tension de frottement sur le fond et : l evaluation de la tension de surface. : hminfrot : minimum water column height used to compute bottom and surface : shear stresses : : hmaxfrot : hauteur maximal de la couche limite de fond (utilisee pour : l evaluation de la tension de frottement sur le fond) : hmaxfrot : maximimum thickness of the bottom boundary layer (used to compute bottom : shear stresses) : : turb_nbeq : nombre d equations du modele de turbulence (0, 1 ou 2) : Parametre du 3D. : turb_nbeq : number of equations in the turbulence model(0, 1 ou 2) : 3D model parameter. : : turb_0eq_option : numero de l option pour le modele a zero equation : 1- coefficients verticaux de viscosite et de diffusion constant : 2- modele de Prandt (longueur de melange + cisaillement) : 3- modele de Quetin (longueur de melange + Richardson) : 4- modele de Pacanovski et Philander (Richardson) : Parametre du 3D. : turb_0eq_option : option number for the zero equation model : 1- constant vertical viscosity and diffusion coefficients : 2- Prandt model (mixing length + shear) : 3- Quetin model (mixing length + Richardson) : 4- Pacanovski and Philander model (Richardson) : 3D model parameter. : : turb_2eq_option : option number for the two-equation model : 3D-model parameter. : 1- k-kl : 2- k-epsilon : 3- k-omega : 4- modele generique : turb_2eq_option : numero de l option pour le modele a deux equations : Parametre du 3D. : 1- k-kl : 2- k-epsilon : 3- k-omega : 4- modele generique : : bgdiff : viscosite et diffusion verticale de fond (background diffusivity) : (en m2/s). Cette valeur est telle que nz >= bgdiff : et kz >= bgdiff : Parametre du 3D. : bgdiff : background vertical viscosity and diffusivity : (in m2/s), with values such that nz >= bgdiff : and kz >= bgdiff : 3D-model parameter. : : kzinit : coefficient de diffusion verticale a l initialisation. : Cette valeur est aussi utilisee dans le cas de modele de : turbulence constant (en m2/s) (turb_0eq_option == 1). : Parametre du 3D seulement. : kzinit : initial vertical diffusivity coefficient. : Coefficient also used thoughout the computation in case of a constant : turbulence model (in m2/s) (turb_0eq_option == 1). For 3D computations only. : : nzinit : coefficient de viscosite verticale a l initialisation. : Cette valeur est aussi utilisee dans le cas de modele de turbulence : constant (en m2/s) (turb_0eq_option == 1). Parametre du 3D seulement. : nzinit : initial vertical viscosity coefficient. : Coefficient also used throughout the computation in case of a constant : turbulence model (in m2/s) (turb_0eq_option == 1). For 3D computations only. : : l_stability : logique positionne a .true. pour avoir une colonne : d eau toujours stable (on melange la partie de la colonne : d eau instable jusqu a avoir la stabilite). : Parametre du 3D. : l_stability : boolean set to .true. in order for the water column to : always be stable (the unstable part of the column is mixed until : stability is reached) : 3D-model parameter. : :----------------------------------------------------------------------------- : &namturb z0bot=0.0001 z0surf=0.5000 l_z0b_hom=.false. file_z0b='../inputs/z0_spat_16nov18_V2.nc' !'/home/datawork-dhysed-s/PROJETS/AMORAD/FORCAGES/NEW_BATHY/z0_init_spatialise_new-red.nc' hminfrot=0.2 hmaxfrot=5.0 turb_nbeq =2 turb_0eq_option =2 turb_2eq_option=2 bgdiff=1.0e-5 nzinit=1.0e-4 kzinit=1.0e-4 l_stability=.false. / : :----------------------------------------------------------------------------- : namfrot :----------------------------------------------------------------------------- : : botstressmax : vitesse critique exprimee en m.s-1 au-dela de laquelle : on augmente, aux points ou cette limite est depassee, : artificiellement le coefficient de frottement. : Parametre du 2D. : botstressmax : critical velocity (in m.s-1) beyond which the shear stress coefficient : is artificially increased in the cells where the critical velocity is exceeded. : 2D model parameter. :----------------------------------------------------------------------------- &namfrot l_incbotstress=.true. ! not used yet botstressmax=3.0 / : :----------------------------------------------------------------------------- : namdiff :----------------------------------------------------------------------------- : : kx : coefficient de diffusion des substances dissoutes, du sel et de la : temperature en m2/s dans la direction W-E. : kx : diffusion coefficient for dissolved matter, salt and : temperature (in m2/s) in the W-E direction. : : ky : coefficient de diffusion des matieres dissoutes, du sel et de la : temperature en m2/s dans la direction S-N. : ky : diffusion coefficient for dissolved matter, salt and : temperature (in m2/s) in the S-N direction. :----------------------------------------------------------------------------- &namdiff kx = 5.0 ky = 5.0 / : :----------------------------------------------------------------------------- : namgridsig : repartition des sigmas ; specifique 3D : namgridsig : sigma levels distribution ; specific to the 3D model :----------------------------------------------------------------------------- : : l_equisig : equisig=.true. donne une equirepartition des niveaux sigma. : equisig=.false. desactive l equirepartition des sigmas et oblige : l utilisateur a donner la repartition des niveaux ci-dessous. : attention ! le nombre de sigmas que l on donne ci-dessous : doit correspondre exactement a kmax du fichier INC/parameters.F90 : Parametre du 3D. Type fortran logical : l_equisig : boolean set to .true. for a regular sigma level distribution. : boolean set to .false. to desactivate the regular distribution of : sigma levels, in which case the user has to define the new : distribution. : look out ! the number of sigma layers prescribed has to be : exactly equal to kmax as defined in INC/parameters.F90 : 3D model parameter. Type fortran logical : : sig : niveaux sigma. Doit etre renseigne si equisig=.false. : Type fortran real(kind=rsh),dimension(kmax) : sig : sigma level vector. Look out! To be defined if equisig=.false. : Type fortran real(kind=rsh),dimension(kmax) : : sig(1)= -0.95,sig(2)= -0.8,sig(3)= -0.65,sig(4)= -0.5, : sig(5)= -0.35,sig(6)= -0.25,sig(7)= -0.15,sig(8)= -0.1, : sig(9)= -0.05,sig(10)= -0.03 : : hc : a la fois profondeur minimale et profondeur intermediaire au- : dessus de laquelle les niveaux sigma generalises sont resserres : hc : either the minimum depth or a shallower depth above which we wish : to have more vertical resolution (generalized sigma) : : b_sig : bottom control parameter for generalized sigmas (bounded by 0 and 1) : : theta_sig : surface control parameter for generalized sigmas (bounded by 0 and 20) : let theta_sig go to zero to get pure sigma coordinates : : grid_angle : angle de rotation de la grille (positif dans le sens des : aiguilles d une montre, le nord est la reference soit 0) : grid_angle : rotation angle of the grid (positive and clockwise from the north) : : l_corio_cst : si .true., la parametre de Coriolis est constant : l_corio_cst : boolean set to .true. for constant Coriolis parameter : :----------------------------------------------------------------------------- : &namgridsig l_equisig=.true. sig(1)= -0.05 hc = 5.0 grid_angle = 0.0 theta_sig=0.000001 b_sig=1 l_corio_cst=.false. / : :----------------------------------------------------------------------------- : namriver : parametres concernant les fleuves. : namriver : river related parameters. :------------------------------------------------------------------------------ : : file_river : fichier decrivant les apports d eau douce (positions, debit...) : Type fortran character(len=lchain) : file_river : file describing fresh water outflows (position, discharge...) : Type fortran character(len=lchain) : : icon : prise en compte ou non des effets de flottabilite : - 0 : pas de calcul de salinite ni de temperature et flottabilite : - 1 : calcul de temperature, de salinite et de flottabilite : (pour le seul calcul de la turbulence). Pas de prise : en compte du gradient de pression interne : - 2 : calcul de temperature, de salinite et de flottabilite : et prise en compte du gradient de pression interne : Type fortran integer. : icon : defines whether buoyancy effects are taken into account : - 0 : no salinity, temperature and buoyancy computation : - 1 : salinity, temperature and buoyancy (for turbulence purpose) : are computed but the internal pressure gradient is not accounted for. : - 2 : salinity, temperature and buoyancy are computed and the internal : pressure gradient is accounted for. : Type fortran integer. : :----------------------------------------------------------------------------- : &namriver file_river = '../inputs/rivers/fleuves_gironde_noCutRiver.dat' icon = 2 / : :----------------------------------------------------------------------------- : namparanum :----------------------------------------------------------------------------- : : dtini : pas de temps initial (dans le cas de regulation du pas : de temps) sinon pas de temps du modele (pas de temps fixe) : en seconde. Il est raisonnable que ce pas de temps soit tel que : : u*dtini/dx ~ 0.2 avec u le courant max rencontre sur la : zone et dx la taille de maille. : Type fortran real(kind=rlg) : dtini : initial time step (in case of an adaptative time step), otherwise : fixed model time step (in s). It is reasonable to choose it : so that: : u*dtini/dx ~ 0.2 where u is the maximum current encountered in the : domain and dx is the grid cell size. : Type fortran real(kind=rlg) : : dtmin : pas de temps minimal en dessous duquel on ne descendra pas : dans la procedure de regulation du pas de temps en seconde. : Type fortran real(kind=rlg) : dtmin : minimum value of the adaptative time step (in s). : Type fortran real(kind=rlg) : : dtmax : pas de temps maximal autorise dans la procedure de : regulation du pas de temps en seconde. Ce pas de temps : doit etre tel que : u*dtmax/dx < 1 avec u le courant max : rencontre sur la zone et dx la taille de maille : : dtmax : maximum value of the adaptative time step (in s). : This time step must be such that : u*dtmax/dx < 1 where u is the maximum current encountered in the : domain and dx the grid cell size. : Type fortran real(kind=rlg) : : l_modele2d : si .true., resolution des equations barotropes quel que : soit le nombre de couches defini par kmax : : l_modele2d : if .true., MARS solves barotropic equations whatever the : number of layers specified by kmax : : cflcrt : nombre de courant maximum autorise dans la procedure de : pas de temps adaptatif. En theorie : ce nombre vaut 1. En pratique le modele est destabilise : des que cflcrt > 0.7. Si ce seuil est franchit le pas de temps : est reduit. : cflcrt : maximum value of the current number used to compute the adaptative : time step. Theoretically equal to 1. Practically, the model is : unstable as soon as cflcrt > 0.7. : The time step is reduced if this threshold is exceeded. : :----------------------------------------------------------------------------- : &namparanum dtini=5.d0 dtmin=5.d0 dtmax=16.d0 l_modele2d=.false. cflcrt=0.6 / : :----------------------------------------------------------------------------- : namthermo :----------------------------------------------------------------------------- : : l_sflx_radlossbot : logique positionne a .true. pour que le fond : ne soit pas une barriere thermique ; ainsi la partie : radiative non encore absorbee dans la colonne d eau est perdue. : Type fortran logical. : l_sflx_radlossbot : boolean set to .true. in order for the bottom to not be : a thermic barrier ; in that case, the radiative part not : absorbed by the water column is lost. : Type fortran logical. : : l_sflx_solarcst : representation schematique des flux de chaleur air-mer : fixes a une constante egale a sflx_solarcst. : Type fortran logical. : l_sflx_solarcst : boolean set to .true. for constant : air-sea heat fluxes (equal to sflx_solarcst). : Type fortran logical. : : sflx_solarcst : flux de chaleur echange a l interface air-mer en : cas d echanges schematiques constants (en W/m2). : Type fortran real(kind=rsh). : sflx_solarcst : heat flux value at the air-sea interface in case of constant : schematic exchanges (in W/m2). : Type fortran real(kind=rsh). : :----------------------------------------------------------------------------- : &namthermo l_sflx_radlossbot=.false. l_sflx_solarcst=.true. sflx_solarcst=0.0 / : :----------------------------------------------------------------------------- : namtraj :----------------------------------------------------------------------------- : : file_trajec : nom du fichier decrivant la position des taches et leurs : caracteristiques : file_trajec : file where the "spots" position and characteristics : are described. : itypetraj : type of definition for patch initial position: : 1: circle, 2: rectangle, 3:netcdf : ndtz : number of time-step within the hydrodynamic model time step : See (Ross and Sharples, 2004, Limnol.Oceanogr.) for sensitivity analysis : To choose depending on the model time resolution, dt/ndtz should be of the order of few seconds : :----------------------------------------------------------------------------- : &namtraj file_trajec='../inputs/TRAJ/traject3d_type3ncdf.dat' itypetraj = 3 ndtz = 50 / : :----------------------------------------------------------------------------- : namdiag : controle des diagnostics. : namdiag : diagnostic control. :----------------------------------------------------------------------------- : : l_diag : logique a positionner a true pour la realisation des diagnostics : implementes dans la routine diagnostic.F90. Type fortran logical. : l_diag : boolean set to .true. in order to compute and write diagnostics implemented in : diagnostic.F90. Type fortran logical. : : date_startdiag : date du debut des diagnostics (chaine de 19 caracteres : JJ/MM/AAAA hh:mm:ss) Type fortran character(len=19). : date_startdiag : diagnostics starting time (19 character string : DD/MM/AAAA hh:mm:ss) Type fortran character(len=19). : : date_enddiag : date de fin des diagnostics (chaine de 19 caracteres : JJ/MM/AAAA hh:mm:ss) Type fortran character(len=19). : date_enddiag : diagnostics finishing time (19 character string : DD/MM/AAAA hh:mm:ss) Type fortran character(len=19). : : file_diag : nom du fichier texte destine a recevoir les series chrono : Type fortran character(len=lchain). : file_diag : ascci file where time series will be written : Type fortran character(len=lchain). : : l_points : logique a positionner a true pour le suivi de points particuliers : l_points : boolean set to .true. for the output of time series at locations : defined in file_point : : file_point : nom du fichier texte contenant la description : d un ensemble de points en lesquels on va sortir les : courants, niveau,salinite, temperature traceur a tous les : pas de temps. : Type fortran character(len=lchain). : file_point : name of the ascii file describing a number of points where : currents, water levels, salinity, : temperature and tracers are written at every time step. : Type fortran character(len=lchain). : :----------------------------------------------------------------------------- : &namdiag l_diag=.false. date_startdiag = '01/01/1996 00:00:00' date_enddiag = '02/01/2001 18:00:01' file_diag='diag.gdg51' l_points=.false. file_point='../inputs/points.dat' / : :------------------------------------------------------------------------------ : path_offline : directory where offline output files (containing conservative : averaged velocities, water depths and sea level) (-Dkey_offline) : and the list of files list_offline.dat are located : : l_bz_offline : to estimate the buoyancy from salinity and temperature fields : (-Dkey_offline) interpolated on time : :------------------------------------------------------------------------------ &namoffline path_offline='./offline_files/' l_bz_offline=.false. / : :----------------------------------------------------------------------------- : namsurge : parametres de conditions limites surcotes + maree : namsurge : parameters to manage obc with surge and tidal signal :----------------------------------------------------------------------------- : : l_tide_surge : logique activant des conditions limites en maree et en surcote : l_tide_surge : logique turning on obc with tidal signal + surge : : file_tide_harmcp_surge : fichier contenant les constantes harmoniques de maree : les plus precises disponibles sur le plateau en France. : Fichier utilise si on active la surcote. : file_tide_harmcp : file of tidal harmonic components, more precise over : Biscay shelf (143 components) : Used when surge signal is turning on. : : file_surge_meteo : fichier de conditions limites cree par le rang precedent avec meteo : file_surge_meteo : obc file issued from a previous rank with meteorology (wind and atmospherical pressure) : : file_surge_nometeo : fichier de conditions limites cree par le rang precedent sans meteo : file_surge_nometeo : obc file issued from a previous rank without any meteorology : :----------------------------------------------------------------------------- : &namsurge l_tide_surge=.false. file_tide_harmcp_surge='../inputs/CSTFRANCE_SHOM2004.nc' file_surge_meteo='../inputs/obc2d_pern.rank2_pern.V8.00_wind' file_surge_nometeo='../inputs/obc2d_pern.rank2_pern.V8.00_nowind' / : :----------------------------------------- : namlagoon : specific with ley_lagoon (used with key_lagoon) : hnivmoy_lagoon : niveau moyen par rapport au zero hydro : nbgrau : nombre de graus (max nbgraumax in parameters) : igrau(1)=172 , igrau(2)=18 , igrau(3)=1 : jgrau(1)=1 , jgrau(2)=12 , jgrau(3)=1 : type_grau(1)=3 ,type_grau(2)=3 ,type_grau(3)=3 : igrau,jgrau : position des graus (autant de valeurs que nbgrau) : type_grau : frontiere Nord (1), Sud (3), Est (2), Ouest(4) (autant de valeurs que nbgrau) : l_open_grau : true si grau ouvert (autant de valeurs que nbgrau) :----------------------------------------- : &namlagoon hnivmoy_lagoon=0.45 nbgrau=1 igrau(1)=172 jgrau(1)=1 type_grau(1)=3 file_obc_ogcm_ts='/home1/datawork/afiandri/MARS_AGRIF/RUN_MARS/THAU/inputs/CLIM/ST_sete_2003_2015_survalQ2.prn' /