!      SUBROUTINE PROF(NHB,LRSTRT,ITAG,LCLAS1)V
!      SUBROUTINE PROF(ITAG,LCLAS1)
      SUBROUTINE PROF_NMM(filename,prefilename,startdate,
     +                    ITAG,INCR)
!
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .
!   SUBROUTINE:  PROF        PROFILE SOUNDINGS
!   PRGRMMR: BLACK           ORG: W/NP22     DATE: 99-04-22
!
! ABSTRACT:  THIS ROUTINE GENERATES THE RAW PROFILE SOUNDING
!            OUTPUT FILES FROM THE FORECAST RESTRT FILE AND
!            AUXILIARY FILES
!
! PROGRAM HISTORY LOG:
!   99-04-22  T BLACK - ORIGINATOR
!   02-07-01  G MANIKIN - FIXED PROBLEM WITH DHCNVC AND DHRAIN
!                          COMPUTATIONS - SEE COMMENTS BELOW
!   03-04-01  M PYLE - BEGAN CONVERTING FOR WRF
!   04-05-26  M PYLE - MADE CHANGES FOR WRF-NMM
!   06-08-01  B ZHOU - ADAPTED FOR SREF WRF-NMM
!   07-08-07  J. Du & B. Zhou -  a new prefilename is defined 
!                for previous forecast file in order to calculate 
!                precip rate during INCHOUR interval
!   08-10-06  J. Du & B. Zhou - changed from IKJ to IJK 
!   12-06-14  M PYLE - non-MPIIO version for WCOSS
!
! USAGE:  CALL PROF FROM PROGRAM POST0
!
!   INPUT ARGUMENT LIST:
!     NHB    - THE UNIT NUMBER FOR READING THE NHB FILE
!     LRSTRT - THE UNIT NUMBER FOR READING THE RESTRT FILE
!     ITAG   - THE FORECAST HOUR WE ARE DEALING WITH
!     LCLAS1 - THE UNIT NUMBER FOR WRITING THE PROFILE DATA
!
!   OUTPUT ARGUMENT LIST:
!     NONE
!
!   SUBPROGRAMS CALLED:
!     UNIQUE:
!
!-----------------------------------------------------------------------
!      use vrbls3d
!      use vrbls2d
!      use soil
!      use masks
       use kinds, only             : i_llong
!
      include 'wrf_io_flags.h'
C      include 'mpif.h'

!     INCLUDE "parmeta"
       INCLUDE "parmsoil"
!-----------------------------------------------------------------------
                             P A R A M E T E R
     & (NSTAT=1500,LCL1ML=15,LCL1SL=52
     &, D608=0.608)
C-----------------------------------------------------------------------
C
C    PARMS FOR HOURLY PROFILER OUTPUT
C      NSTAT - MAX NUMBER OF STATIONS
C      NWORDM - DIMENSION OF OUTPUT ARRAY, MUST BE LARGE ENOUGH
C          TO HOLD ALL VARIABLES
C          (MAX NO MULTI-LAYER VARIABLES*LM + NO OF SINGLE LAYER VARS)
C      LCL1ML - NUMBER OF MULTI-LAYER VARIABLES OUTPUT FOR CLASS 1
C      LCL1SL - NUMBER OF SINGLE LAYER VARIABLES OUTPUT FOR CLASS 1
C
C------------------------------------------------------------------------
                             P A R A M E T E R
     & (ITB=76,JTB=134)
                             P A R A M E T E R
     & (A2=17.2693882,A3=273.16,A4=35.86,PQ0=379.90516,DTR=1.74532925E-2
     &, G=9.81,GI=1./G,RD=287.04,CP=1004.6,CAPA=RD/CP,RHCRIT=0.9999)
C
      PARAMETER (GAMMA=6.5/1000.,ZSL=0.0)
      PARAMETER (TAUCR=RD*GI*290.66,CONST=0.005*G/RD)
      PARAMETER (GORD=G/RD,DP=60.E2)
C------------------------------------------------------------------------
                             R E A L
     & STNLAT(NSTAT),STNLON(NSTAT)
                             I N T E G E R
     & IDSTN(NSTAT),IHINDX(NSTAT),JHINDX(NSTAT)
     &,             IVINDX(NSTAT),JVINDX(NSTAT)

        REAL, ALLOCATABLE:: DETA(:),RDETA(:),AETA(:),UL(:)
     &,RES(:),FIS(:),THS(:),HBOT(:)
     &,CFRACL(:),CFRACM(:),CFRACH(:),SNO(:)
     &,SOILTB(:),SFCEXC(:),SMSTAV(:),SMSTOT(:)
     &,Z0(:),CZEN(:),CZMEAN(:),SR(:)
     &,ACPREC(:),CUPREC(:),ACSNOW(:),ACSNOM(:)
     &,SSROFF(:),BGROFF(:),SFCSHX(:),SFCLHX(:)
     &,SUBSHX(:),SNOPCX(:),ASWIN(:),ASWOUT(:)
     &,ASWTOA(:),ALWIN(:),ALWOUT(:),ALWTOA(:)
     &,TSHLTR(:),TSHLTR_hold(:),QSHLTR(:),PSHLTR(:)
     &,TH10(:),Q10(:),U10(:),V10(:)
     &,TLMIN(:),TLMAX(:)
     &,SMC(:,:),CMC(:),STC(:,:),SH2O(:,:)
     &,VEGFRC(:),POTFLX(:),PSLP(:),PDSL1(:)
     &,EGRID2(:),SM(:),SICE(:)
     &,HBM2(:),FACTR(:)
     &,PTBL(:,:),TTBL(:,:),VEGFRA(:)
     &,T(:,:),Q(:,:),U(:,:),V(:,:),Q2(:,:)
     &,OMGALF(:,:),CWM(:,:),TRAIN(:,:),TCUCN(:,:)
     &,F_RAIN(:,:),F_ICE(:,:),CLDFRA(:,:)
     &,F_RIMEF(:,:)
     &,RSWTT(:,:),RLWTT(:,:),RTOP(:,:)
     &,HTM(:,:),OMGA(:,:)
     &,PRODAT(:),FPACK(:)
     &,STATPR(:),STACPR(:),STAEVP(:)
     &,STAPOT(:),STASHX(:),STASUB(:),STAPCX(:)
     &,STASWI(:),STASWO(:),STALWI(:),STALWO(:)
     &,STALWT(:),STASWT(:),STASNM(:),STASRF(:)
     &,STABRF(:),STASNO(:),DHCNVC(:,:),DHRAIN(:,:)
     &,STADHC(:),STADHR(:),CPRATE(:)
     &,ACPREC0(:),CUPREC0(:),SFCLHX0(:),POTFLX0(:)
     &,SFCSHX0(:),SUBSHX0(:),SNOPCX0(:),ASWIN0(:)
     &,ASWOUT0(:),ALWIN0(:),ALWOUT0(:),ALWTOA0(:)
     &,ASWTOA0(:),ACSNOW0(:),ACSNOM0(:),SSROFF0(:)
     &,BGROFF0(:)
     &,TCUCN0(:,:),TRAIN0(:,:)

        real, allocatable:: DUMSOIL(:),DUMSOIL3(:,:,:)
     &,DUM3D(:,:,:),DUM3D2(:,:,:),DUM3DIKJ(:,:,:)
     &,DUM3D3(:,:,:),DUMMY2(:,:),DUMMY(:,:)
     &,PD(:),PDS(:),GDLAT(:,:),GDLON(:,:)
     &,PMID(:,:),PINT(:,:)
     &,W(:,:),WH(:,:)

	
	real, allocatable:: CROT(:),SROT(:), SLDPTH2(:)
C------------------------------------------------------------------------
C
        integer, allocatable:: IDUM(:,:),LMH(:,:),IDUMMY(:,:)


                             I N T E G E R
     & IDAT(3),IDAT0(3),GDS(200)
       DATA SPVAL/-9999./
C
C------------------------------------------------------------------------
                             L O G I C A L
     & RUN,RESTRT,FRST,PRINT_DIAG
C------------------------------------------------------------------------
                             C H A R A C T E R
     & RSTFIL*98,RESTHR*4,LABEL*32,CISTAT*8,CIDSTN(NSTAT)*8
     &,FNAME*98,ENVAR*98,BLANK*4

C	new stuff
      character(len=31) :: VarName,varin
      character(len=256) :: fileName
      character(len=256) :: prefileName     !Jun Du
      character(len=256) :: newname
      integer :: Status
      character(len=19):: startdate,datestr,datestrold
      character SysDepInfo*80, titlestring*80

	real:: rinc(5)
	integer:: IDATE(8),JDATE(8),IDATENEW(8)

      character*132, allocatable :: datestr_all(:)
      character*132, allocatable :: varname_all(:)
      integer, allocatable       :: domainend_all(:,:)
      integer, allocatable       :: start_block(:)
      integer, allocatable       :: end_block(:)
      integer, allocatable       :: start_byte(:)
      integer, allocatable       :: end_byte(:)
      integer(kind=i_llong), allocatable           :: file_offset(:)
      integer this_offset, this_length

      character*2 cfhr

      real(kind=8) :: timef,btim,curtim

C------------------------------------------------------------------------
      DATA BLANK/'    '/
C------------------------------------------------------------------------
C***
C***  READ IN THE INFORMATION FILE ABOUT THE SOUNDINGS
C***

c	write(6,*) 'filename= ', filename
c	write(6,*) 'startedate= ', startdate

        btim=timef()


	datestr=startdate

      REWIND 19
C
      READ(19)NUMSTA,IDSTN,STNLAT,STNLON
     1,       IHINDX,JHINDX,IVINDX,JVINDX,CIDSTN
      WRITE(6,20)NUMSTA
   20 FORMAT('INIT:  NUMBER OF PROFILE STATIONS ',I5)
	if (ITAG .eq. 0) then
      WRITE(6,30)(IDSTN(N),STNLAT(N)/DTR,STNLON(N)/DTR
     1,               IHINDX(N),JHINDX(N),IVINDX(N),JVINDX(N)
     2,               CIDSTN(N),N=1,NUMSTA)
	endif
   30 FORMAT(2X,I6,2F8.2,4I8,4X,A8)

       FRST=.TRUE.

       if ( frst ) then
         PRINT_DIAG=.true.
         frst = .false.
          status=0
         CALL ext_int_ioinit(SysDepInfo,Status)
          print*,'CALLed ioinit', Status
	  write(6,*) 'filename early in PROF= ', trim(filename)
         CALL ext_int_open_for_read( trim(fileName), 0, 0, " ",
     &  DataHandle, Status)
          print*,'CALLed open for read', Status
       else
           Status = 0
       endif

       if ( Status /= 0 ) then
         print*,'error opening ',fileName, ' Status = ', Status ; stop
       endif

        write(0,*) 'call ext_int_get_dom_ti_char for title'
        call ext_int_get_dom_ti_char(DataHandle,'TITLE',
     &        titlestring,status)
        write(0,*) 'TITLE= ',trim(titlestring)
        write(6,*) 'TITLE= ',trim(titlestring)
        write(0,*) 'status from TITLE get: ', status

C Getting start time

      CALL ext_int_get_dom_ti_char(DataHandle
     1 ,'START_DATE',startdate, status )
        print*,'startdate from read= ',startdate


      jdate=0
      idate=0
      read(startdate,15)iyear,imn,iday,ihrst

        write(0,*) 'read iday: ', iday
         IDATE(1)=iyear
         IDATE(2)=imn
         IDATE(3)=iday
         IDATE(5)=ihrst
        write(0,*) 'defined idate(3): ', idate(3)

 15   format(i4,1x,i2,1x,i2,1x,i2)
      print*,'start yr mo day hr =',iyear,imn,iday,ihrst

! INSERT FROM INITPOST_NMM_BIN
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

C Getting tstart

      tstart=0.
