Custom Query (96 matches)

Passage de tout le code LMDZ au format Fortran 90. La tâche est décomposée en plusieurs tickets:

Inclusion d'un calendrier réaliste dans LMDZ.

Après discussion avec Frédéric et Lionel et interaction avec Phu, il faudrait intervenir à trois niveaux:

1. la lecture des conditions aux limites
2. le calcul de l'angle solaire
3. les sorties

En ce qui concerne les sorties, vu que IOIPSL ne doit pas pouvoir gérer cela (je ne vois pas comment dire à IOIPSL de sortir des moyennes mensuelles avec des mois à longueur variable dans une simulation annuelle), c'est au niveua du script de lancement que sera géré le problème.

Pour la calcul de l'angle solaire, on récupère la routine de calcul utilisée par LMDZ-planète qui sait gérer de "vraies" années.

Pour les conditions aux limites, l'idée c'est d'avoir un compteur de temps dans la dynamique (genre jour julien) et de faire la transformation au calendrier souhaité dans la physique (en utilisant les routines de gestions de calendrier disponibles dans ioipsl) et de lire les différents fichiers de conditions en calculant le numéro du jour dans l'année (comme actuellement mais avec des "vrais" jours).

On devrait donc disposer des calendriers 360jours, 365 jours et grégorien sans trop de travail supplémentaire.

De Pascale:

Je n'ai pas tout suivi mais au fait, avez vous corrigé les niveaux indiqués dans les fichiers histmth. Si ma memoire est bonne l'axe indique des coordonnées pression alors qu'en fait ce sont des niveaux sigma et du coup tout le monde se plante. Les etudiants ne font pas attention et produisent de jolis graphes faux. En praticulier, Charline a du refaire enormement de choses pour sa thèse lorsque j'ai fini par m'en rendre compte. Pas tres agreable à la fin de thèse...

De Frédéric:

On a discuté plusieurs fois de sortir proprement la coordonnées verticale hybride. Mais je ne me souviens plus pourquoi ce n'est pas fait. Il faut en tous cas garder une pseudo-coordonnée pour sortir les choses simplement en mod interactif. Mais il faudrait indiquer que l'axe est en pseudo-pression par exemple

De Lionel : Je crois aussi que les méta-données actuelles ne sont pas bonnes. Premièrement, une convention de base des fichiers NetCDF est que les variables coordonnées, c'est-à-dire les variables dont le nom est identique au nom de leur unique dimension, sont strictement monotones. Ce qui n'est le cas ni avec "Ahyb" ni avec "Bhyb" ("Bhyb" constant nul aux plus hauts niveaux). On se le fait tout de suite signaler par Ferret quand on ouvre un fichier "hist" actuel. Deuxièmement, on est toujours loin de la convention CF, selon laquelle l'en-tête devrait être de ce type :

netcdf histday {
dimensions:
        time = UNLIMITED ; // (30 currently)
        s = 50 ;
        latitude = 73 ;
        longitude = 96 ;
variables:
        float time(time) ;
                time:units = "days since 1980-01-01" ;
                time:calendar = "360d" ;
                time:long_name = "time" ;
        float s(s) ;
                s:positive = "down" ;
                s:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
                s:long_name = "hybrid sigma-pressure coordinate" ;
                s:formula_terms = "ap: ap b: bp ps: psol p0: p0" ;
        float latitude(latitude) ;
                latitude:units = "degrees_north" ;
                latitude:standard_name = "latitude" ;
        float longitude(longitude) ;
                longitude:units = "degrees_east" ;
                longitude:standard_name = "longitude" ;
        float temp(time, s, latitude, longitude) ;
                temp:units = "K" ;
                temp:standard_name = "air_temperature" ;
                temp:long_name = "air temperature" ;
        float ap(s) ;
                ap:title = "Coefficient A pour hybride" ;
                ap:units = "hPa" ;
        float bp(s) ;
                bp:title = "Coefficient B pour hybride" ;
        float psol(time, latitude, longitude) ;
                psol:units = "hPa" ;
                psol:standard_name = "surface_air_pressure" ;
                psol:long_name = "surface pressure" ;
        float p0 ;
                p0:units = "hPa" ;
                p0:long_name = "reference pressure for hybrid coordinate" ;

// global attributes:
                :Conventions = "CF-1.4" ;
                :institution = "IPSL" ;
data:

 time = 360.5, 361.5, ..., 389.5 ;

 s = 0.994072, ..., 0.0001417363 ;

 latitude = 90, 87.5, ..., -90 ;

 longitude = -180, -176.25, ..., 176.25 ;

 p0 = 500 ;

...
}

Une possibilité qui aurait le mérite d'être propre serait de partir d'une définition du fichier conforme à la convention CF comme ci-dessus et d'ajouter pour le confort des variables "alt" et "presnivs" qui dépendent de la coordonnée hybride "s". Dans Ferret par exemple, on pourrait taper :

shade temp[x=@ave, l=1], y[g=temp], alt

pour voir la température en fonction de la pseudo-altitude.