Changeset 5101 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/initfluxsto_p.F
- Timestamp:
- Jul 23, 2024, 8:22:55 AM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/initfluxsto_p.F
r5099 r5101 123 123 dhe=(/ 0,0 /) 124 124 125 call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe,125 CALL flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, 126 126 . 'box',dynu_domain_id) 127 127 128 callhistbeg(trim(infile),iip1, rlong(:,1), jjn, rlat(1,jjb:jje),128 CALL histbeg(trim(infile),iip1, rlong(:,1), jjn, rlat(1,jjb:jje), 129 129 . 1, iip1, 1, jjn, tau0, zjulian, tstep, uhoriid, 130 130 . fileid,dynu_domain_id) … … 156 156 dhe=(/ 0,0 /) 157 157 158 call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe,158 CALL flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, 159 159 . 'box',dynv_domain_id) 160 160 161 callhistbeg('fluxstokev',iip1, rlong(:,1), jjn, rlat(1,jjb:jje),161 CALL histbeg('fluxstokev',iip1, rlong(:,1), jjn, rlat(1,jjb:jje), 162 162 . 1, iip1, 1, jjn,tau0, zjulian, tstep, vhoriid, 163 163 . filevid,dynv_domain_id) … … 167 167 if (mpi_rank==0) then 168 168 169 callhistbeg('defstoke.nc', 1, rl, 1, rl,169 CALL histbeg('defstoke.nc', 1, rl, 1, rl, 170 170 . 1, 1, 1, 1, 171 171 . tau0, zjulian, tstep, dhoriid, filedid) … … 186 186 jjn=jj_nb 187 187 188 callhisthori(fileid, iip1, rlong(:,jjb:jje),jjn,rlat(:,jjb:jje),188 CALL histhori(fileid, iip1, rlong(:,jjb:jje),jjn,rlat(:,jjb:jje), 189 189 . 'scalar','Grille points scalaires', thoriid) 190 190 … … 192 192 C Appel a histvert pour la grille verticale 193 193 C 194 callhistvert(fileid, 'sig_s', 'Niveaux sigma',194 CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 195 195 . 'sigma_level', 196 196 . llm, nivsigs, zvertiid) 197 197 C Pour le fichier V 198 callhistvert(filevid, 'sig_s', 'Niveaux sigma',198 CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 199 199 . 'sigma_level', 200 200 . llm, nivsigs, zvertiid) … … 202 202 if (mpi_rank==0) then 203 203 nivd(1) = 1 204 callhistvert(filedid, 'sig_s', 'Niveaux sigma',204 CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 205 205 . 'sigma_level', 206 206 . 1, nivd, dvertiid) … … 235 235 C Masse 236 236 C 237 callhistdef(fileid, 'masse', 'Masse', 'kg',237 CALL histdef(fileid, 'masse', 'Masse', 'kg', 238 238 . iip1, jjn, thoriid, llm, 1, llm, zvertiid, 239 239 . 32, 'inst(X)', t_ops, t_wrt) … … 241 241 C Pbaru 242 242 C 243 callhistdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',243 CALL histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s', 244 244 . iip1, jjn, uhoriid, llm, 1, llm, zvertiid, 245 245 . 32, 'inst(X)', t_ops, t_wrt) … … 250 250 if (pole_sud) jjn=jj_nb-1 251 251 252 callhistdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',252 CALL histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s', 253 253 . iip1, jjn, vhoriid, llm, 1, llm, zvertiid, 254 254 . 32, 'inst(X)', t_ops, t_wrt) … … 257 257 C 258 258 if (pole_sud) jjn=jj_nb 259 callhistdef(fileid, 'w', 'flx de masse vert', 'kg m/s',259 CALL histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', 260 260 . iip1, jjn, thoriid, llm, 1, llm, zvertiid, 261 261 . 32, 'inst(X)', t_ops, t_wrt) … … 264 264 C Temperature potentielle 265 265 C 266 callhistdef(fileid, 'teta', 'temperature potentielle', '-',266 CALL histdef(fileid, 'teta', 'temperature potentielle', '-', 267 267 . iip1, jjn, thoriid, llm, 1, llm, zvertiid, 268 268 . 32, 'inst(X)', t_ops, t_wrt) … … 272 272 C Geopotentiel 273 273 C 274 callhistdef(fileid, 'phi', 'geopotentiel instantane', '-',274 CALL histdef(fileid, 'phi', 'geopotentiel instantane', '-', 275 275 . iip1, jjn, thoriid, llm, 1, llm, zvertiid, 276 276 . 32, 'inst(X)', t_ops, t_wrt) … … 278 278 C Fin 279 279 C 280 callhistend(fileid)281 callhistend(filevid)282 if (mpi_rank==0) callhistend(filedid)280 CALL histend(fileid) 281 CALL histend(filevid) 282 if (mpi_rank==0) CALL histend(filedid) 283 283 if (ok_sync) then 284 callhistsync(fileid)285 callhistsync(filevid)286 if (mpi_rank==0) callhistsync(filedid)284 CALL histsync(fileid) 285 CALL histsync(filevid) 286 if (mpi_rank==0) CALL histsync(filedid) 287 287 endif 288 288
Note: See TracChangeset
for help on using the changeset viewer.