Changeset 1403 for LMDZ4/trunk/libf/bibio/inithist.F
- Timestamp:
- Jul 1, 2010, 11:02:53 AM (14 years ago)
- Location:
- LMDZ4/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk
- Property svn:mergeinfo changed
-
LMDZ4/trunk/libf/bibio/inithist.F
r1279 r1403 2 2 ! $Id$ 3 3 ! 4 subroutine inithist(infile,day0,anne0,tstep,t_ops,t_wrt,fileid, 5 . filevid) 4 subroutine inithist(day0,anne0,tstep,t_ops,t_wrt) 6 5 7 6 #ifdef CPP_IOIPSL … … 9 8 #endif 10 9 USE infotrac, ONLY : nqtot, ttext 10 use com_io_dyn_mod, only : histid,histvid,histuid, & 11 & dynhist_file,dynhistv_file,dynhistu_file 11 12 12 13 implicit none … … 31 32 C nq: nombre de traceurs 32 33 C 33 C Sortie:34 C fileid: ID du fichier netcdf cree35 C filevid:ID du fichier netcdf pour la grille v36 34 C 37 35 C L. Fairhead, LMD, 03/99 … … 54 52 C Arguments 55 53 C 56 character*(*) infile57 54 integer day0, anne0 58 55 real tstep, t_ops, t_wrt 59 integer fileid, filevid60 56 61 57 #ifdef CPP_IOIPSL … … 83 79 tau0 = itau_dyn 84 80 81 ! ------------------------------------------------------------- 82 ! Creation des 3 fichiers pour les grilles horizontales U,V,Scal 83 ! ------------------------------------------------------------- 84 !Grille U 85 85 do jj = 1, jjp1 86 86 do ii = 1, iip1 … … 90 90 enddo 91 91 92 call histbeg( infile, iip1, rlong(:,1), jjp1, rlat(1,:),92 call histbeg(dynhistu_file, iip1, rlong(:,1), jjp1, rlat(1,:), 93 93 . 1, iip1, 1, jjp1, 94 . tau0, zjulian, tstep, uhoriid, fileid) 95 C 96 C Creation du fichier histoire pour la grille en V (oblige pour l'instant, 97 C IOIPSL ne permet pas de grilles avec des nombres de point differents dans 98 C un meme fichier) 94 . tau0, zjulian, tstep, uhoriid, histuid) 99 95 96 ! Grille V 100 97 do jj = 1, jjm 101 98 do ii = 1, iip1 … … 105 102 enddo 106 103 107 call histbeg( 'dyn_histv.nc', iip1, rlong(:,1), jjm, rlat(1,:),104 call histbeg(dynhistv_file, iip1, rlong(:,1), jjm, rlat(1,:), 108 105 . 1, iip1, 1, jjm, 109 . tau0, zjulian, tstep, vhoriid, filevid) 110 C 111 C Appel a histhori pour rajouter les autres grilles horizontales 112 C 106 . tau0, zjulian, tstep, vhoriid, histvid) 107 108 !Grille Scalaire 113 109 do jj = 1, jjp1 114 110 do ii = 1, iip1 … … 118 114 enddo 119 115 120 call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar', 121 . 'Grille points scalaires', thoriid) 116 call histbeg(dynhist_file, iip1, rlong(:,1), jjp1, rlat(1,:), 117 . 1, iip1, 1, jjp1, 118 . tau0, zjulian, tstep, thoriid, histid) 119 ! ------------------------------------------------------------- 120 C Appel a histvert pour la grille verticale 121 ! ------------------------------------------------------------- 122 call histvert(histid, 'presnivs', 'Niveaux pression','mb', 123 . llm, presnivs/100., zvertiid,'down') 124 call histvert(histvid, 'presnivs', 'Niveaux pression','mb', 125 . llm, presnivs/100., zvertiid,'down') 126 call histvert(histuid, 'presnivs', 'Niveaux pression','mb', 127 . llm, presnivs/100., zvertiid,'down') 122 128 C 123 C Appel a histvert pour la grille verticale 124 C 125 call histvert(fileid, 'sig_s', 'Niveaux sigma','-', 126 . llm, nivsigs, zvertiid) 127 C Pour le fichier V 128 call histvert(filevid, 'sig_s', 'Niveaux sigma','-', 129 . llm, nivsigs, zvertiid) 130 C 129 ! ------------------------------------------------------------- 131 130 C Appels a histdef pour la definition des variables a sauvegarder 131 ! ------------------------------------------------------------- 132 132 C 133 133 C Vents U 134 134 C 135 call histdef( fileid, 'ucov', 'vents u covariants', 'm/s',135 call histdef(histuid, 'u', 'vent u', 'm/s', 136 136 . iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, 137 137 . 32, 'inst(X)', t_ops, t_wrt) … … 139 139 C Vents V 140 140 C 141 call histdef( filevid, 'vcov', 'vents v covariants', 'm/s',141 call histdef(histvid, 'v', 'vent v', 'm/s', 142 142 . iip1, jjm, vhoriid, llm, 1, llm, zvertiid, 143 143 . 32, 'inst(X)', t_ops, t_wrt) … … 146 146 C Temperature potentielle 147 147 C 148 call histdef( fileid, 'teta', 'temperature potentielle', '-',148 call histdef(histid, 'teta', 'temperature potentielle', '-', 149 149 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 150 150 . 32, 'inst(X)', t_ops, t_wrt) … … 152 152 C Geopotentiel 153 153 C 154 call histdef( fileid, 'phi', 'geopotentiel instantane', '-',154 call histdef(histid, 'phi', 'geopotentiel', '-', 155 155 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 156 156 . 32, 'inst(X)', t_ops, t_wrt) … … 158 158 C Traceurs 159 159 C 160 DO iq=1,nqtot 161 call histdef(fileid, ttext(iq), ttext(iq), '-', 162 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 163 . 32, 'inst(X)', t_ops, t_wrt) 164 enddo 165 C 160 ! 161 ! DO iq=1,nqtot 162 ! call histdef(histid, ttext(iq), ttext(iq), '-', 163 ! . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 164 ! . 32, 'inst(X)', t_ops, t_wrt) 165 ! enddo 166 !C 166 167 C Masse 167 168 C 168 call histdef( fileid, 'masse', 'masse', 'kg',169 . iip1, jjp1, thoriid, 1, 1, 1, -99,169 call histdef(histid, 'masse', 'masse', 'kg', 170 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 170 171 . 32, 'inst(X)', t_ops, t_wrt) 171 172 C 172 173 C Pression au sol 173 174 C 174 call histdef( fileid, 'ps', 'pression naturelle au sol', 'Pa',175 call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', 175 176 . iip1, jjp1, thoriid, 1, 1, 1, -99, 176 177 . 32, 'inst(X)', t_ops, t_wrt) 177 178 C 178 C Pressionau sol179 C180 call histdef(fileid, 'phis', 'geopotentiel au sol', '-',181 . iip1, jjp1, thoriid, 1, 1, 1, -99,182 . 32, 'inst(X)', t_ops, t_wrt)183 C179 C Geopotentiel au sol 180 !C 181 ! call histdef(histid, 'phis', 'geopotentiel au sol', '-', 182 ! . iip1, jjp1, thoriid, 1, 1, 1, -99, 183 ! . 32, 'inst(X)', t_ops, t_wrt) 184 !C 184 185 C Fin 185 186 C 186 call histend(fileid) 187 call histend(filevid) 187 call histend(histid) 188 call histend(histuid) 189 call histend(histvid) 188 190 #else 189 191 ! tell the user this routine should be run with ioipsl
Note: See TracChangeset
for help on using the changeset viewer.