Nudging

From GHER

(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
Par nudging on entend " rappel vers des champs de climatologies ".
Par nudging on entend " rappel vers des champs de climatologies ".
-
Par exemple si le modèle présente un shift régulier vers des températures trop élevé, on rajoute un terme à l'équation de la température qui va etre plu ou moins fort selon l'écart du champ de température à un champ de référence donné, le terme supplémentaire à la forme de la force de rappel d'un ressort  ( k(x-x_0) )
+
Par exemple si le modèle présente un shift régulier vers des températures trop élevé, on rajoute un terme à l'équation de la température qui va etre plu ou moins fort selon l'écart du champ de température à un champ de référence donné, le terme supplémentaire à la forme de la force de rappel d'un ressort  ( <math>k(x-x_0)</math> )
Par exemple (SOURII est ici le terme source pour la température..)
Par exemple (SOURII est ici le terme source pour la température..)
Line 17: Line 17:
PRENUD(I,J,K-1) est un masque, soit un fichier 3D, rempli de valeur entre 0 et 1 qui permet de restreindre le nudging par exemple juste en surface ou juste sur un certain détroit.. (donnés en tant que FILEEN dans [[fichier .init|fichier .init]])
PRENUD(I,J,K-1) est un masque, soit un fichier 3D, rempli de valeur entre 0 et 1 qui permet de restreindre le nudging par exemple juste en surface ou juste sur un certain détroit.. (donnés en tant que FILEEN dans [[fichier .init|fichier .init]])
 +
 +
 +
----
 +
 +
'''Paramètres'''
 +
 +
*RNUDG : constante de rappel du ''ressort''
 +
 +
*NNUDG : temps de relaxation  !!!  PAS TROUVE OU C ETAIT APPLIQUE POUR L INSTANT  (Arthur) !!!
 +
 +
*KNU
 +
 +
 +
----
 +
 +
interpollation dans le temps
 +
 +
[[main.f90|main.f90]] appelle getnudg(mjdi)    mjdi temps présent
 +
 +
getnufg interpolle les champs de nudging dans le temps et renvoie les tnud et snud (champ 3D de réf) correspondant
 +
 +
----

Revision as of 14:28, 8 September 2008

Par nudging on entend " rappel vers des champs de climatologies ".

Par exemple si le modèle présente un shift régulier vers des températures trop élevé, on rajoute un terme à l'équation de la température qui va etre plu ou moins fort selon l'écart du champ de température à un champ de référence donné, le terme supplémentaire à la forme de la force de rappel d'un ressort ( k(xx0) )

Par exemple (SOURII est ici le terme source pour la température..)

     DO K=2,KSUPII
       DO I=IINF,ISUP-1
         SOURII(I,K) = ...- RNUDG*( HTII(I,J,K) - HII(I,J)*TNUD(I,J,K-1) )* PRENUD(I,J,K-1)


ici RNUDG est une constante de type constante de ressort données dans le fichier .init ,

HTII(I,J,K) est le champ de température (multiplié par la hauteur HII(i,j),

TNUD(I,J,K-1) est le champ de température de référence qui est interpollé à partir des FILETN (pour la T et FILESN pour la salinité) (donnés dans fichier .init)

PRENUD(I,J,K-1) est un masque, soit un fichier 3D, rempli de valeur entre 0 et 1 qui permet de restreindre le nudging par exemple juste en surface ou juste sur un certain détroit.. (donnés en tant que FILEEN dans fichier .init)



Paramètres

  • RNUDG : constante de rappel du ressort
  • NNUDG : temps de relaxation  !!! PAS TROUVE OU C ETAIT APPLIQUE POUR L INSTANT (Arthur) !!!
  • KNU



interpollation dans le temps

main.f90 appelle getnudg(mjdi) mjdi temps présent

getnufg interpolle les champs de nudging dans le temps et renvoie les tnud et snud (champ 3D de réf) correspondant


Personal tools