cBZhou      call ext_int_get_dom_ti_real(DataHandle,'TSTART',tmp
cBZhou     + ,1,ioutcount,istatus)
cBZhou      tstart=tmp
cBZhou      print*,'status for getting TSTART= ',istatus
cBZhou      IF( abs(istatus-0) .GT. 1)THEN
cBZhou       PRINT*,'you do not have tstart in your WRF output,
cBZhou     + many grid navigation will be read in incorrectly, STOPPING'
cBZhou       STOP
cBZhou      END IF
        print*,'TSTART= ',TSTART

        if(tstart.gt. 1.E-2) THEN
          itstart=nint(tstart)
          rinc(2)=-1.0*tstart
          call w3movdat(rinc,idate,idatenew)
          print *," old start date= ",idate, tstart
          print *," new start date= ",idatenew, tstart
          do ier=1,8
             idate(ier)=idatenew(ier)
          enddo
        endif


!!! added from NAM version
C
C reset imn,iyear,iday,ihrst since they are packed into IDAT which
C is written into the profile output file!
C


          imn=IDATE(2)
          iday=IDATE(3)
          iyear=IDATE(1)
          ihrst=IDATE(5)
C


        call ext_int_get_dom_ti_integer(DataHandle,
     &   'WEST-EAST_GRID_DIMENSION',itmp
     + ,1,ioutcount,istatus)

        write(6,*) 'istatus for west-east dimension: ', istatus
        write(6,*) 'west-east dimension: ', itmp
        IM=itmp-1

        call ext_int_get_dom_ti_integer(DataHandle,
     &   'SOUTH-NORTH_GRID_DIMENSION',itmp
     + ,1,ioutcount,istatus)

        write(6,*) 'south-north dimension: ', itmp
        JM=itmp-1
        call ext_int_get_dom_ti_integer(DataHandle,
     &   'BOTTOM-TOP_GRID_DIMENSION',itmp
     + ,1,ioutcount,istatus)

        write(6,*) 'bottom-top dimension: ', itmp

        LM=itmp-1
        LP1=itmp


        write(6,*) 'LM, LP1: ', LM, LP1
        write(6,*) 'to big allocate block'

!  The end j row is going to be jend_2u for all variables except for V.
	JSTA_2L=1
	JEND_2U=JM
      JS=JSTA_2L
      JE=JEND_2U
      IF (JEND_2U.EQ.JM) THEN
       JEV=JEND_2U+1
      ELSE
       JEV=JEND_2U
      ENDIF

        call ext_int_get_dom_ti_real(DataHandle,'DX',tmp
     + ,1,ioutcount,istatus)
        write(6,*) 'istatus from DX get: ', istatus
        write(6,*) 'tmp: ', tmp
        dxval=nint(tmp*1000.) ! E-grid dlamda in degree
        write(6,*) 'dxval= ', dxval

        call ext_int_get_dom_ti_real(DataHandle,'DY',tmp
     + ,1,ioutcount,istatus)
        dyval=nint(1000.*tmp)
        write(6,*) 'dyval= ', dyval

        call ext_int_get_dom_ti_real(DataHandle,'DT',tmp
     + ,1,ioutcount,istatus)
        DT=tmp
        write(6,*) 'DT= ', DT

        call ext_int_get_dom_ti_real(DataHandle,'CEN_LAT',tmp
     + ,1,ioutcount,istatus)
!        tmp=50.
	TPH0D=tmp
        cenlat=nint(1000.*tmp)
        write(6,*) 'cenlat= ', cenlat

        call ext_int_get_dom_ti_real(DataHandle,'CEN_LON',tmp
     + ,1,ioutcount,istatus)
!        tmp=-111.
	TLM0D=tmp
        cenlon=nint(1000.*tmp)
        write(6,*) 'cenlon= ', cenlon

        call ext_int_get_dom_ti_real(DataHandle,'TRUELAT1',tmp
     + ,1,ioutcount,istatus)
        truelat1=nint(1000.*tmp)
!        write(6,*) 'truelat1= ', truelat1

        call ext_int_get_dom_ti_real(DataHandle,'TRUELAT2',tmp
     + ,1,ioutcount,istatus)
        truelat2=nint(1000.*tmp)
!        write(6,*) 'truelat2= ', truelat2

        call ext_int_get_dom_ti_integer(DataHandle,'MAP_PROJ',itmp
     + ,1,ioutcount,istatus)
        maptype=itmp
        write(6,*) 'maptype is ', maptype

! former parameter statements
        NWORDM=(LCL1ML+1)*LM+2*LCL1SL
        LRECPR=4*(8+9+LCL1ML*LM+LCL1SL)

        write(6,*) 'NWORDM, LRECPR: ', NWORDM, LRECPR
! former parameter statements

       ALLOCATE( DETA(LM),RDETA(LM),AETA(LM),UL(2*LM)
     &,RES(NUMSTA),FIS(NUMSTA),THS(NUMSTA),HBOT(NUMSTA)
     &,CFRACL(NUMSTA),CFRACM(NUMSTA),CFRACH(NUMSTA),SNO(NUMSTA)
     &,SOILTB(NUMSTA),SFCEXC(NUMSTA),SMSTAV(NUMSTA),SMSTOT(NUMSTA)
     &,Z0(NUMSTA),CZEN(NUMSTA),CZMEAN(NUMSTA),SR(NUMSTA)
     &,ACPREC(NUMSTA),CUPREC(NUMSTA),ACSNOW(NUMSTA),ACSNOM(NUMSTA)
     &,SSROFF(NUMSTA),BGROFF(NUMSTA),SFCSHX(NUMSTA),SFCLHX(NUMSTA)
     &,SUBSHX(NUMSTA),SNOPCX(NUMSTA),ASWIN(NUMSTA),ASWOUT(NUMSTA)
     &,ASWTOA(NUMSTA),ALWIN(NUMSTA),ALWOUT(NUMSTA),ALWTOA(NUMSTA)
     &,TSHLTR(NUMSTA),TSHLTR_hold(NUMSTA),QSHLTR(NUMSTA),PSHLTR(NUMSTA)
     &,TH10(NUMSTA),Q10(NUMSTA),U10(NUMSTA),V10(NUMSTA)
     &,TLMIN(NUMSTA),TLMAX(NUMSTA)
     &,SMC(NUMSTA,NSOIL),CMC(NUMSTA),STC(NUMSTA,NSOIL)
     &,SH2O(NUMSTA,NSOIL),SLDPTH2(NSOIL)
     &,VEGFRC(NUMSTA),POTFLX(NUMSTA),PSLP(NUMSTA),PDSL1(NUMSTA)
     &,EGRID2(NUMSTA),SM(NUMSTA),SICE(NUMSTA)
     &,HBM2(NUMSTA),FACTR(NUMSTA)
     &,PTBL(ITB,JTB),TTBL(JTB,ITB),VEGFRA(NUMSTA)
     &,T(NUMSTA,LM),Q(NUMSTA,LM),U(NUMSTA,LM),V(NUMSTA,LM)
     &,Q2(NUMSTA,LM)
     &,OMGALF(NUMSTA,LM),CWM(NUMSTA,LM),TRAIN(NUMSTA,LM)
     &,F_RAIN(NUMSTA,LM),F_ICE(NUMSTA,LM),CLDFRA(NUMSTA,LM)
     &,TCUCN(NUMSTA,LM),F_RIMEF(NUMSTA,LM)
     &,RSWTT(NUMSTA,LM),RLWTT(NUMSTA,LM)
     &,RTOP(NUMSTA,LM)
     &,HTM(NUMSTA,LM),OMGA(NUMSTA,LM)
     &,PRODAT(NWORDM),FPACK(NWORDM)
     &,STATPR(NUMSTA),STACPR(NUMSTA),STAEVP(NUMSTA)
     &,STAPOT(NUMSTA),STASHX(NUMSTA),STASUB(NUMSTA),STAPCX(NUMSTA)
     &,STASWI(NUMSTA),STASWO(NUMSTA),STALWI(NUMSTA),STALWO(NUMSTA)
     &,STALWT(NUMSTA),STASWT(NUMSTA),STASNM(NUMSTA),STASRF(NUMSTA)
     &,STABRF(NUMSTA),STASNO(NUMSTA),DHCNVC(LM,NUMSTA)
     &,DHRAIN(LM,NUMSTA)
!?     &,STADHC(LM),STADHR(LM),CPRATE(LM)
     &,STADHC(LM),STADHR(LM),CPRATE(NUMSTA)
     &,ACPREC0(NUMSTA),CUPREC0(NUMSTA),SFCLHX0(NUMSTA),POTFLX0(NUMSTA)
     &,SFCSHX0(NUMSTA),SUBSHX0(NUMSTA),SNOPCX0(NUMSTA),ASWIN0(NUMSTA)
     &,ASWOUT0(NUMSTA),ALWIN0(NUMSTA),ALWOUT0(NUMSTA),ALWTOA0(NUMSTA)
     &,ASWTOA0(NUMSTA),ACSNOW0(NUMSTA),ACSNOM0(NUMSTA),SSROFF0(NUMSTA)
     &,BGROFF0(NUMSTA)
     &,TCUCN0(NUMSTA,LM),TRAIN0(NUMSTA,LM))

       ALLOCATE ( DUMSOIL(NSOIL),DUMSOIL3(IM,NSOIL,JM)
     &,DUM3D(IM,JM,LM),DUM3D2(IM,JM,LM+1),DUM3DIKJ(IM,LM,JM)
     &,DUM3D3(IM,JM,LM+1),DUMMY2(IM,JM),DUMMY(IM,JM)
     &,PD(NUMSTA),PDS(NUMSTA),GDLAT(IM,JM),GDLON(IM,JM)
     &,PMID(NUMSTA,LM),PINT(NUMSTA,LM+1)
     &,W(NUMSTA,LM+1),WH(NUMSTA,LM) )

      ALLOCATE(IDUM(IM,JM),LMH(IM,JM),IDUMMY(IM,JM))




      VarName='TOYVAR'
        write(6,*) 'call getVariableBikj for : ', VarName
        write(6,*) 'size(DUM3D): ', size(DUM3D,dim=1),
     &                              size(DUM3D,dim=2),
     &                              size(DUM3D,dim=3)
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)


      VarName='LU_INDEX'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,     &
     +  IM,1,JM,1,IM,JS,JE,1)


!      VarName='LMH'
!       write(6,*) 'call getIVariable for : ', VarName
!      call getIVariable(fileName,DateStr,DataHandle,VarName,IDUMMY,     &
!        IM,1,JM,1,IM,JS,JE,1)

!      VarName='LMV'
!       write(6,*) 'call getIVariable for : ', VarName
!      call getIVariable(fileName,DateStr,DataHandle,VarName,IDUMMY,     &
!        IM,1,JM,1,IM,JS,JE,1)

      VarName='HBM2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='SM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

!        if (ierr /= 0) then
!          print*,"Error reading ", VarName,"Assigned missing values"
!          SM=SPVAL
!        else
          DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    SM(N)=DUMMY2(I,J)
          END DO
!        endif

      VarName='SICE'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

!        if (ierr /= 0) then
!          print*,"Error reading ", VarName,"Assigned missing values"
!          SICE=SPVAL
!        else
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   SICE(N)=DUMMY2(I,J)
          END DO
!        endif

