{{{
#!html
INCA TECHNICAL DOCUMENTATION
}}}
----
----
[[PageOutline]]
= INCA5 =
INCA5 is the version developped for CMIP6.
== Several Reference versions ==
* '''INCA5.1.0''' --> first version of INCA5 = INCA4 + Nitrates
* '''INCA5.1.1''' --> now outputs are managed by XIOS library
= INCA 4 =
INCA4 is the 39 level version of the INCA model. It's fit with version 5 of the model LMDZ (LMDZ5)
== Several Reference versions ==
* '''INCA4.1.0''' --> first version of INCA4 (running on titane (CCRT))
* '''INCA4.1.1''' --> We fixed some bugs, in particular in optaer_5wv and optaer_2bd (there was an error on the number of vertical levels). This version can run on Curie (TGCC) and is compatible with the version rev2019 og LMDZORINCA_v5.
* '''INCA4.1.2''' --> same version with a new parameter file (inca.def) to manage some parameters in offline.
* '''INCA4.1.3''' --> version with modification to fit with LMDZ radiation in inca radlwsw_inca (version use with IPSLCM5CHT_v5) -- update 3/12/2013
= Pre-processor =
With the INCA model you can choose between several chemistry configuration (AER, NMHC_AER or GES). For each one we give a pre-process file that will be interpret by a pre-processor. These files are store in INCA/INP/ directory.
Its documentation is available [http://www-lsceinca.cea.fr/pdf/Doc_inca_1.pdf here].
= Technicals Notes on INCA =
== How change the molar mass in pre-processing ==
For example if you want give the carbon molar mass to the CO2, you need to modify the INP/*def file of your chemistry configuration like this
{{{
CO2OC -> CO2
become
CO2OC -> CO2@12
}}}
= Input Files =
== NMHC_AER ==
|| File Name || Resolution || dimensions ||
||o3clim.nc || [lonxlat-lon+2]x47x12 || vector pressure time ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||so4.nc || [lonxlat-lon+2]x19x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
||rhv.dat || [lonxlat-lon+2] lignes || ||
||wth.dat || [lonxlat-lon+2] lignes || ||
||cly.dat || [lonxlat-lon+2] lignes || ||
||rhvEC.txt || 51520 lignes || ||
||wthEC.txt || 51520 lignes || ||
||clyEC.txt || 51520 lignes || ||
== AER ==
|| File Name || Resolution || dimensions ||
||oxydants.nc || lonx[lat+1]x19x12 || x y presnivs time_counter ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
||rhv.dat || [lonxlat-lon+2] lignes || ||
||wth.dat || [lonxlat-lon+2] lignes || ||
||cly.dat || [lonxlat-lon+2] lignes || ||
||rhvEC.txt || 51520 lignes || ||
||wthEC.txt || 51520 lignes || ||
||clyEC.txt || 51520 lignes || ||
== GES ==
|| File Name || Resolution || dimensions ||
||oxydants.nc || lonx[lat+1]x19x12 || x y presnivs time_counter ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
||andres.nc || [lon+1]x[lat+1]x12 || LON LAT TIME ||
||casa_m.nc || [lon+1]x[lat+1]x12 || LON LAT TIME ||
||casa_h.nc || [lon+1]x[lat+1]x8*nb_day_in_mth || LON LAT TIME ||
||taka.nc || [lon+1]x[lat+1]x12 || LON LAT TIME ||
== NMHC (old configuration) ==
|| File Name || Resolution || dimensions ||
||o3clim.nc || [lonxlat-lon+2]x47x12 || vector pressure time ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||so4.nc || [lonxlat-lon+2]x19x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
== CH4_AER (old configuration) ==
|| File Name || Resolution || dimensions ||
||o3clim.nc || [lonxlat-lon+2]x47x12 || vector pressure time ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||so4.nc || [lonxlat-lon+2]x19x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
||rhv.dat || [lonxlat-lon+2] lignes || ||
||wth.dat || [lonxlat-lon+2] lignes || ||
||cly.dat || [lonxlat-lon+2] lignes || ||
||rhvEC.txt || 51520 lignes || ||
||wthEC.txt || 51520 lignes || ||
||clyEC.txt || 51520 lignes || ||
== CH4 (old configuration) ==
|| File Name || Resolution || dimensions ||
||o3clim.nc || [lonxlat-lon+2]x47x12 || vector pressure time ||
||npp.nc || [lonxlat-lon+2]x12 || vector time ||
||aircraft.nc || [lonxlat-lon+2]x26x12 || vector lev time ||
||so4.nc || [lonxlat-lon+2]x19x12 || vector lev time ||
||sflx.nc || [lonxlat-lon+2]x12 || vector time ||
||phototable.dat || || ||
||landuse.nc || [lonxlat-lon+2]x11 || vectore type ||
== Description du fichier sflx en fonction des configurations chimiques ==
=== les espèces ===
Les flux de surfaces des espèces sont contenus dans le fichier sflx.nc lu par la routine sflx_inti.F90. [[BR]]
Suivant la configuration que l'on utilise ce fichier ne doit pas contenir les même données. Voici une liste
rappelant l'ensemble de ces espèces. [[BR]]
|| '''CH4''' || '''NMHC''' ||'''GES'''||'''AER''' || '''AERONLY''' ||
||flx_n2o ||flx_co2 ||flx_n2o ||flx_co2tot ||flx_bc ||
||flx_ch4 ||flx_n2o ||flx_ch4 ||flx_co2oce ||flx_pom ||
||flx_sf6 ||flx_mcf ||flx_mcf ||flx_co2bio ||flx_bbbc ||
||flx_co2 ||flx_nox ||flx_co ||flx_bc ||flx_bbpom ||
||flx_kr85 ||flx_h2 || ||flx_pom ||flx_so2 ||
||flx_f11 ||flx_co || ||flx_bbbc ||flx_so4 ||
||flx_f12 ||flx_cobbg || ||flx_bbpom ||flx_h2s ||
||flx_f113 ||flx_ch4 || ||flx_so2 ||conc_dms ||
||flx_mcf ||flx_ch3oh || ||flx_so4 ||fractnat_aer ||
||flx_ccl4 ||flx_c2h5oh || ||flx_h2s ||fractnat_bc ||
||flx_ch3cl ||flx_c2h6 || ||conc_dms ||fractnat_pom ||
||flx_ch3br ||flx_c3h8 || ||fractnat_aer ||fractnat_so4 ||
||flx_f22 ||flx_akan || ||fractnat_bc ||fractnat_dust ||
||flx_f134 ||flx_c2h4 || ||fractnat_pom ||fractnat_ss ||
||flx_f141 ||flx_c3h6 || ||fractnat_so4 || ||
||flx_f142 ||flx_c2h2 || ||fractnat_dust || ||
||flx_h1301 ||flx_alken || ||fractnat_ss || ||
||flx_h1211 ||flx_arom || || || ||
||flx_co ||flx_ch2o || || || ||
||flx_co_sec ||flx_ch3cho || || || ||
||flx_nox ||flx_ch3coch3 || || || ||
||flx_h2 ||flx_mek || || || ||
|| ||flx_mvk || || || ||
|| ||flx_ch3cooh || || || ||
|| ||flx_isoprene || || || ||
|| ||flx_terpenes || || || ||
Si vous travaillez avec une configuration utilisant les aérosols (CH4_AER ou NMHC_AER) en plus de la chimie vous devez
ajouter la colonne AER à la colonne de votre configuration d'origine. Ex : CH4 + AER (et non pas AERONLY) [[BR]]
==== Nouveaux fichiers flx (septembre 2018) ====
||= flx_NAT =||= flx_ANT =||= flx_MISC =||= flx_GHG =||= flx_BBG =||
|| flx_n2o_nat || flx_ch4_ant || conc_dms || flx_sf6_ghg || flx_ch4_bbg ||
|| flx_ch4_nat || flx_ch3oh_ant || conc_nh4 || flx_co2_ghg || flx_ch3oh_bbg ||
|| flx_ch3oh_nat || flx_c2h5oh_ant || fractnat_aer || flx_n2o_ghg || flx_c2h5oh_bbg ||
|| flx_c2h5oh_nat || flx_c2h6_ant || fractnat_bc || flx_kr85_ghg || flx_c2h6_bbg ||
|| flx_c2h6_nat || flx_c3h8_ant || fractnat_pom || flx_f11_ghg || flx_c3h8_bbg ||
|| flx_c3h8_nat || flx_alkan_ant || fractnat_so4 || flx_f12_ghg || flx_alkan_bbg ||
|| flx_alkan_nat || flx_c2h4_ant || fractnat_dust || flx_f113_ghg || flx_c2h4_bbg ||
|| flx_c2h4_nat || flx_c3h6_ant || fractnat_ss || flx_h1301_ghg || flx_c3h6_bbg ||
|| flx_c3h6_nat || flx_c2h2_ant || || flx_h1211_ghg || flx_c2h2_bbg ||
|| flx_c2h2_nat || flx_alken_ant || || flx_ccl4_ghg || flx_alken_bbg ||
|| flx_alken_nat || flx_arom_ant || || flx_c2cl4_ghg || flx_arom_bbg ||
|| flx_arom_nat || flx_ch2o_ant || || flx_ch3cl_ghg || flx_ch2o_bbg ||
|| flx_ch2o_nat || flx_ch3cho_ant || || flx_ch3br_ghg || flx_ch3cho_bbg ||
|| flx_ch3cho_nat || flx_ch3coch3_ant || || flx_mcf_ghg || flx_ch3coch3_bbg ||
|| flx_ch3coch3_nat || flx_mek_ant || || flx_f22_ghg || flx_mek_bbg ||
|| flx_mek_nat || flx_mvk_ant || || flx_f134_ghg || flx_mvk_bbg ||
|| flx_mvk_nat || flx_ch3cooh_ant || || flx_f141_ghg || flx_ch3cooh_bbg ||
|| flx_ch3cooh_nat || flx_co_ant || || flx_f142_ghg || flx_co_bbg ||
|| flx_isoprene_nat || flx_h2_ant || || || flx_h2_bbg ||
|| flx_terpenes_nat || flx_so2_ant || || || flx_so2_bbg ||
|| flx_co_nat || flx_nh3_ant || || || flx_h2s_bbg ||
|| flx_h2_nat || flx_hono_ant || || || flx_ocs_bbg ||
|| flx_so2_nat || flx_nox_ant || || || flx_nh3_bbg ||
|| flx_h2s_nat || flx_bc_ant || || || flx_hono_bbg ||
|| flx_nh3_nat || flx_pom_ant || || || flx_nox_bbg ||
|| flx_hono_nat || || || || flx_bc_bbg ||
|| flx_ocs_nat || || || || flx_pom_bbg ||
|| flx_nox_nat || || || || ||
[[BR]]
[[BR]]
=== Entête des fichiers suivant la résolution du modèle ===
-d 96x72x19
{{{
ncdump -h sflx.nc >>
dimensions :
vector = 6818;
time = UNLIMITED; //(12 currently)
}}}
-d 96x71x19
{{{
ncdump -h sflx.nc >>
dimensions :
vector = 6722;
time = UNLIMITED; //(12 currently)
}}}
= Les Output =
== Fichiers d'ouput netcdf ==
Le code INCA crée différents fichiers d'output :
* inca_inst = fichier avec une valeur journalière instantanée pour chaque variable (fichier désactivé par défaut)
* inca_avgr = fichier avec une valeur journalière moyennée pour chaque variable (fichier sauvegardé par défaut)
* forcage = différents diagnostiques propres aux aérosols
== Comment choisir les variables des fichiers de sorties ==
Le modèle INCA permet de sortir deux types de fichiers netcdf : [[BR]]
* inca_avgr.nc avec des moyennes journalières
* inca_inst.nc avec des données instantanées
Les variables contenues dans ces fichiers sont définies dans le fichier def en entrée du pré-processeur (INP/inca_AER.def, INP/inca_NMHC_AER.def ...)
Il contient une rubrique Outputs du type :
{{{
Outputs
Write frequency = 1d
Start File Number = 001
Density = 4
Retention time = 1y
Species
Inst
All
Endlst
Avgr
All
Endlst
Endent
(...)
Endent
}}}
Cette rubrique contient plusieurs catégories chacune représentant un type de variable :
Liste des catégories :
Species - Group members - Surface flux - Deposition Velocity - Photorates - Reaction rates - Washout rates - External forcing - temperature - water vapor - surface pressure - production - loss.
||Nom du champs dans inca || Signification ||Dimension ||Unité ||catégorie ||Nom dans le fichier de sortie ||
||mmr ||Rapport de mélange des traceurs ||3D ||VMR ||species || ||
||nas ||Rapport de mélange des membres d'un groupe ||3D ||MMR ||group members || ||
||eflux ||Flux de surface (émission) ||2D ||Kg/m2/s ||surface flux ||Emi_ ||
||dvel ||Vitesse de dépôt ||2D ||Cm/s ||deposition velocity||Dep_ ||
||temp ||temperature ||3D ||K ||temperature || ||
||h2o ||Rapport de mélange de la vapeur d'eau ||3D ||MMR ||water vapor || ||
||ps ||Pression de surface ||2D ||Pa ||surface pressure || ||
||jrates ||Taux de photolyse ||3D ||cm-3 s-1 ||photorates ||j001... ||
||reaction_rates ||Taux de réaction ||3D ||cm-3 s-1 ||reaction rates ||k001... ||
||hrates ||Taux de lessivage ||3D ||s-1 ||washout rates ||w001... ||
||extfrc ||Forçage extérieur ||3D ||cm-3 s-1 ||external forcing ||ext001... ||
||pox ||Taux de production ||3D ||cm-3 s-1 ||production || ||
||lox ||Taux de destruction ||3D ||cm-3 s-1 ||loss || ||
||dflux ||Flux de surface (dépôt) ||2D ||kg/m2/s ||deposition flux || ||
Pour chacune des catégories vous pouvez choisir soit All, soit indiquer la liste des variables qui vous interessent
{{{
Surface Flux
Inst
All
Endlst
Avgr
Pb210
CO,NO,NO2,SO2,DMS
CIDUSTM, CIN
CSSSM, CSSO4M, CSN
AIBCM, AIPOMM, AIN
ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN
SSSSM, SSN
Endlst
Endent
}}}
Vous pouvez également indiquer que vous souhaitez garder cette catégorie vide
{{{
Surface Flux
Inst
Endlst
Avgr
All
Endlst
Endent
}}}
== Choisir un fichier de sortie netcdf ==
Par défaut nous ne proposons pas de conserver le fichier inca_inst.nc Si vous voulez l'ajouter il faut l'indiquer
dans la carte inca.card (modipsl/config/LMDZINCA_v2/EXP_.../COMP/inca.card) [[BR]]
{{{
[OutputFiles]
List= (inca_avgr.nc, ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc, NONE) \
(inca_inst.nc, ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_inst.nc, NONE)
}}}
De même si vous voulez ne pas conserver le fichier inca_avgr il faut l'enlever de cette liste.
== Fichiers d'ouput texte ==
Les sorties texte de INCA (print *) sont stockées dans le fichier inca.out que vous retrouverez dans le répertoire IGCM_OUT/.../!JobName/CHM/Debug/ ou IGCM_OUT/../!JobName/DEBUG/