Merge pull request #14434 from mcgratta/master

FDS Source: Remove initialization of pointers to NULL()

Author: mcgratta

Source/ccib.f90

@@ -83,9 +83,9 @@ MODULE CC_SCALARS
 
 TYPE(WALL_TYPE),           POINTER :: WC
 TYPE(CFACE_TYPE),          POINTER :: CFA
-TYPE(CC_CUTCELL_TYPE),     POINTER :: CC=>NULL()
-TYPE(CC_CUTFACE_TYPE),     POINTER :: CF=>NULL()
-TYPE(CC_RCFACE_TYPE),      POINTER :: RCF=>NULL()
+TYPE(CC_CUTCELL_TYPE),     POINTER :: CC
+TYPE(CC_CUTFACE_TYPE),     POINTER :: CF
+TYPE(CC_RCFACE_TYPE),      POINTER :: RCF
 TYPE(BOUNDARY_COORD_TYPE), POINTER :: BC
 TYPE(EXTERNAL_WALL_TYPE),  POINTER :: EWC
 TYPE(BOUNDARY_PROP1_TYPE), POINTER :: B1
@@ -390,8 +390,8 @@ SUBROUTINE CUTFACE_VELOCITIES(NM,UU,VV,WW,CUTFACES)
 INTEGER, INTENT(IN) :: NM
 REAL(EB), POINTER, DIMENSION(:,:,:), INTENT(INOUT) :: UU,VV,WW
 LOGICAL, INTENT(IN) :: CUTFACES
-TYPE(CC_CUTFACE_TYPE), POINTER :: CF=>NULL()
-TYPE(CFACE_TYPE), POINTER :: CFA=>NULL()
+TYPE(CC_CUTFACE_TYPE), POINTER :: CF
+TYPE(CFACE_TYPE), POINTER :: CFA
 TYPE(BOUNDARY_COORD_TYPE), POINTER :: BC
 TYPE(BOUNDARY_PROP1_TYPE), POINTER :: B1
 INTEGER :: ICF,ICFA,JCF,I,J,K,X1AXIS
@@ -1697,9 +1697,9 @@ SUBROUTINE GET_FH_FROM_PRHS_AND_BCS(NM,DT,CYL_FCT,UNKH,NUNKH,IPZ,F_H)
 
 ! Local variables:
 INTEGER :: I,J,K,IROW,IW,IIG,JJG,KKG,IOR,ICFACE,IFACE,JFACE,ILH,JLH,KLH,IRC
-TYPE (WALL_TYPE), POINTER :: WC=>NULL()
-TYPE (EXTERNAL_WALL_TYPE), POINTER :: EWC=>NULL()
-TYPE(CFACE_TYPE), POINTER :: CFA=>NULL()
+TYPE (WALL_TYPE), POINTER :: WC
+TYPE (EXTERNAL_WALL_TYPE), POINTER :: EWC
+TYPE(CFACE_TYPE), POINTER :: CFA
 REAL(EB) :: IDX, AF, VAL, BCV
 
 ! Dummy Assignment.
@@ -5565,8 +5565,8 @@ SUBROUTINE CC_COMPUTE_VISCOSITY(DT,NM)
 
 ! Local Variables:
 INTEGER :: I,J,K
-REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP=>NULL(),UU=>NULL(),VV=>NULL(),WW=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP,UU,VV,WW
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP
 REAL(EB) :: NU_EDDY,DELTA,A_IJ(3,3),DUDX,DUDY,DUDZ,DVDX,DVDY,DVDZ,DWDX,DWDY,DWDZ
 
 REAL(EB) :: TNOW
@@ -5789,7 +5789,7 @@ SUBROUTINE ROTATED_CUBE_VELOCITY_FLUX(NM,TLEVEL)
 ! Local Variables:
 REAL(EB) :: XGLOB(2,1), XLOC(2,1), FGLOB(2,1), FLOC(2,1), X_I, Y_J, NU
 INTEGER  :: I,J,K
-REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP
 REAL(EB) :: SIN_T, COS_T, COS_KX, SIN_KX, COS_KY, SIN_KY
 
 CALL POINT_TO_MESH(NM)
@@ -6602,7 +6602,7 @@ SUBROUTINE CC_VELOCITY_CUTFACES(APPLY_TO_ESTIMATED_VARIABLES)
 INTEGER  :: NM,ICF,I,J,K,X1AXIS
 REAL(EB) :: AREATOT, A_CART
 REAL(EB), POINTER, DIMENSION(:,:,:) :: UP,VP,WP
-TYPE(CC_CUTFACE_TYPE), POINTER :: CF=>NULL()
+TYPE(CC_CUTFACE_TYPE), POINTER :: CF
 
 ! Mesh Loop:
 MESH_LOOP : DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
@@ -7510,7 +7510,7 @@ SUBROUTINE SET_EXIMADVFLX_3D(NM,UU,VV,WW)
 
 ! Local Variables:
 INTEGER :: N,I,J,K,X1AXIS,IFACE
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 ! Loop on scalars:
 SPECIES_LOOP : DO N=1,N_TOTAL_SCALARS
 
@@ -7555,9 +7555,9 @@ SUBROUTINE SET_EXIMRHOZZLIM_3D(NM,N)
 
 ! Local Variables:
 INTEGER :: I,J,K,X1AXIS
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: FX_ZZ=>NULL(),FY_ZZ=>NULL(),FZ_ZZ=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: FX_ZZ,FY_ZZ,FZ_ZZ
 INTEGER :: IFACE,IW
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 
 FX_ZZ=>SWORK1
 FY_ZZ=>SWORK2
@@ -7614,10 +7614,10 @@ SUBROUTINE SET_EXIMRHOHSLIM_3D(NM)
 
 ! Local Variables:
 INTEGER :: I,J,K,X1AXIS
-REAL(EB), POINTER, DIMENSION(:,:,:) :: FX_H_S=>NULL(),FY_H_S=>NULL(),FZ_H_S=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:) :: FX_H_S,FY_H_S,FZ_H_S
 REAL(EB) :: H_S,ZZ_GET(1:N_TRACKED_SPECIES),TMP_F_GAS,VELC2
 INTEGER :: IFACE,IW
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 
 FX_H_S=>WORK2
 FY_H_S=>WORK3
@@ -7685,7 +7685,7 @@ SUBROUTINE SET_EXIMDIFFLX_3D(NM,RHO_D_DZDX,RHO_D_DZDY,RHO_D_DZDZ)
 
 ! Local Variables:
 INTEGER :: N,I,J,K,X1AXIS,IFACE
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 
 ! First, set all diffusive fluxes to zero on CC_SOLID faces:
 ! IAXIS:
@@ -8462,7 +8462,7 @@ SUBROUTINE CCSPECIES_ADVECTION
 REAL(EB), PARAMETER :: SGNFCT=1._EB
 INTEGER :: IOR, ICFA
 TYPE(BOUNDARY_PROP1_TYPE), POINTER :: B1
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: RGF=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: RGF
 
 U_DOT_DEL_RHO_Z_VOL=>WORK7
 U_DOT_DEL_RHO_Z_VOL=0._EB
@@ -8995,7 +8995,7 @@ SUBROUTINE CCENTHALPY_ADVECTION
 REAL(EB) :: RHO_H_S_PV(-2:1), VELC, VELC2, FN_H_S, TMP_F_GAS
 INTEGER  :: IOR, ICFA
 TYPE(BOUNDARY_PROP1_TYPE), POINTER :: B1
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 LOGICAL :: DO_LO, DO_HI
 
 ! IAXIS faces:
@@ -10676,8 +10676,8 @@ SUBROUTINE CC_DIFFUSIVE_MASS_FLUXES(NM)
 
 ! Local Variables:
 INTEGER :: N,I,J,K,X1AXIS,ISIDE,IFACE,ICC,JCC,ICF
-REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:) :: RHOP
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP
 REAL(EB) :: D_Z_N(0:I_MAX_TEMP),CCM1,CCP1,IDX,DIFF_FACE,D_Z_TEMP(-1:0),MUV(-1:0),MU_DNSV(-1:0), &
             RHOPV(-1:0),TMPV(-1:0),ZZPV(-1:0),X1F,PRFCT
 REAL(EB), ALLOCATABLE, DIMENSION(:) :: ZZ_GET,RHO_D_DZDN_GET
@@ -11522,10 +11522,10 @@ SUBROUTINE GET_EXPLICIT_ADVDIFFVECTOR_SCALAR_3D(N)
 INTEGER :: X1AXIS,IFACE,IND(LOW_IND:HIGH_IND),IND_LOC(LOW_IND:HIGH_IND),ICF
 INTEGER :: LOCROW_1,LOCROW_2,ILOC,IROW,ICC,JCC,ISIDE,IW
 REAL(EB):: AF,VELC,RHO_Z_PV(-2:1),RHOPV(-2:1),FCT,ZZ_GET_N,FN_ZZ
-REAL(EB), POINTER, DIMENSION(:,:,:)  :: RHOP=>NULL(),UP=>NULL(),VP=>NULL(),WP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:)  :: UU=>NULL(),VV=>NULL(),WW=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:):: ZZP=>NULL()
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:)  :: RHOP,UP,VP,WP
+REAL(EB), POINTER, DIMENSION(:,:,:)  :: UU,VV,WW
+REAL(EB), POINTER, DIMENSION(:,:,:,:):: ZZP
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 LOGICAL :: DO_LO,DO_HI
 INTEGER :: IIG,JJG,KKG,IOR
 REAL(EB) :: UN
@@ -11984,8 +11984,8 @@ SUBROUTINE CC_CHECK_MASS_DENSITY
 INTEGER, PARAMETER :: MAX_SURR_CELLS=20
 INTEGER :: NM, NCELL, ICC, JCC, I, J, K, IFC, IFACE, IFC1, JFC1, ICC1, JCC1, IC, LOHI, ILH, X1AXIS, &
            II, JJ, KK, IIF, JJF, KKF, IRC, N
-REAL(EB), POINTER, DIMENSION(:,:,:)   :: DELTA_RHO=>NULL(),DELTA_RHO_ZZ=>NULL(),RHOP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: RHO_ZZ=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:)   :: DELTA_RHO,DELTA_RHO_ZZ,RHOP
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: RHO_ZZ
 REAL(EB) :: MASS_C, MASS_N(1:MAX_SURR_CELLS), RHO_CELL(1:MAX_SURR_CELLS), VOL(1:MAX_SURR_CELLS), &
             RHO_CUT, SIGN_FACTOR, SUM_MASS_N, SUM_RHO_ZZ, CONST, PRFCT, CC_RHOP, CC_RHO_ZZP, RHO_ZZ_MIN, &
             RHO_ZZ_MAX, RHO_ZZ_CUT, RHO_ZZ_TEST
@@ -12648,8 +12648,8 @@ SUBROUTINE GET_RHOZZ_CC_3D
 
 ! Local Variables:
 INTEGER :: NM,N,I,J,K,ICC,JCC
-REAL(EB), POINTER, DIMENSION(:,:,:)   :: RHOP=>NULL(),UP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:)   :: RHOP,UP
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP
 REAL(EB) :: VOLTOT
 INTEGER :: NMX
 
@@ -12842,7 +12842,7 @@ SUBROUTINE PUT_RHOZZVECTOR_SCALAR_3D(N)
 
 ! Local Variables:
 INTEGER :: NM,I,J,K,IROW_LOC,ICC,JCC
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP
 
 ! Loop meshes:
 MESH_LOOP : DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