! HTM always 1 for sigma coord

          DO L = 1,LM
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            HTM(N,L) = 1.0
           END DO 
          END DO

      VarName='PD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   pd(N)=DUMMY2(I,J)
         ENDDO
	
	write(6,*) 'PD(20) (a): ', PD(20)

       curtim=timef()
       write(0,*) 'elapsed to PD read: ',1.e-3*(curtim-btim)




      VarName='FIS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   FIS(N)=DUMMY2(I,J)
          END DO

      VarName='RES'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,     &
     +  IM,1,JM,1,IM,JS,JE,1)

          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   RES(N)=DUMMY2(I,J)
          END DO


      VarName='T'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            T(N,L) = DUM3D(I,J,LL)
           END DO 
          END DO

      VarName='Q'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            Q(N,L) = DUM3D(I,J,LL)
           END DO 
          END DO

      VarName='U'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IVINDX(N)	
	    J=JVINDX(N)
            U(N,L) = DUM3D(I,J,LL)
           END DO 
          END DO

      VarName='V'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IVINDX(N)	
       	    J=JVINDX(N)
            V(N,L) = DUM3D(I,J,LL)
           END DO 
          END DO

       curtim=timef()
       write(0,*) 'elapsed past V read: ',1.e-3*(curtim-btim)

        varname='DX_NMM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

       varname='DETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='AETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='ETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='DETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='AETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

        varname='ETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

        varname='DY_NMM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

        varname='DLMD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

        varname='DPHD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

       varname='PDTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  1,1,1,1,1,1,1,1)
          PDTOP=DUMMY2(1,1)

	write(6,*) 'PDTOP: ', PDTOP

        varname='PT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'PT ', DUMMY2(1,1)
        PT=DUMMY2(1,1)

        varname='PBLH'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='MIXHT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='USTAR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

!used?
      varname='Z0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   Z0(N)=DUMMY2(I,J)
          END DO


      varname='THS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   THS(N)=DUMMY2(I,J)
          END DO

      VarName='QS'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='TWBS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='QWBS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='PREC' ! instantaneous precip rate?
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='APREC' ! instantaneous precip rate?
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ACPREC' ! accum total precip
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

           DO N=1,NUMSTA
	     I=IHINDX(N)	
	     J=JHINDX(N)
	     ACPREC(N)=DUMMY2(I,J)
	     if (IDSTN(N) .eq. 724050) then
	       write(6,*) 'ACPREC at DCA: ', ACPREC(N)
	     endif
           END DO

      varname='CUPREC' ! accum cumulus precip
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   CUPREC(N)=DUMMY2(I,J)
          END DO 

      VarName='SNO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SNO(N) = dummy ( i, j )
          END DO

      VarName='SI'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='CLDEFI'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TH10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   TH10(N)=DUMMY2(I,J)
          END DO

      varname='Q10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   Q10(N)=DUMMY2(I,J)
          END DO

      varname='PSHLTR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   PSHLTR(N)=DUMMY2(I,J)
          END DO

      varname='TSHLTR'   !potential Temp'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   TSHLTR_hold(N)=DUMMY2(I,J)
          END DO

      varname='QSHLTR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
CHC CONVERT FROM MIXING RATIO TO SPECIFIC HUMIDITY
	   QSHLTR(N) = DUMMY2(I,J) / (1.0 + DUMMY2(I,J))
          END DO

      VarName='Q2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            Q2(N,L) = DUM3D (I,J,LL)
           END DO
          END DO

      varname='AKHS_OUT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='AKMS_OUT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ALBASE'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ALBEDO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

!
!
! Very confusing story ...
!
! Retrieve hbot => It is named CNVBOT in the model and
! with HBOTS (shallow conv) and HBOTD (deep conv) represent
! the 3 sets of convective cloud base/top arrays tied to the frequency
! that history files are written.
!
! IN THE *MODEL*, array HBOT is similar to CNVBOT but is
! used in radiation and is tied to the frequency of radiation updates.
!
! For historical reasons model array CNVBOT is renamed HBOT
! and manipulated throughout the sounding post.
!
      VarName='CNVBOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           HBOT(N)= float(LM)-dummy2(i,j)+1.0
          END DO

      VarName='CNVTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CZEN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   CZEN(N)=DUMMY2(I,J)
          END DO

      varname='CZMEAN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   CZMEAN(N)=DUMMY2(I,J)
          END DO

      varname='EPSR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='GLAT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='GLON'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TSK'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='MXSNAL'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='RADOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='SIGT4'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TGROUND'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)



      varname='CWM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            CWM(N,L) = DUM3D(I,J,LL)
           END DO 
          END DO

      varname='F_ICE'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
            I=IHINDX(N)
            J=JHINDX(N)
            F_ICE(N,L) = DUM3D(I,J,LL)
           END DO
          END DO

      varname='F_RAIN'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
            I=IHINDX(N)
            J=JHINDX(N)
            F_RAIN(N,L) = DUM3D(I,J,LL)
           END DO
          END DO

      varname='F_RIMEF'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
          DO L = 1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
            I=IHINDX(N)
            J=JHINDX(N)
            F_RIMEF(N,L) = DUM3D(I,J,LL)
           END DO
          END DO

       curtim=timef()
       write(0,*) 'elapsed past F_RIMEF read: ',1.e-3*(curtim-btim)

!      varname='CLDFRA'
!          DO L = 1,LM
!           LL=LM-L+1
!           DO N=1,NUMSTA
!            I=IHINDX(N)
!            J=JHINDX(N)
!            CLDFRA(N,L) = DUM3D(I,J,LL)
!           END DO
!          END DO

      varname='SR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           SR(N)=DUMMY2(I,J)
          ENDDO

      varname='CFRACH'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           CFRACH(N)=DUMMY2(I,J)
          ENDDO

      varname='CFRACL'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           CFRACL(N)=DUMMY2(I,J)
          ENDDO

      varname='CFRACM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           CFRACM(N)=DUMMY2(I,J)
          ENDDO

      varname='ISLOPE'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='DZSOIL'
      call getVariableB(fileName,DateStr,DataHandle,VarName,SLDPTH2,
     +    1,1,1,NSOIL,1,1,1,NSOIL)

      VarName='SLDPTH'
      call getVariableB(fileName,DateStr,DataHandle,VarName,SLDPTH2,
     +     1,1,1,NSOIL,1,1,1,NSOIL)

      VarName='CMC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           CMC(N) = DUMMY2(i,j)
          END DO

        varname='GRNFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='PCTSNO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='SOILTB'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
  	   SOILTB(N)=DUMMY2(I,J)
          ENDDO

      varname='VEGFRC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
  	   VEGFRC(N)=DUMMY2(I,J)
          ENDDO

      VarName='SH2O'
        write(6,*) 'call getVariableB for : ', VarName
        write(6,*) 'know that NSOIL is: ', NSOIL
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)
          DO L = 1,NSOIL
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            SH2O(N,L)=DUM3D(I,J,NSOIL-L+1)
           ENDDO
          ENDDO

      VarName='SMC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)
          DO L = 1,NSOIL
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            SMC(N,L) = DUM3D(i,j,nsoil-l+1)
           END DO
          END DO

      VarName='STC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)
          DO L = 1,NSOIL
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            STC(N,L) = DUM3D(i,j,nsoil-l+1)
           END DO
          END DO

      VarName='PINT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D2,      &
     + IM,1,JM,LM+1,IM,JS,JE,LM+1)

          DO L=1,LP1
           DO N=1,NUMSTA 
            LL=LP1-L+1
	    I=IHINDX(N)	
	    J=JHINDX(N)
             PINT (N,L) = dum3d2 (I,J,LL) 
            END DO
           END DO 

        DO L=2,LP1
         DO N=1,NUMSTA
           PMID(N,L-1)=(PINT(N,L-1)+PINT(N,L))*0.5
           if (n.eq.1 .or. n.eq.100) then
             print *,'PMID:',n,l,PMID(N,L-1),PINT(N,L),PINT(N,L-1)
           endif
         end do
        end do

!!!!! CONSTRAIN Q TO A PARTICULAR RH VALUE, FOLLOWING CALRH OF WRFPOST

        do N=1,NUMSTA
         do L=1,LM

           QC= (PQ0/PMID(N,L)) *EXP(A2*(T(N,L)-A3)/(T(N,L)-A4))
           RH=Q(N,L)/QC

           IF (RH .gt. RHCRIT) THEN
           Q(N,L)=0.999*RHCRIT*QC
           ENDIF

         enddo
        enddo

!!!!! END RH CONSTRAIN

      VarName='W'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D2,      &
     + IM,1,JM,LM+1,IM,JS,JE,LM+1)
          DO L = 1,LP1
           LL=LP1-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
            W(N,L)=DUM3D2(I,J,LL)
          END DO
         END DO

      DO L = 1,LM
       DO N=1,NUMSTA
	I=IHINDX(N)	
	J=JHINDX(N)
        WH(N,L) = (W(N,L)+W(N,L+1))*0.5
       END DO
      END DO


      VarName='ACFRCV'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ACFRST'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SSROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SSROFF(N) = DUMMY(i,j)
          ENDDO

c reading UNDERGROUND RUNOFF
      VarName='BGROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           BGROFF(N) = dummy ( i, j )
          END DO

      VarName='RLWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RLWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ALWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWIN(N)=dummy(i,j)
          END DO

      VarName='ALWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWOUT(N)=dummy(i,j)
          END DO

      VarName='ALWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWTOA(N)=dummy(i,j)
          END DO

      VarName='RSWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RSWINC'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RSWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ASWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ASWIN(N)=dummy(i,j)
          END DO

      VarName='ASWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
      	   J=JHINDX(N)
           ASWOUT(N)=dummy(i,j)
          END DO

      VarName='ASWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ASWTOA(N)=dummy(i,j)
          END DO

      VarName='SFCSHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SFCSHX(N)=dummy(i,j)
          END DO

      VarName='SFCLHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SFCLHX(N)=dummy(i,j)
	   if (IDSTN(N) .eq. 724050) then
	     write(6,*) 'DCA value for SFCLHX: ', SFCLHX(N)
	   endif
          END DO

      VarName='SUBSHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SUBSHX(N)=dummy(i,j)
          END DO

      VarName='SNOPCX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SNOPCX(N)=dummy(i,j)
          END DO

      VarName='SFCUVX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='POTEVP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='POTFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           POTFLX(N)=dummy(i,j)
          END DO

      VarName='TLMIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           TLMIN(N)=dummy(i,j)
          END DO

      VarName='TLMAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           TLMAX(N)=dummy(i,j)
          END DO

! big insert

      VarName='T02_MIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='T02_MAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RH02_MIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RH02_MAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='RLWTT'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    RLWTT(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO

      varname='RSWTT'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    RSWTT(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO

      varname='TCUCN'
          TCUCN=SPVAL

      varname='TRAIN'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    TRAIN(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO
!

      VarName='NCFRCV'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='NCFRST'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='NPHS0'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName, 
     +  NPHS,1,1,1,1,1,1,1,1)
      write(6,*) 'NPHS= ', NPHS

      VarName='NPREC'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName, 
     +    IDUMMY(1,1),1,1,1,1,1,1,1,1)
          NPREC=IDUMMY(1,1)
      write(6,*) 'NPREC= ', NPREC

      VarName='NCLOD'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,  
     +    IDUMMY(1,1),1,1,1,1,1,1,1,1)
          NCLOD=IDUMMY(1,1)
      write(6,*) 'NCLOD= ', NCLOD

      VarName='NHEAT'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NHEAT,  
     +    1,1,1,1,1,1,1,1)
      write(6,*) 'NHEAT= ', NHEAT

      VarName='NRDLW'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NRDLW, 
     +    1,1,1,1,1,1,1,1)
      write(6,*) 'NRDLW= ', NRDLW

      VarName='NRDSW'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NRDSW, 
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'NRDSW= ', NRDSW

      VarName='NSRFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NSRFC, 
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'NSRFC= ', NSRFC

      VarName='AVRAIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          AVRAIN=DUMMY2(1,1)

      VarName='AVCNVC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          AVCNVC=DUMMY2(1,1)
!
      VarName='ACUTIM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ACUTIM=DUMMY2(1,1)
      write(6,*) 'ACUTIM= ', ACUTIM

      VarName='ARDLW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ARDLW=DUMMY2(1,1)
      write(6,*) 'ARDLW= ', ARDLW
                                                                                             
      VarName='ARDSW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ARDSW=DUMMY2(1,1)
        write(6,*) 'ARDSW= ', ARDSW
           
      VarName='ASRFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ASRFC=DUMMY2(1,1)

      VarName='APHTIM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          APHTIM=DUMMY2(1,1)

      VarName='MAX10MW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAX10U'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAX10V'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXUPDR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXDNDR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXHLCY'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXDBZ'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='LANDMASK'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='PSFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='TH2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

