Changes between Version 2 and Version 3 of Ticket #21


Ignore:
Timestamp:
Aug 5, 2009, 12:59:28 PM (15 years ago)
Author:
lguez
Comment:

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #21 – Description

    v2 v3  
    1 "create_etat0_limit" est actuellement incompatible avec "gcm" en ce qui concerne la gestion du temps. Le problème vient de la définition de la variable "day_ini". "day_ini" est évalué dans "gcm" au premier jour de la simulation courante et ne change plus au cours de la simulation. Sauf erreur, avant les modifications de Laurent, "day_ini" était le nombre de jours depuis le 1er janvier de "annee_ref". Mais Laurent a semblé considérer que "day_ini" était le nombre de jours depuis la date de référence "(annee_ref, day_ref)". Cela se voit à la ligne 222 de "leapfrog.F" :
     1"create_etat0_limit" est actuellement incompatible avec "gcm" en ce qui concerne la définition de la variable "day_ini". "day_ini" est évalué dans "gcm" au premier jour de la simulation courante et ne change plus au cours de la simulation. Sauf erreur, avant les modifications de Laurent, "day_ini" était le nombre de jours depuis le 1er janvier de "annee_ref". Mais Laurent a semblé considérer que "day_ini" était le nombre de jours depuis la date de référence "(annee_ref, day_ref)". Cela se voit à la ligne 222 de "leapfrog.F" :
    22{{{
    3       jD_cur = jD_ref + (day_ini - 1) + int (itau * dtvr / daysec)
     3jD_cur = jD_ref + (day_ini - 1) + int (itau * dtvr / daysec)
    44}}}
    55Le problème est que "etat0_netcdf" est toujours écrit selon l'ancienne définition de "day_ini". Si bien que si on crée par exemple un état initial avec "create_etat0_limit" avec :
    … …  
    99}}}
    1010calendrier de 360 jours,
    11 on démarre dans "gcm" avec un jour courant au 1er novembre 1981. Et les erreurs qui vont avec ça dans "physiq" pour la position sur l'orbite, etc.
     11on démarre dans "gcm" avec :
     12{{{
     13day_ini = 331
     14}}}
     15d'où un jour courant au 1er novembre 1981. Et les erreurs qui vont avec ça dans "physiq" pour la position sur l'orbite, etc.
    1216
    1317L'ancienne définition de "day_ini" semble aussi avoir cours dans "pbl_surface" et "cpl_init" :