- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/initdynav_loc.f90
r5158 r5159 18 18 USE lmdz_comgeom 19 19 20 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 21 USE lmdz_paramet 20 22 IMPLICIT NONE 21 23 22 ! 24 23 25 ! Routine d'initialisation des ecritures des fichiers histoires LMDZ 24 26 ! au format IOIPSL. Initialisation du fichier histoire moyenne. 25 ! 27 26 28 ! Appels succesifs des routines: histbeg 27 29 ! histhori … … 29 31 ! histdef 30 32 ! histend 31 ! 33 32 34 ! Entree: 33 ! 35 34 36 ! day0,anne0: date de reference 35 37 ! tstep : frequence d'ecriture 36 38 ! t_ops: frequence de l'operation pour IOIPSL 37 39 ! t_wrt: frequence d'ecriture sur le fichier 38 ! 40 39 41 ! Sortie: 40 42 ! fileid: ID du fichier netcdf cree 41 ! 43 42 44 ! L. Fairhead, LMD, 03/99 43 ! 45 44 46 ! ===================================================================== 45 ! 47 46 48 ! Declarations 47 INCLUDE "dimensions.h" 48 INCLUDE "paramet.h" 49 50 49 51 50 52 ! Arguments 51 ! 53 52 54 INTEGER(kind = 4) day0, anne0 53 55 REAL :: tstep, t_ops, t_wrt … … 55 57 ! This routine needs IOIPSL 56 58 ! Variables locales 57 ! 59 58 60 INTEGER :: tau0 59 61 REAL :: zjulian … … 82 84 IF (adjust) return 83 85 84 ! 86 85 87 ! Initialisations 86 ! 88 87 89 pi = 4. * atan (1.) 88 ! 90 89 91 ! Appel a histbeg: creation du fichier netcdf et initialisations diverses 90 92 ! … … 192 194 193 195 194 ! 196 195 197 ! Appel a histvert pour la grille verticale 196 ! 198 197 199 CALL histvert(histaveid, 'presnivs', 'Niveaux Pression& 198 200 & approximatifs', 'mb', llm, presnivs / 100., zvertiid, 'down') … … 202 204 & approximatifs', 'mb', llm, presnivs / 100., zvertiidu, 'down') 203 205 204 ! 206 205 207 ! Appels a histdef pour la definition des variables a sauvegarder 206 ! 208 207 209 ! Vents U 208 ! 210 209 211 jjn = jj_nb 210 212 CALL histdef(histuaveid, 'u', 'vent u moyen ', & … … 212 214 32, 'ave(X)', t_ops, t_wrt) 213 215 214 ! 216 215 217 ! Vents V 216 ! 218 217 219 IF (pole_sud) jjn = jj_nb - 1 218 220 CALL histdef(histvaveid, 'v', 'vent v moyen', & … … 220 222 32, 'ave(X)', t_ops, t_wrt) 221 223 222 ! 224 223 225 ! Temperature 224 ! 226 225 227 jjn = jj_nb 226 228 CALL histdef(histaveid, 'temp', 'temperature moyenne', 'K', & 227 229 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 228 230 32, 'ave(X)', t_ops, t_wrt) 229 ! 231 230 232 ! Temperature potentielle 231 ! 233 232 234 CALL histdef(histaveid, 'theta', 'temperature potentielle', 'K', & 233 235 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & … … 235 237 236 238 237 ! 239 238 240 ! Geopotentiel 239 ! 241 240 242 CALL histdef(histaveid, 'phi', 'geopotentiel moyen', '-', & 241 243 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 242 244 32, 'ave(X)', t_ops, t_wrt) 243 ! 245 244 246 ! Traceurs 245 ! 247 246 248 ! DO iq=1,nqtot 247 249 ! CALL histdef(histaveid, tracers(iq)%name, … … 250 252 ! . 32, 'ave(X)', t_ops, t_wrt) 251 253 ! enddo 252 ! 254 253 255 ! Masse 254 ! 256 255 257 CALL histdef(histaveid, 'masse', 'masse moyenne', 'kg', & 256 258 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 257 259 32, 'ave(X)', t_ops, t_wrt) 258 ! 260 259 261 ! Pression au sol 260 ! 262 261 263 CALL histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', & 262 264 iip1, jjn, thoriid, 1, 1, 1, -99, & 263 265 32, 'ave(X)', t_ops, t_wrt) 264 ! 266 265 267 ! Geopotentiel au sol 266 ! 268 267 269 ! CALL histdef(histaveid, 'phis', 'geopotentiel au sol', '-', 268 270 ! . iip1, jjn, thoriid, 1, 1, 1, -99, 269 271 ! . 32, 'ave(X)', t_ops, t_wrt) 270 ! 272 271 273 ! Fin 272 ! 274 273 275 CALL histend(histaveid) 274 276 CALL histend(histuaveid)
Note: See TracChangeset
for help on using the changeset viewer.