c
c reading 10 m wind
      VarName='U10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           U10(N) = dummy ( i, j )
          END DO

      VarName='V10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           V10(N) = dummy ( i, j )
          END DO
  
      VarName='SMSTAV'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SMSTAV (N) = dummy ( i, j )
          ENDDO

      VarName='SMSTOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SMSTOT (N) = dummy ( i, j )
          ENDDO

      VarName='SFROFF'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
!
      VarName='UDROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='IVGTYP'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY       &
     +  ,IM,1,JM,1,IM,JS,JE,1)

      VarName='ISLTYP'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY       &
     +  ,IM,1,JM,1,IM,JS,JE,1)

      VarName='VEGFRA'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
    	   J=JHINDX(N)
  	   VEGFRA(N)=DUMMY(I,J)
          ENDDO

      VarName='SFCEVP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='GRDFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SFCEXC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SFCEXC (N) = dummy ( i, j )
          ENDDO

      VarName='ACSNOW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ACSNOW (N) = dummy ( i, j )
          END DO
 
      VarName='ACSNOM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ACSNOM (N) = dummy ( i, j )
          END DO

      VarName='SNOW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='CANWAT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SST'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='WEASD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='THZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='QZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='UZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='VZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='QSFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='HTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='HBOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

!      VarName='CPRATE'
!          DO N=1,NUMSTA
!           I=IHINDX(N)
!           J=JHINDX(N)
!           CPRATE(N) = dummy ( i, j )
!          END DO

!!!! does closing the datahandle help?
!       close(iunit)
	call ext_int_ioclose(DataHandle,Status)
        print *,DataHandle,Status

	write(6,*) 'PD(20) (b): ', PD(20)


!!!! DONE GETTING

!      DO L=1,LM+1
      DO L=1,LM
      DO N=1,NUMSTA
	I=IHINDX(N)	
	J=JHINDX(N)
	if (N .eq. 60 .and. L .eq. 1) then
	write(6,*) 'G, RD, D608: ', G, RD, D608
	endif

!            IF(ABS(T(N,L)).GT.1.0E-3)
!     &        OMGA(N,L) = -W(N,L)*PINT(N,L)*G/
!     &                 (RD*T(N,L)*(1.+D608*Q(N,L)))

        WMID=(W(N,L)+W(N,L+1))/2.
            IF(ABS(T(N,L)).GT.1.0E-3)
     &        OMGA(N,L) = -WMID*PMID(N,L)*G/
     &                 (RD*T(N,L)*(1.+D608*Q(N,L)))

!	if (mod(N,60) .eq. 0 .and. mod(L,10) .eq. 0) then
!	write(6,*) 'N, L, W, PINT, T, Q, OMGA: ', N,L,W(N,L),PINT(N,L),
!     &  T(N,L),Q(N,L),OMGA(N,L)
!	endif
	
      END DO
      END DO

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! DONE INSERT

C------------------------------------------------------------------------
C***
C***  READ QUANTITIES NEEDED FROM THE NHB FILE
C***

!	write(6,*) 'set LMH to : ', LM
!	write(6,*) 'IM,jm: ', Im,jm
	do J=1,JM
	 do I=1,IM
	   LMH(I,J)=LM
	 enddo
	enddo

C
c     DO L=1,LM
c       READ(NHB)
c     ENDDO
C
c     READ(NHB)DY,CPGFV,EN,ENT,R,PT,TDDAMP
c    1,        F4D,F4Q,EF4T,DETA,RDETA,AETA


	write(6,*) 'TLM0D, TPH0D: ', TLM0D, TPH0D

C
C	GET ROTATION ANGLES FOR WIND
C
	
!	write(6,*) 'hardwired DT: '
!	DT=18

      NTSPH=INT(3600./DT+0.50)

	write(6,*) 'rot angles defined, and ntsph= ', ntsph
C
C------------------------------------------------------------------------
C***
c       DO L = 1, LM
c       DO N=1,NUMSTA
!	  Q2(N,L)=0.
c       ENDDO
c       ENDDO

c       DO L=1,LM
c       DO N=1,NUMSTA
!	  TRAIN(N,L)=-9999.
!	  TCUCN(N,L)=-9999.
c       ENDDO
c     ENDDO
C
!      DO N=1,NUMSTA
!	CZEN(N)=-9999.
!      ENDDO
C
c     DO N=1,NUMSTA
c       I=IHINDX(N)	
c       J=JHINDX(N)
c       HBOT(N)=-9999.
c     ENDDO
C
c     DO N=1,NUMSTA
c      PSLP(N)=-9999.
c     ENDDO
c
c     DO N=1,NUMSTA
!	SOILTB(N)=-9999.
!	SMSTOT(N)=-9999.
!	SFCEXC(N)=-9999.
c     ENDDO
C
c     DO N=1,NUMSTA
!	CZMEAN(N)=-9999.
c       U00(N)=-9999.
!       SR(N)=-9999.
c     ENDDO
C
c     DO N=1,NUMSTA
!       SFCSHX(N)=-9999.
!       SFCLHX(N)=-9999.
!       SUBSHX(N)=-9999.
!       SNOPCX(N)=-9999.
c     ENDDO
C
c     DO N=1,NUMSTA
!        ASWIN(N)=-9999.
!        ASWOUT(N)=-9999.
!        ASWTOA(N)=-9999.
!        ALWOUT(N)=-9999.
!        ALWTOA(N)=-9999.
c     ENDDO

!      DO N=1,NUMSTA
!        TH10(N)=-9999.
!        Q10(N)=-9999.
!      ENDDO

C
C------------------------------------------------------------------------
!     DO N=1,NUMSTA
!       POTFLX(N)=-9999.
!       TLMIN(N)=-9999.
!       TLMAX(N)=-9999.
!     ENDDO
C------------------------------------------------------------------------
C***
C***  READ RADIATIVE TEMPERATURE TENDENCIES
C***
!      DO L=1,LM
!        DO N=1,NUMSTA
!         RSWTT(N,L)=-9999.
!         RLWTT(N,L)=-9999.
!        ENDDO
!      ENDDO
C
c     CLOSE(LRSTRT)
C------------------------------------------------------------------------
C***
C***  THE FORECAST HOUR
C***  
c     IFHR=NTSD/NTSPH
	write(6,*) 'past block of code'
cZhou	IFHR=ITAG
cZhou	write(6,*) 'IFHR: ', IFHR
C------------------------------------------------------------------------
      IF(ITAG.GT.0)THEN
        
C***
C***  GENERATE THE NAME OF THE PRECEDING RESTRT FILE
C***
        ITAG0=ITAG-INCR
        

	write(6,*) 'PD(20) (c): ', PD(20)

        RINC(1)=0.
        RINC(2)=float(ITAG0)
        RINC(3)=0.
        RINC(4)=0.
        RINC(5)=0.

	write(6,*) 'RINC(2): ', rinc(2)

cZhou        call w3movdat(rinc,idate,jdate)
cZhou        write(DateStrold,301) JDATE(1),JDATE(2),JDATE(3),JDATE(5)
cZhou 301    format(i4,'-',i2.2,'-',i2.2,'_',i2.2,':00:00')
cZhou	write(6,*) 'filename later in PROF: ', filename, '_END'
cZhou        len=lnblnk_(filename)
cZhou	write(6,*) 'LEN= ', LEN
cZhou        write(6,*) 'carried over part: ', filename(1:len-19)
cZhou        filename=filename(1:len-19)//DateStrold

cZhou        write(6,*) 'old filename is ', trim(filename)
cZhou        write(6,*) 'date for old file is: ', datestrold
     
         if(ITAG0.lt.10) then
          write(cfhr,'(1a,i1)') '0', ITAG0
         else
          write(cfhr,'(i2)')  ITAG0
         end if

         len=lnblnk(filename)          !call C code which has been compiled to object file


         filename=prefilename(1:len)
c        filename=filename(1:len-2)//cfhr
         write(*,*) 'Previous file=', trim(filename)

	write(6,*) 'PD(20) (ca): ', PD(20)
	
        CALL ext_int_ioinit(SysDepInfo,Status)
	write(6,*) 'PD(20) (caa): ', PD(20)
         CALL ext_int_open_for_read( trim(fileName), 0, 0, " ",
     &  DataHandle, Status)
	write(6,*) 'PD(20) (cab): ', PD(20)

       if ( Status /= 0 ) then
         print*,'error opening ',fileName, ' Status = ', Status 
	print*, 'quitting at open of EARLIER file'
	STOP
       endif
	write(6,*) 'PD(20) (cac): ', PD(20)

       curtim=timef()
       write(0,*) 'elapsed past full read first file : ',
     &        1.e-3*(curtim-btim)

C***  READ THE PREVIOUS RESTRT FILE
C***
C

!!!! INSERT FULL READ FROM ABOVE

	write(6,*) 'FROM EARLIER FILE'
        call ext_int_get_dom_ti_real(DataHandle,'DX',tmp
     + ,1,ioutcount,istatus)
        dxval=nint(tmp*1000.) ! E-grid dlamda in degree

        call ext_int_get_dom_ti_real(DataHandle,'DY',tmp
     + ,1,ioutcount,istatus)
        dyval=nint(1000.*tmp)

        call ext_int_get_dom_ti_real(DataHandle,'DT',tmp
     + ,1,ioutcount,istatus)
        DT=tmp

        call ext_int_get_dom_ti_real(DataHandle,'CEN_LAT',tmp
     + ,1,ioutcount,istatus)
!        tmp=50.0
	TPH0D=tmp
        cenlat=nint(1000.*tmp)

        call ext_int_get_dom_ti_real(DataHandle,'CEN_LON',tmp
     + ,1,ioutcount,istatus)
!        tmp=-111.0
	TLM0D=tmp
        cenlon=nint(1000.*tmp)

        call ext_int_get_dom_ti_real(DataHandle,'TRUELAT1',tmp
     + ,1,ioutcount,istatus)
        truelat1=nint(1000.*tmp)

        call ext_int_get_dom_ti_real(DataHandle,'TRUELAT2',tmp
     + ,1,ioutcount,istatus)
        truelat2=nint(1000.*tmp)

        call ext_int_get_dom_ti_integer(DataHandle,'MAP_PROJ',itmp
     + ,1,ioutcount,istatus)
        maptype=itmp


! dup block from above

      VarName='TOYVAR'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      VarName='LU_INDEX'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,     &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='HBM2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='SM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='SICE'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='PD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='FIS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='RES'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,     &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='T'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      VarName='Q'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      VarName='U'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      VarName='V'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

        varname='DX_NMM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

       varname='DETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='AETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='ETA1'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA1,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='DETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

       varname='AETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

        varname='ETA2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,ETA2,       &
     +  LM,1,1,1,LM,1,1,1)

        varname='DY_NMM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

        varname='DLMD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

        varname='DPHD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM0D,       &
     +  1,1,1,1,1,1,1,1)

       varname='PDTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  1,1,1,1,1,1,1,1)
          PDTOP=DUMMY2(1,1)

	write(6,*) 'PDTOP: ', PDTOP

        varname='PT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'PT ', DUMMY2(1,1)
        PT=DUMMY2(1,1)

        varname='PBLH'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='MIXHT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='USTAR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='Z0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='THS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='QS'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='TWBS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='QWBS'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='PREC' ! instantaneous precip rate?
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='APREC' ! instantaneous precip rate?
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ACPREC' ! accum total precip
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

           DO N=1,NUMSTA
	     I=IHINDX(N)	
	     J=JHINDX(N)
	     ACPREC0(N)=DUMMY2(I,J)
	     if (IDSTN(N) .eq. 724050) then
	       write(6,*) 'ACPREC0 at DCA: ', ACPREC0(N)
	     endif
           END DO

      varname='CUPREC' ! accum cumulus precip
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
	   CUPREC0(N)=DUMMY2(I,J)
          END DO 

      VarName='SNO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SI'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='CLDEFI'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TH10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='Q10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='PSHLTR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TSHLTR'   !potential Temp'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='QSHLTR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='Q2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      varname='AKHS_OUT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='AKMS_OUT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ALBASE'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ALBEDO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

