Changeset 1357 for LMDZ4/branches/LMDZ4V5.0-dev/libf/bibio/initdynav.F
- Timestamp:
- Apr 14, 2010, 4:03:19 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/bibio/initdynav.F
r1279 r1357 2 2 ! $Id$ 3 3 ! 4 subroutine initdynav(infile,day0,anne0,tstep,t_ops,t_wrt 5 . ,fileid) 4 subroutine initdynav(day0,anne0,tstep,t_ops,t_wrt) 6 5 7 6 #ifdef CPP_IOIPSL … … 9 8 #endif 10 9 USE infotrac, ONLY : nqtot, ttext 11 10 use com_io_dyn_mod, only : histaveid,histvaveid,histuaveid, & 11 & dynhistave_file,dynhistvave_file,dynhistuave_file 12 12 implicit none 13 13 … … 30 30 C t_wrt: frequence d'ecriture sur le fichier 31 31 C 32 C Sortie:33 C fileid: ID du fichier netcdf cree34 32 C 35 33 C L. Fairhead, LMD, 03/99 … … 52 50 C Arguments 53 51 C 54 character*(*) infile55 52 integer day0, anne0 56 53 real tstep, t_ops, t_wrt 57 integer fileid58 54 59 55 #ifdef CPP_IOIPSL … … 61 57 C Variables locales 62 58 C 63 integer thoriid, zvertiid64 59 integer tau0 65 60 real zjulian 66 61 integer iq 67 62 real rlong(iip1,jjp1), rlat(iip1,jjp1) 63 integer uhoriid, vhoriid, thoriid, zvertiid 68 64 integer ii,jj 69 65 integer zan, dayref … … 88 84 enddo 89 85 90 call histbeg(infile, iip1, rlong(:,1), jjp1, rlat(1,:), 86 ! Creation de 3 fichiers pour les differentes grilles horizontales 87 ! Restriction de IOIPSL: seulement 2 coordonnees dans le meme fichier 88 ! Grille Scalaire 89 call histbeg(dynhistave_file, iip1, rlong(:,1), jjp1, rlat(1,:), 91 90 . 1, iip1, 1, jjp1, 92 . tau0, zjulian, tstep, thoriid, fileid) 93 91 . tau0, zjulian, tstep, thoriid,histaveid) 92 93 C Creation du fichier histoire pour les grilles en V et U (oblige pour l'instant, 94 C IOIPSL ne permet pas de grilles avec des nombres de point differents dans 95 C un meme fichier) 96 ! Grille V 97 do jj = 1, jjm 98 do ii = 1, iip1 99 rlong(ii,jj) = rlonv(ii) * 180. / pi 100 rlat(ii,jj) = rlatv(jj) * 180. / pi 101 enddo 102 enddo 103 104 call histbeg(dynhistvave_file, iip1, rlong(:,1), jjm, rlat(1,:), 105 . 1, iip1, 1, jjm, 106 . tau0, zjulian, tstep, vhoriid,histvaveid) 107 ! Grille U 108 do jj = 1, jjp1 109 do ii = 1, iip1 110 rlong(ii,jj) = rlonu(ii) * 180. / pi 111 rlat(ii,jj) = rlatu(jj) * 180. / pi 112 enddo 113 enddo 114 115 call histbeg(dynhistuave_file, iip1, rlong(:,1),jjp1, rlat(1,:), 116 . 1, iip1, 1, jjp1, 117 . tau0, zjulian, tstep, uhoriid,histuaveid) 94 118 C 95 119 C Appel a histvert pour la grille verticale 96 120 C 97 call histvert(fileid, 'sigss', 'Niveaux sigma','Pa', 98 . llm, nivsigs, zvertiid) 121 call histvert(histaveid,'presnivs','Niveaux Pression 122 & approximatifs','mb',llm, presnivs/100., zvertiid,'down') 123 call histvert(histuaveid,'presnivs','Niveaux Pression 124 & approximatifs','mb',llm, presnivs/100., zvertiid,'down') 125 call histvert(histvaveid,'presnivs','Niveaux Pression 126 & approximatifs','mb',llm, presnivs/100., zvertiid,'down') 99 127 C 100 128 C Appels a histdef pour la definition des variables a sauvegarder … … 102 130 C Vents U 103 131 C 104 write(6,*)'inithistave',tstep 105 call histdef(fileid, 'u', 'vents u scalaires moyennes', 106 . 'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 107 . 32, 'ave(X)', t_ops, t_wrt) 108 109 C 132 ! write(6,*)'inithistave',tstep 133 call histdef(histuaveid, 'u', 'vent u moyen ', 134 . 'm/s', iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, 135 . 32, 'ave(X)', t_ops, t_wrt) 136 110 137 C Vents V 111 138 C 112 call histdef( fileid, 'v', 'vents v scalaires moyennes',113 . 'm/s', iip1, jj p1, thoriid, llm, 1, llm, zvertiid,139 call histdef(histvaveid, 'v', 'vent v moyen', 140 . 'm/s', iip1, jjm, vhoriid, llm, 1, llm, zvertiid, 114 141 . 32, 'ave(X)', t_ops, t_wrt) 115 142 … … 117 144 C Temperature 118 145 C 119 call histdef( fileid, 'temp', 'temperature moyennee', 'K',146 call histdef(histaveid, 'temp', 'temperature moyenne', 'K', 120 147 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 121 148 . 32, 'ave(X)', t_ops, t_wrt) … … 123 150 C Temperature potentielle 124 151 C 125 call histdef(fileid, 'theta', 'temperature potentielle', 'K', 126 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 127 . 32, 'ave(X)', t_ops, t_wrt) 128 129 152 call histdef(histaveid, 'theta', 'temperature potentielle', 'K', 153 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 154 . 32, 'ave(X)', t_ops, t_wrt) 130 155 C 131 156 C Geopotentiel 132 157 C 133 call histdef( fileid, 'phi', 'geopotentiel moyenne', '-',158 call histdef(histaveid, 'phi', 'geopotentiel moyen', '-', 134 159 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 135 160 . 32, 'ave(X)', t_ops, t_wrt) … … 137 162 C Traceurs 138 163 C 139 DO iq=1,nqtot140 call histdef(fileid, ttext(iq), ttext(iq), '-',141 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid,142 . 32, 'ave(X)', t_ops, t_wrt)143 enddo164 ! DO iq=1,nqtot 165 ! call histdef(histaveid, ttext(iq), ttext(iq), '-', 166 ! . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 167 ! . 32, 'ave(X)', t_ops, t_wrt) 168 ! enddo 144 169 C 145 170 C Masse 146 171 C 147 call histdef(fileid, 'masse', 'masse', 'kg', 172 call histdef(histaveid, 'masse', 'masse', 'kg', 173 . iip1, jjp1, thoriid, llm, 1, llm, zvertiid, 174 . 32, 'ave(X)', t_ops, t_wrt) 175 C 176 C Pression au sol 177 C 178 call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', 148 179 . iip1, jjp1, thoriid, 1, 1, 1, -99, 149 180 . 32, 'ave(X)', t_ops, t_wrt) 150 181 C 151 C Pression au sol 152 C 153 call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', 154 . iip1, jjp1, thoriid, 1, 1, 1, -99, 155 . 32, 'ave(X)', t_ops, t_wrt) 156 C 157 C Pression au sol 158 C 159 call histdef(fileid, 'phis', 'geopotentiel au sol', '-', 160 . iip1, jjp1, thoriid, 1, 1, 1, -99, 161 . 32, 'ave(X)', t_ops, t_wrt) 162 C 182 C Geopotentiel au sol 183 C 184 ! call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', 185 ! . iip1, jjp1, thoriid, 1, 1, 1, -99, 186 ! . 32, 'ave(X)', t_ops, t_wrt) 187 !C 163 188 C Fin 164 189 C 165 call histend(fileid) 190 call histend(histaveid) 191 call histend(histuaveid) 192 call histend(histvaveid) 166 193 #else 167 194 ! tell the user this routine should be run with ioipsl
Note: See TracChangeset
for help on using the changeset viewer.