[1632] | 1 | ! $Id: initdynav_p.F 1279 2009-12-10 09:02:56Z fairhead $ |
---|
[5099] | 2 | |
---|
[5118] | 3 | SUBROUTINE initdynav_loc(day0, anne0, tstep, t_ops, t_wrt) |
---|
[1632] | 4 | |
---|
[5105] | 5 | ! This routine needs IOIPSL |
---|
[5118] | 6 | USE IOIPSL |
---|
| 7 | USE parallel_lmdz |
---|
| 8 | USE lmdz_write_field |
---|
| 9 | USE misc_mod |
---|
| 10 | ! USE infotrac |
---|
| 11 | USE com_io_dyn_mod, ONLY: histaveid, histvaveid, histuaveid, & |
---|
| 12 | dynhistave_file, dynhistvave_file, dynhistuave_file |
---|
| 13 | USE comconst_mod, ONLY: pi |
---|
| 14 | USE comvert_mod, ONLY: presnivs |
---|
| 15 | USE temps_mod, ONLY: itau_dyn |
---|
| 16 | USE lmdz_description, ONLY: descript |
---|
| 17 | USE lmdz_iniprint, ONLY: lunout, prt_level |
---|
[5136] | 18 | USE lmdz_comgeom |
---|
[1632] | 19 | |
---|
[5159] | 20 | USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm |
---|
| 21 | USE lmdz_paramet |
---|
[5118] | 22 | IMPLICIT NONE |
---|
[1632] | 23 | |
---|
[5159] | 24 | |
---|
[5105] | 25 | ! Routine d'initialisation des ecritures des fichiers histoires LMDZ |
---|
| 26 | ! au format IOIPSL. Initialisation du fichier histoire moyenne. |
---|
[5159] | 27 | |
---|
[5105] | 28 | ! Appels succesifs des routines: histbeg |
---|
| 29 | ! histhori |
---|
| 30 | ! histver |
---|
| 31 | ! histdef |
---|
| 32 | ! histend |
---|
[5159] | 33 | |
---|
[5105] | 34 | ! Entree: |
---|
[5159] | 35 | |
---|
[5105] | 36 | ! day0,anne0: date de reference |
---|
| 37 | ! tstep : frequence d'ecriture |
---|
| 38 | ! t_ops: frequence de l'operation pour IOIPSL |
---|
| 39 | ! t_wrt: frequence d'ecriture sur le fichier |
---|
[5159] | 40 | |
---|
[5105] | 41 | ! Sortie: |
---|
| 42 | ! fileid: ID du fichier netcdf cree |
---|
[5159] | 43 | |
---|
[5105] | 44 | ! L. Fairhead, LMD, 03/99 |
---|
[5159] | 45 | |
---|
[5105] | 46 | ! ===================================================================== |
---|
[5159] | 47 | |
---|
[5105] | 48 | ! Declarations |
---|
[1632] | 49 | |
---|
[5159] | 50 | |
---|
| 51 | |
---|
[5105] | 52 | ! Arguments |
---|
[5159] | 53 | |
---|
[5118] | 54 | INTEGER(kind = 4) day0, anne0 |
---|
[5116] | 55 | REAL :: tstep, t_ops, t_wrt |
---|
[1632] | 56 | |
---|
[5105] | 57 | ! This routine needs IOIPSL |
---|
| 58 | ! Variables locales |
---|
[5159] | 59 | |
---|
[5116] | 60 | INTEGER :: tau0 |
---|
| 61 | REAL :: zjulian |
---|
| 62 | INTEGER :: iq |
---|
[5118] | 63 | REAL :: rlong(iip1, jjp1), rlat(iip1, jjp1) |
---|
[5116] | 64 | INTEGER :: uhoriid, vhoriid, thoriid |
---|
[5118] | 65 | INTEGER :: zvertiid, zvertiidv, zvertiidu |
---|
| 66 | INTEGER :: ii, jj |
---|
[5116] | 67 | INTEGER :: zan, dayref |
---|
[5118] | 68 | INTEGER :: jjb, jje, jjn |
---|
[1632] | 69 | |
---|
[5105] | 70 | ! definition du domaine d'ecriture pour le rebuild |
---|
[5103] | 71 | |
---|
[5118] | 72 | INTEGER, DIMENSION(2) :: ddid |
---|
| 73 | INTEGER, DIMENSION(2) :: dsg |
---|
| 74 | INTEGER, DIMENSION(2) :: dsl |
---|
| 75 | INTEGER, DIMENSION(2) :: dpf |
---|
| 76 | INTEGER, DIMENSION(2) :: dpl |
---|
| 77 | INTEGER, DIMENSION(2) :: dhs |
---|
| 78 | INTEGER, DIMENSION(2) :: dhe |
---|
[5103] | 79 | |
---|
[5105] | 80 | INTEGER :: dynhistave_domain_id |
---|
| 81 | INTEGER :: dynhistvave_domain_id |
---|
| 82 | INTEGER :: dynhistuave_domain_id |
---|
[1632] | 83 | |
---|
[5117] | 84 | IF (adjust) return |
---|
[1632] | 85 | |
---|
[5159] | 86 | |
---|
[5105] | 87 | ! Initialisations |
---|
[5159] | 88 | |
---|
[5105] | 89 | pi = 4. * atan (1.) |
---|
[5159] | 90 | |
---|
[5105] | 91 | ! Appel a histbeg: creation du fichier netcdf et initialisations diverses |
---|
| 92 | ! |
---|
[5103] | 93 | |
---|
[5105] | 94 | zan = anne0 |
---|
| 95 | dayref = day0 |
---|
| 96 | CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) |
---|
| 97 | tau0 = itau_dyn |
---|
[1632] | 98 | |
---|
[5158] | 99 | DO jj = 1, jjp1 |
---|
| 100 | DO ii = 1, iip1 |
---|
[5118] | 101 | rlong(ii, jj) = rlonv(ii) * 180. / pi |
---|
| 102 | rlat(ii, jj) = rlatu(jj) * 180. / pi |
---|
[5105] | 103 | enddo |
---|
| 104 | enddo |
---|
[1632] | 105 | |
---|
| 106 | |
---|
[5105] | 107 | ! Creation de 3 fichiers pour les differentes grilles horizontales |
---|
| 108 | ! Restriction de IOIPSL: seulement 2 coordonnees dans le meme fichier |
---|
| 109 | ! Grille Scalaire |
---|
[1632] | 110 | |
---|
[5118] | 111 | jjb = jj_begin |
---|
| 112 | jje = jj_end |
---|
| 113 | jjn = jj_nb |
---|
[1632] | 114 | |
---|
[5118] | 115 | ddid = (/ 1, 2 /) |
---|
| 116 | dsg = (/ iip1, jjp1 /) |
---|
| 117 | dsl = (/ iip1, jjn /) |
---|
| 118 | dpf = (/ 1, jjb /) |
---|
| 119 | dpl = (/ iip1, jje /) |
---|
| 120 | dhs = (/ 0, 0 /) |
---|
| 121 | dhe = (/ 0, 0 /) |
---|
[1632] | 122 | |
---|
[5118] | 123 | CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, & |
---|
| 124 | 'box', dynhistave_domain_id) |
---|
[5103] | 125 | |
---|
[5118] | 126 | CALL histbeg(dynhistave_file, iip1, rlong(:, 1), jjn, & |
---|
| 127 | rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, & |
---|
| 128 | zjulian, tstep, thoriid, & |
---|
| 129 | histaveid, dynhistave_domain_id) |
---|
[1632] | 130 | |
---|
| 131 | |
---|
[5105] | 132 | ! Creation du fichier histoire pour les grilles en V et U (oblige pour l'instant, |
---|
| 133 | ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans |
---|
| 134 | ! un meme fichier) |
---|
| 135 | ! Grille V |
---|
[5103] | 136 | |
---|
[5118] | 137 | jjb = jj_begin |
---|
| 138 | jje = jj_end |
---|
| 139 | jjn = jj_nb |
---|
| 140 | IF (pole_sud) jjn = jjn - 1 |
---|
| 141 | IF (pole_sud) jje = jje - 1 |
---|
[1632] | 142 | |
---|
[5158] | 143 | DO jj = jjb, jje |
---|
| 144 | DO ii = 1, iip1 |
---|
[5118] | 145 | rlong(ii, jj) = rlonv(ii) * 180. / pi |
---|
| 146 | rlat(ii, jj) = rlatv(jj) * 180. / pi |
---|
[5105] | 147 | enddo |
---|
| 148 | enddo |
---|
[1632] | 149 | |
---|
[5118] | 150 | ddid = (/ 1, 2 /) |
---|
| 151 | dsg = (/ iip1, jjm /) |
---|
| 152 | dsl = (/ iip1, jjn /) |
---|
| 153 | dpf = (/ 1, jjb /) |
---|
| 154 | dpl = (/ iip1, jje /) |
---|
| 155 | dhs = (/ 0, 0 /) |
---|
| 156 | dhe = (/ 0, 0 /) |
---|
[1632] | 157 | |
---|
[5118] | 158 | CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, & |
---|
| 159 | 'box', dynhistvave_domain_id) |
---|
[1632] | 160 | |
---|
[5118] | 161 | CALL histbeg(dynhistvave_file, iip1, rlong(:, 1), jjn, & |
---|
| 162 | rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, & |
---|
| 163 | zjulian, tstep, vhoriid, & |
---|
| 164 | histvaveid, dynhistvave_domain_id) |
---|
[5103] | 165 | |
---|
[5105] | 166 | ! Grille U |
---|
[2475] | 167 | |
---|
[5158] | 168 | DO jj = 1, jjp1 |
---|
| 169 | DO ii = 1, iip1 |
---|
[5118] | 170 | rlong(ii, jj) = rlonu(ii) * 180. / pi |
---|
| 171 | rlat(ii, jj) = rlatu(jj) * 180. / pi |
---|
[5105] | 172 | enddo |
---|
| 173 | enddo |
---|
[1632] | 174 | |
---|
[5118] | 175 | jjb = jj_begin |
---|
| 176 | jje = jj_end |
---|
| 177 | jjn = jj_nb |
---|
[1632] | 178 | |
---|
[5118] | 179 | ddid = (/ 1, 2 /) |
---|
| 180 | dsg = (/ iip1, jjp1 /) |
---|
| 181 | dsl = (/ iip1, jjn /) |
---|
| 182 | dpf = (/ 1, jjb /) |
---|
| 183 | dpl = (/ iip1, jje /) |
---|
| 184 | dhs = (/ 0, 0 /) |
---|
| 185 | dhe = (/ 0, 0 /) |
---|
[1632] | 186 | |
---|
[5118] | 187 | CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, & |
---|
| 188 | 'box', dynhistuave_domain_id) |
---|
[5103] | 189 | |
---|
[5118] | 190 | CALL histbeg(dynhistuave_file, iip1, rlong(:, 1), jjn, & |
---|
| 191 | rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, & |
---|
| 192 | zjulian, tstep, uhoriid, & |
---|
| 193 | histuaveid, dynhistuave_domain_id) |
---|
[5103] | 194 | |
---|
| 195 | |
---|
[5159] | 196 | |
---|
[5105] | 197 | ! Appel a histvert pour la grille verticale |
---|
[5159] | 198 | |
---|
[5118] | 199 | CALL histvert(histaveid, 'presnivs', 'Niveaux Pression& |
---|
| 200 | & approximatifs', 'mb', llm, presnivs / 100., zvertiid, 'down') |
---|
| 201 | CALL histvert(histuaveid, 'presnivs', 'Niveaux Pression& |
---|
| 202 | & approximatifs', 'mb', llm, presnivs / 100., zvertiidv, 'down') |
---|
| 203 | CALL histvert(histvaveid, 'presnivs', 'Niveaux Pression& |
---|
| 204 | & approximatifs', 'mb', llm, presnivs / 100., zvertiidu, 'down') |
---|
[1632] | 205 | |
---|
[5159] | 206 | |
---|
[5105] | 207 | ! Appels a histdef pour la definition des variables a sauvegarder |
---|
[5159] | 208 | |
---|
[5105] | 209 | ! Vents U |
---|
[5159] | 210 | |
---|
[5118] | 211 | jjn = jj_nb |
---|
[5105] | 212 | CALL histdef(histuaveid, 'u', 'vent u moyen ', & |
---|
[5118] | 213 | 'm/s', iip1, jjn, uhoriid, llm, 1, llm, zvertiidu, & |
---|
| 214 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[1632] | 215 | |
---|
[5159] | 216 | |
---|
[5105] | 217 | ! Vents V |
---|
[5159] | 218 | |
---|
[5118] | 219 | IF (pole_sud) jjn = jj_nb - 1 |
---|
[5105] | 220 | CALL histdef(histvaveid, 'v', 'vent v moyen', & |
---|
[5118] | 221 | 'm/s', iip1, jjn, vhoriid, llm, 1, llm, zvertiidv, & |
---|
| 222 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[1632] | 223 | |
---|
[5159] | 224 | |
---|
[5105] | 225 | ! Temperature |
---|
[5159] | 226 | |
---|
[5118] | 227 | jjn = jj_nb |
---|
[5105] | 228 | CALL histdef(histaveid, 'temp', 'temperature moyenne', 'K', & |
---|
[5118] | 229 | iip1, jjn, thoriid, llm, 1, llm, zvertiid, & |
---|
| 230 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[5159] | 231 | |
---|
[5105] | 232 | ! Temperature potentielle |
---|
[5159] | 233 | |
---|
[5105] | 234 | CALL histdef(histaveid, 'theta', 'temperature potentielle', 'K', & |
---|
[5118] | 235 | iip1, jjn, thoriid, llm, 1, llm, zvertiid, & |
---|
| 236 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[1632] | 237 | |
---|
[5105] | 238 | |
---|
[5159] | 239 | |
---|
[5105] | 240 | ! Geopotentiel |
---|
[5159] | 241 | |
---|
[5105] | 242 | CALL histdef(histaveid, 'phi', 'geopotentiel moyen', '-', & |
---|
[5118] | 243 | iip1, jjn, thoriid, llm, 1, llm, zvertiid, & |
---|
| 244 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[5159] | 245 | |
---|
[5105] | 246 | ! Traceurs |
---|
[5159] | 247 | |
---|
[5105] | 248 | ! DO iq=1,nqtot |
---|
| 249 | ! CALL histdef(histaveid, tracers(iq)%name, |
---|
| 250 | ! . tracers(iq)%longName, '-', |
---|
| 251 | ! . iip1, jjn, thoriid, llm, 1, llm, zvertiid, |
---|
| 252 | ! . 32, 'ave(X)', t_ops, t_wrt) |
---|
| 253 | ! enddo |
---|
[5159] | 254 | |
---|
[5105] | 255 | ! Masse |
---|
[5159] | 256 | |
---|
[5105] | 257 | CALL histdef(histaveid, 'masse', 'masse moyenne', 'kg', & |
---|
[5118] | 258 | iip1, jjn, thoriid, llm, 1, llm, zvertiid, & |
---|
| 259 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[5159] | 260 | |
---|
[5105] | 261 | ! Pression au sol |
---|
[5159] | 262 | |
---|
[5105] | 263 | CALL histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', & |
---|
[5118] | 264 | iip1, jjn, thoriid, 1, 1, 1, -99, & |
---|
| 265 | 32, 'ave(X)', t_ops, t_wrt) |
---|
[5159] | 266 | |
---|
[5105] | 267 | ! Geopotentiel au sol |
---|
[5159] | 268 | |
---|
[5105] | 269 | ! CALL histdef(histaveid, 'phis', 'geopotentiel au sol', '-', |
---|
| 270 | ! . iip1, jjn, thoriid, 1, 1, 1, -99, |
---|
| 271 | ! . 32, 'ave(X)', t_ops, t_wrt) |
---|
[5159] | 272 | |
---|
[5105] | 273 | ! Fin |
---|
[5159] | 274 | |
---|
[5105] | 275 | CALL histend(histaveid) |
---|
| 276 | CALL histend(histuaveid) |
---|
| 277 | CALL histend(histvaveid) |
---|
[5116] | 278 | END SUBROUTINE initdynav_loc |
---|