!
!
! Very confusing story ...
!
! Retrieve hbot => It is named CNVBOT in the model and
! with HBOTS (shallow conv) and HBOTD (deep conv) represent
! the 3 sets of convective cloud base/top arrays tied to the frequency
! that history files are written.
!
! IN THE *MODEL*, array HBOT is similar to CNVBOT but is
! used in radiation and is tied to the frequency of radiation updates.
!
! For historical reasons model array CNVBOT is renamed HBOT
! and manipulated throughout the sounding post.
!
      VarName='CNVBOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='CNVTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CZEN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CZMEAN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='EPSR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='GLAT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='GLON'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TSK'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='MXSNAL'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='RADOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='SIGT4'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='TGROUND'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CWM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      varname='F_ICE'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      varname='F_RAIN'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      varname='F_RIMEF'
        write(6,*) 'call getVariableBikj for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)

      varname='SR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CFRACH'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CFRACL'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='CFRACM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='ISLOPE'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='DZSOIL'
      call getVariableB(fileName,DateStr,DataHandle,VarName,SLDPTH2,
     +    1,1,1,NSOIL,1,1,1,NSOIL)

      VarName='SLDPTH'
      call getVariableB(fileName,DateStr,DataHandle,VarName,SLDPTH2,
     +     1,1,1,NSOIL,1,1,1,NSOIL)

      VarName='CMC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='GRNFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

        varname='PCTSNO'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='SOILTB'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='VEGFRC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SH2O'
        write(6,*) 'call getVariableB for : ', VarName
        write(6,*) 'know that NSOIL is: ', NSOIL
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)

      VarName='SMC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)

      VarName='STC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableBikj(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,NSOIL)

      VarName='PINT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D2,      &
     + IM,1,JM,LM+1,IM,JS,JE,LM+1)

      VarName='W'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D2,      &
     + IM,1,JM,LM+1,IM,JS,JE,LM+1)

      VarName='ACFRCV'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ACFRST'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SSROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SSROFF0(N) = DUMMY(i,j)
          ENDDO

c reading UNDERGROUND RUNOFF
      VarName='BGROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           BGROFF0(N) = dummy ( i, j )
          END DO

      VarName='RLWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RLWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ALWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWIN0(N)=dummy(i,j)
          END DO

      VarName='ALWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWOUT0(N)=dummy(i,j)
          END DO

      VarName='ALWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ALWTOA0(N)=dummy(i,j)
          END DO

      VarName='RSWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RSWINC'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RSWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ASWIN'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ASWIN0(N)=dummy(i,j)
          END DO

      VarName='ASWOUT'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
      	   J=JHINDX(N)
           ASWOUT0(N)=dummy(i,j)
          END DO

      VarName='ASWTOA'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ASWTOA0(N)=dummy(i,j)
          END DO

      VarName='SFCSHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SFCSHX0(N)=dummy(i,j)
          END DO

      VarName='SFCLHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SFCLHX0(N)=dummy(i,j)
	   if (IDSTN(N) .eq. 724050) then
	     write(6,*) 'DCA value for SFCLHX0: ', SFCLHX0(N)
	   endif
          END DO

      VarName='SUBSHX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SUBSHX0(N)=dummy(i,j)
          END DO

      VarName='SNOPCX'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           SNOPCX0(N)=dummy(i,j)
          END DO

      VarName='SFCUVX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='POTEVP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='POTFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           POTFLX0(N)=dummy(i,j)
          END DO

      VarName='TLMIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           TLMIN(N)=dummy(i,j)
          END DO

      VarName='TLMAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
           I=IHINDX(N)
           J=JHINDX(N)
           TLMAX(N)=dummy(i,j)
          END DO

! big insert

      VarName='T02_MIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='T02_MAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RH02_MIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='RH02_MAX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      varname='RLWTT'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    RLWTT(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO

      varname='RSWTT'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    RSWTT(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO

      varname='TCUCN'
          TCUCN0=SPVAL

      varname='TRAIN'
       write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUM3D,      &
     + IM,1,JM,LM,IM,JS,JE,LM)
	  DO L=1,LM
           LL=LM-L+1
           DO N=1,NUMSTA
	    I=IHINDX(N)	
	    J=JHINDX(N)
	    TRAIN0(N,L)=DUM3D(I,J,LL)
	   ENDDO	
	  ENDDO
!

      VarName='NCFRCV'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='NCFRST'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY,      &
     +  IM,1,JM,1,IM,JS,JE,1)

      VarName='NPHS0'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName, 
     +  NPHS,1,1,1,1,1,1,1,1)
      write(6,*) 'NPHS= ', NPHS

      VarName='NPREC'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName, 
     +    IDUMMY(1,1),1,1,1,1,1,1,1,1)
          NPREC=IDUMMY(1,1)
      write(6,*) 'NPREC= ', NPREC

      VarName='NCLOD'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,  
     +    IDUMMY(1,1),1,1,1,1,1,1,1,1)
          NCLOD=IDUMMY(1,1)
      write(6,*) 'NCLOD= ', NCLOD

      VarName='NHEAT'
        write(6,*) 'call getIVariable for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NHEAT,  
     +    1,1,1,1,1,1,1,1)
      write(6,*) 'NHEAT= ', NHEAT

      VarName='NRDLW'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NRDLW, 
     +    1,1,1,1,1,1,1,1)
      write(6,*) 'NRDLW= ', NRDLW

      VarName='NRDSW'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NRDSW, 
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'NRDSW= ', NRDSW

      VarName='NSRFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,NSRFC, 
     +  1,1,1,1,1,1,1,1)
        write(6,*) 'NSRFC= ', NSRFC

      VarName='AVRAIN'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          AVRAIN0=DUMMY2(1,1)

      VarName='AVCNVC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          AVCNVC0=DUMMY2(1,1)
!
      VarName='ACUTIM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ACUTIM0=DUMMY2(1,1)

      VarName='ARDLW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ARDLW=DUMMY2(1,1)
      write(6,*) 'ARDLW= ', ARDLW
                                                                                             
      VarName='ARDSW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ARDSW=DUMMY2(1,1)
        write(6,*) 'ARDSW= ', ARDSW
           
      VarName='ASRFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          ASRFC=DUMMY2(1,1)

      VarName='APHTIM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY2(1,1), 
     +  1,1,1,1,1,1,1,1)
          APHTIM0=DUMMY2(1,1)

      VarName='MAX10MW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAX10U'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAX10V'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXUPDR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXDNDR'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXHLCY'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='MAXDBZ'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='LANDMASK'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='PSFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='TH2'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

c
c reading 10 m wind
      VarName='U10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='V10'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
  
      VarName='SMSTAV'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SMSTOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SFROFF'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
!
      VarName='UDROFF'
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='IVGTYP'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY       &
     +  ,IM,1,JM,1,IM,JS,JE,1)

      VarName='ISLTYP'
        write(6,*) 'call getVariableB for : ', VarName
      call getIVariableB(fileName,DateStr,DataHandle,VarName,IDUMMY       &
     +  ,IM,1,JM,1,IM,JS,JE,1)

      VarName='VEGFRA'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SFCEVP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='GRDFLX'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SFCEXC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='ACSNOW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ACSNOW0 (N) = dummy ( i, j )
          END DO
 
      VarName='ACSNOM'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)
          DO N=1,NUMSTA
	   I=IHINDX(N)	
	   J=JHINDX(N)
           ACSNOM0 (N) = dummy ( i, j )
          END DO

      VarName='SNOW'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='CANWAT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='SST'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='WEASD'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='THZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='QZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='UZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='VZ0'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='QSFC'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='HTOP'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

      VarName='HBOT'
        write(6,*) 'call getVariableB for : ', VarName
      call getVariableB(fileName,DateStr,DataHandle,VarName,DUMMY,      &
     + IM,1,JM,1,IM,JS,JE,1)

       curtim=timef()
       write(0,*) 'elapsed past full read second file : ',
     + 1.e-3*(curtim-btim)

! end dup block from above

	write(6,*) 'PD(20) (d): ', PD(20)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!! END INSERT FULL READ

c       DO N=1,NUMSTA
c         POTFLX0(N)=-9999.
c       ENDDO
C

      ENDIF

C
Cmp 	IDEALLY, WON'T NEED MANY MODS BELOW THIS POINT
C
C------------------------------------------------------------------------
C***
C***  ALL THE DATA IS NOW IN.
C***  CALCULATE CLOUD FRACTION AND CLOUD WATER/ICE ID NUMBER.
C***
C------------------------------------------------------------------------
      UTIM=1.
      US=1.
      CCLIMIT=1.E-3
      CLIMIT =1.E-20
C-----------------------------------------------------------------------
!$OMP parallel do 
      DO N=1,NUMSTA
        PDSL1(N)=PD(N)*RES(N)
      ENDDO
C----------------------------------------------------------------------
C----------------------------------------------------------------------
C***
C***  BEGIN THE PROFILE POSTING CODE.
C***
C----------------------------------------------------------------------
C----------------------------------------------------------------------
C***
C***  USE ZERO IN ACCUMULATION ARRAYS AT APPROPRIATE TIMES
C***
       IF(ITAG .eq. 0) THEN
C
C set TLMAX,TLMIN= Lowest Layer Temp if forecast hour (ITAG)=0
C
         DO N=1,NUMSTA
           TLMIN(N)=T(N,LM)
           TLMAX(N)=T(N,LM)
         ENDDO
	write(6,*) 'here (2)...should only see once'
C
C what would appropriate if test be here?
C
        DO N=1,NUMSTA
C
C*** ZERO ACCUMLATION ARRAYS.
C
          STATPR(N)=0.
          STACPR(N)=0.
          STAEVP(N)=0.
          STAPOT(N)=0.
          STASHX(N)=0.
          STASUB(N)=0.
          STAPCX(N)=0.
          STASWI(N)=0.
          STASWO(N)=0.
          STALWI(N)=0.
          STALWO(N)=0.
          STALWT(N)=0.
          STASWT(N)=0.
          STASNM(N)=0.
          STASNO(N)=0.
          STASRF(N)=0.
          STABRF(N)=0.
          DO L=1,LM
            DHCNVC(L,N)=0.
            DHRAIN(L,N)=0.
          ENDDO
        ENDDO
C
        GO TO 300
       ENDIF
C---------------------------------------------------------------------
C***
C***  WE MUST CHECK TO SEE IF WE ARE 1 HOUR AFTER ANY OF THE 
C***  ACCUMULATION BUCKETS HAVE BEEN EMPTIED.  IF WE ARE AT SUCH A 
C***  TIME THEN WE NEED TO SET TO ZERO THE VARIABLES USED TO HOLD
C***  THE PRECEDING HOUR'S VALUES.
C***
C---------------------------------------------------------------------
C
C
C 	At this point, accumulation buckets are a foreign concept in
C	the WRF model.  
C
C

c     TIME=(NTSD-1)*DT

	write(6,*) 'hardwired NPHS'

	if (DT .eq. 18) then
	NPHS=10
	elseif (DT .eq. 20) then
	NPHS=20
	elseif (DT .ge. 25 .and. DT .le. 28) then 
        NPHS=6
	elseif (DT .ge. 8.9 .and. DT .le. 11) then 
        NPHS=12
	elseif (DT .le. 8.9 ) then 
        NPHS=8
	endif

        IFCST=3600*ITAG
	TIME=IFCST
	NTSD=(TIME/DT)+1

	write(6,*) 'TIME, DT, NTSD: ', TIME, DT, NTSD

      RESET0=TIME-(NTSD/NPREC)*NPREC*DT
      RESET1=(NPHS-1)*DT+3600.

