Le changement de bathymetrie Orca nous oblige a recalculer les fichiers de poids (Sophie)
 
 

Le calcul des restarts Oasis:

1. Restart oceaniques:

Routine: /home/dasp/maisonna/Aux_oasis/Forcage/Predicate/rw_Rlevitus_Roasis.f90
On part des climatologies Levitus (valeurs interpolees partout sur ls continents) que l'on lit-ecrit au format netcdf. Valeurs du mois de janvier pour les champs de temperature. La glace est calculee comme suit: 1 si la temperature est superieure a -0.77, 0 sinon. L'albedo est de 0.065 sur les points mer et de 0.67 sur les points glace. On obtient le restart oce2oas.nc.
 

2. Restart atmospheriques:

1ere partie: 
Dans /home/dasp/maisonna/Aux_oasis/Forcage/Oasis, on interpole les champs de temperature et d'albedo du restart oceanique Oasis precedemment calcule (un champ a la fois, puis concatenation des deux resultats dans le fichier flxocean).
2eme partie:
Routine: /home/dasp/maisonna/Aux_oasis/Forcage/Predicate/rw_Rarpege_Roasis.f90
On lit les climatologies arpeges construites par t2 (lecture Vairmer) , et on y ajoute les champs de temperature et albedo que l'on vient de creer(lecture binaire). On obtient le restart atm2oas.nc.
 

3 Problemes resolus:

• L'examen des champs de vents avant et apres interpolation fait apparaitre un defaut, notamment au large des cotes du Chili . Le probleme venait du fait que notre grille at42 (et son masque) avait ete definie dans le sens inverse (Nord-Sud)  de la convention Oasis (Sud-Nord).
• L'operation SUBGRID ne se fait pas (debordement de tableau). En fait, nos deux SUBGRID utilisent le deuxieme fichier de poids MOZAIC alors que la variable jpnfs de parameter.h est a 1. Du coup, SUBGRID lit le fichier mozaic et mets le resultat dans le deuxieme element du tableau de subgrid ... qui a ete dimensionne a 1. Solution: on change les numero des operations mozaic: celui qui est utilise par SUBGRID porte le numero 1 .
• Le fichier de poids du runoff avait ete constitue sans tenir compte des points de recouvrement de la grille ORCA (i=181 et i=182). Il y avait donc deux points de deversement en trop dans l'Antarctique (meme probleme pour le fichier de poids du calving).
• Le runoff de la mer Caspienne se deversait sur un point de grille de tres faible surface, risquant de generer ponctuellement des baisses importantes de salinite. 5 points sont desormais en embouchure de la Volga.
• Probleme legerement different poutr le Danube: on augmente egalement le nombre de points embouchures et on les deplace en Mer Noire . En effet, nous avions choisi de deverser le bassin versant du Danube en Mer Egee, ce qui inversait le sens du flux d'eau en sortie de la Mer Noire (Dardannelles).: la mer Noire est un bassin d'evaporation, mais cette evaporation est compensee (et surpassee) par les runoffs. Si ces runoffs sont envoyes directement en Mediterranee, E-P est positif sur le point de deversement de la mer Noire en Mediterranee. Au contraire, il est negatif si le Danube est correctement envoye dans la mer Noire.