My Project
Data Types | Functions/Subroutines | Variables
mod_par Module Reference

Data Types

interface  aprint
 
interface  pprint
 

Functions/Subroutines

subroutine aprint_vec (IUNIT, VARP, VART, NOW, ILOC, MSG)
 
subroutine pprint_vec (IUNIT, VARP, VART, NOW, ILOC, MSG)
 
subroutine aprint_arr (IUNIT, VARP, VART, NOW, ILOC, K1, K2, MSG)
 
subroutine pprint_arr (IUNIT, VARP, VART, NOW, ILOC, K1, K2, MSG)
 

Variables

integer, dimension(:), pointer el_pid
 
integer, dimension(:), pointer elid
 
integer, dimension(:), pointer nlid
 
integer, dimension(:), pointer elid_x
 
integer, dimension(:), pointer nlid_x
 
integer, dimension(:), pointer egid
 
integer, dimension(:), pointer ngid
 
integer, dimension(:), pointer egid_x
 
integer, dimension(:), pointer ngid_x
 
integer nhe
 
integer, dimension(:), pointer he_lst
 
integer, dimension(:), pointer he_own
 
integer nbn
 
integer mx_mlt
 
integer, dimension(:), pointer bn_lst
 
integer, dimension(:), pointer bn_loc
 
integer, dimension(:), pointer bn_mlt
 
integer, dimension(:,:), pointer bn_ney
 
integer, dimension(:), pointer nde_id
 
integer nhn
 
integer, dimension(:), pointer hn_lst
 
integer, dimension(:), pointer hn_own
 

Function/Subroutine Documentation

◆ aprint_arr()

subroutine mod_par::aprint_arr ( integer, intent(in)  IUNIT,
real(sp), dimension(:,:), intent(in), allocatable, target  VARP,
character(len=*), intent(in)  VART,
type(time), intent(in)  NOW,
integer, intent(in)  ILOC,
integer, intent(in)  K1,
integer, intent(in), optional  K2,
character(len=*), intent(in), optional  MSG 
)

Definition at line 288 of file mod_par.f90.

288  USE mod_time
289  IMPLICIT NONE
290  TYPE(TIME), INTENT(IN) :: NOW
291  INTEGER, INTENT(IN) :: IUNIT,ILOC,K1
292  INTEGER, INTENT(IN),OPTIONAL :: K2
293  REAL(SP), ALLOCATABLE, INTENT(IN),TARGET :: VARP(:,:)
294  CHARACTER(LEN=*), INTENT(IN) :: VART
295  CHARACTER(LEN=*), INTENT(IN),OPTIONAL :: MSG
296  REAL(SP), POINTER:: VARP_O(:,:)
297 
298  varp_o => varp
299 
300  IF (PRESENT(k2)) THEN
301  IF(PRESENT(msg))THEN
302  CALL pprint_arr(iunit,varp_o,vart,now,iloc,k1,k2,msg)
303  ELSE
304  CALL pprint_arr(iunit,varp_o,vart,now,iloc,k1,k2)
305  END IF
306  ELSE
307  IF(PRESENT(msg))THEN
308  CALL pprint_arr(iunit,varp_o,vart,now,iloc,k1,k1,msg)
309  ELSE
310  CALL pprint_arr(iunit,varp_o,vart,now,iloc,k1)
311  END IF
312  END IF
mod_time
Definition: mod_time.f90:40
Here is the call graph for this function:

◆ aprint_vec()

subroutine mod_par::aprint_vec ( integer, intent(in)  IUNIT,
real(sp), dimension(:), intent(in), allocatable, target  VARP,
character(len=*), intent(in)  VART,
type(time), intent(in)  NOW,
integer, intent(in)  ILOC,
character(len=*), intent(in), optional  MSG 
)

Definition at line 145 of file mod_par.f90.

145  USE mod_time
146  IMPLICIT NONE
147  TYPE(TIME), INTENT(IN) :: NOW
148  INTEGER, INTENT(IN) :: IUNIT,ILOC
149  REAL(SP), ALLOCATABLE, INTENT(IN),TARGET :: VARP(:)
150  CHARACTER(LEN=*), INTENT(IN) :: VART
151  CHARACTER(LEN=*), INTENT(IN),OPTIONAL :: MSG
152  REAL(SP), POINTER:: VARP_O(:)
153 
154  varp_o => varp
155 
156  IF(PRESENT(msg))THEN
157  CALL pprint_vec(iunit,varp_o,vart,now,iloc,msg)
158  ELSE
159  CALL pprint_vec(iunit,varp_o,vart,now,iloc)
160  END IF
161 
mod_time
Definition: mod_time.f90:40
Here is the call graph for this function:

◆ pprint_arr()

subroutine mod_par::pprint_arr ( integer, intent(in)  IUNIT,
real(sp), dimension(:,:), intent(in), pointer  VARP,
character(len=*), intent(in)  VART,
type(time), intent(in)  NOW,
integer, intent(in)  ILOC,
integer, intent(in)  K1,
integer, intent(in), optional  K2,
character(len=*), intent(in), optional  MSG 
)

Definition at line 317 of file mod_par.f90.

317  USE control
318  USE lims
319 
320  IMPLICIT NONE
321  TYPE(TIME), INTENT(IN) :: NOW
322  INTEGER, INTENT(IN) :: IUNIT,ILOC,K1
323  INTEGER, INTENT(IN),OPTIONAL :: K2
324  REAL(SP), POINTER, INTENT(IN) :: VARP(:,:)
325  CHARACTER(LEN=*), INTENT(IN) :: VART
326  CHARACTER(LEN=*), INTENT(IN),OPTIONAL :: MSG
327 
328  CHARACTER(LEN=80),parameter :: VAR_E = "element"
329  CHARACTER(LEN=80),parameter :: VAR_N = "node"
330 
331  CHARACTER(LEN=100) :: STRNG
332  CHARACTER(LEN=20) :: short
333 
334  INTEGER :: I,J,K,PROCMAX,II,IBND,Kopt,IERR
335  LOGICAL :: PRINT_PROC
336 
337  !==============================================================================|
338 
339  !------------------------------------------------------------------------------|
340  ! Process Iunit for Errors |
341  !------------------------------------------------------------------------------|
342  IF(iunit /= ipt .and. (iunit < 300 .or. iunit > 7000) )THEN
343  CALL fatal_error('ERROR IN PPRINT',&
344  & 'FILE UNIT < 300 AND UNIT > 7000 ARE RESERVED FOR FVCOM I/O',&
345  & 'PLEASE INCREASE IUNIT TO 300+')
346  END IF
347 
348  !------------------------------------------------------------------------------|
349  ! Process Vartype for Errors |
350  !------------------------------------------------------------------------------|
351  IF(vart /= var_e .AND. vart /= var_n)THEN
352  CALL fatal_error('VART IN PPRINT NOT CORRECT :'//trim(vart),&
353  & 'SHOULD BE "'//trim(var_e)//'" or "'//trim(var_n)//'"')
354  END IF
355 
356  !------------------------------------------------------------------------------|
357  ! Process optional sigma level range |
358  !------------------------------------------------------------------------------|
359 
360  IF(PRESENT(k2)) THEN
361  kopt=k2
362  ELSE
363  kopt=k1
364  END IF
365 
366  !------------------------------------------------------------------------------|
367  ! Process string output |
368  !------------------------------------------------------------------------------|
369  IF(PRESENT(msg)) strng=trim(msg)//"; IINT"
370 
371 
372  IF(abs(iint) .lt. 1000) THEN
373  WRITE(short,'(I5)') iint
374  ELSE IF(abs(iint) .lt. 1000000) THEN
375  WRITE(short,'(I8)') iint
376  ELSE
377  WRITE(short,*) iint
378  END IF
379 
380  IF(use_real_world_time) THEN
381  strng = trim(strng)//trim(short)//", Date/Time:"&
382  &//trim(write_datetime(now,3,timezone))//"; ILOC "
383  ELSE
384  strng = trim(strng)//trim(short)//", Time(s):"
385  WRITE(short,'(f16.8)') seconds(now)
386  strng = trim(strng)//trim(short)//"; ILOC="
387  END IF
388 
389  WRITE(short,'(I8)') iloc
390  strng = trim(strng)//trim(short)//"; VALUES="
391  !------------------------------------------------------------------------------|
392  ! Single Processor Case |
393  !------------------------------------------------------------------------------|
394  IF(nprocs == 1)THEN
395  WRITE(iunit,*) trim(strng)
396  WRITE(iunit,*) (varp(iloc,k),k=k1,kopt)
397  END IF
398 
399  !------------------------------------------------------------------------------|
400  ! Multi Processor Case with Element Based Variable (u,v,ww, etc) |
401  ! Transform to Local Element ID with "ELID" |
402  !------------------------------------------------------------------------------|
403 
404  IF(nprocs /= 1 .AND. vart == var_e .AND. elid(iloc) /= 0)THEN
405 
406  WRITE(iunit,*) trim(strng)
407  WRITE(iunit,*) (varp(elid(iloc),k),k=k1,kopt)
408 
409  END IF
410 
411  !------------------------------------------------------------------------------|
412  ! Multi Processor Case with Node Based Variable (s1,t1,rho1,e1, etc) |
413  ! Transform to Local Node ID with "NLID" |
414  ! If Node is Interprocessor Boundary Node, Choose Processor with Highest |
415  ! ID Number to Write Values to File |
416  !------------------------------------------------------------------------------|
417 
418  IF(nprocs /= 1 .AND. vart == var_n .AND. nlid(iloc) > 0)THEN
419 
420  print_proc = .true.
421  IF(nde_id(nlid(iloc)) == 1)THEN !!BOUNDARY NODE
422 
423  DO ii=1,nbn
424  IF(bn_lst(ii) == iloc) ibnd = ii
425  END DO
426 
427  procmax = 10000
428  DO j=1,nprocs
429  IF(bn_ney(ibnd,j)==1) THEN
430  IF(j < procmax) procmax = j
431  END IF
432  END DO
433 
434  IF(procmax /= myid) print_proc = .false. !!NOT RESPONSIBLE FOR OUTPUT
435  END IF
436 
437  IF(print_proc)THEN
438  WRITE(iunit,*) trim(strng)
439  WRITE(iunit,*) (varp(nlid(iloc),k),k=k1,kopt)
440  END IF
441 
442  END IF
443 
444 
445  RETURN
mod_time::abs
Definition: mod_time.f90:64
lims::myid
integer myid
Definition: mod_main.f90:67
lims::nprocs
integer, target nprocs
Definition: mod_main.f90:72
control::iint
integer(itime) iint
Definition: mod_main.f90:850
control::timezone
character(len=80) timezone
Definition: mod_main.f90:126
control::use_real_world_time
logical use_real_world_time
Definition: mod_main.f90:131
mod_time::seconds
real(dp) function seconds(MJD)
Definition: mod_time.f90:742
mod_utils::fatal_error
subroutine fatal_error(ER1, ER2, ER3, ER4)
Definition: mod_utils.f90:230
control::ipt
integer ipt
Definition: mod_main.f90:922
mod_time::write_datetime
character(len=80) function write_datetime(mjdin, prec, TZONE)
Definition: mod_time.f90:682
Here is the call graph for this function:
Here is the caller graph for this function:

◆ pprint_vec()

subroutine mod_par::pprint_vec ( integer, intent(in)  IUNIT,
real(sp), dimension(:), intent(in), pointer  VARP,
character(len=*), intent(in)  VART,
type(time), intent(in)  NOW,
integer, intent(in)  ILOC,
character(len=*), intent(in), optional  MSG 
)

Definition at line 166 of file mod_par.f90.

166  USE control
167  USE lims
168 
169  IMPLICIT NONE
170  TYPE(TIME), INTENT(IN) :: NOW
171  INTEGER, INTENT(IN) :: IUNIT,ILOC
172  REAL(SP), POINTER, INTENT(IN) :: VARP(:)
173  CHARACTER(LEN=*), INTENT(IN) :: VART
174  CHARACTER(LEN=*), INTENT(IN),OPTIONAL :: MSG
175 
176  CHARACTER(LEN=80),parameter :: VAR_E = "element"
177  CHARACTER(LEN=80),parameter :: VAR_N = "node"
178 
179  CHARACTER(LEN=100) :: STRNG
180  CHARACTER(LEN=20) :: short
181 
182  INTEGER :: I,J,K,PROCMAX,II,IBND,Kopt,IERR
183  LOGICAL :: PRINT_PROC
184 
185  !==============================================================================|
186 
187  !------------------------------------------------------------------------------|
188  ! Process Iunit for Errors |
189  !------------------------------------------------------------------------------|
190  IF(iunit /= ipt .and. (iunit < 300 .or. iunit > 7000) )THEN
191  CALL fatal_error('ERROR IN PPRINT',&
192  & 'FILE UNIT < 300 AND UNIT > 7000 ARE RESERVED FOR FVCOM I/O',&
193  & 'PLEASE INCREASE IUNIT TO 300+')
194  END IF
195 
196  !------------------------------------------------------------------------------|
197  ! Process Vartype for Errors |
198  !------------------------------------------------------------------------------|
199  IF(vart /= var_e .AND. vart /= var_n)THEN
200  CALL fatal_error('VART IN PPRINT NOT CORRECT :'//trim(vart),&
201  & 'SHOULD BE "'//trim(var_e)//'" or "'//trim(var_n)//'"')
202  END IF
203 
204 
205  !------------------------------------------------------------------------------|
206  ! Process string output |
207  !------------------------------------------------------------------------------|
208  IF(PRESENT(msg)) strng=trim(msg)//"; IINT"
209 
210 
211  IF(abs(iint) .lt. 1000) THEN
212  WRITE(short,'(I5)') iint
213  ELSE IF(abs(iint) .lt. 1000000) THEN
214  WRITE(short,'(I8)') iint
215  ELSE
216  WRITE(short,*) iint
217  END IF
218 
219  IF(use_real_world_time) THEN
220  strng = trim(strng)//trim(short)//", Date/Time:"&
221  &//trim(write_datetime(now,3,timezone))//"; ILOC= "
222  ELSE
223  strng = trim(strng)//trim(short)//", Time(s):"
224  WRITE(short,'(f16.8)') seconds(now)
225  strng = trim(strng)//trim(short)//"; ILOC="
226  END IF
227 
228  WRITE(short,'(I8)') iloc
229  strng = trim(strng)//trim(short)//"; VALUE="
230 
231  !------------------------------------------------------------------------------|
232  ! Single Processor Case |
233  !------------------------------------------------------------------------------|
234  IF(nprocs == 1)THEN
235  WRITE(iunit,*) trim(strng),varp(iloc)
236  END IF
237 
238  !------------------------------------------------------------------------------|
239  ! Multi Processor Case with Element Based Variable (u,v,ww, etc) |
240  ! Transform to Local Element ID with "ELID" |
241  !------------------------------------------------------------------------------|
242 
243  IF(nprocs /= 1 .AND. vart == var_e .AND. elid(iloc) /= 0)THEN
244 
245  WRITE(iunit,*) trim(strng),varp(elid(iloc))
246 
247  END IF
248 
249  !------------------------------------------------------------------------------|
250  ! Multi Processor Case with Node Based Variable (s1,t1,rho1,e1, etc) |
251  ! Transform to Local Node ID with "NLID" |
252  ! If Node is Interprocessor Boundary Node, Choose Processor with Highest |
253  ! ID Number to Write Values to File |
254  !------------------------------------------------------------------------------|
255 
256  IF(nprocs /= 1 .AND. vart == var_n .AND. nlid(iloc) > 0)THEN
257 
258  print_proc = .true.
259  IF(nde_id(nlid(iloc)) == 1)THEN !!BOUNDARY NODE
260 
261  DO ii=1,nbn
262  IF(bn_lst(ii) == iloc) ibnd = ii
263  END DO
264 
265  procmax = 10000
266  DO j=1,nprocs
267  IF(bn_ney(ibnd,j)==1) THEN
268  IF(j < procmax) procmax = j
269  END IF
270  END DO
271 
272  IF(procmax /= myid) print_proc = .false. !!NOT RESPONSIBLE FOR OUTPUT
273  END IF
274 
275  IF(print_proc)THEN
276  WRITE(iunit,*) trim(strng),varp(nlid(iloc))
277  END IF
278 
279  END IF
280 
281 
282 
283  RETURN
mod_time::abs
Definition: mod_time.f90:64
lims::myid
integer myid
Definition: mod_main.f90:67
lims::nprocs
integer, target nprocs
Definition: mod_main.f90:72
control::iint
integer(itime) iint
Definition: mod_main.f90:850
control::timezone
character(len=80) timezone
Definition: mod_main.f90:126
control::use_real_world_time
logical use_real_world_time
Definition: mod_main.f90:131
mod_time::seconds
real(dp) function seconds(MJD)
Definition: mod_time.f90:742
mod_utils::fatal_error
subroutine fatal_error(ER1, ER2, ER3, ER4)
Definition: mod_utils.f90:230
control::ipt
integer ipt
Definition: mod_main.f90:922
mod_time::write_datetime
character(len=80) function write_datetime(mjdin, prec, TZONE)
Definition: mod_time.f90:682
Here is the call graph for this function:
Here is the caller graph for this function:

Variable Documentation

◆ bn_loc

integer, dimension(:), pointer mod_par::bn_loc

Definition at line 78 of file mod_par.f90.

78  INTEGER, POINTER :: BN_LOC(:) !!LOCAL IDENTITY OF BOUNDARY NODES

◆ bn_lst

integer, dimension(:), pointer mod_par::bn_lst

Definition at line 77 of file mod_par.f90.

77  INTEGER, POINTER :: BN_LST(:) !!GLOBAL IDENTITY OF BOUNDARY NODES

◆ bn_mlt

integer, dimension(:), pointer mod_par::bn_mlt

Definition at line 79 of file mod_par.f90.

79  INTEGER, POINTER :: BN_MLT(:) !!MULTIPLICITY OF BOUNDARY NODES

◆ bn_ney

integer, dimension(:,:), pointer mod_par::bn_ney

Definition at line 80 of file mod_par.f90.

80  INTEGER, POINTER :: BN_NEY(:,:) !!NODE OWNER LIST

◆ egid

integer, dimension(:), pointer mod_par::egid

Definition at line 60 of file mod_par.f90.

60  INTEGER, POINTER :: EGID(:) !!GLOBAL ID OF LOCAL ELEMENT

◆ egid_x

integer, dimension(:), pointer mod_par::egid_x

Definition at line 62 of file mod_par.f90.

62  INTEGER, POINTER :: EGID_X(:) !!GLOBAL ID OF LOCAL ELEMENT

◆ el_pid

integer, dimension(:), pointer mod_par::el_pid

Definition at line 51 of file mod_par.f90.

51  INTEGER, POINTER :: EL_PID(:) !!PROCESSOR OWNER OF GLOBAL ELEMENT

◆ elid

integer, dimension(:), pointer mod_par::elid

Definition at line 53 of file mod_par.f90.

53  INTEGER, POINTER :: ELID(:) !!LOCAL VALUE OF GLOBAL ELEMENT

◆ elid_x

integer, dimension(:), pointer mod_par::elid_x

Definition at line 55 of file mod_par.f90.

55  INTEGER, POINTER :: ELID_X(:) !!LOCAL VALUE OF GLOBAL ELEMENT INCLUDING HALOS

◆ he_lst

integer, dimension(:), pointer mod_par::he_lst

Definition at line 69 of file mod_par.f90.

69  INTEGER, POINTER :: HE_LST(:) !!GLOBAL IDENTITIES OF HALO ELEMENTS

◆ he_own

integer, dimension(:), pointer mod_par::he_own

Definition at line 70 of file mod_par.f90.

70  INTEGER, POINTER :: HE_OWN(:) !!OWNER OF HALO ELEMENTS

◆ hn_lst

integer, dimension(:), pointer mod_par::hn_lst

Definition at line 87 of file mod_par.f90.

87  INTEGER, POINTER :: HN_LST(:) !!LIST OF HALO NODES

◆ hn_own

integer, dimension(:), pointer mod_par::hn_own

Definition at line 88 of file mod_par.f90.

88  INTEGER, POINTER :: HN_OWN(:) !!PRIMARY OWNER OF HALO NODES

◆ mx_mlt

integer mod_par::mx_mlt

Definition at line 76 of file mod_par.f90.

76  INTEGER :: MX_MLT !!MAX MULTIPLICITY OF BOUNDARY NODES

◆ nbn

integer mod_par::nbn

Definition at line 75 of file mod_par.f90.

75  INTEGER :: NBN !!NUMBER OF BOUNDARY NODES

◆ nde_id

integer, dimension(:), pointer mod_par::nde_id

Definition at line 81 of file mod_par.f90.

81  INTEGER, POINTER :: NDE_ID(:) !! = 0 IF INTERNAL, 1 IF ON INTERNAL BOUNDARY

◆ ngid

integer, dimension(:), pointer mod_par::ngid

Definition at line 61 of file mod_par.f90.

61  INTEGER, POINTER :: NGID(:) !!GLOBAL ID OF LOCAL NODE

◆ ngid_x

integer, dimension(:), pointer mod_par::ngid_x

Definition at line 63 of file mod_par.f90.

63  INTEGER, POINTER :: NGID_X(:) !!GLOBAL ID OF LOCAL NODE

◆ nhe

integer mod_par::nhe

Definition at line 68 of file mod_par.f90.

68  INTEGER :: NHE !!NUMBER OF HALO ELEMENTS

◆ nhn

integer mod_par::nhn

Definition at line 86 of file mod_par.f90.

86  INTEGER :: NHN !!NUMBER OF HALO NODES

◆ nlid

integer, dimension(:), pointer mod_par::nlid

Definition at line 54 of file mod_par.f90.

54  INTEGER, POINTER :: NLID(:) !!LOCAL VALUE OF GLOBAL NODE

◆ nlid_x

integer, dimension(:), pointer mod_par::nlid_x

Definition at line 56 of file mod_par.f90.

56  INTEGER, POINTER :: NLID_X(:) !!LOCAL VALUE OF GLOBAL NODE INCLUDING HALOS
Generated by   doxygen 1.8.15