C
!?      IF(MOD(NTSD,NPREC).GE.NPHS.AND.RESET0.LE.RESET1)THEN

	write(6,*) 'RESET0, RESET1 for PREC: ', RESET0, RESET1
	write(6,*) 'NTSD, NPREC, NPHS: ', NTSD, NPREC, NPHS
	write(6,*) 'RESET0 < RESET1 for precip? : ', RESET0, RESET1
	write(6,*) 'MOD < NPHS for precip? : ', mod(NTSD,NPREC),NPHS

      IF(MOD(NTSD,NPREC).GE.NPHS.AND.RESET0.LE.RESET1)THEN
	write(6,*) 'settin STATPR, etc. to zero'
        DO N=1,NUMSTA
          STATPR(N)=0.
          STACPR(N)=0.
          STASNM(N)=0.
          STASNO(N)=0.
          STASRF(N)=0.
          STABRF(N)=0.
        ENDDO
      ELSE
	write(6,*) 'NOT resetting PRECIP accumulators'
        DO N=1,NUMSTA
          STATPR(N)=ACPREC0(N)*1.E3
          STACPR(N)=CUPREC0(N)*1.E3
!          STASNM(N)=ACSNOM0(N)*1.E3
!          STASNO(N)=ACSNOW0(N)*1.E3
!          STASRF(N)=SSROFF0(N)*1.E3
!          STABRF(N)=BGROFF0(N)*1.E3
          STASNM(N)=ACSNOM0(N)
          STASNO(N)=ACSNOW0(N)
          STASRF(N)=SSROFF0(N)
          STABRF(N)=BGROFF0(N)
        ENDDO
      ENDIF          
C
      RESET0=TIME-(NTSD/NRDSW)*NRDSW*DT
	
	write(6,*) 'TIME, DT: ', TIME, DT
	write(6,*) 'RESET0 < RESET1? : ', RESET0, RESET1
	write(6,*) 'MOD < NPHS?: ', MOD(NTSD,NRDSW),NPHS
	
      IF( MOD(NTSD,NRDSW) .GE. NPHS .AND. RESET0 .LE. RESET1)THEN
!      IF( MOD(NTSD,NRDSW) .LE. NPHS .AND. RESET0 .LE. RESET1)THEN
	write(6,*) 'resetting SW accumulators'
       DO N=1,NUMSTA
         STASWI(N)=0.
         STASWO(N)=0.
         STASWT(N)=0.
       ENDDO
      ELSE
	write(6,*) 'NOT resetting SW accumulators'
        DO N=1,NUMSTA
          STASWI(N)=ASWIN0(N)
          STASWO(N)=ASWOUT0(N)
          STASWT(N)=ASWTOA0(N)
        ENDDO
      ENDIF

!	write(6,*) 'STASWI(90): ', STASWI(90)
C
      RESET0=TIME-(NTSD/NRDLW)*NRDLW*DT
!	write(6,*) 'TIME, DT: ', TIME, DT
!	write(6,*) 'RESET0, RESET1 : ', RESET0, RESET1
!	write(6,*) 'NTSD, NRDLW, NPHS: ', NTSD, NRDLW, NPHS
      IF(MOD(NTSD,NRDLW).GE.NPHS.AND.RESET0.LE.RESET1)THEN
!      IF(MOD(NTSD,NRDLW).LE.NPHS.AND.RESET0.LE.RESET1)THEN
	write(6,*) 'resetting STALWI'
       DO N=1,NUMSTA
         STALWI(N)=0.
         STALWO(N)=0.
         STALWT(N)=0.
       ENDDO
      ELSE
	write(6,*) 'setting STALWI to ALWIN0'
        DO N=1,NUMSTA
          STALWI(N)=ALWIN0(N)
          STALWO(N)=ALWOUT0(N)
          STALWT(N)=-ALWTOA0(N)
        ENDDO
      ENDIF
C
      RESET0=TIME-(NTSD/NSRFC)*NSRFC*DT
      IF(MOD(NTSD,NSRFC).GE.NPHS.AND.RESET0.LE.RESET1)THEN
!      IF(MOD(NTSD,NSRFC).LE.NPHS.AND.RESET0.LE.RESET1)THEN
       DO N=1,NUMSTA
         STAEVP(N)=0.
         STAPOT(N)=0.
         STASHX(N)=0.
         STASUB(N)=0.
         STAPCX(N)=0.
       ENDDO
      ELSE
        DO N=1,NUMSTA
          STAEVP(N)=SFCLHX0(N)
          STAPOT(N)=POTFLX0(N)
          STASHX(N)=SFCSHX0(N)
          STASUB(N)=SUBSHX0(N)
          STAPCX(N)=SNOPCX0(N)
        ENDDO
      ENDIF
C
      RESET0=TIME-(NTSD/NHEAT)*NHEAT*DT
      IF(MOD(NTSD,NHEAT).GE.NCNVC.AND.RESET0.LE.RESET1)THEN
       DO N=1,NUMSTA
         DO L=1,LM
           DHCNVC(L,N)=0.
           DHRAIN(L,N)=0.
         ENDDO
       ENDDO
      ELSE
       DO N=1,NUMSTA
         DO L=1,LM
            DHCNVC(L,N)=TCUCN0(N,L)
            DHRAIN(L,N)=TRAIN0(N,L)
         ENDDO
       ENDDO
      ENDIF
 
C------------------------------------------------------------------
  300 CONTINUE
C------------------------------------------------------------------
C
C***  FOR ROTATION OF WINDS FROM E-GRID TO GEODETIC ORIENTATION
C***  WE NEED THE TWO QUANTITIES BELOW.
C
      SINPH0=SIN(TPH0D*DTR)
      COSPH0=COS(TPH0D*DTR)
C
C***  INITIAL CALCULATIONS/PREPARATIONS.  WE LOAD SEVERAL
C***  ARRAYS WITH PROFILE VARIABLES.
C
!$OMP parallel do
      DO N=1,NUMSTA
        IF(CZMEAN(N).GT.0.)THEN
          FACTR(N)=CZEN(N)/CZMEAN(N)
        ELSE
          FACTR(N)=0.
        ENDIF
      ENDDO
C
C***  ADJUST SHORTAVE TENDENCIES TO ACCOUNT FOR CHANGE OF SOLAR POSITION
C***  BETWEEN CALLS TO RADIATION
C
!$OMP parallel do
      DO L=1,LM
        DO N=1,NUMSTA
          RSWTT(N,L)=RSWTT(N,L)*FACTR(N)
        ENDDO
      ENDDO
C
C***  COMPUTE RTOP
C
!$OMP parallel do
      DO L=1,LM
        DO N=1,NUMSTA
!          APEL=PT+AETA(L)*PDSL1(N)
          APEL=PMID(N,L)
          RTOP(N,L)=RD*T(N,L)*(1.+0.608*Q(N,L))/APEL
        ENDDO
      ENDDO

	write(6,*) 'PD(20) (dd): ', PD(20)
C
C***  PDS IS SURFACE PRESSURE.
C
!$OMP parallel do 
	DO N=1,NUMSTA
        PDS(N)=PD(N)+PDTOP+PT
	if (PDS(N) .lt. 50000.) then
	write(6,*) 'N,PD,PDTOP,PDS: ', N,PD(N),PDTOP,PDS(N)
	endif
	ENDDO
C
C***  EGRID2 IS THE SURFACE TEMPERATURE.
C
!$OMP parallel do 
      DO N=1,NUMSTA
	TSHLTR(N)= TSHLTR_hold(N)*(PDS(N)*1.E-5)**CAPA  !convert potential to real temperature
        EGRID2(N)= THS(N)*(PDS(N)*1.E-5)**CAPA

	if (EGRID2(N) .lt. 200.) then
	write(6,*) 'THETA, PDS, CAPA, TSK: ', N,THS(N),
     &           PDS(N),CAPA,EGRID2(N)
	endif

        IF(ACPREC(N).LT.0.)ACPREC(N)=0.
        IF(CUPREC(N).LT.0.)CUPREC(N)=0.

!!! constrain surface RH

           QC=(PQ0/PSHLTR(N))*EXP(A2*(TSHLTR(N)-A3)/(TSHLTR(N)-A4))
           RH=QSHLTR(N)/QC
           IF (RH .gt. RHCRIT) THEN
c          write(6,*) 'reducing surface RH from: ', RH, ' at N: ', N
           QSHLTR(N)=0.999*RHCRIT*QC
	   ENDIF

      ENDDO
C
C** Compute PSLP using NMC reduction
C
      DO N=1,NUMSTA
       PSFC = PDS(N)
       ZSFC = FIS(N)*GI
       PSLP(N) = PSFC
C
C    COMPUTE LAYER TAU (VIRTUAL TEMP*RD/G).
       TVRT = T(N,LM)*(1.0+D608*Q(N,LM))
       TAU  = TVRT*RD*GI
C
C    COMPUTE TAU AT THE GROUND (Z=ZSFC) AND SEA LEVEL (Z=0)
C    ASSUMING A CONSTANT LAPSE RATE OF GAMMA=6.5DEG/KM.
       TVRSFC = TVRT + (ZSFC- ZSL)*GAMMA
       TAUSFC = TVRSFC*RD*GI
       TVRSL  = TVRT + (ZSFC- ZSL)*GAMMA
       TAUSL  = TVRSL*RD*GI
C
C    IF NEED BE APPLY SHEULL CORRECTION.
       IF ((TAUSL.GT.TAUCR).AND.(TAUSFC.LE.TAUCR)) THEN
          TAUSL=TAUCR
       ELSEIF ((TAUSL.GT.TAUCR).AND.(TAUSFC.GT.TAUCR)) THEN
          TAUSL = TAUCR-CONST*(TAUSFC-TAUCR)**2
       ENDIF
C
C    COMPUTE MEAN TAU.
       TAUAVG = 0.5*(TAUSL+TAUSFC)
C
C    COMPUTE SEA LEVEL PRESSURE.
       IF (FIS(N).GT.1.0) PSLP(N) = PSFC*EXP(ZSFC/TAUAVG)
c      print *,n,idstn(n),pslp(n),tvrt
       
      ENDDO
C
C***  SET CYCLE, DATE, AND FORECAST TIME.
C
c      IHR  =NTSD/NTSPH+0.5

!MAYBE      read(datestr,15)iyear,imn,iday,ihrst

         IDATE(2)=imn
         IDATE(3)=iday
         IDATE(1)=iyear
         IDATE(5)=ihrst

	IDAT(3)=IDATE(1)
	IDAT(1)=IDATE(2)
	IDAT(2)=IDATE(3)

      IYR  =IDAT(3)
      IMNTH=IDAT(1)
      IDAY =IDAT(2)
c      IFCST=(NTSD-1)*DT

!      IFCST=NTSPH*ITAG
      IFCST=3600*ITAG
	IHR=ITAG
	write(6,*) 'IFCST: ', IFCST
C
      WRITE(6,*)' POST PROFILE FOR ',
     1                       IYR,IMNTH,IDAY,IHR
	write(6,*) 'IHRST= ', IHRST
C
C***  SET RTSPH,RTSCU,RTSRA TO 1. OVER THE NUMBER OF TIMES THE
C***  VARIOUS PHYSICS ROUTINES HAVE BEEN
C***  CALLED SINCE LAST OUTPUT OF PROFILER DATA.  NECESSARY FOR
C***  CORRECT AVERAGING OF VARIABLES.
C
      IF(APHTIM.GT.0.)THEN
        RTSPH=1./(APHTIM-APHTIM0)
      ELSE
        RTSPH=1.
      ENDIF
	write(6,*) 'APHTIM, RTSPH: ', APHTIM, APHTIM-APHTIM0, RTSPH
C
C  Counter AVRAIN is reset every 6 or 3-h in WRF-NMM, so for stable
c  precip heating rate use counter for convection (ACUTIM)
c  since in the WRF-NMM microphysics (for stable precip physics)
c  is called at the same time as the convection.
C
      IF(ACUTIM.GT.0.)THEN
        RTSCU=1./(ACUTIM-ACUTIM0)
        RTSRA=1./(ACUTIM-ACUTIM0)
      ELSE
        RTSCU=1.
        RTSRA=1.
      ENDIF
	write(6,*) 'ACUTIM, RTSCU: ', ACUTIM, ACUTIM-ACUTIM0, RTSCU