@@ -12972,9 +12972,9 @@ SUBROUTINE GET_ADVDIFFVECTOR_SCALAR_3D(N)
 INTEGER :: X1AXIS,IFACE,IND(LOW_IND:HIGH_IND),IND_LOC(LOW_IND:HIGH_IND),ICF,IND1,IND2,IOR
 INTEGER :: LOCROW_1,LOCROW_2,ILOC,IROW,ICC,JCC,IW
 REAL(EB):: AF,VELC,RHO_Z,FN_ZZ,FCT
-REAL(EB), POINTER, DIMENSION(:,:,:)  :: RHOP=>NULL(),UP=>NULL(),VP=>NULL(),WP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:)::  ZZP=>NULL()
-TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:)  :: RHOP,UP,VP,WP
+REAL(EB), POINTER, DIMENSION(:,:,:,:)::  ZZP
+TYPE(CC_REGFACEZ_TYPE),  POINTER, DIMENSION(:) :: REGFACE_Z
 LOGICAL :: DO_LO,DO_HI
 
 ! This routine computes RHS due to boundary conditions prescribed in immersed solids
@@ -13214,8 +13214,8 @@ SUBROUTINE CC_VELOCITY_BC(T,NM,APPLY_TO_ESTIMATED_VARIABLES,DO_IBEDGES)
 INTEGER :: I,J,K,II,JJ,KK,ICF,X1AXIS,IE,IEC,I_SGN,ICD,ICD_SGN,IIF,JJF,KKF,FAXIS,ICDO,ICDO_SGN,IS
 REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: UVW_EP
 REAL(EB), ALLOCATABLE, DIMENSION(:,:,:,:) :: DUVW_EP
-REAL(EB), POINTER, DIMENSION(:,:,:) :: UU=>NULL(),VV=>NULL(),WW=>NULL(),RHOP=>NULL()
-REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:) :: UU,VV,WW,RHOP
+REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZP
 REAL(EB) :: NU,MU_FACE,RHO_FACE,DXN_STRM_UB,SLIP_FACTOR,SRGH,U_TAU,Y_PLUS,TNOW,MU_EP,MU_DUIDXJ_USE(2),DUIDXJ_USE(2),&
             DUIDXJ(-2:2),MU_DUIDXJ(-2:2),DXX(2),DF, DE, UE, UF, UB, A_CART, AREATOT, VEL_CART, PRFCT
 TYPE(CC_EDGE_TYPE), POINTER :: CC_EDGE
@@ -13530,7 +13530,7 @@ SUBROUTINE CC_EDGE_TAU_OMG
 LOGICAL :: ALTERED_GRADIENT(-2:2),CORNER_EDGE
 REAL(EB):: XB_IB,OMEV_EP(-2:2),TAUV_EP(-2:2),DWDY,DVDZ,DUDZ,DWDX,DUDY,DVDX, &
            DUIDXJ_EP(-2:2),MU_DUIDXJ_EP(-2:2),EC_B(-2:2),EC_EP(-2:2),DEL_EP,VEL_GAS,CEP,CB,U1_T,U2_T
-TYPE(SURFACE_TYPE), POINTER :: SF=>NULL()
+TYPE(SURFACE_TYPE), POINTER :: SF
 REAL(EB) :: VLG(-2:2),NUV(-2:2),RGH(-2:2),UTA(-2:2)
 VLG(-2:2)=0._EB; NUV(-2:2)=0._EB; RGH(-2:2)=0._EB; UTA(-2:2)=0._EB
 
@@ -15436,7 +15436,7 @@ SUBROUTINE CC_CHECK_DIVERGENCE(T,DT,PREDVEL)
 
 REAL(EB):: PRFCT, DIV, RES, DIVVOL, DIV_JCC, VOL, DPCC, DIV2,TLOC,DTLOC
 
-REAL(EB), POINTER, DIMENSION(:,:,:)  :: UP=>NULL(), VP=>NULL(), WP=>NULL(), DP=>NULL()
+REAL(EB), POINTER, DIMENSION(:,:,:)  :: UP, VP, WP, DP
 REAL(EB), ALLOCATABLE, DIMENSION(:)  :: RESMAXV, RESVOLMX
 REAL(EB), ALLOCATABLE, DIMENSION(:,:):: DIVMNX, DIVVOLMNX, VOLMNX
 INTEGER,  ALLOCATABLE, DIMENSION(:,:):: IJKRM, RESICJCMX
