Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/inithist_loc.F90
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/inithist_loc.F90
r5158 r5159 17 17 USE lmdz_comgeom 18 18 19 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 20 USE lmdz_paramet 19 21 IMPLICIT NONE 20 22 21 ! 23 22 24 ! Routine d'initialisation des ecritures des fichiers histoires LMDZ 23 25 ! au format IOIPSL 24 ! 26 25 27 ! Appels succesifs des routines: histbeg 26 28 ! histhori … … 28 30 ! histdef 29 31 ! histend 30 ! 32 31 33 ! Entree: 32 ! 34 33 35 ! day0,anne0: date de reference 34 36 ! tstep: duree du pas de temps en seconde … … 36 38 ! t_wrt: frequence d'ecriture sur le fichier 37 39 ! nq: nombre de traceurs 38 ! 39 ! 40 41 40 42 ! L. Fairhead, LMD, 03/99 41 ! 43 42 44 ! ===================================================================== 43 ! 45 44 46 ! Declarations 45 INCLUDE "dimensions.h" 46 INCLUDE "paramet.h" 47 48 47 49 48 50 ! Arguments 49 ! 51 50 52 INTEGER :: day0, anne0 51 53 REAL :: tstep, t_ops, t_wrt … … 53 55 ! This routine needs IOIPSL 54 56 ! Variables locales 55 ! 57 56 58 INTEGER :: tau0 57 59 REAL :: zjulian … … 80 82 IF (adjust) return 81 83 82 ! 84 83 85 ! Initialisations 84 ! 86 85 87 pi = 4. * atan (1.) 86 ! 88 87 89 ! Appel a histbeg: creation du fichier netcdf et initialisations diverses 88 90 ! … … 200 202 llm, presnivs / 100., zvertiidu, 'down') 201 203 202 ! 204 203 205 ! ------------------------------------------------------------- 204 206 ! Appels a histdef pour la definition des variables a sauvegarder 205 207 ! ------------------------------------------------------------- 206 ! 208 207 209 ! Vents U 208 ! 210 209 211 jjn = jj_nb 210 212 CALL histdef(histuid, 'u', 'vent u', & … … 212 214 32, 'inst(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(histvid, 'v', 'vent v', & … … 220 222 32, 'inst(X)', t_ops, t_wrt) 221 223 222 ! 224 223 225 ! Temperature 224 ! 226 225 227 jjn = jj_nb 226 228 CALL histdef(histid, 'temp', 'temperature', 'K', & 227 229 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 228 230 32, 'inst(X)', t_ops, t_wrt) 229 ! 231 230 232 ! Temperature potentielle 231 ! 233 232 234 CALL histdef(histid, '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(histid, 'phi', 'geopotentiel', '-', & 241 243 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 242 244 32, 'inst(X)', t_ops, t_wrt) 243 ! 245 244 246 ! Traceurs 245 ! 247 246 248 ! DO iq=1,nqtot 247 249 ! CALL histdef(histid, tracers(iq)%name, … … 250 252 ! . 32, 'inst(X)', t_ops, t_wrt) 251 253 ! enddo 252 ! 254 253 255 ! Masse 254 ! 256 255 257 CALL histdef(histid, 'masse', 'masse', 'kg', & 256 258 iip1, jjn, thoriid, llm, 1, llm, zvertiid, & 257 259 32, 'inst(X)', t_ops, t_wrt) 258 ! 260 259 261 ! Pression au sol 260 ! 262 261 263 CALL histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', & 262 264 iip1, jjn, thoriid, 1, 1, 1, -99, & 263 265 32, 'inst(X)', t_ops, t_wrt) 264 ! 266 265 267 ! Geopotentiel au sol 266 ! 268 267 269 ! CALL histdef(histid, 'phis', 'geopotentiel au sol', '-', 268 270 ! . iip1, jjn, thoriid, 1, 1, 1, -99, 269 271 ! . 32, 'inst(X)', t_ops, t_wrt) 270 ! 272 271 273 ! Fin 272 ! 274 273 275 CALL histend(histid) 274 276 CALL histend(histuid)
Note: See TracChangeset
for help on using the changeset viewer.