C
c     IF(AVRAIN.GT.0.)THEN
c       RTSRA=1./(AVRAIN-AVRAIN0)
c     ELSE
c       RTSRA=1.
c     ENDIF
c      write(6,*) 'AVRAIN, RTSRA: ', AVRAIN, AVRAIN-AVRAIN0, RTSRA
C
C--------------------------------------------------------------------------
C--------------------------------------------------------------------------
C***
C***  OUTPUT PROFILE DATA.  THE FOLLOWING LOOP IS OVER ALL PROFILE SITES.
C***
C--------------------------------------------------------------------------
	LCLAS1=79
       curtim=timef()
       write(0,*) 'elapsed to begin preparing for write : ',
     + 1.e-3*(curtim-btim)

      OPEN(UNIT=LCLAS1,ACCESS='DIRECT',RECL=LRECPR,IOSTAT=IER)
C--------------------------------------------------------------------------
	write(6,*) 'RECORD LENGTH = ', LRECPR

      DO 1000 N=1,NUMSTA
C
C***  ZERO OUTPUT ARRAY.
C
      DO K=1,NWORDM
        PRODAT(K)=0.
        FPACK(K) =0.
      ENDDO
C
C***  CONSTRUCT HEADER FOR CURRENT PROFILE SITE.  THE HEADER CONTAINS
C***  THE FOLLOWING INFORMATION:  PACKED CYCLE-DATE, FORECAST TIME,
C***  INTEGER STATION ID, STATION LATITUDE, STATION LONGITUDE, STATION
C***  ELEVATION, NUMBER OF VERTICAL LEVELS IN PROFILE, NUMBER OF MULTI-
C***  LEVEL PARAMETERS, NUMBER OF SINGLE LEVEL PARAMETERS, TOTAL LENGTH
C***  (IN WORDS) OF MULTI- AND SINGLE LEVEL DATA, PROFILE CLASS FLAG,
C***  AND A DUMMY WORD FOR FUTURE USE.
C
	if (N .eq. 20) write(6,*) 'PD(20) (e): ', PD(20)

      IH=IHINDX(N)
      JH=JHINDX(N)
      LMHK     = LMH(IH,JH)
      NWORD2   = 2*LMHK
      NWORD3   = 3*LMHK
      NWORD4   = 4*LMHK
      NWORD5   = 5*LMHK
      NWORD6   = 6*LMHK
      NWORD7   = 7*LMHK
      NWORD8   = 8*LMHK
      NWORD9   = 9*LMHK
      NWORD10  = 10*LMHK
      NWORD11  = 11*LMHK
      NWORD12  = 12*LMHK
      NWORD13  = 13*LMHK
      NWORD14  = 14*LMHK
      NWORD15  = 15*LMHK
      ISTAT    = IDSTN(N)
      CISTAT   = CIDSTN(N)
C
      FPACK(1) = STNLAT(N)/DTR
      FPACK(2) = -STNLON(N)/DTR
!CHANGEDFORNMM      FPACK(2) = STNLON(N)/DTR
      IF(FPACK(2).LT.-180.)FPACK(2)=FPACK(2)+360.
      FPACK(3) = FIS(N)*GI
      FPACK(4) = FLOAT(LMHK)
      FPACK(5) = LCL1ML
      FPACK(6) = LCL1SL
      FPACK(7) = 9+FPACK(5)*FPACK(4)+FPACK(6)
      FPACK(8) = 999.
      FPACK(9) = 999.

C
C***  WIND ROTATION SINES AND COSINES
C

	write(6,*) 'STNLON, TLM0D, COSPH0, SINPH0: ', N,STNLON(N), 
     &                                       TLM0D, COSPH0, SINPH0

      DLM    = STNLON(N)+TLM0D*DTR
      XX     = COSPH0*COS(STNLAT(N))*COS(DLM)
     1        +SINPH0*SIN(STNLAT(N))
      YY     = -COS(STNLAT(N))*SIN(DLM)
      TLON   = ATAN(YY/XX)
      ALPHA  = ASIN(SINPH0*SIN(TLON)/COS(STNLAT(N)))
      SINALP = SIN(ALPHA)
      COSALP = COS(ALPHA)

!      SINALP = SROT(N)
!      COSALP = CROT(N)
C
C------------------------------------------------------------------
C***  EXTRACT PRESSURE AND TEMPERATURE PROFILES.
C***  EXTRACT/ROTATE U AND V WIND COMPONENT PROFILES.
C***  EXTRACT SPECIFIC HUMIDITY AND TEMPERATURE TENDENCY.
C***  EXTRACT CLOUD WATER, HEATING DUE TO CONVECTION, LARGE
C***  SCALE RAIN, SHORT WAVE RADIATION, LONG WAVE RADIATION,
C***  AND CLOUD FRACTION.
C------------------------------------------------------------------
C
      DO LV=1,LMHK
        LVL=LMHK-LV+1
!        PRODAT(LVL)      = PDSL1(N)*AETA(LV)+PT
        PRODAT(LVL)      = PMID(N,LV)
        PRODAT(LMHK+LVL) = T(N,LV)

C***  ROTATE WINDS
C
        UT     = U(N,LV)
        VT     = V(N,LV)

        PRODAT(NWORD2+LVL) = UT*COSALP+VT*SINALP
        PRODAT(NWORD3+LVL) = VT*COSALP-UT*SINALP
C
        PRODAT(NWORD4+LVL) = Q(N,LV)
C
        IF(RTOP(N,LV).GT.1.E-12) THEN
           PRODAT(NWORD5+LVL) = OMGA(N,LV)
	ENDIF

c       if(N.eq.1 .or. N.eq.100) then
c         write(*,'(a5,i5,f8.0,3f8.2)')'TEST:',N,PRODAT(LVL),
c    &        PRODAT(LMHK+LVL),PRODAT(NWORD2+LVL),
c    &        PRODAT(NWORD3+LVL)
c       end if


!
!--- Separate posting of cloud water (NWORD6), rain (NWORD14), and
!    ice (cloud ice + snow, NWORD13).  The ice and rain fields are
!    *APPENDED* at the end of the section dealing with vertical
!    profile fields  (Ferrier/Manikin 11/30/04)
!
        QLIQ = (1.-F_ice(N,LV))*CWM(N,LV)
        PRODAT(NWORD13+LVL) = F_rain(N,LV)*QLIQ
        PRODAT(NWORD6+LVL) = QLIQ-PRODAT(NWORD13+LVL)
        PRODAT(NWORD14+LVL) = F_ice(N,LV)*CWM(N,LV)
        PRODAT(NWORD7+LVL) = TCUCN(N,LV)
        PRODAT(NWORD8+LVL) = TRAIN(N,LV)
        PRODAT(NWORD9+LVL) = RSWTT(N,LV)
        PRODAT(NWORD10+LVL)= RLWTT(N,LV)
        PRODAT(NWORD11+LVL)= CLDFRA(N,LV)*100.
C
        IF(LV.EQ.1)THEN
          PRODAT(NWORD12+LVL)=Q2(N,LV)
        ELSE
          PRODAT(NWORD12+LVL)=(Q2(N,LV)+Q2(N,LV-1))*0.5
        ENDIF
      ENDDO

C
C***  MODIFY ACCUMLATIONS SO AS TO REPRESENT ACCUMULATED
C***  CHANGE SINCE LAST PROFILE OUTPUT TIME.
C
CGSM  MODIFIED CODE TO ACCOUNT FOR DHCNVC AND DHRAIN BEING
C       COMPUTED FROM TOP DOWN WHILE PRODAT IS FILLED FROM
C       BOTTOM UP 
C
C Suspect this is no longer necessary in WRF-NMM
C
      DO LL=1,LMHK
        LVL=LMHK-LL+1
         STADHC(LL) = PRODAT(NWORD7+LL) - DHCNVC(LVL,N)
         STADHR(LL) = PRODAT(NWORD8+LL) - DHRAIN(LVL,N)
C
         DHCNVC(LVL,N) = PRODAT(NWORD7+LL)
         DHRAIN(LVL,N) = PRODAT(NWORD8+LL)
C
       IF(MOD(NTSD,NHEAT).LT.NCNVC)THEN
          DHCNVC(LVL,N) = 0.
          DHRAIN(LVL,N) = 0.
       ENDIF
      ENDDO
C
C***  EXTRACT SINGLE LEVEL DATA.   EGRID2 IS SURFACE TEMPERATURE.
C
      PRODAT(NWORD15+1)  = PSLP  (N)
      PRODAT(NWORD15+2)  = PDS   (N)
      PRODAT(NWORD15+3)  = EGRID2(N)
      PRODAT(NWORD15+4)  = TLMIN (N)
      PRODAT(NWORD15+5)  = TLMAX (N)
c     print *,'tlmax, tlmin ',n,tlmax(n),tlmin(n)
      PRODAT(NWORD15+6)  = SMSTAV(N)*100.
      PRODAT(NWORD15+7)  = ACPREC(N)*1000.
      PRODAT(NWORD15+8)  = CUPREC(N)*1000.
      PRODAT(NWORD15+27) = Z0    (N)
C
      STAPRX=PRODAT(NWORD15+7)-STATPR(N)
	if (IDSTN(N) .eq. 724050) then
	write(6,*) 'STAPRX at DCA: ', STAPRX
	endif
      STACRX=PRODAT(NWORD15+8)-STACPR(N)

!	if (STAPRX .gt. 0) then
!	write(6,*) '1hr precip: ',  N,STAPRX
!	endif
C
C***  ROTATE WINDS
C
      UT     = U10(N)
      VT     = V10(N)
      PRODAT(NWORD15+28) = UT*COSALP+VT*SINALP
      PRODAT(NWORD15+29) = VT*COSALP-UT*SINALP
C
      PRODAT(NWORD15+30) = TH10  (N)
      PRODAT(NWORD15+31) = Q10   (N)
      PRODAT(NWORD15+32) = TSHLTR(N)
      PRODAT(NWORD15+33) = QSHLTR(N)
      PRODAT(NWORD15+34) = SFCEXC(N)
      PRODAT(NWORD15+35) = VEGFRC(N)*100.
      PRODAT(NWORD15+36) = CMC   (N)*1000.
      PRODAT(NWORD15+37) = SMC   (N,1)
      PRODAT(NWORD15+38) = SMC   (N,2)
      PRODAT(NWORD15+39) = SMC   (N,3)
      PRODAT(NWORD15+40) = SMC   (N,4)
      PRODAT(NWORD15+41) = STC   (N,1)
      PRODAT(NWORD15+42) = STC   (N,2)
      PRODAT(NWORD15+43) = STC   (N,3)
      PRODAT(NWORD15+44) = STC   (N,4)
      PRODAT(NWORD15+45) = SM    (N) + SICE(N)
      PRODAT(NWORD15+46) = CFRACL(N)*100.
      PRODAT(NWORD15+47) = CFRACM(N)*100.
      PRODAT(NWORD15+48) = CFRACH(N)*100.
      PRODAT(NWORD15+49) = SR    (N)*100.
      PRODAT(NWORD15+50) = NINT(HBOT(N))
      PRODAT(NWORD15+51) = CPRATE(N)
      PRODAT(NWORD15+52) = F_RIMEF(N,LMHK)
C
      PRODAT(NWORD15+9)   = SFCLHX(N)
      PRODAT(NWORD15+10)  = POTFLX(N)
      PRODAT(NWORD15+11)  = SFCSHX(N)
      PRODAT(NWORD15+12)  = SUBSHX(N)
      PRODAT(NWORD15+13)  = SNOPCX(N)
      PRODAT(NWORD15+14)  = ASWIN (N)
      PRODAT(NWORD15+15)  = ASWOUT(N)
      PRODAT(NWORD15+16)  = ALWIN (N)
      PRODAT(NWORD15+17)  = ALWOUT(N)
      PRODAT(NWORD15+18)  =-ALWTOA(N)
      PRODAT(NWORD15+19)  = ASWTOA(N)
