probes::myprobe Interface Reference

Public Member Functions
subroutine	myprobe_vec (PROBE, VEC)
subroutine	myprobe_arr (PROBE, ARR)

Detailed Description

Definition at line 114 of file mod_probe.f90.

Member Function/Subroutine Documentation

myprobe_arr()

subroutine probes::myprobe::myprobe_arr	(	type(probe_obj)	PROBE,
		real(sp), dimension(:,:), allocatable, target	ARR
	)

Definition at line 359 of file mod_probe.f90.

     USE all_vars, only:xmc,ymc,h1,h,xm,ym,h,lon,lat,lonc,latc
     IMPLICIT NONE
     TYPE(PROBE_OBJ) :: PROBE
     REAL(SP), ALLOCATABLE, TARGET :: ARR(:,:)
     CHARACTER(LEN=80):: cstr1,cstr2,cstr3
     INTEGER :: I,J,IBND,PROCMAX
     !============ENSURE SIGMA RANGE IS TENABLE=============================

     if(dbg_set(dbg_sbr)) write(ipt,*) "START: MYPROBE_ARR"

     IF(probe%K_ONE > ubound(arr,2) .or. probe%K_TWO > ubound(arr,2) )THEN
        CALL fatal_error('ERROR IN PROBE SETUP: PROBE LEVEL RANGE NOT CORRECT FOR VARIABLE: '//trim(probe_variable),&
             & 'MAKE SURE PROBE LEVELS ARE LESS THAN OR EQUAL TO THE NUMBER OF MODEL LEVELS') 
     END IF

     IF(probe%K_ONE < 1 .or. probe%K_TWO <1)THEN
        CALL fatal_error('ERROR IN PROBE SETUP: PROBE LEVEL RANGE NOT CORRECT FOR VARIABLE: '//trim(probe_variable),&
             & 'MAKE SURE PROBE LEVELS ARE GREATER THAN OR EQUAL TO ONE') 
     END IF
     
     
     IF(probe_levels(1) >  probe_levels(2) )THEN
        CALL fatal_error&
             &('ERROR IN PROBE SETUP: PROBE LEVEL RANGE NOT CORRECT FOR VARIABLE: '//trim(probe_variable),&
             & 'THE PROBE LEVEL INTERVAL MUST SPECIFY A VALID RANGE a:b') 
     END IF
     
     
     IF(probe%D_LOC_GL  <1)THEN 
        write(cstr1,'(i8)') probe%D_LOC_GL
        write(cstr2,'(i8)') ngl
        
        CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
             & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
     END IF
     
     IF( ubound(arr,1) == nt)THEN

        IF(probe%D_LOC_GL > ngl)THEN 
           
           write(cstr1,'(i8)') probe%D_LOC_GL
           write(cstr2,'(i8)') ngl
           
           CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
                & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
        END IF

        IF(elid(probe%D_LOC_GL) /= 0) THEN
           
           probe%D_LOC = elid(probe%D_LOC_GL)
           probe%MINE=.true.
           
           ! METERS
           probe%XLOC = xmc(probe%D_LOC)
           probe%YLOC = ymc(probe%D_LOC)
           ! SPHERICAL
           probe%LONLOC = lonc(probe%D_LOC)
           probe%LATLOC = latc(probe%D_LOC)
           ! DEPTH 
           probe%DPTH = h1(probe%D_LOC) 
           
           probe%VEC => arr(probe%D_LOC,probe%K_ONE:probe%K_TWO)
        END IF


     ELSE IF ( ubound(arr,1) == mt)THEN

        IF(probe%D_LOC_GL > mgl)THEN 
           write(cstr1,'(i8)') probe%D_LOC_GL
           write(cstr2,'(i8)') mgl
           
           CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
                & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
           
        END IF

        IF(nlid(probe%D_LOC_GL) == 0) RETURN

        IF(nlid(probe%D_LOC_GL) > 0) THEN
           
           IF(nde_id(nlid(probe%D_LOC_GL)) == 1)THEN   !!BOUNDARY NODE
              DO j=1,nbn
                 IF(bn_lst(j) == probe%D_LOC_GL) ibnd = j
              END DO
              !----Choose Processor of Lowest ID to be responsible for node
              procmax = 10000
              DO j=1,nprocs
                 IF(bn_ney(ibnd,j)==1) THEN
                    IF(j < procmax) procmax = j
                 END IF
              END DO
              
              IF(procmax == myid) THEN
                 probe%MINE = .true. !!NOT RESPONSIBLE FOR TIME SERIES
                 probe%D_LOC = nlid(probe%D_LOC_GL)

                 ! METERS
                 probe%XLOC = xm(probe%D_LOC)
                 probe%YLOC = ym(probe%D_LOC)
                 ! SPHERICAL
                 probe%LONLOC = lon(probe%D_LOC)
                 probe%LATLOC = lat(probe%D_LOC)
                 ! DEPTH 
                 probe%DPTH = h(probe%D_LOC) 

                 probe%VEC => arr(probe%D_LOC,probe%K_ONE:probe%K_TWO)
              END IF
           ELSE
             probe%MINE = .true. !!NOT RESPONSIBLE FOR TIME SERIES
             probe%D_LOC = nlid(probe%D_LOC_GL)

             ! METERS
             probe%XLOC = xm(probe%D_LOC)
             probe%YLOC = ym(probe%D_LOC)
             ! SPHERICAL
             probe%LONLOC = lon(probe%D_LOC)
             probe%LATLOC = lat(probe%D_LOC)
             ! DEPTH 
             probe%DPTH = h(probe%D_LOC) 

             probe%VEC => arr(probe%D_LOC,probe%K_ONE:probe%K_TWO)
       
           END IF
        END IF
        

     ELSE
        CALL fatal_error('MYPROBE: INVALID VARIABLE SIZE?'&
             &,'Variable:'//trim(probe%VAR))
     END IF

     if(dbg_set(dbg_sbr)) write(ipt,*) "END: MYPROBE_ARR"

myprobe_vec()

subroutine probes::myprobe::myprobe_vec	(	type(probe_obj)	PROBE,
		real(sp), dimension(:), allocatable, target	VEC
	)

Definition at line 239 of file mod_probe.f90.

     USE all_vars, only:xmc,ymc,h1,h,xm,ym,h,lon,lat,lonc,latc
     IMPLICIT NONE
     TYPE(PROBE_OBJ) :: PROBE
     REAL(SP), ALLOCATABLE, TARGET :: VEC(:)
     CHARACTER(LEN=80):: cstr1,cstr2,cstr3
     INTEGER :: I,J,IBND,PROCMAX
     !============SEE IF DATA POINT IS IN THE GLOBAL DOMAIN=================
     
     if(dbg_set(dbg_sbr)) write(ipt,*) "START: MYPROBE_VEC"
     
     IF (probe%K_ONE /= -1 .or. probe%K_TWO /= -1) CALL fatal_error&
          &("ERROR IN PROBE SETUP: PROBE_LEVELS should not be set for vector variables",&
          & "Do not specify it in the PROBE Namelist object for "//trim(probe_variable))

     IF(probe%D_LOC_GL  <1)THEN 
        write(cstr1,'(i8)') probe%D_LOC_GL
        write(cstr2,'(i8)') ngl
        
        CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
             & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
     END IF
     
     IF( ubound(vec,1) == nt)THEN

        IF(probe%D_LOC_GL > ngl)THEN 
           
           write(cstr1,'(i8)') probe%D_LOC_GL
           write(cstr2,'(i8)') ngl
           
           CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
                & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
        END IF

        IF(elid(probe%D_LOC_GL) /= 0) THEN
           
           probe%D_LOC = elid(probe%D_LOC_GL)
           probe%MINE=.true.
           
           ! METERS
           probe%XLOC = xmc(probe%D_LOC)
           probe%YLOC = ymc(probe%D_LOC)
           ! SPHERICAL
           probe%LONLOC = lonc(probe%D_LOC)
           probe%LATLOC = latc(probe%D_LOC)
           ! DEPTH 
           probe%DPTH = h1(probe%D_LOC) 
           
           probe%SCL => vec(probe%D_LOC)
        END IF
        
     ELSE IF ( ubound(vec,1) == mt)THEN
        
        IF(probe%D_LOC_GL > mgl)THEN 
           write(cstr1,'(i8)') probe%D_LOC_GL
           write(cstr2,'(i8)') mgl
           
           CALL fatal_error('ERROR IN PROBE SETUP: DATA LOCATION'//trim(cstr1)//' FOR TIME SERIES FILE: '//trim(probe_variable),&
                & 'IS NOT IN GLOBAL DOMAIN: 1 --> '//trim(cstr2))
           
        END IF

        IF(nlid(probe%D_LOC_GL) == 0) RETURN

        IF(nlid(probe%D_LOC_GL) > 0) THEN
           
           IF(nde_id(nlid(probe%D_LOC_GL)) == 1)THEN   !!BOUNDARY NODE
              DO j=1,nbn
                 IF(bn_lst(j) == probe%D_LOC_GL) ibnd = j
              END DO
              !----Choose Processor of Lowest ID to be responsible for node
              procmax = 10000
              DO j=1,nprocs
                 IF(bn_ney(ibnd,j)==1) THEN
                    IF(j < procmax) procmax = j
                 END IF
              END DO
              
              IF(procmax == myid) THEN
                 probe%MINE = .true. !!NOT RESPONSIBLE FOR TIME SERIES
                 probe%D_LOC = nlid(probe%D_LOC_GL)
                 
                 ! METERS
                 probe%XLOC = xm(probe%D_LOC)
                 probe%YLOC = ym(probe%D_LOC)
                 ! SPHERICAL
                 probe%LONLOC = lon(probe%D_LOC)
                 probe%LATLOC = lat(probe%D_LOC)
                 ! DEPTH 
                 probe%DPTH = h(probe%D_LOC) 

                 probe%SCL => vec(probe%D_LOC)
              END IF
            ELSE
              probe%MINE = .true. !!NOT RESPONSIBLE FOR TIME SERIES
              probe%D_LOC = nlid(probe%D_LOC_GL)
                 
              ! METERS
              probe%XLOC = xm(probe%D_LOC)
              probe%YLOC = ym(probe%D_LOC)
              ! SPHERICAL
              probe%LONLOC = lon(probe%D_LOC)
              probe%LATLOC = lat(probe%D_LOC)
              ! DEPTH 
              probe%DPTH = h(probe%D_LOC) 

              probe%SCL => vec(probe%D_LOC)
        
           END IF
        END IF
        

     ELSE
        CALL fatal_error('MYPROBE: INVALID VARIABLE SIZE (Not equal MT or NT?)'&
             &,'Variable:'//trim(probe%VAR))
     END IF
     if(dbg_set(dbg_sbr)) write(ipt,*) "END: MYPROBE_VEC"

---

The documentation for this interface was generated from the following file:

• FVCOMg01/mod_probe.f90