!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE A3DMAX0 ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE a3dmax0(a,m1,m2,i1,i2,n1,n2,j1,j2,l1,l2,k1,k2, & 322,1
amax,amin)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Find the maximum and minimum of a 3-D array, a, in a specified
! subdomain.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 10/01/91.
!
! MODIFICATION HISTORY:
!
! 5/06/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! a A 3-D array whose max. and min. values are sought here.
!
! m1,m2 i-index of array, a.
! i1,i2 i-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! n1,n2 j-index of array, a.
! j1,j2 j-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! l1,l2 k-index of array, a.
! k1,k2 k-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
!
! OUTPUT:
!
! amax The maximum value of an array, a, in a specified
! subdomain.
! amin The minimum value of an array, a, in a specified
! subdomain.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: m1,n1,l1,m2,n2,l2
REAL :: a(m1:m2,n1:n2,l1:l2)
REAL :: amin, amax
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i1,i2,j1,j2,k1,k2,i,j,k
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'globcst.inc'
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
amax=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
amax = MAX(amax, a(i,j,k))
END DO
END DO
END DO
amin=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
amin = MIN(amin, a(i,j,k))
END DO
END DO
END DO
IF (mp_opt > 0) CALL mpmax0
(amax,amin)
RETURN
END SUBROUTINE a3dmax0
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE A3DMAX0LCL ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE a3dmax0lcl(a,m1,m2,i1,i2,n1,n2,j1,j2,l1,l2,k1,k2, & 57
amax,amin)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Find the maximum and minimum of a 3-D array, a, in a specified
! subdomain but only on the local processor when running message
! passing version.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 10/01/91.
!
! MODIFICATION HISTORY:
!
! 5/06/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! a A 3-D array whose max. and min. values are sought here.
!
! m1,m2 i-index of array, a.
! i1,i2 i-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! n1,n2 j-index of array, a.
! j1,j2 j-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! l1,l2 k-index of array, a.
! k1,k2 k-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
!
! OUTPUT:
!
! amax The maximum value of an array, a, in a specified
! subdomain.
! amin The minimum value of an array, a, in a specified
! subdomain.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: m1,n1,l1,m2,n2,l2
REAL :: a(m1:m2,n1:n2,l1:l2)
REAL :: amin, amax
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i1,i2,j1,j2,k1,k2,i,j,k
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
amax=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
amax = MAX(amax, a(i,j,k))
END DO
END DO
END DO
amin=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
amin = MIN(amin, a(i,j,k))
END DO
END DO
END DO
RETURN
END SUBROUTINE a3dmax0lcl
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE A3DMAX ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE a3dmax(a,m1,m2,i1,i2,n1,n2,j1,j2,l1,l2,k1,k2, & 42,1
amax,amin,imax,jmax,kmax, imin,jmin,kmin)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Find the maximum, minimum and the index locations of an array, a,
! in a specified subdomain.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 10/01/91.
!
! MODIFICATION HISTORY:
!
! 5/06/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! a A 3-D array whose max. and min. values are sought here.
!
! m1,m2 i-index of array, a.
! i1,i2 i-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! n1,n2 j-index of array, a.
! j1,j2 j-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! l1,l2 k-index of array, a.
! k1,k2 k-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
!
! OUTPUT:
!
! amax The maximum value of an array, a, in a specified
! subdomain.
! amin The minimum value of an array, a, in a specified
! subdomain.
! imax,jmax,kmax The index location of the maximum value (i,j,k)
! imin,jmin,kmin The index location of the minimum value (i,j,k)
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: m1,n1,l1,m2,n2,l2
REAL :: a(m1:m2,n1:n2,l1:l2)
REAL :: amin, amax
INTEGER :: imax,jmax,kmax, imin,jmin,kmin
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i1,i2,j1,j2,k1,k2,i,j,k
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'globcst.inc'
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
imax = i1
jmax = j1
kmax = k1
amax=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
IF( a(i,j,k) > amax ) THEN
amax = a(i,j,k)
imax = i
jmax = j
kmax = k
END IF
END DO
END DO
END DO
imin = i1
jmin = j1
kmin = k1
amin=a(i1,j1,k1)
DO k=k1,k2
DO j=j1,j2
DO i=i1,i2
IF( a(i,j,k) < amin ) THEN
amin = a(i,j,k)
imin = i
jmin = j
kmin = k
END IF
END DO
END DO
END DO
IF (mp_opt > 0) CALL mpmax(amax,amin,m2,n2,l2,imax,jmax,kmax,imin,jmin,kmin)
RETURN
END SUBROUTINE a3dmax
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE WRIGAR ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE wrigar(a,nx0,nx1,ny0,ny1,nz0,nz1,i0,i1,j0,j1,k0,k1, & 101,2
title,tkoff,mode)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Orchestrate the formatted printing of select 2-d slices from 3-D
! arrays. A subdomain (i0:i1,j0:j1,k0:k1) is selected for printing.
! A maximum width of 32 numbers can be printed in a horizontal
! direction. A portion of the array table may be truncated. In the
! case when the arrays are large, the data are printed for every
! other data point.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 12/10/89.
!
! MODIFICATION HISTORY:
!
! 5/20/92 (M. Xue)
! Added full documentation.
!
! 9/20/93 (A. Sathye)
! Fixed problem with late nx,ny.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT :
!
! a A 3-D array whose content will be printed.
!
! nx0,nx1 i-index of array a.
! ny0,ny1 j-index of array a.
! nz0,nz1 k-index of array a.
!
! i0,i1 i-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
! j0,j1 j-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
! k0,k1 k-index defining a subdomain of array, a, where
! the max. and min. of the array is sought
!
! title The name of the field (character string).
! tkoff A value to be subtracted from array, a (for printing).
!
! Mode A print control option for the selection of the slice
! orientation.
! = 0, all slices will be printed.
! = 1, x-y slices will be printed.
! = 2, x-z slices will be printed.
! = 3, y-z slices will be printed.
!
! OUTPUT:
!
! None.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: nx0,nx1,ny0,ny1,nz0,nz1
INTEGER :: i0,i1,j0,j1,k0,k1
REAL :: a(nx0:nx1,ny0:ny1,nz0:nz1)
CHARACTER (LEN=*) :: title
REAL :: tkoff
INTEGER :: mode
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i,j,k,modd,nap,nexp
REAL :: armin,armax,scale,factor
INTEGER :: nplan, mstep
PARAMETER(nplan=1000,mstep=10)
INTEGER :: iap(nplan)
INTEGER :: is,js,ks,ndigit,nunit,leng,nxmi,nymi,nzmi
PARAMETER(ndigit=2 ,nunit=6,leng=132)
PARAMETER(nxmi=33,nymi=5,nzmi=20)
REAL :: eps
PARAMETER(eps=1.0E-30)
INTEGER :: ixp(9),iyp(9),izp(9)
SAVE ixp,iyp,izp
DATA ixp /nxmi,-1,-1,-1,-1,-1,-1,-1,-1/
DATA iyp /-1,-1,-1,-1,-1,-1,-1,-1,-1/
DATA izp /-1,-1,-1,-1,-1,-1,-1,-1,-1/
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
DO i = 1, nplan
iap(i) = i
END DO
!
!-----------------------------------------------------------------------
!
! Calculate the max and min for scaling:
!
!-----------------------------------------------------------------------
!
!
armin=a(i0,j0,k0)
armax=armin
DO k=k0,k1
DO j=j0,j1
DO i=i0,i1
armin=MIN( armax=MAX( END DO
END DO
END DO
!
!-----------------------------------------------------------------------
!
! Determine the scaling factor:
!
!-----------------------------------------------------------------------
!
!
WRITE(nunit,1001) title,armin,armax
IF(armax == armin .OR. mode == 4) RETURN
armax=armax-tkoff
armin=armin-tkoff
scale=1.0E-30+MAX(ABS(armax),ABS(armin))
nexp=INT(ALOG10(scale))-ndigit
IF(nint(scale*10.**(-nexp)) == 10**(-nexp+1)) nexp=nexp+1
IF(ALOG10(scale) <= 0.) nexp=nexp-1
factor=10.**(-nexp)
IF(k0 == k1) THEN
modd=1
ELSE IF(j0 == j1) THEN
modd=2
ELSE IF(i0 == i1) THEN
modd=3
ELSE
modd=mode
END IF
is = MAX(1, nint( (i1-i0+1)/32.0 ))
js = MAX(1, nint( (j1-j0+1)/32.0 ))
ks = 1
IF(modd == 0) THEN
CALL outarr(a,nx0,nx1,i0,i1,is,ixp,ny0,ny1,j0,j1,js,iyp, &
nz0,nz1,k0,k1,ks,izp,nplan,leng,title,nunit, &
scale,nexp,factor,ndigit,tkoff,modd)
ELSE
IF(modd == 1) THEN
nap=k1-k0+1
ELSE IF(modd == 2) THEN
nap=j1-j0+1
ELSE
nap=i1-i0+1
END IF
CALL outarr(a,nx0,nx1,i0,i1,is,iap,ny0,ny1,j0,j1,js, &
iap, &
nz0,nz1,k0,k1,ks,iap,nap,leng,title,nunit, &
scale,nexp,factor,ndigit,tkoff,modd)
END IF
1001 FORMAT(/t2,a8,t30,'Min=',e14.7,' Max=',e14.7)
RETURN
END SUBROUTINE wrigar
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE OUTARR ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE outarr (a,l0,l1,i0,i1,is,ixp,m0,m1,j0,j1,js,iyp, & 2
n0,n1,k0,k1,ks,isp,nplan,leng,title,nchan,scale,nexp, &
factor,ndigit,tkoff,mode)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Print a formatted table of a 2-d array. This subroutine is called
! by WRIGRA.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 12/10/89.
!
! MODIFICATION HISTORY:
!
! 5/06/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
IMPLICIT NONE
INTEGER :: l0,l1,m0,m1,n0,n1
REAL :: a(l0:l1,m0:m1,n0:n1)
REAL :: tkoff, factor, scale
INTEGER :: mode, ndigit
CHARACTER (LEN=*) :: title
CHARACTER (LEN=9) :: cmod(3)
CHARACTER (LEN=2) :: numero(39)
CHARACTER (LEN=20) :: forma1
CHARACTER (LEN=10) :: forma2
CHARACTER (LEN=13) :: forma3
INTEGER :: nplan
INTEGER :: ixp(nplan),iyp(nplan),isp(nplan)
INTEGER :: leng,lengl,i0,i1,j0,j1,k0,k1,is,js,ks
INTEGER :: nchan,nexp,lfield,line,ndec,ndig
INTEGER :: i11,i12,li,imax,kp,j11,j12,mi,jmax,jp,ip
INTEGER :: i,j,k
SAVE cmod,forma1,numero,forma2,forma3
DATA cmod/'plane is=','plane iy=','plane ix='/
DATA numero/'01','02','03','04','05','06','07','08','09' &
,'10','11','12','13','14','15','16','17','18','19' &
,'20','21','22','23','24','25','26','27','28','29' &
,'30','31','32','33','34','35','36','37','38','39'/
DATA forma1/'(1x,i2,''*'',T6,31I04)'/
DATA forma2/'(t6,31i04)'/
DATA forma3/'(i4,04e30.23)'/
!
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
lengl=leng-6
lfield=ndigit+2
line=lengl/lfield
forma1(18:19)=numero(ndigit+2)
forma1(15:16)=numero(line)
forma2(8:9)=numero(ndigit+2)
forma2(5:6)=numero(line)
ndec=MAX(1,ndigit+2-7)
ndig=ndec+7
forma3(8:9)=numero(ndig)
forma3(11:12)=numero(ndec)
forma3(5:6)=numero(line)
IF(mode == 0 .OR. mode == 1) THEN
i11=i0+MIN(i1-i0,(line-1)*is)
imax=MIN(lengl,(lfield*((i1-i0)/is+1)))
li=i0+(line-1)*is
i12=MIN(i1,li)
DO kp=1,nplan
! IF(isp(kp).ge.n0) THEN
k=k0-1+isp(kp)
WRITE(nchan,10) title,tkoff,cmod(1),k,nexp
WRITE(nchan,forma2) (i,i=i0,i11,is)
WRITE(nchan,4) ('*',i=2,imax)
DO j=j1,j0,-js
IF(ndigit < 15) THEN
WRITE(nchan,forma1) j,( nint((a(i,j,k)-tkoff)*factor) &
,i=i0,i12,is)
ELSE
WRITE(nchan,forma3) j,(a(i,j,k),i=i0,i12,is)
END IF
END DO
WRITE(nchan,4) ('*',i=2,imax)
! ENDIF
END DO
END IF
IF(mode == 0 .OR. mode == 2) THEN
i11=i0+MIN(i1-i0,(line-1)*is)
imax=MIN(lengl,(lfield*((i1-i0)/is+1)))
li=i0+(line-1)*is
i12=MIN(i1,li)
DO jp=1,nplan
! IF(iyp(jp).ge.j0) THEN
j=j0-1+iyp(jp)
WRITE(nchan,10) title,tkoff,cmod(2),j,nexp
WRITE(nchan,forma2) (i,i=i0,i11,is)
WRITE(nchan,4) ('*',i=2,imax)
DO k=k1,k0,-ks
IF(ndigit < 15) THEN
WRITE(nchan,forma1) k,(nint((a(i,j,k)-tkoff)*factor) &
,i=i0,i12,is)
ELSE
WRITE(nchan,forma3) k,(a(i,j,k),i=i0,i12,is)
END IF
END DO
WRITE(nchan,4) ('*',i=2,imax)
! ENDIF
END DO
END IF
IF(mode == 0 .OR. mode == 3) THEN
j11=j0+MIN(j1-j0,(line-1)*js)
jmax=MIN(lengl,(lfield*((j1-j0)/js+1)))
mi=j0+(line-1)*js
j12=MIN(j1,mi)
DO ip=1,nplan
! IF(ixp(ip).ge.i0) THEN
i=i0-1+ixp(ip)
WRITE(nchan,10) title,tkoff,cmod(3),i,nexp
WRITE(nchan,forma2) (j,j=j11,j0,-js)
WRITE(nchan,4) ('*',j=2,jmax)
DO k=k1,k0,-ks
IF(ndigit < 15) THEN
WRITE(nchan,forma1) k,(nint((a(i,j,k)-tkoff)*factor) &
,j=j12,j0,-js)
ELSE
WRITE(nchan,forma3) k,(a(i,j,k),j=j12,j0,-js)
END IF
END DO
WRITE(nchan,4) ('*',j=2,jmax)
! ENDIF
END DO
END IF
RETURN
! 2 FORMAT(t30,'Units of 10**',i3)
! 3 FORMAT(t6,31I4)
4 FORMAT(t7,125A1)
! 5 FORMAT(1X,i2,'*',t6,31I4)
10 FORMAT(//t7,a,'(add:',e10.3,')',3X,a,i4,5X,'Units of 10**' &
,i3,/)
END SUBROUTINE outarr
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTBASFN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gtbasfn(fnkey,dirname,ldirnam,hdmpfmt,mgrid,nestgrd, & 10,1
basdmpfn,lbasdmpf )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return a unique name for the grid and base state array dump file.
! The naming convention of the history data dump is:
!
! fnkey.bingrdbas.aa unformatted binary data set.
! fnkey.ascgrdbas.aa formatted ASCII data set.
! fnkey.hdfgrdbas.aa HDF data set
! fnkey.pakgrdbas.aa packed binary data set
!
! where fnkey is a string for name construction and aa is a two
! digit number appended to the data set name in case a data file
! called fnkey.hdfgrdbas.(aa-1) already exists on the disk.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 11/20/91
!
! MODIFICATION HISTORY:
!
! 5/05/92 (M. Xue)
! Added full documentation.
!
! 7/30/92 (M. Xue)
! Added option for packed data set.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! fnkey A character string used to name the files
! dirname A file directory name to be attached to the file name.
! ldirnam Length of the directory string name
! hdmpfmt Parameter specifying the format of output data set.
! mgrid The grid number
! nestgrd Flag for nested grid run.
!
!-----------------------------------------------------------------------
!
! OUTPUT:
!
! basdmpfn The name of the data dump file.
! lbasdmpf The length of character string basdmpfn.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: fnkey ! A character string used to name the
! files
CHARACTER (LEN=* ) :: dirname ! A string giving the file directory
! name
INTEGER :: ldirnam ! The length of the directory name
INTEGER :: hdmpfmt ! Parameter specifying the format of
! output data set.
INTEGER :: mgrid ! The grid number
INTEGER :: nestgrd ! Flag for nested grid run.
CHARACTER(LEN=*) :: basdmpfn ! Name of grid & base state data file.
INTEGER :: lbasdmpf ! Length of character string basdmpfn.
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: lfnkey
CHARACTER (LEN=80) :: temchar
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
lfnkey = LEN( fnkey)
IF( hdmpfmt == 1.OR.hdmpfmt == 0) THEN ! Unformatted binary dump
basdmpfn = fnkey(1:lfnkey)//'.bingrdbas'
ELSE IF( hdmpfmt == 2 ) THEN ! Formatted ASCII dump
basdmpfn = fnkey(1:lfnkey)//'.ascgrdbas'
ELSE IF( hdmpfmt == 3 ) THEN ! HDF data dump
basdmpfn = fnkey(1:lfnkey)//'.hdfgrdbas'
ELSE IF( hdmpfmt == 4 ) THEN ! Packed binary dump
basdmpfn = fnkey(1:lfnkey)//'.pakgrdbas'
ELSE IF( hdmpfmt == 5 ) THEN ! Savi3D data dump
!-----------------------------------------------------------------------
! For Savi3D data dump, the grid and base state information is
! always written together with the other fields.
!-----------------------------------------------------------------------
ELSE IF( hdmpfmt == 6 ) THEN ! Binary with skipping
basdmpfn = fnkey(1:lfnkey)//'.bn2grdbas'
ELSE IF( hdmpfmt == 7 ) THEN ! NetCDF format
basdmpfn = fnkey(1:lfnkey)//'.netgrdbas'
ELSE IF( hdmpfmt == 8 ) THEN ! Packed NetCDF format
basdmpfn = fnkey(1:lfnkey)//'.npkgrdbas'
ELSE IF( hdmpfmt == 9 ) THEN ! GrADS data dump
!-----------------------------------------------------------------------
! For GrADS data dump, the grid and base state information is
! always written together with the other fields.
!-----------------------------------------------------------------------
ELSE IF( hdmpfmt == 10 ) THEN ! GRIB format
basdmpfn = fnkey(1:lfnkey)//'.grbgrdbas'
END IF
lbasdmpf = 10 + LEN(fnkey)
IF(nestgrd == 1) THEN
WRITE(basdmpfn((lbasdmpf+1):(lbasdmpf+4)),'(a,i2.2)')'.g',mgrid
lbasdmpf = lbasdmpf + 4
END IF
IF (mp_opt > 0) THEN
temchar = basdmpfn
WRITE(basdmpfn,'(a,a,2i2.2)') trim(temchar),'_',loc_x,loc_y
lbasdmpf = lbasdmpf + 5
END IF
IF( dirname /= ' ' ) THEN
temchar = basdmpfn
basdmpfn = dirname(1:ldirnam)//'/'//temchar
lbasdmpf = lbasdmpf + ldirnam + 1
END IF
CALL fnversn(basdmpfn, lbasdmpf)
RETURN
END SUBROUTINE gtbasfn
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTDMPFN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gtdmpfn(fnkey,dirname,ldirnam,curtim,hdmpfmt, & 17,2
mgrid,nestgrd, hdmpfn, ldmpf )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return a unique name for the history data dump at time 'curtim'.
! The naming convention of the history data dump is:
!
! fnkey.hdfnnnnnn.aa for hdf data set
! fnkey.binnnnnnn.aa for unformatted binary data set.
!
! where fnkey is a string for name construction, nnnnnn indicates
! the time of the data set in hour/minute/second format, and aa is
! a two digit number appended to the data set name in case a data
! file called fnkey.hdfnnnnnn.(aa-1) already exists on the disk.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 11/20/91
!
! MODIFICATION HISTORY:
!
! 5/05/92 (M. Xue)
! Added full documentation.
!
! 7/30/92 (M. Xue)
! Added option for packed data set.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! fnkey A character string for the name of this run.
! dirname A file directory name to be attached to the file name.
! ldirnam Length of the directory string name
! curtim The model time in seconds.
! hdmpfmt Parameter specifying the format of output data set.
! mgrid The grid number
! nestgrd Flag for nested grid run.
!
!-----------------------------------------------------------------------
!
! OUTPUT:
!
! hdmpfn The name of the history data dump file.
! ldmpf The length of character string hdmpfn.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: fnkey ! A character string used to name the
! files
CHARACTER (LEN=* ) :: dirname ! A string giving the file directory
! name.
INTEGER :: ldirnam ! The length of the directory name.
REAL :: curtim ! Current model time.
INTEGER :: hdmpfmt ! Parameter specifying the format of
! output data set.
INTEGER :: mgrid ! The grid number.
INTEGER :: nestgrd ! Flag for nested grid run.
CHARACTER(LEN=*) :: hdmpfn ! Name of the history data dump file.
INTEGER :: ldmpf ! Length of character string hdmpfn.
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
CHARACTER (LEN=80) :: timsnd
CHARACTER (LEN=80) :: temchar
INTEGER :: tmstrln
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
CALL cvttsnd( curtim, timsnd, tmstrln )
IF( hdmpfmt == 1.OR.hdmpfmt == 0) THEN ! Unformatted binary dump
WRITE(hdmpfn,'(a,a)') fnkey, '.bin'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 2 ) THEN ! Formatted ASCII dump
WRITE(hdmpfn,'(a,a)') fnkey, '.asc'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 3 ) THEN ! HDF data dump
WRITE(hdmpfn,'(a,a)') fnkey, '.hdf'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 4 ) THEN ! Packed binary dump
WRITE(hdmpfn,'(a,a)') fnkey, '.pak'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 5 ) THEN ! Data dump for Savi3D
WRITE(hdmpfn,'(a,a)') fnkey, '.svi'
ELSE IF( hdmpfmt == 6 ) THEN ! Binary with skipping
WRITE(hdmpfn,'(a,a)') fnkey, '.bn2'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 7 ) THEN ! NetCDF format
WRITE(hdmpfn,'(a,a)') fnkey, '.net'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 8 ) THEN ! Packed NetCDF format
WRITE(hdmpfn,'(a,a)') fnkey, '.npk'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 9 ) THEN ! Data dump for GrADS
WRITE(hdmpfn,'(a,a)') fnkey, '.gad'
ELSE IF( hdmpfmt == 10 ) THEN ! Data dump for GrADS
WRITE(hdmpfn,'(a,a)') fnkey, '.grb'//timsnd(1:tmstrln)
ELSE IF( hdmpfmt == 11 ) THEN ! Data dump for Vis5D
WRITE(hdmpfn,'(a,a)') fnkey, '.v5d'//timsnd(1:tmstrln)
END IF
IF( hdmpfmt == 5 .OR. hdmpfmt == 9 ) THEN
ldmpf = LEN(fnkey) + 4
ELSE
ldmpf = LEN(fnkey) + 4 + tmstrln
END IF
IF(nestgrd == 1) THEN
WRITE(hdmpfn((ldmpf+1):(ldmpf+4)), '(a,i2.2)') '.g',mgrid
ldmpf = ldmpf + 4
END IF
IF (mp_opt > 0) THEN
temchar = hdmpfn
WRITE(hdmpfn,'(a,a,2i2.2)') trim(temchar),'_',loc_x,loc_y
ldmpf = ldmpf + 5
END IF
IF( dirname /= ' ' ) THEN
temchar = hdmpfn
hdmpfn = dirname(1:ldirnam)//'/'//temchar
ldmpf = ldmpf + ldirnam + 1
END IF
CALL fnversn(hdmpfn, ldmpf)
RETURN
END SUBROUTINE gtdmpfn
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTRSTFN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gtrstfn(fnkey,dirname,ldirnam,curtim, & 1,2
mgrid,nestgrd, rstoutf, lrstof )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return a unique name for the restart data dump at time 'curtim'.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 03/12/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! fnkey A character string used to name the files
! dirname A file directory name to be attached to the file name.
! ldirnam Length of the directory string name
! curtim The model time in seconds.
! mgrid The grid number
! nestgrd Flag for nested grid run.
!
!-----------------------------------------------------------------------
!
! OUTPUT:
!
! rstoutf The name of the restart dump file.
! lrstof The length of the file name.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: fnkey ! A character string used to name the
! files
CHARACTER (LEN=* ) :: dirname ! A string giving the file directory
! name.
INTEGER :: ldirnam ! The length of the directory name.
REAL :: curtim ! Current model time.
INTEGER :: mgrid ! The grid number.
INTEGER :: nestgrd ! Flag for nested grid run.
CHARACTER (LEN=*) :: rstoutf ! Name of the history data dump file.
INTEGER :: lrstof ! Length of character string hdmpfn.
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
CHARACTER (LEN=80) :: timsnd
CHARACTER (LEN=80) :: temchar
INTEGER :: tmstrln
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!
!-----------------------------------------------------------------------
!
! Construct the file name:
!
! The format of the restart file names is
!
! fnkey.rstnnnnnn
!
! where fnkey are characters used to name the files.
! nnnnnn is a 6 digit integer number indicating the time of
! the data in seconds.
!
!-----------------------------------------------------------------------
!
CALL cvttsnd( curtim, timsnd, tmstrln )
WRITE(rstoutf,'(a,a)') fnkey, '.rst'//timsnd(1:tmstrln)
lrstof = LEN(fnkey) + 4 + tmstrln
IF( nestgrd == 1 ) THEN
!-----------------------------------------------------------------------
!
! Attach the grid number to the file name
!
!-----------------------------------------------------------------------
WRITE(rstoutf((lrstof+1):(lrstof+4)), '(a,i2.2)') '.g',mgrid
lrstof = lrstof + 4
END IF
IF (mp_opt > 0) THEN
temchar = rstoutf
WRITE(rstoutf,'(a,a,2i2.2)') trim(temchar),'_',loc_x,loc_y
lrstof = lrstof + 5
END IF
IF( dirname /= ' ' ) THEN
temchar = rstoutf(1:lrstof)
rstoutf = dirname(1:ldirnam)//'/'//temchar
lrstof = lrstof + ldirnam + 1
END IF
!
!-----------------------------------------------------------------------
!
! Append a version number to the file name if the named file already
! exists.
!
!-----------------------------------------------------------------------
!
CALL fnversn(rstoutf, lrstof)
RETURN
END SUBROUTINE gtrstfn
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTLOGFN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gtlogfn( fnkey, mgrid, nestgrd, logfn, llogfn ) 2,1
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return a unique name for the input log file. The name convention
! of the log file is:
!
! fnkey.log.aa
!
! where fnkey is a character string for naming the files and aa is
! a two digit number appended to the data set name in case the file
! fnkey.log.(aa-1) already exists.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 3/17/1991.
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! fnkey A character string used to name the files
! mgrid The grid number
! nestgrd Flag for nested grid run.
!
! OUTPUT:
!
! logfn Log file filename.
! llogfn The length of character string logfn.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: fnkey ! A character string used to name the
! files
INTEGER :: mgrid ! The grid number
INTEGER :: nestgrd ! Flag for nested grid run.
CHARACTER (LEN=80 ) :: logfn ! Log file filename.
INTEGER :: llogfn ! The length of character string logfn.
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
WRITE(logfn,'(a,a)') fnkey, '.log'
llogfn = 4 + LEN(fnkey)
IF( nestgrd == 1 ) THEN
WRITE(logfn((llogfn+1):(llogfn+4)), '(a,i2.2)') '.g',mgrid
llogfn = llogfn + 4
END IF
CALL fnversn(logfn, llogfn )
RETURN
END SUBROUTINE gtlogfn
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE FNVERSN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE fnversn( filename, fnlen ) 40,1
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Append the file version number to a file name if the named
! file already exists.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 3/17/1991.
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! filename name of the file on input
! fnlen file name length on input
!
! OUTPUT:
!
! filename name of the file on output
! fnlen file name length on output
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: filename ! name of the file on input/output
CHARACTER (LEN=132) :: CHAR
INTEGER :: fnlen ! file name length on input/output
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
LOGICAL :: iexist, iexist1, iexist2
INTEGER :: nnn, fnlen_old, fnlen_tem
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
CHAR = filename
fnlen_old = fnlen
fnlen_tem = fnlen
nnn = 0
200 INQUIRE(FILE=CHAR(1:fnlen_tem), EXIST=iexist)
INQUIRE(FILE=CHAR(1:fnlen_tem)//'.Z' , EXIST=iexist1)
INQUIRE(FILE=CHAR(1:fnlen_tem)//'.gz', EXIST=iexist2)
IF(iexist.OR.iexist1.OR.iexist2)THEN
nnn = nnn+1
IF( nnn > 99) THEN
WRITE(6,'(/1x,a,/1x,a/)') &
'An alternative name could not be found for ', &
CHAR(1:fnlen_old),' Job stopped in FNVERSN.'
CALL arpsstop(' ',1)
END IF
WRITE(CHAR((fnlen_old+1):(fnlen_old+3)),'(a,i2.2)')'.',nnn
fnlen_tem = fnlen_old + 3
GO TO 200
END IF
fnlen = fnlen_tem
filename(1:fnlen) = CHAR(1:fnlen_tem)
RETURN
END SUBROUTINE fnversn
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE STRLNTH ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE strlnth( string, length ) 176
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return the length of the non-blank part of a character string.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 11/20/91
!
! MODIFICATION HISTORY:
!
! 5/05/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! string A character string
! length The declared length of the character string 'string'.
!
! OUTPUT:
!
! length The length of the non-blank part of the string.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: string ! A character string for the name of
! this run.
INTEGER :: length ! The length of the non-blank part
! of a string.
INTEGER :: i
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
IF (length == 0) length = LEN(string)
DO i = length,1,-1
IF(string(i:i) /= ' ') EXIT
END DO
! 200 CONTINUE
length = MAX(1,i)
RETURN
END SUBROUTINE strlnth
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE STRMIN ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE strmin( string, length )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Minimize a string length by removing consecutive blank spaces.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 1/15/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! string A character string
! length The declared length of the character string 'string'.
!
! OUTPUT:
!
! length The length of string with consecutive blank spaces
! removed.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: string ! A character string for the name of
! this run.
INTEGER :: length ! The length of the non-blank part
! of a string.
CHARACTER (LEN=1) :: str_1
CHARACTER (LEN=256) :: str
INTEGER :: i,len_old
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
IF( length > 256) THEN
PRINT*,'Work string defined in STRMIN was too small.'
PRINT*,'The output from this subroutine may not be correct.'
length=256
END IF
len_old = length
length = 1
str = string
DO i = 2,len_old
str_1 = str(i-1:i-1)
IF(.NOT.(str(i:i) == ' '.AND. &
(str_1 == ' '.OR.str_1 == '('.OR.str_1 == '='))) THEN
length=length+1
string(length:length)=str(i:i)
END IF
END DO
! 200 CONTINUE
RETURN
END SUBROUTINE strmin
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE COMLNTH ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE comlnth( string, length )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return the length of the non-blank part of a character string.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 11/20/91
!
! MODIFICATION HISTORY:
!
! 5/05/92 (M. Xue)
! Added full documentation.
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! string A character string
! length The declared length of the character string 'string'.
!
! OUTPUT:
!
! length The length of the non-blank part of the string.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: string ! A character string for the name of
! this run.
INTEGER :: length ! The length of the non-blank part
! of a string.
INTEGER :: i
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
DO i = 1, 80
IF(string(i:i) == ' ') THEN
IF (string(i+1:i+1) == ' ') EXIT
END IF
END DO
length = MAX(1,i-1)
RETURN
END SUBROUTINE comlnth
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE CVTTIM ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE cvttim(timsnd, timhms) 1
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Convert time, given in seconds, into a 6 character string
! containing time in the hour/minute/second format.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 11/02/92
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! timsnd Time in seconds
!
! OUTPUT:
!
! timhms string contain time in hour:minute:second format
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
REAL :: timsnd ! Time in seconds
CHARACTER (LEN=6 ) :: timhms ! string contain time in
! hour:minute:second format
INTEGER :: h,m,s
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
h = INT(timsnd/3600.0 )
m = INT((timsnd-h*3600.0)/60.0)
s = nint(timsnd-h*3600.0-m*60.0)
IF( s == 60) THEN
m = m+1
s = 0
END IF
WRITE(timhms,'(3i2.2)') h,m,s
RETURN
END SUBROUTINE cvttim
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE CVTTSND ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE cvttsnd(time, timsnd, tmstrln) 18,2
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Convert time given in second into a character string
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 07/17/2000
!
! MODIFICATION HISTORY:
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! time Time in seconds
!
! OUTPUT:
!
! timsnd Time string in seconds
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
REAL :: time ! Time in seconds
CHARACTER (LEN=*) :: timsnd ! Time string in seconds
INTEGER :: tmstrln ! Length of time string
INTEGER :: itime
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
itime = nint(time)
IF ( itime < 999999 ) THEN
tmstrln=6
IF( LEN(timsnd) < tmstrln) GO TO 10
WRITE(timsnd,'(i6.6)') itime
ELSE IF ( itime < 9999999 ) THEN
tmstrln=7
IF( LEN(timsnd) < tmstrln) GO TO 10
WRITE(timsnd,'(i7.7)') itime
ELSE IF ( itime < 99999999 ) THEN
tmstrln=8
IF( LEN(timsnd) < tmstrln) GO TO 10
WRITE(timsnd,'(i8.8)') itime
ELSE IF ( itime < 999999999 ) THEN
tmstrln=9
IF( LEN(timsnd) < tmstrln) GO TO 10
WRITE(timsnd,'(i9.9)') itime
ELSE
WRITE (6,'(a/a,i4/a,e16.8,a)') &
'WARNING: The time is too large to fit in 9 characters', &
' tmstrln = ',tmstrln, &
' time = ',time, ' seconds.'
CALL arpsstop('arpsstop called from CVTTSND time too large ',1)
END IF
RETURN
10 CONTINUE
WRITE(6,'(a/a/a)') &
'String timsnd passed into CVTTSND not long enough.', &
'Need to be at least ',tmstrln,' charactere long.', &
'Job stopped in CVTTSND.'
CALL arpsstop('arpsstop called from CVTTSND string insufficient &
& length',1)
RETURN
END SUBROUTINE cvttsnd
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTLFNKEY ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gtlfnkey( runname, lfnkey ) 19
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Find out the number of characters to be used to construct file
! names.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 03/15/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
! 12/12/1996 (Yuhe Liu)
! Removed the restrict of 6 characters to runname.
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
CHARACTER (LEN=* ) :: runname ! Input
INTEGER :: lfnkey ! Output
INTEGER :: lenstr, firstb, firstc
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
lenstr = LEN( runname )
firstb = INDEX( runname, ' ')
firstc = INDEX( runname, ',')
IF( firstb == 0) firstb = lenstr+1
IF( firstc == 0) firstc = lenstr+1
lfnkey = MAX(1, MIN( lenstr, firstb-1, firstc-1 ) )
RETURN
END SUBROUTINE gtlfnkey
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GETUNIT ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE getunit( nunit ) 93,1
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Find a free FORTRAN I/O unit from a list and return that unit.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 4/21/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! OUTPUT:
!
! nunit A free fortran I/O unit number.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: nunit
LOGICAL :: used
INTEGER :: list(30), nfree
SAVE list, nfree
DATA list /10,11,12,13,14,15,16,17,18,19, &
20,21,22,23,24,25,26,27,28,29, &
30,31,32,33,34,35,36,37,38,39/
DATA nfree /30/
INCLUDE 'mp.inc' ! Message passing parameters.
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
5 IF( nfree < 1) THEN
IF (myproc == 0) WRITE(6,'(1x,a,a)') &
'No more unit number is available from the list. ', &
'Job stopped in GETUNIT.'
CALL arpsstop('arpsstop called from GETINIT out of file numbers',1)
END IF
nunit = list(nfree)
nfree = nfree-1
INQUIRE( UNIT=nunit, OPENED=used)
IF( used ) GO TO 5
IF (myproc == 0) WRITE(6,'(1x,a,i3,a)') 'Fortran I/O unit ',nunit, &
' picked from the free list.'
RETURN
!
!
!##################################################################
!##################################################################
!###### ######
!###### ENTRY RETUNIT ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
ENTRY retunit( nunit )
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Return a freed unit to the list.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 4/21/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! nunit A freed fortran I/O unit number to be returned.
!
!-----------------------------------------------------------------------
!
INQUIRE( UNIT=nunit, OPENED=used)
IF( used ) RETURN
nfree = nfree + 1
IF( nfree <= 30 ) list( nfree ) = nunit
nfree = MIN( nfree, 30)
IF (myproc == 0 ) WRITE(6,'(1x,a,i3,a)') 'Fortran I/O unit ',nunit, &
' returned to the free list.'
RETURN
END SUBROUTINE getunit
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE SECTHRZ ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE secthrz(nx,ny,nz,s,z,ss1) 2
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Interpolate 3-D data to a given horizontal level.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue & Hao Jin
! 12/18/92.
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
! INPUT :
!
! nx Number of grid points in the x-direction (east/west)
! ny Number of grid points in the y-direction (north/south)
! nz Number of grid points in the vertical
!
! s 3-dimensional array of data to contour
! s1 2-dimensional array of data to contour
!
! x x coordinate of grid points in physical/comp. space (m)
! y y coordinate of grid points in physical/comp. space (m)
! z z coordinate of grid points in physical space (m)
!
! z01 value of x for first i grid point to plot
!
! OUTPUT:
! ss1 interpolated 3-D data to a given horizontal level
!
!-----------------------------------------------------------------------
!
! Parameters of output
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: nx,ny,nz ! Number of grid points in 3 directions
REAL :: s(nx,ny,nz) ! 3-dimensional array of data to contour
REAL :: z(nx,ny,nz) ! z coordinate of grid points
! in physical space (m)
REAL :: ss1(nx,ny) ! interpolated 3-D data to a
! given horizontal level
!
!-----------------------------------------------------------------------
!
! Common blocks for plotting control parameters
!
!-----------------------------------------------------------------------
!
REAL :: z01 ! the given height of the slice
COMMON /sliceh/z01
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i,j,k
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'globcst.inc'
!
!-----------------------------------------------------------------------
!
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!
!-----------------------------------------------------------------------
!
! Find index for interpolation
!
!-----------------------------------------------------------------------
!
DO i=1,nx-1
DO j=1,ny-1
IF(z01 <= z(i,j,1)) GO TO 11
IF(z01 >= z(i,j,nz-1)) GO TO 12
DO k=2,nz-2
IF(z01 >= z(i,j,k).AND.z01 < z(i,j,k+1)) GO TO 15
END DO
11 k=1
GO TO 15
12 k=nz-1
GO TO 15
15 ss1(i,j)=s(i,j,k)+(s(i,j,k+1)-s(i,j,k))* &
(z01-z(i,j,k))/(z(i,j,k+1)-z(i,j,k))
!-----------------------------------------------------------------------
!
! If the data point is below the ground level, set the
! data value to the missing value.
!
!-----------------------------------------------------------------------
IF( z01 < z(i,j,2) ) ss1(i,j) = -9999.0
END DO
END DO
RETURN
END SUBROUTINE secthrz
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE SECTVRT ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE sectvrt(nx,ny,nz,s,x,y,z,ss1,zs1,n,xp,yp) 1
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Interpolate a 3-D data to 2-d vectical plane.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue & Hao Jin
! 12/18/92.
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
! INPUT :
!
! nx Number of grid points in the x-direction (east/west)
! ny Number of grid points in the y-direction (north/south)
! nz Number of grid points in the vertical
!
! s 3-dimensional array of data to contour
! s1 2-dimensional array of data to contour
!
! x x coordinate of grid points in physical/comp. space (m)
! y y coordinate of grid points in physical/comp. space (m)
! z z coordinate of grid points in computational space (m)
!
! OUTPUT:
! ss1 interpolated a 3-D data to 2-d vectical plane
!
!
!-----------------------------------------------------------------------
!
! Parameters of output
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: nx,ny,nz ! Number of grid points in 3 directions
INTEGER :: n
REAL :: s(nx,ny,nz) ! 3-dimensional array of data to contour
REAL :: x(nx,ny,nz) ! x coordinate of grid points
! in physical/comp. space (m)
REAL :: y(nx,ny,nz) ! y coordinate of grid points
! in physical/comp. space (m)
REAL :: z(nx,ny,nz) ! z coordinate of grid points
! in computational space (m)
REAL :: ss1(n,nz)
REAL :: zs1(n,nz)
REAL :: xp(n)
REAL :: yp(n)
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: i,k
INTEGER :: is,js
REAL :: s1,s2,s3,s4,sgrid,xs1,ys1
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
INCLUDE 'globcst.inc'
INCLUDE 'grid.inc' ! Grid & map parameters.
!
!-----------------------------------------------------------------------
!
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!
!-----------------------------------------------------------------------
!
! Intepolate the date to the given point.
!
!-----------------------------------------------------------------------
!
DO k=1,nz-1
DO i=1,n
xs1=xp(i)
ys1=yp(i)
is = MAX(1, MIN(nx-2, INT( (xs1-x(1,2,2))/dx+1 ) ))
js = MAX(1, MIN(ny-2, INT( (ys1-y(2,1,2))/dy+1 ) ))
s1 = (xs1-x(is ,js ,k))*(ys1-y(is ,js ,k))
s2 =-(xs1-x(is+1,js ,k))*(ys1-y(is+1,js ,k))
s3 = (xs1-x(is+1,js+1,k))*(ys1-y(is+1,js+1,k))
s4 =-(xs1-x(is ,js+1,k))*(ys1-y(is ,js+1,k))
sgrid = (x(is+1,js,k)-x(is,js,k))*(y(is,js+1,k)-y(is,js,k))
zs1(i,k) =(z(is ,js ,k)*s3+z(is+1,js ,k)*s4 &
+z(is+1,js+1,k)*s1+z(is ,js+1,k)*s2)/sgrid
ss1(i,k) =(s(is ,js ,k)*s3+s(is+1,js ,k)*s4 &
+s(is+1,js+1,k)*s1+s(is ,js+1,k)*s2)/sgrid
END DO
END DO
RETURN
END SUBROUTINE sectvrt
!
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE REFLEC ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE reflec(nx,ny,nz, rhobar, qr, qs, qh, reflc ) 3
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Compute the radar reflectivity factor following Kessler (1969).
! Here, arg=Z (mm**6/m**3), and dBz = 10log10 (arg).
!
!-----------------------------------------------------------------------
!
! AUTHOR: K. Droegemeier and M.Xue
! 4/19/93
!
! MODIFICATION HISTORY:
!
! 9/10/94 (Weygandt & Y. Lu)
! Cleaned up documentation.
!
! 12/6/95 (J. Zong and M. Xue)
! Added qs and qh to the argument list of this subroutine to
! facilitate inclusion of the contributions of qs and qh to reflec-
! tivity. A relation between radar reflectivity factor and snow
! content is adopted from Rogers and Yau (1989) and extended to
! represent the effects of snow and graupel/hail on the
! reflectivity. globcst.inc is included to pass the value of ice.
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! nx Number of grid points in the x-direction (east/west)
! ny Number of grid points in the y-direction (north/south)
! nz Number of grid points in the vertical
!
! rhobar Base state density (kg/m**3)
! qr Rainwater mixing ratio (kg/kg)
! qs Snow mixing ratio (kg/kg)
! qh Hail mixing ratio (kg/kg)
!
! OUTPUT:
!
! reflc Radar reflectivity factor.
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
! Include files:
!
!-----------------------------------------------------------------------
!
!
INCLUDE 'globcst.inc'
!
!-----------------------------------------------------------------------
!
INTEGER :: nx,ny,nz
REAL :: rhobar(nx,ny,nz) ! Base state air density (kg/m**3)
REAL :: qr (nx,ny,nz) ! Rain water mixing ratio (kg/kg)
REAL :: qs (nx,ny,nz) ! Snow mixing ratio (kg/kg)
REAL :: qh (nx,ny,nz) ! Hail mixing ratio (kg/kg)
REAL :: reflc (nx,ny,nz) ! Radar reflectivity (dBZ)
!
!-----------------------------------------------------------------------
!
! Misc. local variables
!
!-----------------------------------------------------------------------
!
INTEGER :: i,j,k
REAL :: arg,svnfrth
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!
!-----------------------------------------------------------------------
!
! Compute the radar reflectivity factor following Kessler (1969).
! Here, arg=Z (mm**6/m**3), and dBz = 10log10 (arg).
!
!-----------------------------------------------------------------------
!
svnfrth = 7./4.
DO k = 1,nz-1
DO j = 1,ny-1
DO i = 1,nx-1
arg = 17300.0*( rhobar(i,j,k)*1000.0 &
* MAX(0.0,qr(i,j,k)) )**svnfrth
IF (ice == 1) THEN
arg = arg + 38000.0*( rhobar(i,j,k)*1000.0 &
* MAX(0.0,qs(i,j,k)+qh(i,j,k)) )**2.2
END IF
reflc(i,j,k) = 10.0*ALOG10( MAX(arg,1.0) )
END DO
END DO
END DO
RETURN
END SUBROUTINE reflec
!########################################################################
!########################################################################
!######### #########
!######### SUBROUTINE reflec_ferrier #########
!######### #########
!######### Developed by #########
!######### Center for Analysis and Prediction of Storms #########
!######### University of Oklahoma #########
!######### #########
!########################################################################
!########################################################################
SUBROUTINE reflec_ferrier(nx,ny,nz, rho, qr, qs, qh, t, rff)
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! This subroutine estimates logarithmic radar reflectivity using
! equations customized for use with the ARPS Lin-Tao microphysics
! package.
!
! The equations (after some algebra to optimize code) are as follows:
!
! dBZ = 10*LOG10(Ze)
!
! Ze = Zer + Zes + Zeh (contributions from rain, snow, and hail).
!
! k * 720 1.75
! Zer = --------------------------- * (rho * qr)
! 1.75 0.75 1.75
! pi * N0r * rhor
!
! Zes = Zesnegf for "dry" snow (T < 0 C), or Zes = Zesposf for "wet"
! snow (T > 0 C)
!
! 2 0.25
! k * 720 * |K| * rhos
! ice 1.75
! Zesnegf = --------------------------------------- * (rho * qs)
! 1.75 2 0.75 2
! pi * |K| * N0s * rhoice
! water
!
! k * 720 1.75
! Zesposf = --------------------------- * (rho * qs)
! 1.75 0.75 1.75
! pi * N0s * rhos
!
! / k * 720 \ 0.95 1.6625
! Zeh = | --------------------------- | * (rho * qh)
! \ 1.75 0.75 1.75 /
! \ pi * N0h * rhoh /
!
!-----------------------------------------------------------------------
!
! REFERENCES:
!
! Jahn, D., D. Weber, E. Kemp, and H. Neeman, 2000: Evidence of
! convective-induced turbulence outside the immediate storm region:
! Part III. CAPS report submitted to AlliedSignal/Honeywell, 37pp.
!
! Ferrier, B. S., W.-K. Tao, and J. Simpson, 1995: A double-moment
! multiple-phase four-class bulk ice scheme. Part II: Simulations
! of convective storms in different large-scale environments and
! comparisons with other bulk parameterizations. J. Atmos. Sci.,
! 45, 3846-3879.
!
! McCumber, M., W.-K. Tao, J. Simpson, R. Penc, and S.-T. Soong, 1991:
! Comparison of ice-phase microphysical parameterization schemes using
! numerical simulations of tropical convection. J. Appl. Meteor.,
! 30, 985-1004.
!
! Smith, P. L., 1984: Equivalent radar reflectivity factors for snow
! and ice particle. J. Climate Appl. Meteor., 23, 1258-1260.
!
! Smith, P. L., Jr., C. G. Myers, and H. D. Orville, 1975: Radar
! reflectivity factor calculations in numerical cloud models using
! bulk parameterization of precipitation. J. Appl. Meteor., 14,
! 1156-1165.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Henry Neeman, Spring 2000.
!
! MODIFICATION HISTORY:
!
! Eric Kemp, 8 October 2001
! Reformatted code for ARPS Fortran 90 standard.
!
! Ming Xue, 16 Oct. 2001
! Changed ni,nj,nk to nx,ny,nz. Change loop bounds. Changed the order
! of argument list. Removed IF( ice == 0 ) check.
!
!-----------------------------------------------------------------------
!
! Force explicit declarations.
!
!-----------------------------------------------------------------------
IMPLICIT NONE
!-----------------------------------------------------------------------
! Include files.
!-----------------------------------------------------------------------
INCLUDE 'globcst.inc'
!-----------------------------------------------------------------------
! Declare arguments.
!-----------------------------------------------------------------------
INTEGER, INTENT(IN) :: nx,ny,nz ! Dimensions of grid
REAL, INTENT(IN) :: rho(nx,ny,nz) ! Air density (kg m**-3)
REAL, INTENT(IN) :: qr(nx,ny,nz) ! Rain mixing ratio (kg kg**-1)
REAL, INTENT(IN) :: qs(nx,ny,nz) ! Snow mixing ratio (kg kg**-1)
REAL, INTENT(IN) :: qh(nx,ny,nz) ! Hail mixing ratio (kg kg**-1)
REAL, INTENT(IN) :: t(nx,ny,nz) ! Temperature (K)
REAL, INTENT(OUT) :: rff(nx,ny,nz) ! Reflectivity (dBZ)
!-----------------------------------------------------------------------
! Declare local parameters.
!-----------------------------------------------------------------------
REAL,PARAMETER :: ki2 = 0.176 ! Dielectric factor for ice if other
! than melted drop diameters are used.
! REAL,PARAMETER :: ki2 = 0.208 ! Dielectric factor for ice melted drop
! ! diameters are used.
REAL,PARAMETER :: kw2=0.93 ! Dielectric factor for water.
REAL,PARAMETER :: degKtoC=273.15 ! Conversion factor from degrees K to
! degrees C
REAL,PARAMETER :: m3todBZ=1.0E+18 ! Conversion factor from m**3 to
! mm**6 m**-3.
REAL,PARAMETER :: pi=3.1415926 ! Pi.
REAL,PARAMETER :: pipowf=7.0/4.0 ! Power to which pi is raised.
REAL,PARAMETER :: N0r=8.0E+06 ! Intercept parameter in 1/(m^4) for rain.
REAL,PARAMETER :: N0s=3.0E+06 ! Intercept parameter in 1/(m^4) for snow.
REAL,PARAMETER :: N0h=4.0E+04 ! Intercept parameter in 1/(m^4) for hail.
REAL,PARAMETER :: N0xpowf=3.0/4.0 ! Power to which N0r,N0s & N0h are
! raised.
REAL,PARAMETER :: rhoxpowf=7.0/4.0 ! Powers to which rhor, rhos & rhoh
! are raised.
REAL,PARAMETER :: approxpow=0.95 ! Approximation power for hail
! integral.
REAL,PARAMETER :: rqrpowf=7.0/4.0 ! Power to which product rho * qr
! is raised.
REAL,PARAMETER :: rqsnpowf=7.0/4.0 ! Power to which product rho * qs
! is raised (wet snow).
REAL,PARAMETER :: rqsppowf=7.0/4.0 ! Power to which product rho * qs
! is raised (dry snow).
REAL,PARAMETER :: rqhpowf=(7.0/4.0)*approxpow ! Power to which product
! rho * qh is raised.
REAL,PARAMETER :: rhoi=917 ! Density of ice (kg m**-3)
REAL,PARAMETER :: rhor=1000 ! Density of rain (kg m**-3)
REAL,PARAMETER :: rhos=100 ! Density of snow (kg m**-3)
REAL,PARAMETER :: rhoh=913 ! Density of hail (kg m**-3)
REAL,PARAMETER :: rhoipowf=2.0 ! Power to which rhoi is raised.
REAL,PARAMETER :: rhorpowf=7.0/4.0 ! Power to which rhor is raised.
REAL,PARAMETER :: rhospowf=1.0/4.0 ! Power to which rhos is raised.
REAL,PARAMETER :: rhohpowf=7.0/4.0 ! Power to which rhoh is raised.
REAL,PARAMETER :: Zefact=720.0 ! Multiplier for Ze components.
REAL,PARAMETER :: lg10mul=10.0 ! Log10 multiplier
!-----------------------------------------------------------------------
! Declare local variables.
!-----------------------------------------------------------------------
REAL :: rcomp,scomp,hcomp,sumcomp
INTEGER :: i,j,k
REAL :: Zerf,Zesnegf,Zesposf,Zehf
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!-----------------------------------------------------------------------
! First gather all the constants together. (These are treated as
! variables because Fortran 90 does not allow real exponents when
! calculating parameters).
!-----------------------------------------------------------------------
Zerf = (m3todBZ * Zefact) / &
((pi ** pipowf) * (N0r ** N0xpowf) * &
(rhor ** rhoxpowf))
Zesnegf = ((m3todBZ * Zefact * Ki2 * (rhos ** rhospowf)) / &
((pi ** pipowf) * Kw2 * (N0s ** N0xpowf) * &
(rhoi ** rhoipowf)))
Zesposf = ((m3todBZ * Zefact) / &
((pi ** pipowf) * (N0s ** N0xpowf) * &
(rhos ** rhoxpowf)))
Zehf = (((m3todBZ * Zefact) / &
((pi ** pipowf) * (N0h ** N0xpowf) * &
(rhoh ** rhoxpowf))) ** approxpow)
!-----------------------------------------------------------------------
! Now loop through the grid points.
!-----------------------------------------------------------------------
DO k = 1,nz-1
DO j = 1,ny-1
DO i = 1,nx-1
!-----------------------------------------------------------------------
! Check for bad air density value.
!-----------------------------------------------------------------------
IF (rho(i,j,k) <= 0.0) THEN
rff(i,j,k) = 0.0
ELSE
!-----------------------------------------------------------------------
! Calculate reflectivity contribution from rain.
!-----------------------------------------------------------------------
IF (qr(i,j,k) <= 0.0) THEN
rcomp = 0.0
ELSE
rcomp = Zerf * ((rho(i,j,k) * qr(i,j,k)) ** rqrpowf)
END IF
!-----------------------------------------------------------------------
! Calculate reflectivity contribution from snow (dry or wet).
!-----------------------------------------------------------------------
IF (qs(i,j,k) <= 0.0) THEN
scomp = 0.0
ELSE IF (t(i,j,k) <= degKtoC) THEN
scomp = Zesnegf * ((rho(i,j,k) * qs(i,j,k)) ** rqsnpowf)
ELSE
scomp = Zesposf * ((rho(i,j,k) * qs(i,j,k)) ** rqsppowf)
END IF
!-----------------------------------------------------------------------
! Calculate reflectivity contribution from hail.
!-----------------------------------------------------------------------
IF (qh(i,j,k) <= 0.0) THEN
hcomp = 0.0
ELSE
hcomp = Zehf * ((rho(i,j,k) * qh(i,j,k)) ** rqhpowf)
END IF
!-----------------------------------------------------------------------
! Now add the contributions and convert to logarithmic reflectivity
! factor dBZ.
!-----------------------------------------------------------------------
sumcomp = rcomp + scomp + hcomp
rff(i,j,k) = lg10mul * LOG10(MAX(sumcomp,1.0))
END IF ! IF (rho(i,j,k) <= 0.0) ... ELSE ...
END DO ! DO i
END DO ! DO j
END DO ! DO k
RETURN
END SUBROUTINE reflec_ferrier
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE SETCORNERLL ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE setcornerll( nx,ny, x, y ) 8,6
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Set the latitudes and longitudes set all corner points.
!
! Before calling this subroutine, the map projection should have
! been set up.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Yuhe Liu
! 09/30/1997
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! nx Number of grid points for the model
! grid in the east-west direction.
! ny Number of grid points for the model
! grid in the north-south direction.
!
! x Analysis grid points in the e-w direction
! (in grid units)
! y Analysis grid points in the n-s direction
! (in grid units)
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: nx ! Number of model grid points
! in the east-west direction.
INTEGER :: ny ! Number of model grid points
! in the north-south direction
REAL :: x(nx) ! 2-D model grid points east-west
! direction (model grid units)
REAL :: y(ny) ! 2-D model grid points north-south
! direction (model grid units)
!
!-----------------------------------------------------------------------
!
INCLUDE 'globcst.inc'
INCLUDE 'grid.inc' ! Grid & map parameters.
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
REAL :: tema1, tema2, temb1, temb2
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
tema1 = 0.5*(x(1)+x(2))
tema2 = 0.5*(x(nx-1)+x(nx)) + dx ! The nx'th scalar point
temb1 = 0.5*(y(1)+y(2))
temb2 = 0.5*(y(ny-1)+y(ny)) + dy ! The ny'th scalar point
CALL xytoll(1,1, tema1,temb1, swlats,swlons) ! for scalar grid
CALL xytoll(1,1, tema2,temb2, nelats,nelons) ! for scalar grid
CALL xytoll(1,1, x(1), temb1, swlatu,swlonu) ! for u-wind grid
CALL xytoll(1,1, x(nx),temb2, nelatu,nelonu) ! for u-wind grid
CALL xytoll(1,1, tema1,y(1), swlatv,swlonv) ! for v-wind grid
CALL xytoll(1,1, tema2,y(ny), nelatv,nelonv) ! for v-wind grid
RETURN
END SUBROUTINE setcornerll
SUBROUTINE wrtvar(nx,ny,nz, var, varnam, time, runname,dirname) 5,5
!
!-----------------------------------------------------------------------
!
! Purpose:
! To write an array 'var' into a file in binary format.
!
! Author: Ming Xue
!
! Input:
!
! nx, ny, nz Dimensions of input array 'var'.
! var Input array to be written out.
! varnam String of length 6 to designate the input array
! time The time of the input data array in seconds.
!
! Output:
!
! A disk file containing array 'var'.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: nx,ny,nz
REAL :: var(nx,ny,nz)
CHARACTER (LEN=6) :: varnam
CHARACTER (LEN=*) :: runname
CHARACTER (LEN=*) :: dirname
CHARACTER (LEN=80) :: timsnd
CHARACTER (LEN=132) :: vfnam
CHARACTER (LEN=132) :: temchar
INTEGER :: lvfnam, ierr, istat, nunit
INTEGER :: tmstrln
INTEGER :: lrunnam,ldirnam
REAL :: time
INTEGER :: i,j,k
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
CALL cvttsnd( time, timsnd, tmstrln )
CALL gtlfnkey( runname, lrunnam )
vfnam = runname(1:lrunnam)//'.'//varnam//timsnd(1:tmstrln)
lvfnam = len_trim(vfnam)
ldirnam = LEN_trim(dirname)
IF( ldirnam == 0 ) THEN
dirname = '.'
ldirnam = 1
END IF
PRINT*,'ldirnam,dirname=',ldirnam,dirname
IF( dirname /= ' ' ) THEN
temchar = vfnam
vfnam = dirname(1:ldirnam)//'/'//temchar
lvfnam = lvfnam + ldirnam + 1
END IF
PRINT*,'lvfnam,vfnam=',lvfnam,vfnam
PRINT*,'lvfnam,vfnam=',lvfnam,vfnam(1:lvfnam)
CALL getunit( nunit)
CALL asnctl ('NEWLOCAL', 1, ierr)
CALL asnfile(vfnam, '-F f77 -N ieee', ierr)
OPEN(UNIT=nunit,FILE=trim(vfnam(1:lvfnam)),STATUS='unknown', &
FORM='unformatted',IOSTAT= istat )
WRITE(nunit) varnam
WRITE(nunit) nx,ny, nz
DO k=1,nz
WRITE(nunit) ((var(i,j,k),i=1,nx),j=1,ny)
END DO
CLOSE(UNIT=nunit)
CALL retunit(nunit)
RETURN
END SUBROUTINE wrtvar
SUBROUTINE readvar(nx,ny,nz, varnam, time, var, runname),5
!-----------------------------------------------------------------------
!
! Purpose:
! To read in array 'var' from a file.
!
! Author: Ming Xue
!
! Modifications:
! 2/8/1998 (M.Xue)
!
! Added check on the arrays bounds read in from data.
!
! Input:
!
! nx, ny, nz Dimensions of array 'var' to be read in.
! varnam String of length 6 to designate the array
! time The time of the array in seconds.
!
! Output:
!
! var Array readin from a disk file.
!
!-----------------------------------------------------------------------
IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: nx,ny,nz
CHARACTER (LEN=6) :: varnam_in
CHARACTER (LEN=*) :: varnam
CHARACTER (LEN=*) :: runname
REAL :: time
REAL :: var(nx,ny,nz)
CHARACTER (LEN=80) :: timsnd
CHARACTER (LEN=132) :: vfnam
INTEGER :: ierr, istat, nunit, lvfnam, lvar
INTEGER :: tmstrln
INTEGER :: lrunnam,ldirnam
INTEGER :: nx_in,ny_in,nz_in
INTEGER :: i,j,k
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
CALL cvttsnd( time, timsnd, tmstrln )
CALL gtlfnkey( runname, lrunnam )
lvar = LEN(varnam)
vfnam = runname(1:lrunnam)//'.'//varnam(1:lvar)//timsnd(1:tmstrln)
CALL getunit( nunit)
CALL asnctl ('NEWLOCAL', 1, ierr)
CALL asnfile(vfnam, '-F f77 -N ieee', ierr)
WRITE(6,'(''READING ARRAY '',A,'' FROM FILE '',A)') &
varnam(1:lvar),vfnam
OPEN(UNIT=nunit,FILE=trim(vfnam),STATUS='old', &
FORM='unformatted',IOSTAT= istat )
READ(nunit) varnam_in
READ(nunit) nx_in,ny_in,nz_in
IF(nx_in /= nx .OR. ny_in /= ny .OR. nz_in /= nz) THEN
WRITE(6,'(a,/a,a,/a,3I5,/a,3I5)') &
'Warning in subroutine READVAR: Dimensions of data file ', &
vfnam,' do not agree with the expected dimensions.', &
'nx, ny and nz in the data are ',nx_in,ny_in,nz_in, &
'nx, ny and nz expected are ',nx,ny,nz
END IF
DO k=1,nz
READ (nunit) ((var(i,j,k),i=1,nx),j=1,ny)
END DO
CLOSE(UNIT=nunit)
CALL retunit(nunit)
WRITE(6,'(''ARRAY '',A,'' READ FROM FILE '',A)') &
varnam_in,vfnam
RETURN
END SUBROUTINE readvar
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE WRTVAR1 ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE wrtvar1(nx,ny,nz,var, varid,varname,varunits, & 13,5
time,runname,dirname, STATUS)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
! To write an array 'var' into a file in binary format.
!
! AUTHOR: Ming Xue
!
! MODIFICATION HISTORY:
!
! 1998/03/17 (R. Carpenter)
! Introduced I/O status variable and array name and units.
!
! INPUT:
!
! nx, ny, nz Dimensions of input array 'var'.
! var Array to be written out.
! varid String of length 6 (padded with _ as necessary) to
! designate the input array. (E.g., 'w_____')
! varname String describing the field (e.g., 'Vertical Velocity')
! varunits String describing the units (e.g., 'm/s')
! time The model time in seconds.
! runname Run name
! dirname Output directory (use '.' for current directory)
!
! OUTPUT:
!
! status Exit status (0=okay, 1=warning, 2=error)
!
! I/O:
!
! An unformatted binary file named dirname/runname.varid{time} is created:
! nx, ny, nz
! var
! varname
! varunits
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
! Misc. local variables:
!
!-----------------------------------------------------------------------
!
INTEGER :: nx,ny,nz, STATUS
REAL :: var(nx,ny,nz)
CHARACTER (LEN=*) :: varunits
CHARACTER (LEN=*) :: varname
CHARACTER (LEN=40) :: varunits1
CHARACTER (LEN=40) :: varname1
CHARACTER (LEN=6) :: varid
CHARACTER (LEN=132) :: tmpdirnam
CHARACTER (LEN=132) :: vfnam
INTEGER :: ierr, IOSTAT, nunit
REAL :: time
CHARACTER (LEN=*) :: dirname
CHARACTER (LEN=*) :: runname
CHARACTER (LEN=80):: timestring
INTEGER :: lfnkey
INTEGER :: ltimestring, lvfnam
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
IOSTAT = 0
STATUS = 0
CALL gtlfnkey( runname, lfnkey )
tmpdirnam = dirname
IF ( LEN(tmpdirnam) == 0 .OR. tmpdirnam == ' ' ) THEN
tmpdirnam = '.'
END IF
WRITE (vfnam,'(4A,A6)') trim(tmpdirnam),'/',runname(1:lfnkey),'.',varid
CALL cvttsnd( time, timestring, ltimestring )
lvfnam = len_trim(vfnam)
vfnam(lvfnam+1:lvfnam+ltimestring) = timestring(1:ltimestring)
PRINT *, 'WRTVAR1: Writing array to ', trim(vfnam)
CALL getunit( nunit)
CALL asnctl ('NEWLOCAL', 1, ierr)
CALL asnfile(vfnam, '-F f77 -N ieee', ierr)
OPEN(UNIT=nunit,FILE=vfnam,STATUS='unknown',FORM='unformatted', &
ERR=9000, IOSTAT=IOSTAT)
WRITE(nunit, ERR=9000, IOSTAT=IOSTAT) nx,ny,nz
WRITE(nunit, ERR=9000, IOSTAT=IOSTAT) var
varname1 = varname
varunits1= varunits
WRITE(nunit, ERR=9000, IOSTAT=IOSTAT) varname1
WRITE(nunit, ERR=9000, IOSTAT=IOSTAT) varunits1
CLOSE(UNIT=nunit)
CALL retunit(nunit)
RETURN ! Normal return
! I/O error handling
! Note that IOSTAT < 0 should not occur in this subroutine.
9000 CONTINUE ! I/O errors
CLOSE(UNIT=nunit)
CALL retunit(nunit)
IF (IOSTAT < 0) THEN
STATUS = 2
PRINT *, &
'WRTVAR1: I/O ERRORS OCCURRED ', &
'(possible end of record or file): ', &
IOSTAT, STATUS
ELSE IF (IOSTAT > 0) THEN
STATUS = 2
PRINT *, 'WRTVAR1: UNRECOVERABLE I/O ERRORS OCCURRED: ', &
IOSTAT, STATUS
END IF
RETURN
END SUBROUTINE wrtvar1
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE READVAR1 ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
SUBROUTINE readvar1(nx,ny,nz,var, varid,varname,varunits, &,5
time,runname,dirname, STATUS)
!-----------------------------------------------------------------------
!
! PURPOSE:
! To read in array 'var' from a file.
!
! AUTHOR: Ming Xue
!
! MODIFICATION HISTORY:
!
! 2/8/1998 (M.Xue)
! Added check on the arrays bounds read in from data.
!
! 1998/03/17 (R. Carpenter)
! Introduced I/O status variable and array name and units.
!
! INPUT:
!
! nx, ny, nz Dimensions of input array 'var'.
! varid String of length 6 (padded with _ as necessary) to
! designate the input array. (E.g., 'w_____')
! time The model time in seconds.
! runname Run name
! dirname Input directory (use '.' for current directory)
!
! OUTPUT:
!
! var Array to be written out.
! varname String describing the field (e.g., 'Vertical Velocity')
! varunits String describing the units (e.g., 'm/s')
! status Exit status (0=okay, 1=warning, 2=error)
!
! I/O:
!
! An unformatted binary file named dirname/runname.varid{time} is read:
! nx, ny, nz
! var
! varname
! varunits
!
!-----------------------------------------------------------------------
IMPLICIT NONE
INTEGER :: nx,ny,nz, STATUS
CHARACTER (LEN=6) :: varid
CHARACTER (LEN=*) :: varname, varunits
CHARACTER (LEN=40) :: varname1, varunits1
CHARACTER (LEN=*) :: runname, dirname
REAL :: time
REAL :: var(nx,ny,nz)
CHARACTER (LEN=132) :: tmpdirnam
CHARACTER (LEN=132) :: vfnam
CHARACTER (LEN=80 ) :: timestring
INTEGER :: ierr, IOSTAT, nunit
INTEGER :: nx_in,ny_in,nz_in, lfnkey
INTEGER :: lvfnam, ltimestring
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
STATUS = 0
CALL gtlfnkey( runname, lfnkey )
tmpdirnam = dirname
IF ( LEN(tmpdirnam) == 0 .OR. tmpdirnam == ' ' ) THEN
tmpdirnam = '.'
END IF
WRITE (vfnam,'(4A,A6)') trim(tmpdirnam),'/',runname(1:lfnkey),'.',varid
CALL cvttsnd( time, timestring, ltimestring )
lvfnam = len_trim(vfnam)
vfnam(lvfnam+1:lvfnam+ltimestring) = timestring(1:ltimestring)
PRINT *, 'READVAR1: Reading array from ',trim(vfnam)
CALL getunit( nunit)
CALL asnctl ('NEWLOCAL', 1, ierr)
CALL asnfile(vfnam, '-F f77 -N ieee', ierr)
OPEN(UNIT=nunit,FILE=trim(vfnam),STATUS='old', FORM='unformatted', &
ERR=9000, IOSTAT=IOSTAT)
READ(nunit, ERR=9000, END=9000, IOSTAT=IOSTAT) nx_in,ny_in,nz_in
IF(nx_in /= nx .OR. ny_in /= ny .OR. nz_in /= nz) THEN
WRITE(6,'(a,/a,a,/a,3I5,/a,3I5)') &
'Warning in subroutine READVAR1: Dimensions of data file ', &
vfnam,' do not agree with the expected dimensions.', &
'nx, ny and nz in the data are ',nx_in,ny_in,nz_in, &
'nx, ny and nz expected are ',nx,ny,nz
END IF
READ(nunit, ERR=9000, END=9000, IOSTAT=IOSTAT) var
READ(nunit, ERR=9000,END=9000,IOSTAT=IOSTAT) varname1
READ(nunit, ERR=9000,END=9000,IOSTAT=IOSTAT) varunits1
varname = varname1
varunits = varunits1
CLOSE(UNIT=nunit)
CALL retunit(nunit)
RETURN ! Normal return
! I/O error handling
! Note that warning (status=1) (e.g., end of record or file) is implemented as
! error (status=2) because of ambiguities in IOSTAT.
9000 CONTINUE ! I/O errors
CLOSE(UNIT=nunit)
CALL retunit(nunit)
IF (IOSTAT < 0) THEN
STATUS = 2
PRINT *, 'READVAR1: I/O ERRORS OCCURRED ', &
'(possible end of record or file): ',IOSTAT, STATUS
ELSE IF (IOSTAT > 0) THEN
STATUS = 2
PRINT *, 'READVAR1: UNRECOVERABLE I/O ERRORS OCCURRED: ', &
IOSTAT,STATUS
END IF
RETURN
END SUBROUTINE readvar1
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GET_INPUT_FILE_NAMES ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE get_input_file_names(hinfmt, & 5,2
grdbasfn,hisfile,nhisfile)
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 7/17/2000.
!
! MODIFICATION HISTORY:
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! hinfmt
!
! OUTPUT:
!
! grdbasfn,hisfile,nhisfile
!
!-----------------------------------------------------------------------
!
! Variable Declarations.
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
!
INTEGER :: hinfmt ! input data format
INTEGER :: hdmpinopt
CHARACTER (LEN=132) :: hdmpftrailer
CHARACTER (LEN=132) :: hdmpfheader
REAL :: tintv_dmpin, tbgn_dmpin, tend_dmpin
INTEGER :: nhisfile_max,nhisfile
PARAMETER(nhisfile_max=200)
CHARACTER (LEN=132) :: grdbasfn
CHARACTER (LEN=132) :: hisfile(nhisfile_max)
NAMELIST /history_data/ hinfmt, nhisfile, grdbasfn, hisfile, &
hdmpinopt,hdmpfheader,hdmpftrailer, &
tintv_dmpin, tbgn_dmpin, tend_dmpin
INTEGER :: lengbf,nf,length
!
!-----------------------------------------------------------------------
!
! Get the names of the input data files.
!
!-----------------------------------------------------------------------
!
READ(5,history_data,ERR=100)
WRITE(6,'(a)')'Namelist history_data was successfully read.'
WRITE(6,'(a,i3)') 'Input hdmpinopt=', hdmpinopt
WRITE(6,'(a,i3)') 'Input hinfmt =', hinfmt
WRITE(6,'(a,i3)') 'Input nhisfile =', nhisfile
IF( hdmpinopt == 1 ) THEN
CALL gthinfns(hdmpfheader,hdmpftrailer,hinfmt, &
tintv_dmpin, tbgn_dmpin, tend_dmpin, &
grdbasfn,hisfile,nhisfile)
END IF
IF( nhisfile > nhisfile_max ) THEN
WRITE(6,'(/a,a,I5,/a/)') &
'The number of history files to be processed', &
' exceeded the maximum number ',nhisfile_max, &
'Please increase the size of array hisfile.'
END IF
lengbf =len_trim( grdbasfn )
WRITE(6,'(/a,a)')'The grid base-state file name is ', &
grdbasfn(1:lengbf)
DO nf=1,nhisfile
WRITE(6,'(a,i3,a,a)') &
'History file No.',nf,' is ',trim(hisfile(nf))
END DO
!
GO TO 10
100 WRITE(6,'(a,a)') 'Error reading NAMELIST file. ', &
'Job stopped in GET_INPUT_FILE_NAMES.'
CALL arpsstop('arpsstop called from iGET_INPUT_FILE_NAMES error'// &
'reading namelist',1)
10 CONTINUE
RETURN
END SUBROUTINE get_input_file_names
!
!##################################################################
!##################################################################
!###### ######
!###### SUBROUTINE GTHINFNS ######
!###### ######
!###### Developed by ######
!###### Center for Analysis and Prediction of Storms ######
!###### University of Oklahoma ######
!###### ######
!##################################################################
!##################################################################
!
SUBROUTINE gthinfns(hdmpfheader,hdmpftrailer,hinfmt, & 1,3
tintv_dmpin, tbgn_dmpin, tend_dmpin, &
grdbasfn,hisfile,nhisfile)
!
!-----------------------------------------------------------------------
!
! PURPOSE:
!
! Retrun a list of history file names given the start, end times
! and time interval, as well as the name header and trailing string.
!
!-----------------------------------------------------------------------
!
! AUTHOR: Ming Xue
! 4/7/2000
!
! MODIFICATION HISTORY:
!
!-----------------------------------------------------------------------
!
! INPUT:
!
! hdmpfheader,hdmpftrailer,hinfmt
! tintv_dmpin,tbgn_dmpin,tend_dmpin
!
! OUTPUT:
!
! grdbasfn,hisfile,nhisfile
!
!-----------------------------------------------------------------------
!
! Variable Declarations:
!
!-----------------------------------------------------------------------
!
IMPLICIT NONE
INTEGER :: hinfmt
CHARACTER (LEN=132) :: hdmpfheader,hdmpftrailer
REAL :: tintv_dmpin, tbgn_dmpin, tend_dmpin
INTEGER :: nhisfile_max,nhisfile
CHARACTER (LEN=132) :: grdbasfn
PARAMETER(nhisfile_max=200)
CHARACTER (LEN=132) :: hisfile(nhisfile_max)
REAL :: time
CHARACTER (LEN=80) :: timsnd
INTEGER :: tmstrln
INTEGER :: lheader,ltrailer
INTEGER :: length,i
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
! Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
lheader = len_trim(hdmpfheader)
IF( hinfmt == 1.OR.hinfmt == 0) THEN
grdbasfn = hdmpfheader(1:lheader)//'.bingrdbas'
ELSE IF( hinfmt == 2 ) THEN ! Formatted ASCII dump
grdbasfn = hdmpfheader(1:lheader)//'.ascgrdbas'
ELSE IF( hinfmt == 3 ) THEN ! HDF data dump
grdbasfn = hdmpfheader(1:lheader)//'.hdfgrdbas'
ELSE IF( hinfmt == 4 ) THEN ! Packed binary dump
grdbasfn = hdmpfheader(1:lheader)//'.pakgrdbas'
ELSE IF( hinfmt == 5 ) THEN ! Savi3D data dump
!-----------------------------------------------------------------------
! For Savi3D data dump, the grid and base state information is
! always written together with the other fields.
!-----------------------------------------------------------------------
ELSE IF( hinfmt == 6 ) THEN ! Binary with skipping
grdbasfn = hdmpfheader(1:lheader)//'.bn2grdbas'
ELSE IF( hinfmt == 7 ) THEN ! NetCDF format
grdbasfn = hdmpfheader(1:lheader)//'.netgrdbas'
ELSE IF( hinfmt == 8 ) THEN ! Packed NetCDF format
grdbasfn = hdmpfheader(1:lheader)//'.npkgrdbas'
ELSE IF( hinfmt == 9 ) THEN ! GrADS data dump
!-----------------------------------------------------------------------
! For GrADS data dump, the grid and base state information is
! always written together with the other fields.
!-----------------------------------------------------------------------
ELSE IF( hinfmt == 10 ) THEN ! GRIB format
grdbasfn = hdmpfheader(1:lheader)//'.grbgrdbas'
END IF
ltrailer= len_trim(hdmpftrailer)
IF( ltrailer > 0 ) THEN
length = LEN( grdbasfn)
CALL strlnth( grdbasfn, length )
grdbasfn(length+1:length+ltrailer)=hdmpftrailer(1:ltrailer)
END IF
time = tbgn_dmpin
nhisfile = 0
DO i=1,100000000
IF( time > tend_dmpin + 0.01*tintv_dmpin ) GO TO 105
nhisfile = nhisfile + 1
IF( nhisfile > nhisfile_max) THEN
WRITE(6,'(1x,a,i3,/1x,a,/1x,a,/1x,a)') &
'The number of history files to be processed exceeded ', &
nhisfile_max, &
' please reduce the number of files to be processed', &
'in a single job, or edit the program and re-set parameter', &
' nhisfile_max to a larger value. '
CALL arpsstop('arpsstop called from GTHINFNS number of files'// &
'exceeded',1)
END IF
CALL cvttsnd( time , timsnd, tmstrln )
hisfile(i)=' '
IF( hinfmt == 1.OR.hinfmt == 0) THEN ! Unformatted binary dump
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.bin'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 2 ) THEN ! Formatted ASCII dump
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.asc'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 3 ) THEN ! HDF data dump
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.hdf'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 4 ) THEN ! Packed binary dump
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.pak'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 5 ) THEN ! Data dump for Savi3D
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.svi'
ELSE IF( hinfmt == 6 ) THEN ! Binary with skipping
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.bn2'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 7 ) THEN ! NetCDF format
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.net'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 8 ) THEN ! Packed NetCDF format
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.npk'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 9 ) THEN ! Data dump for GrADS
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.gad'
ELSE IF( hinfmt == 10 ) THEN ! Data dump for GrADS
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.grb'//timsnd(1:tmstrln)
ELSE IF( hinfmt == 11 ) THEN ! Data dump for Vis5D
hisfile(i)(1:lheader+4+tmstrln) &
=hdmpfheader(1:lheader)//'.v5d'//timsnd(1:tmstrln)
END IF
IF( ltrailer > 0 ) THEN
length = len_trim( hisfile(i) )
hisfile(i)(length+1:length+ltrailer) &
=hdmpftrailer(1:ltrailer)
END IF
time = time + tintv_dmpin
END DO
105 CONTINUE
RETURN
END SUBROUTINE gthinfns