@@ -15887,7 +15887,7 @@ SUBROUTINE GET_CRTCFCC_INT_STENCILS
 INTEGER :: IS,I_SGN,ICD,ICD_SGN,IIF,JJF,KKF,FAXIS,IEC,IE,SKIP_FCT,IEP,JEP,KEP,INDS(1:2,IAXIS:KAXIS),AX,ICEDG,JCEDG,LOHI
 REAL(EB):: DXX(2),AREA_CF,XB_IB,DEL_EP,DEL_IBEDGE,XYZ1(IAXIS:KAXIS),XYZ2(IAXIS:KAXIS)
 
-TYPE(CC_CUTEDGE_TYPE), POINTER :: CE=>NULL()
+TYPE(CC_CUTEDGE_TYPE), POINTER :: CE
 
 REAL(EB) CPUTIME,CPUTIME_START,CPUTIME_START_LOOP
 CHARACTER(100) :: MSEGS_FILE
@@ -18775,7 +18775,7 @@ END SUBROUTINE SET_CC_MATVEC_DATA
 SUBROUTINE GET_RCEDGE_FACE_LIST
 
 INTEGER :: NM,ICF,JCF,I,J,K,X1AXIS,EAXIS,LOHI_AXIS,IRC,ISIDE,HILO
-TYPE(MESH_TYPE), POINTER :: M=>NULL()
+TYPE(MESH_TYPE), POINTER :: M
 
 MESH_LOOP : DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
    M=>MESHES(NM)
@@ -20739,7 +20739,7 @@ SUBROUTINE GET_GASPHASE_REGRCFACES_DATA
 LOGICAL :: FLGIN
 INTEGER, PARAMETER :: OZPOS=0, ICPOS=1, JCPOS=2, IFPOS=3
 INTEGER :: IERR
-TYPE(MESH_TYPE), POINTER :: M=>NULL()
+TYPE(MESH_TYPE), POINTER :: M
 
 ! Mesh loop:
 MAIN_MESH_LOOP : DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
@@ -21927,7 +21927,7 @@ SUBROUTINE GET_LINKED_FACE_INDEXES_F
            OFACE(3),OLO_UNKZ,OHI_UNKZ,OICF,OJCF,IECE,JECE,ILINK,ICC,JCC,IW,COUNT
 REAL(EB):: ACRT,CCVOL_THRES
 LOGICAL :: ALL_FLG,CC_LINKED
-TYPE(MESH_TYPE), POINTER :: M=>NULL()
+TYPE(MESH_TYPE), POINTER :: M
 INTEGER :: ILOC,SIZE_FACE
 INTEGER, ALLOCATABLE, DIMENSION(:,:) :: FACE_LIST,FACELAUX
 REAL(EB), ALLOCATABLE,DIMENSION(:)   :: FACE_AREA,FACEARAUX

Source/ctrl.f90

@@ -23,7 +23,7 @@ SUBROUTINE UPDATE_CONTROLS(T,DT,CTRL_STOP_STATUS,RUN_START)
 INTEGER :: NC,N
 LOGICAL :: CTRL_STOP_STATUS
 LOGICAL, INTENT(IN) :: RUN_START
-TYPE(DEVICE_TYPE), POINTER :: DV=>NULL()
+TYPE(DEVICE_TYPE), POINTER :: DV
 
 CTRL_STOP_STATUS = .FALSE.
 

Source/data.f90

@@ -2258,7 +2258,7 @@ SUBROUTINE THERMO_TABLE_GAS(I_TMP,CP,SPEC_INDEX,SS_INDEX,RCON,G_F)
 REAL(EB) :: TE,H,S,TB(4),A(7,3),BC(2,3)
 INTEGER :: B,BANDS
 CHARACTER(LABEL_LENGTH):: POLY
