Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature/SatDiagnUpdates (netCDF Satellite Diagnostics) #1134

Closed
1 change: 1 addition & 0 deletions AUTHORS.txt
Original file line number Diff line number Diff line change
Expand Up @@ -245,6 +245,7 @@ Viral Shah
Jingyuan Shao
Lu Shen
Tomas Sherwen
Joshua Shutter
Sam Silva
Nicole Smith-Downey
Anne Laerke Soerensen
Expand Down
20 changes: 10 additions & 10 deletions GeosCore/diag51_mod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -158,9 +158,9 @@ MODULE DIAG51_MOD
INTEGER :: IU_ND51

! Arrays
INTEGER, ALLOCATABLE :: GOOD(:)
INTEGER, ALLOCATABLE :: GOOD_CT(:)
INTEGER, ALLOCATABLE :: COUNT_CHEM3D(:,:,:)
REAL(fp), ALLOCATABLE :: GOOD(:)
REAL(fp), ALLOCATABLE :: GOOD_CT(:)
REAL(fp), ALLOCATABLE :: COUNT_CHEM3D(:,:,:)
REAL(fp),ALLOCATABLE :: Q(:,:,:,:)

CONTAINS
Expand Down Expand Up @@ -295,7 +295,7 @@ SUBROUTINE GET_LOCAL_TIME( Input_Opt, State_Grid )
! GOOD indicates which boxes have local times between HR1 and HR2
IF ( LT >= Input_Opt%ND51_HR1 .and. &
LT <= Input_Opt%ND51_HR2 ) THEN
GOOD(I) = 1
GOOD(I) = 1e+0_fp
ENDIF
ENDDO

Expand Down Expand Up @@ -552,7 +552,7 @@ SUBROUTINE ACCUMULATE_DIAG51( Input_Opt, State_Chm, &
! Accumulate data
Q(X,Y,K,W) = Q(X,Y,K,W) + &
( State_Chm%Species(I,J,L,id_OH) * GOOD(X) ) * &
( State_Met%AIRDEN(I,J,L) * CONV_OH )
( State_Met%AIRDEN(I,J,L) * CONV_OH * 1.0e+3_fp) !1.0e+3_fp added by JDS for correct unit conversion (1/20/22)


ELSE IF ( N == 502 .and. IS_NOy ) THEN
Expand Down Expand Up @@ -1082,7 +1082,7 @@ SUBROUTINE ACCUMULATE_DIAG51( Input_Opt, State_Chm, &
ENDDO
!$OMP END PARALLEL DO

GOOD(:) = 0
GOOD(:) = 0e+0_fp

! Free pointers
Spc => NULL()
Expand Down Expand Up @@ -1944,22 +1944,22 @@ SUBROUTINE INIT_DIAG51( Input_Opt, State_Grid, RC )
! Array denoting where LT is between HR1 and HR2
ALLOCATE( GOOD( State_Grid%NX ), STAT=AS )
IF ( AS /= 0 ) CALL ALLOC_ERR( 'GOOD' )
GOOD = 0
GOOD = 0e+0_fp

! Counter of "good" times per day at each grid box
ALLOCATE( GOOD_CT( ND51_NI ), STAT=AS )
IF ( AS /= 0 ) CALL ALLOC_ERR( 'GOOD_CT' )
GOOD_CT = 0
GOOD_CT = 0e+0_fp

! Accumulating array
ALLOCATE( Q( ND51_NI, ND51_NJ, ND51_NL, Input_Opt%N_ND51),STAT=AS)
IF ( AS /= 0 ) CALL ALLOC_ERR( 'Q' )
Q = 0d0
Q = 0e+0_fp

! Accumulating array
ALLOCATE( COUNT_CHEM3D( ND51_NI, ND51_NJ, ND51_NL ), STAT=AS )
IF ( AS /= 0 ) CALL ALLOC_ERR( 'COUNT_CHEM3D' )
COUNT_CHEM3D = 0
COUNT_CHEM3D = 0e+0_fp

END SUBROUTINE INIT_DIAG51
!EOC
Expand Down
Loading