Atmospherical forcing

Scientific description

Formulations Scientific part. see rapport under MARS website

CPP keys

tab

Parameterization

NetCDF input file

1. Specify the name of the file you want to read in paracom.txt (file_meteo)

2. In paraspec.txt choose the variables you want to read and their names in the netCDF file:

   l_sflx_rpa=.false.
   l_sflx_rwind=.false.
   l_sflx_rsat=.false.
   l_sflx_rrh=.false.
   l_sflx_rsolar=.false.
   l_sflx_rir=.false.
   l_sflx_rrain=.false.
   l_sflx_rfcc=.false.
   l_sflx_rath2m=.true.
   l_sflx_flux=.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='swhf'    ! solar radiation flux (short wave heat flux)
   name_sflx_ir   ='lwhf'    ! 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 precipitation
   

To read the fluxes in a flux form (and to not use the use of bulk formulae to estimate the turbulent fluxes) set l_sflx_rflux to .true. as well as l_sflx_rwind, l_sflx_rsolar, l_sflx_rrain.

Frequent encountered problems

1. Dans sflxread.F90 au premier pas de temps, on lit chaque variable et regarde les min max pour voir si c’est dans les valeurs sont réalistes (en particulier pour arrêter le code s’il y a un pb d’unité, ou si la couverture nuageuse par exemple d ‘est pas comprise en 0 et 1). Parfois certains codes atmosphériques donnent de faibles valeurs négatives. Dans ce cas, tl faut changer les bornes des tests dans sflxread.F90

2. Vérifier les signes, unités des variables atmosphériques contenues dans le fichier ainsi que leur contenu: champs cumulés / instantannés.

   Le code MARS attend:

   vent (m/s)
   flux solaire ... (W.m-2) positive quand cela chauffe l'ocean donc downward
   Evaporation Upward (positive when the ocean evaporates), (kg/m2 = mm/s)
   pluie Downward (positive when it is raining), (kg/m2 = mm/s)
   pression (Pa)
   

   Parfois il faut changer le signe du champ pluie voire du flux infra-rouge dans sflxread.F90. Il faut aussi modifier sflxread.F90 selon si les flux sont cumulés ou instantannés dans le fichier.

3. La date d’origine des temps.

   L’origine des temps de MARS est 01-01-1900 00:00:00. Si l’origine des temps est différente dasn le fichier météo, il faut recaler la date avec la ligne suivante dans sflxread.F90:

   torigin=tool_datosec('01/01/1900 00:00:00')-tool_datosec('01/01/0000 00:00:00')
   

Note

Les points 2 et 3 ne peuvent se produire que si on part des champs météorologiques diffusés par météofrance. Quand on utilise les champs remaniés par notre outil chainop, aucun souci sur les signes, unités et champs cumulés, ni sur la date d’origine des temps.