!      PRODAT(NWORD15+20)  = ACSNOW(N)*1000.
      PRODAT(NWORD15+20)  = ACSNOW(N)
!      PRODAT(NWORD15+21)  = SMSTOT(N)*1000.
      PRODAT(NWORD15+21)  = SMSTOT(N)
!      PRODAT(NWORD15+22)  = SNO   (N)*1000.
      PRODAT(NWORD15+22)  = SNO   (N)
!      PRODAT(NWORD15+23)  = ACSNOM(N)*1000.
      PRODAT(NWORD15+23)  = ACSNOM(N)
!      PRODAT(NWORD15+24)  = SSROFF(N)*1000.
      PRODAT(NWORD15+24)  = SSROFF(N)
!      PRODAT(NWORD15+25)  = BGROFF(N)*1000.
      PRODAT(NWORD15+25)  = BGROFF(N)
      PRODAT(NWORD15+26)  = SOILTB(N)
C
C***  ACCUMULATED CHANGE SINCE LAST PROFILE OUTPUT TIME.
C
      PSFCEVP  = PRODAT(NWORD15+9 ) - STAEVP(N)
      PPOTEVP  = PRODAT(NWORD15+10) - STAPOT(N)
      PSFCSHX  = PRODAT(NWORD15+11) - STASHX(N)
      PSFCSUB  = PRODAT(NWORD15+12) - STASUB(N)
      PSNOPCX  = PRODAT(NWORD15+13) - STAPCX(N)
      PRSWIN   = PRODAT(NWORD15+14) - STASWI(N)
!       write(6,*) 'N, SUBSHX(N),STASUB(N),PSFCSUB ', N,
!    &     SUBSHX(N), STASUB(N), PSFCSUB
      PRSWOUT  = PRODAT(NWORD15+15) - STASWO(N)
      PRLWIN   = PRODAT(NWORD15+16) - STALWI(N)
!	if (N .eq. 90) then
!	write(6,*) 'N, PRODAT(NWORD15+16) , STALWI(N), PRLWIN: ', N,
!     &     PRODAT(NWORD15+16), STALWI(N), PRLWIN
!	endif
      PRLWOUT  = PRODAT(NWORD15+17) - STALWO(N)
      PRLWTOA  = PRODAT(NWORD15+18) - STALWT(N)
      PRSWTOA  = PRODAT(NWORD15+19) - STASWT(N)
      PACSNOW  = PRODAT(NWORD15+20) - STASNO(N)
      PACSNOM  = PRODAT(NWORD15+23) - STASNM(N)
      PSSROFF  = PRODAT(NWORD15+24) - STASRF(N)
      PBGROFF  = PRODAT(NWORD15+25) - STABRF(N)
!	write(6,*) 'N, PACSNOW: ', N, PACSNOW
C***
C***  TRANSFER STATION PROFILE DATA TO "PACKED" OUTPUT ARRAY.
C***
      NN   = 0
      NLEN = FPACK(7)
C	write(6,*) 'NWORD15+41,NWORD15+32 ', NWORD15+41,NWORD15+32
C	write(6,*) 'SOIL TEMP ', PRODAT(NWORD13+41)
C        write(6,*) 'SHELT TEMP ', PRODAT(NWORD13+32) 
C
      DO NL = 10,NLEN
        NN = NL-9
        FPACK(NL) = PRODAT(NN)
!       if(NN.eq.NWORD15+22) then
!         print *,'fpack ',nl,nn,nword15+22,fpack(nl),prodat(nn)
!       endif
      ENDDO
C
C***  REPLACE ACCUMULATED QUANTITIES WITH ACCUMULATION
C***  SINCE LAST PROFILE OUTPUT TIME.
C
      DO LL = 1,LMHK
        FPACK(9+NWORD7+LL) = STADHC(LL)*RTSCU
        FPACK(9+NWORD8+LL) = STADHR(LL)*RTSRA
!        FPACK(9+NWORD7+LL) = -9999.
!        FPACK(9+NWORD8+LL) = -9999.
      ENDDO
C
      FPACK(9+NWORD15+7)  = STAPRX
!	write(6,*) 'precip written to FPACK element: ', 9+NWORD13+7
      FPACK(9+NWORD15+8)  = STACRX
      FPACK(9+NWORD15+9)  = PSFCEVP * RTSPH
      FPACK(9+NWORD15+10) = PPOTEVP * RTSPH

!	if (IDSTN(N) .eq. 724050) then
!	write(6,*) 'PPOTEVP, RTSPH: ', PPOTEVP, RTSPH
!	write(6,*) 'packing PPOTEVP*RTSPH= ', FPACK(9+NWORD13+10)
!	endif

      FPACK(9+NWORD15+11) = PSFCSHX * RTSPH
      FPACK(9+NWORD15+12) = PSFCSUB * RTSPH
!       write(6,*) 'N, SUBSHX(N),STASUB(N),PSFCSUB,FPACK ', N,
!    &     SUBSHX(N), STASUB(N), PSFCSUB, FPACK(9+NWORD15+12)
!      FPACK(9+NWORD15+12) =-PSFCSUB * RTSPH
      FPACK(9+NWORD15+13) = PSNOPCX * RTSPH
      FPACK(9+NWORD15+14) = PRSWIN  * RTSPH
	if (N .eq. 90) then
	write(6,*) 'N, RTSPH, SWRD: ', N, RTSPH, FPACK(9+NWORD13+14)
	endif

!	if (mod(N,NUMSTA/5) .eq. 0) then
!	write(6,*) 'N, RTSPH, SWRD: ', N, RTSPH, FPACK(9+NWORD13+14)
!	endif
      FPACK(9+NWORD15+15) = PRSWOUT * RTSPH
      FPACK(9+NWORD15+16) = PRLWIN  * RTSPH
      FPACK(9+NWORD15+17) = PRLWOUT * RTSPH
      FPACK(9+NWORD15+18) = PRLWTOA * RTSPH
      FPACK(9+NWORD15+19) = PRSWTOA * RTSPH
      FPACK(9+NWORD15+20) = PACSNOW
      FPACK(9+NWORD15+23) = PACSNOM
      FPACK(9+NWORD15+24) = PSSROFF
      FPACK(9+NWORD15+25) = PBGROFF
C
!      IF(RESTRT)THEN
      IF(ITAG .eq. 0)THEN
        DO LL = 1,LMHK
          FPACK(9+NWORD7+LL) = 0.
          FPACK(9+NWORD8+LL) = 0.
        ENDDO
C
        FPACK(9+NWORD15+7)  = 0.
        FPACK(9+NWORD15+8)  = 0.
        FPACK(9+NWORD15+9)  = 0.
        FPACK(9+NWORD15+10) = 0.
        FPACK(9+NWORD15+11) = 0.
        FPACK(9+NWORD15+12) = 0.
        FPACK(9+NWORD15+13) = 0.
        FPACK(9+NWORD15+14) = 0.
        FPACK(9+NWORD15+15) = 0.
        FPACK(9+NWORD15+16) = 0.
        FPACK(9+NWORD15+17) = 0.
        FPACK(9+NWORD15+18) = 0.
        FPACK(9+NWORD15+19) = 0.
        FPACK(9+NWORD15+20) = 0.
        FPACK(9+NWORD15+23) = 0.
        FPACK(9+NWORD15+24) = 0.
        FPACK(9+NWORD15+25) = 0.
      ENDIF
C---------------------------------------------------------------------
C***
C***  WRITE PROFILE DATA
C***
      
cZhou      NREC=(IFHR/INCR)*NUMSTA+N
        NREC=N
      write(6,*) 'NREC,IHRST,IDAT: ', NREC,IHRST,IDAT,IFCST

!        if (mod(NREC,10) .eq. 0) then
        write(6,*) 'IDAT, IFCST, ISTAT, CISTAT: ',
     &  IDAT, IFCST, ISTAT, CISTAT
!        endif

      WRITE(LCLAS1,REC=NREC)IHRST,IDAT,IFCST,ISTAT,CISTAT
     1,                    (FPACK(NL),NL=1,NLEN)
      if(nrec.le.2.or.nrec.ge.1376) then
        write(6,*)'IHRST,IDAT,IFCST,ISTAT,CISTAT,NLEN ',
     1     IHRST,IDAT,IFCST,ISTAT,CISTAT,NLEN
      endif

c     if(NREC.eq.1 .or. NREC.eq.100) then
c        write(6,*)'IHRST,IDAT,IFCST,ISTAT,CISTAT ',
c    1     IHRST,IDAT,IFCST,ISTAT,CISTAT

c       write(*,'(a5,10f8.0)') 'F1-10:',(FPACK(k),k=1,10)
c      do kk=1,LMHK
c       write(*,'(a5,f8.0,3f8.3)') 'FPACK',FPACK(9+kk),
c    &    FPACK(9+LMHK+kk),FPACK(9+2*LMHK+kk),
c    &    FPACK(9+3*LMHK+kk)
c      end do
c     end if
C---------------------------------------------------------------------
 1000 CONTINUE
       curtim=timef()
       write(0,*) 'elapsed to end of all writes : ',
     + 1.e-3*(curtim-btim)
      CLOSE(LCLAS1)
C
C***  END OF PROFILE SITE LOOP
C
C***  END PROFILE POSTING CODE.
C---------------------------------------------------------------------

  979  continue
        DEALLOCATE( DETA,RDETA,AETA,UL
     &,RES,FIS,THS,HBOT
     &,CFRACL,CFRACM,CFRACH,SNO
     &,SOILTB,SFCEXC,SMSTAV,SMSTOT
     &,Z0,CZEN,CZMEAN,SR
     &,ACPREC,CUPREC,ACSNOW,ACSNOM
     &,SSROFF,BGROFF,SFCSHX,SFCLHX
     &,SUBSHX,SNOPCX,ASWIN,ASWOUT
     &,ASWTOA,ALWIN,ALWOUT,ALWTOA
     &,TSHLTR,TSHLTR_hold,QSHLTR,PSHLTR
     &,TH10,Q10,U10,V10
     &,TLMIN,TLMAX
     &,SMC,CMC,STC,SH2O
     &,VEGFRC,POTFLX,PSLP,PDSL1
     &,EGRID2,SM,SICE
     &,HBM2,FACTR
     &,PTBL,TTBL,VEGFRA
     &,T,Q,U,V,Q2
     &,OMGALF,CWM,TRAIN,TCUCN
     &,F_RAIN,F_ICE,CLDFRA
     &,RSWTT,RLWTT,RTOP
     &,HTM,OMGA
     &,PRODAT,FPACK
     &,STATPR,STACPR,STAEVP
     &,STAPOT,STASHX,STASUB,STAPCX
     &,STASWI,STASWO,STALWI,STALWO
     &,STALWT,STASWT,STASNM,STASRF
     &,STABRF,STASNO,DHCNVC,DHRAIN
     &,STADHC,STADHR,CPRATE
     &,ACPREC0,CUPREC0,SFCLHX0,POTFLX0
     &,SFCSHX0,SUBSHX0,SNOPCX0,ASWIN0
     &,ASWOUT0,ALWIN0,ALWOUT0,ALWTOA0
     &,ASWTOA0,ACSNOW0,ACSNOM0,SSROFF0
     &,BGROFF0,SLDPTH2)
   
      print *," deallocate 1"
      print *, "PROF_NMM DONE!!!"

!      call mpi_file_close(iunit,ierr)
!      call mpi_finalize(mpi_comm_world, ierr)

      RETURN
      END

!
!       function written early Dec. 1999 by M. Pyle to support  workstation
!       Eta for users with etime but not timef functionality (like  certain
!mp     HPs)  Designed to duplicate timef (elapsed time in milliseconds)
!
        function timef()
        real et(2)
        real*8 timef
        timef=etime(et)
        timef=timef*1.e3
        end

        function rtc()
        real et(2)
        real*8 rtc
        rtc=etime(et)
        rtc=rtc*1.e3
        end