-TYPE (THERMO_DATA_TYPE), POINTER :: TD=>NULL()
+TYPE (THERMO_DATA_TYPE), POINTER :: TD
 
 G_F = 0._EB !kJ/mol
 IF (SPEC_INDEX==-1 .OR. CONSTANT_SPECIFIC_HEAT_RATIO) THEN
@@ -2348,7 +2348,7 @@ SUBROUTINE THERMO_TABLE_LIQUID(I_TMP,CP,H_V,T_REF,T_MELT,T_BOIL,SPEC_INDEX,FUEL,
 REAL(EB) :: TE
 INTEGER :: B
 LOGICAL,INTENT(OUT) :: FUEL
-TYPE (THERMO_DATA_TYPE), POINTER :: TD=>NULL()
+TYPE (THERMO_DATA_TYPE), POINTER :: TD
 
 !Default to water properteis
 
@@ -2412,7 +2412,7 @@ SUBROUTINE GAS_PROPS(GAS_INDEX,SIGMA,EPSOK,PR_GAS,MW,FORMULA,LISTED,ATOM_COUNTS,
 CHARACTER(FORMULA_LENGTH), INTENT(INOUT) :: FORMULA
 REAL(EB) :: SIGMAIN,EPSOKIN,MWIN,ATOM_COUNTS(118),H_FIN,PR_GASIN
 CHARACTER(FORMULA_LENGTH) :: FORMULAIN
-TYPE(THERMO_DATA_TYPE), POINTER :: TD=>NULL()
+TYPE(THERMO_DATA_TYPE), POINTER :: TD
 
 SIGMAIN = SIGMA
 EPSOKIN = EPSOK
@@ -2512,7 +2512,7 @@ SUBROUTINE CALC_GAS_PROPS(J,N,D_TMP,MU_TMP,K_TMP,CP_TMP,G_F_TMP)
 INTEGER :: I(1)
 REAL(EB), INTENT(OUT) :: D_TMP,MU_TMP,K_TMP,CP_TMP,G_F_TMP
 REAL(EB) :: SIGMA2,EPSIJ,AMW,TSTAR,OMEGA,TMP
-TYPE(SPECIES_TYPE),POINTER :: SS=>NULL(),SS0=>NULL()
+TYPE(SPECIES_TYPE),POINTER :: SS,SS0
 
 SS=>SPECIES(N)
 
@@ -2661,7 +2661,7 @@ SUBROUTINE CALC_MIX_PROPS(J,D_TMP,MU_TMP,K_TMP,CP_TMP,EPSK,SIG,D_USER,MU_USER,K_
 REAL(EB), INTENT(IN) :: EPSK,SIG,D_USER,MU_USER,K_USER,MW,CP_USER,PR_USER
 REAL(EB), INTENT(INOUT) :: D_TMP,MU_TMP,K_TMP,CP_TMP
 REAL(EB) :: SIGMA2,EPSIJ,AMW,TSTAR,OMEGA,TMP,PR_TMP
-TYPE(SPECIES_TYPE),POINTER :: SS0=>NULL()
+TYPE(SPECIES_TYPE),POINTER :: SS0
 
 TMP = REAL(J,EB)
 IF (SIG > 0._EB .AND. EPSK > 0._EB) THEN
@@ -3100,7 +3100,7 @@ SUBROUTINE SHUTDOWN_ATOM(ATOM_COUNT,NR,ERR)
 REAL(EB), INTENT(IN) :: ATOM_COUNT(118),ERR
 REAL(EB) :: REACTION_BALANCE(118)
 LOGICAL :: EX
-TYPE (REACTION_TYPE), POINTER :: RN=>NULL()
+TYPE (REACTION_TYPE), POINTER :: RN
 
 IF (MY_RANK==0) THEN
    CALL WRITE_ATOM_BALANCE_DIAGNOSTICS(LU_ERR)