63 CHARACTER(CHRLEN),
PARAMETER,
PRIVATE :: &
64 RCSIdentString =
'$RCSfile: RFLU_ModENSIGHTUtils.F90,v $ $Revision: 1.4 $'
115 LOGICAL,
INTENT(IN) :: emptypartflag
116 TYPE(t_region
),
POINTER :: pregion
122 CHARACTER(80) :: dummystring
123 INTEGER :: errorflag,
i,icl,ifl,ipatch,ivg,ivl,nbvert,nbvertest
124 INTEGER,
DIMENSION(:),
ALLOCATABLE :: vlist
126 TYPE(t_grid),
POINTER :: pgrid
127 TYPE(t_patch),
POINTER :: ppatch
133 global => pregion%global
136 'RFLU_ModENSIGHTUtils.F90')
138 pgrid => pregion%grid
140 global%postPartNumber = global%postPartNumber + 1
147 WRITE(if_ens_geometry) dummystring
148 WRITE(if_ens_geometry) global%postPartNumber
150 WRITE(dummystring,
'(A,I5.5)')
'VOL_',pregion%iRegionGlobal
151 WRITE(if_ens_geometry) dummystring
153 IF ( emptypartflag .EQV. .false. )
THEN
159 dummystring =
'coordinates'
160 WRITE(if_ens_geometry) dummystring
161 WRITE(if_ens_geometry) pgrid%nVertTot
162 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(XCOORD,ivg),KIND=SPREAL), &
163 ivg=1,pgrid%nverttot)
164 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(YCOORD,ivg),KIND=SPREAL), &
165 ivg=1,pgrid%nverttot)
166 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(ZCOORD,ivg),KIND=SPREAL), &
167 ivg=1,pgrid%nverttot)
177 IF ( pgrid%nTets > 0 )
THEN
178 dummystring =
'tetra4'
179 WRITE(if_ens_geometry) dummystring
180 WRITE(if_ens_geometry) pgrid%nTets
181 WRITE(if_ens_geometry) ((pgrid%tet2v(ivl,icl),ivl=1,4),icl=1,pgrid%nTets)
184 IF ( pgrid%nTetsTot > pgrid%nTets )
THEN
185 dummystring =
'g_tetra4'
186 WRITE(if_ens_geometry) dummystring
187 WRITE(if_ens_geometry) pgrid%nTetsTot-pgrid%nTets
188 WRITE(if_ens_geometry) ((pgrid%tet2v(ivl,icl),ivl=1,4), &
189 icl=pgrid%nTets+1,pgrid%nTetsTot)
196 IF ( pgrid%nHexs > 0 )
THEN
197 dummystring =
'hexa8'
198 WRITE(if_ens_geometry) dummystring
199 WRITE(if_ens_geometry) pgrid%nHexs
200 WRITE(if_ens_geometry) ((pgrid%hex2v(ivl,icl),ivl=1,8),icl=1,pgrid%nHexs)
203 IF ( pgrid%nHexsTot > pgrid%nHexs )
THEN
204 dummystring =
'g_hexa8'
205 WRITE(if_ens_geometry) dummystring
206 WRITE(if_ens_geometry) pgrid%nHexsTot-pgrid%nHexs
207 WRITE(if_ens_geometry) ((pgrid%hex2v(ivl,icl),ivl=1,8), &
208 icl=pgrid%nHexs+1,pgrid%nHexsTot)
215 IF ( pgrid%nPris > 0 )
THEN
216 dummystring =
'penta6'
217 WRITE(if_ens_geometry) dummystring
218 WRITE(if_ens_geometry) pgrid%nPris
219 WRITE(if_ens_geometry) ((pgrid%pri2v(ivl,icl),ivl=1,6),icl=1,pgrid%nPris)
222 IF ( pgrid%nPrisTot > pgrid%nPris )
THEN
223 dummystring =
'g_penta6'
224 WRITE(if_ens_geometry) dummystring
225 WRITE(if_ens_geometry) pgrid%nPrisTot-pgrid%nPris
226 WRITE(if_ens_geometry) ((pgrid%pri2v(ivl,icl),ivl=1,6) , &
227 icl=pgrid%nPris+1,pgrid%nPrisTot)
234 IF ( pgrid%nPyrs > 0 )
THEN
235 dummystring =
'pyramid5'
236 WRITE(if_ens_geometry) dummystring
237 WRITE(if_ens_geometry) pgrid%nPyrs
238 WRITE(if_ens_geometry) ((pgrid%pyr2v(ivl,icl),ivl=1,5),icl=1,pgrid%nPyrs)
241 IF ( pgrid%nPyrsTot > pgrid%nPyrs )
THEN
242 dummystring =
'g_pyramid5'
243 WRITE(if_ens_geometry) dummystring
244 WRITE(if_ens_geometry) pgrid%nPyrsTot-pgrid%nPyrs
245 WRITE(if_ens_geometry) ((pgrid%pyr2v(ivl,icl),ivl=1,5), &
246 icl=pgrid%nPyrs+1,pgrid%nPyrsTot)
254 DO ipatch = 1,pgrid%nPatches
255 ppatch => pregion%patches(ipatch)
265 IF ( ppatch%nBTris > 0 )
THEN
266 global%postPartNumber = global%postPartNumber + 1
269 WRITE(if_ens_geometry) dummystring
270 WRITE(if_ens_geometry) global%postPartNumber
271 WRITE(dummystring,
'(A,I3.3,A,I5.5)')
'PAT_',ipatch,
'_TRI-A_', &
272 pregion%iRegionGlobal
273 WRITE(if_ens_geometry) dummystring
275 IF ( emptypartflag .EQV. .false. )
THEN
279 nbvertest = 3*ppatch%nBTris
281 ALLOCATE(vlist(nbvertest),stat=errorflag)
282 global%error = errorflag
283 IF ( global%error /= err_none )
THEN
284 CALL
errorstop(global,err_allocate,__line__,
'vList')
288 3,ppatch%nBTris,vlist,nbvertest,nbvert)
292 dummystring =
'coordinates'
293 WRITE(if_ens_geometry) dummystring
294 WRITE(if_ens_geometry) nbvert
295 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(XCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
296 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(YCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
297 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(ZCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
301 dummystring =
'tria3'
302 WRITE(if_ens_geometry) dummystring
303 WRITE(if_ens_geometry) ppatch%nBTris
306 ppatch%bTri2v(1:3,1:ppatch%nBTris),nbvert, &
309 WRITE(if_ens_geometry) ((ppatch%bTri2v(ivl,ifl),ivl=1,3), &
313 ppatch%bTri2v(1:3,1:ppatch%nBTris),nbvert, &
318 DEALLOCATE(vlist,stat=errorflag)
319 global%error = errorflag
320 IF ( global%error /= err_none )
THEN
321 CALL
errorstop(global,err_deallocate,__line__,
'vList')
330 IF ( ppatch%nBTrisTot > ppatch%nBTris )
THEN
331 global%postPartNumber = global%postPartNumber + 1
334 WRITE(if_ens_geometry) dummystring
335 WRITE(if_ens_geometry) global%postPartNumber
336 WRITE(dummystring,
'(A,I3.3,A,I5.5)')
'PAT_',ipatch,
'_TRI-V_', &
337 pregion%iRegionGlobal
338 WRITE(if_ens_geometry) dummystring
340 IF ( emptypartflag .EQV. .false. )
THEN
344 nbvertest = 3*(ppatch%nBTrisTot-ppatch%nBTris)
346 ALLOCATE(vlist(nbvertest),stat=errorflag)
347 global%error = errorflag
348 IF ( global%error /= err_none )
THEN
349 CALL
errorstop(global,err_allocate,__line__,
'vList')
353 ppatch%bTri2v(1:3,ppatch%nBTris+1:ppatch%nBTrisTot), &
354 3,ppatch%nBTrisTot-ppatch%nBTris,vlist,nbvertest,nbvert)
358 dummystring =
'coordinates'
359 WRITE(if_ens_geometry) dummystring
360 WRITE(if_ens_geometry) nbvert
361 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(XCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
362 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(YCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
363 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(ZCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
367 dummystring =
'g_tria3'
368 WRITE(if_ens_geometry) dummystring
369 WRITE(if_ens_geometry) ppatch%nBTrisTot-ppatch%nBTris
372 ppatch%bTri2v(1:3,ppatch%nBTris+1:ppatch%nBTrisTot),nbvert, &
375 WRITE(if_ens_geometry) ((ppatch%bTri2v(ivl,ifl),ivl=1,3), &
376 ifl=ppatch%nBTris+1,ppatch%nBTrisTot)
379 ppatch%bTri2v(1:3,ppatch%nBTris+1:ppatch%nBTrisTot),nbvert, &
384 DEALLOCATE(vlist,stat=errorflag)
385 global%error = errorflag
386 IF ( global%error /= err_none )
THEN
387 CALL
errorstop(global,err_deallocate,__line__,
'vList')
400 IF ( ppatch%nBQuads > 0 )
THEN
401 global%postPartNumber = global%postPartNumber + 1
404 WRITE(if_ens_geometry) dummystring
405 WRITE(if_ens_geometry) global%postPartNumber
406 WRITE(dummystring,
'(A,I3.3,A,I5.5)')
'PAT_',ipatch,
'_QUAD-A_', &
407 pregion%iRegionGlobal
408 WRITE(if_ens_geometry) dummystring
410 IF ( emptypartflag .EQV. .false. )
THEN
414 nbvertest = 4*ppatch%nBQuads
416 ALLOCATE(vlist(nbvertest),stat=errorflag)
417 global%error = errorflag
418 IF ( global%error /= err_none )
THEN
419 CALL
errorstop(global,err_allocate,__line__,
'vList')
423 ppatch%bQuad2v(1:4,1:ppatch%nBQuads), &
424 4,ppatch%nBQuads,vlist,nbvertest,nbvert)
428 dummystring =
'coordinates'
429 WRITE(if_ens_geometry) dummystring
430 WRITE(if_ens_geometry) nbvert
431 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(XCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
432 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(YCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
433 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(ZCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
437 dummystring =
'quad4'
438 WRITE(if_ens_geometry) dummystring
439 WRITE(if_ens_geometry) ppatch%nBQuads
442 ppatch%bQuad2v(1:4,1:ppatch%nBQuads),nbvert, &
445 WRITE(if_ens_geometry) ((ppatch%bQuad2v(ivl,ifl),ivl=1,4) , &
446 ifl=1,ppatch%nBQuads)
449 ppatch%bQuad2v(1:4,1:ppatch%nBQuads),nbvert, &
454 DEALLOCATE(vlist,stat=errorflag)
455 global%error = errorflag
456 IF ( global%error /= err_none )
THEN
457 CALL
errorstop(global,err_deallocate,__line__,
'vList')
466 IF ( ppatch%nBQuadsTot > ppatch%nBQuads )
THEN
467 global%postPartNumber = global%postPartNumber + 1
470 WRITE(if_ens_geometry) dummystring
471 WRITE(if_ens_geometry) global%postPartNumber
472 WRITE(dummystring,
'(A,I3.3,A,I5.5)')
'PAT_',ipatch,
'_QUAD-V_', &
473 pregion%iRegionGlobal
474 WRITE(if_ens_geometry) dummystring
476 IF ( emptypartflag .EQV. .false. )
THEN
480 nbvertest = 4*(ppatch%nBQuadsTot-ppatch%nBQuads)
482 ALLOCATE(vlist(nbvertest),stat=errorflag)
483 global%error = errorflag
484 IF ( global%error /= err_none )
THEN
485 CALL
errorstop(global,err_allocate,__line__,
'vList')
489 ppatch%bQuad2v(1:4,ppatch%nBQuads+1:ppatch%nBQuadsTot), &
490 4,ppatch%nBQuadsTot-ppatch%nBQuads,vlist,nbvertest,nbvert)
494 dummystring =
'coordinates'
495 WRITE(if_ens_geometry) dummystring
496 WRITE(if_ens_geometry) nbvert
497 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(XCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
498 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(YCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
499 WRITE(if_ens_geometry) (
REAL(pGrid%xyz(ZCOORD,vList(ivl)), &
KIND=SPREAL),ivl=1,nbvert)
503 dummystring =
'g_quad4'
504 WRITE(if_ens_geometry) dummystring
505 WRITE(if_ens_geometry) ppatch%nBQuadsTot-ppatch%nBQuads
508 ppatch%bQuad2v(1:4,ppatch%nBQuads+1:ppatch%nBQuadsTot), &
509 nbvert,vlist(1:nbvert))
511 WRITE(if_ens_geometry) ((ppatch%bQuad2v(ivl,ifl),ivl=1,4), &
512 ifl=ppatch%nBQuads+1,ppatch%nBQuadsTot)
515 ppatch%bQuad2v(1:4,ppatch%nBQuads+1:ppatch%nBQuadsTot), &
516 nbvert,vlist(1:nbvert))
520 DEALLOCATE(vlist,stat=errorflag)
521 global%error = errorflag
522 IF ( global%error /= err_none )
THEN
523 CALL
errorstop(global,err_deallocate,__line__,
'vList')
572 LOGICAL,
INTENT(IN) :: emptypartflag
573 INTEGER,
INTENT(IN) :: ifile
574 REAL(RFREAL),
DIMENSION(:),
POINTER :: var
575 TYPE(t_region
),
POINTER :: pregion
581 CHARACTER(80) :: dummystring
582 INTEGER :: errorflag,icl,ifl,ipatch,offs
584 TYPE(t_grid),
POINTER :: pgrid
585 TYPE(t_patch),
POINTER :: ppatch
591 global => pregion%global
594 'RFLU_ModENSIGHTUtils.F90')
596 pgrid => pregion%grid
598 global%postPartNumber = global%postPartNumber + 1
605 WRITE(ifile) dummystring
606 WRITE(ifile) global%postPartNumber
608 IF ( emptypartflag .EQV. .false. )
THEN
614 IF ( pgrid%nTets > 0 )
THEN
615 dummystring =
'tetra4'
616 WRITE(ifile) dummystring
617 WRITE(ifile) (
REAL(var(pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
621 IF ( pgrid%nTetsTot > pgrid%nTets )
THEN
622 dummystring =
'g_tetra4'
623 WRITE(ifile) dummystring
624 WRITE(ifile) (
REAL(var(pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
625 icl=pgrid%ntets+1,pgrid%ntetstot)
632 IF ( pgrid%nHexs > 0 )
THEN
633 dummystring =
'hexa8'
634 WRITE(ifile) dummystring
635 WRITE(ifile) (
REAL(var(pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
639 IF ( pgrid%nHexsTot > pgrid%nHexs )
THEN
640 dummystring =
'g_hexa8'
641 WRITE(ifile) dummystring
642 WRITE(ifile) (
REAL(var(pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
643 icl=pgrid%nhexs+1,pgrid%nhexstot)
650 IF ( pgrid%nPris > 0 )
THEN
651 dummystring =
'penta6'
652 WRITE(ifile) dummystring
653 WRITE(ifile) (
REAL(var(pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
657 IF ( pgrid%nPrisTot > pgrid%nPris )
THEN
658 dummystring =
'g_penta6'
659 WRITE(ifile) dummystring
660 WRITE(ifile) (
REAL(var(pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
661 icl=pgrid%npris+1,pgrid%npristot)
668 IF ( pgrid%nPyrs > 0 )
THEN
669 dummystring =
'pyramid5'
670 WRITE(ifile) dummystring
671 WRITE(ifile) (
REAL(var(pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
675 IF ( pgrid%nPyrsTot > pgrid%nPyrs )
THEN
676 dummystring =
'g_pyramid5'
677 WRITE(ifile) dummystring
678 WRITE(ifile) (
REAL(var(pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
679 icl=pgrid%npyrs+1,pgrid%npyrstot)
687 DO ipatch = 1,pgrid%nPatches
688 ppatch => pregion%patches(ipatch)
698 IF ( ppatch%nBTris > 0 )
THEN
699 global%postPartNumber = global%postPartNumber + 1
702 WRITE(ifile) dummystring
703 WRITE(ifile) global%postPartNumber
705 IF ( emptypartflag .EQV. .false. )
THEN
706 dummystring =
'tria3'
707 WRITE(ifile) dummystring
708 WRITE(ifile) (
REAL(var(pPatch%bf2c(ifl)),KIND=SPREAL), &
717 IF ( ppatch%nBTrisTot > ppatch%nBTris )
THEN
718 global%postPartNumber = global%postPartNumber + 1
721 WRITE(ifile) dummystring
722 WRITE(ifile) global%postPartNumber
724 offs = ppatch%nBTris + ppatch%nBQuads
726 IF ( emptypartflag .EQV. .false. )
THEN
727 dummystring =
'g_tria3'
728 WRITE(ifile) dummystring
729 WRITE(ifile) (
REAL(var(pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
730 ifl=1,ppatch%nbtristot-ppatch%nbtris)
742 IF ( ppatch%nBQuads > 0 )
THEN
743 global%postPartNumber = global%postPartNumber + 1
746 WRITE(ifile) dummystring
747 WRITE(ifile) global%postPartNumber
751 IF ( emptypartflag .EQV. .false. )
THEN
752 dummystring =
'quad4'
753 WRITE(ifile) dummystring
754 WRITE(ifile) (
REAL(var(pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
755 ifl=1,ppatch%nbquads)
763 IF ( ppatch%nBQuadsTot > ppatch%nBQuads )
THEN
764 global%postPartNumber = global%postPartNumber + 1
767 WRITE(ifile) dummystring
768 WRITE(ifile) global%postPartNumber
770 offs = ppatch%nBTrisTot + ppatch%nBQuads
772 IF ( emptypartflag .EQV. .false. )
THEN
773 dummystring =
'g_quad4'
774 WRITE(ifile) dummystring
775 WRITE(ifile) (
REAL(var(pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
776 ifl=1,ppatch%nbquadstot-ppatch%nbquads)
824 LOGICAL,
INTENT(IN) :: emptypartflag
825 INTEGER,
INTENT(IN) :: ifile
826 REAL(RFREAL),
DIMENSION(:,:),
POINTER :: var
827 TYPE(t_region
),
POINTER :: pregion
833 CHARACTER(80) :: dummystring
834 INTEGER :: errorflag,icl,ifl,ipatch,offs
836 TYPE(t_grid),
POINTER :: pgrid
837 TYPE(t_patch),
POINTER :: ppatch
843 global => pregion%global
846 'RFLU_ModENSIGHTUtils.F90')
848 pgrid => pregion%grid
850 global%postPartNumber = global%postPartNumber + 1
857 WRITE(ifile) dummystring
858 WRITE(ifile) global%postPartNumber
860 IF ( emptypartflag .EQV. .false. )
THEN
866 IF ( pgrid%nTets > 0 )
THEN
867 dummystring =
'tetra4'
868 WRITE(ifile) dummystring
869 WRITE(ifile) (
REAL(var(1,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
871 WRITE(ifile) (
REAL(var(2,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
873 WRITE(ifile) (
REAL(var(3,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
877 IF ( pgrid%nTetsTot > pgrid%nTets )
THEN
878 dummystring =
'g_tetra4'
879 WRITE(ifile) dummystring
880 WRITE(ifile) (
REAL(var(1,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
881 icl=pgrid%ntets+1,pgrid%ntetstot)
882 WRITE(ifile) (
REAL(var(2,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
883 icl=pgrid%ntets+1,pgrid%ntetstot)
884 WRITE(ifile) (
REAL(var(3,pGrid%tet2CellGlob(icl)),KIND=SPREAL), &
885 icl=pgrid%ntets+1,pgrid%ntetstot)
892 IF ( pgrid%nHexs > 0 )
THEN
893 dummystring =
'hexa8'
894 WRITE(ifile) dummystring
895 WRITE(ifile) (
REAL(var(1,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
897 WRITE(ifile) (
REAL(var(2,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
899 WRITE(ifile) (
REAL(var(3,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
903 IF ( pgrid%nHexsTot > pgrid%nHexs )
THEN
904 dummystring =
'g_hexa8'
905 WRITE(ifile) dummystring
906 WRITE(ifile) (
REAL(var(1,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
907 icl=pgrid%nhexs+1,pgrid%nhexstot)
908 WRITE(ifile) (
REAL(var(2,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
909 icl=pgrid%nhexs+1,pgrid%nhexstot)
910 WRITE(ifile) (
REAL(var(3,pGrid%hex2CellGlob(icl)),KIND=SPREAL), &
911 icl=pgrid%nhexs+1,pgrid%nhexstot)
918 IF ( pgrid%nPris > 0 )
THEN
919 dummystring =
'penta6'
920 WRITE(ifile) dummystring
921 WRITE(ifile) (
REAL(var(1,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
923 WRITE(ifile) (
REAL(var(2,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
925 WRITE(ifile) (
REAL(var(3,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
929 IF ( pgrid%nPrisTot > pgrid%nPris )
THEN
930 dummystring =
'g_penta6'
931 WRITE(ifile) dummystring
932 WRITE(ifile) (
REAL(var(1,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
933 icl=pgrid%npris+1,pgrid%npristot)
934 WRITE(ifile) (
REAL(var(2,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
935 icl=pgrid%npris+1,pgrid%npristot)
936 WRITE(ifile) (
REAL(var(3,pGrid%pri2CellGlob(icl)),KIND=SPREAL), &
937 icl=pgrid%npris+1,pgrid%npristot)
944 IF ( pgrid%nPyrs > 0 )
THEN
945 dummystring =
'pyramid5'
946 WRITE(ifile) dummystring
947 WRITE(ifile) (
REAL(var(1,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
949 WRITE(ifile) (
REAL(var(2,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
951 WRITE(ifile) (
REAL(var(3,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
955 IF ( pgrid%nPyrsTot > pgrid%nPyrs )
THEN
956 dummystring =
'g_pyramid5'
957 WRITE(ifile) dummystring
958 WRITE(ifile) (
REAL(var(1,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
959 icl=pgrid%npyrs+1,pgrid%npyrstot)
960 WRITE(ifile) (
REAL(var(2,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
961 icl=pgrid%npyrs+1,pgrid%npyrstot)
962 WRITE(ifile) (
REAL(var(3,pGrid%pyr2CellGlob(icl)),KIND=SPREAL), &
963 icl=pgrid%npyrs+1,pgrid%npyrstot)
971 DO ipatch = 1,pgrid%nPatches
972 ppatch => pregion%patches(ipatch)
982 IF ( ppatch%nBTris > 0 )
THEN
983 global%postPartNumber = global%postPartNumber + 1
986 WRITE(ifile) dummystring
987 WRITE(ifile) global%postPartNumber
989 IF ( emptypartflag .EQV. .false. )
THEN
990 dummystring =
'tria3'
991 WRITE(ifile) dummystring
992 WRITE(ifile) (
REAL(var(1,pPatch%bf2c(ifl)),KIND=SPREAL), &
994 WRITE(ifile) (
REAL(var(2,pPatch%bf2c(ifl)),KIND=SPREAL), &
996 WRITE(ifile) (
REAL(var(3,pPatch%bf2c(ifl)),KIND=SPREAL), &
1005 IF ( ppatch%nBTrisTot > ppatch%nBTris )
THEN
1006 global%postPartNumber = global%postPartNumber + 1
1008 dummystring =
'part'
1009 WRITE(ifile) dummystring
1010 WRITE(ifile) global%postPartNumber
1012 offs = ppatch%nBTris + ppatch%nBQuads
1014 IF ( emptypartflag .EQV. .false. )
THEN
1015 dummystring =
'g_tria3'
1016 WRITE(ifile) dummystring
1017 WRITE(ifile) (
REAL(var(1,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1018 ifl=1,ppatch%nbtristot-ppatch%nbtris)
1019 WRITE(ifile) (
REAL(var(2,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1020 ifl=1,ppatch%nbtristot-ppatch%nbtris)
1021 WRITE(ifile) (
REAL(var(3,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1022 ifl=1,ppatch%nbtristot-ppatch%nbtris)
1034 IF ( ppatch%nBQuads > 0 )
THEN
1035 global%postPartNumber = global%postPartNumber + 1
1037 dummystring =
'part'
1038 WRITE(ifile) dummystring
1039 WRITE(ifile) global%postPartNumber
1041 offs = ppatch%nBTris
1043 IF ( emptypartflag .EQV. .false. )
THEN
1044 dummystring =
'quad4'
1045 WRITE(ifile) dummystring
1046 WRITE(ifile) (
REAL(var(1,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1047 ifl=1,ppatch%nbquads)
1048 WRITE(ifile) (
REAL(var(2,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1049 ifl=1,ppatch%nbquads)
1050 WRITE(ifile) (
REAL(var(3,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1051 ifl=1,ppatch%nbquads)
1059 IF ( ppatch%nBQuadsTot > ppatch%nBQuads )
THEN
1060 global%postPartNumber = global%postPartNumber + 1
1062 dummystring =
'part'
1063 WRITE(ifile) dummystring
1064 WRITE(ifile) global%postPartNumber
1066 offs = ppatch%nBTrisTot + ppatch%nBQuads
1068 IF ( emptypartflag .EQV. .false. )
THEN
1069 dummystring =
'g_quad4'
1070 WRITE(ifile) dummystring
1071 WRITE(ifile) (
REAL(var(1,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1072 ifl=1,ppatch%nbquadstot-ppatch%nbquads)
1073 WRITE(ifile) (
REAL(var(2,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1074 ifl=1,ppatch%nbquadstot-ppatch%nbquads)
1075 WRITE(ifile) (
REAL(var(3,pPatch%bf2c(ifl+offs)),KIND=SPREAL), &
1076 ifl=1,ppatch%nbquadstot-ppatch%nbquads)
1124 subroutine, public rflu_buildconnvertlist(global, connList, connListDim1, connListDim2, vList, vListDimMax, vListDim)
subroutine registerfunction(global, funName, fileName)
subroutine, public rflu_denumberlist(global, listDim1, listDim2, list, keyDim, key)
subroutine, public rflu_ens_writevector(pRegion, var, iFile, emptyPartFlag)
subroutine, public rflu_ens_writegrid(pRegion, emptyPartFlag)
subroutine, public rflu_renumberlist(global, listDim1, listDim2, list, keyDim, key)
subroutine, public rflu_ens_writescalar(pRegion, var, iFile, emptyPartFlag)
subroutine errorstop(global, errorCode, errorLine, addMessage)
subroutine deregisterfunction(global)