!-------------------------------------------------------------------------------------------------- ! DOXYGEN !> @brief GENERATE SEASONAL OZONE DISTRIBUTION !> @details COMPUTES THE SEASONAL CLIMATOLOGY OF OZONE USING 81-LAYER DATA FROM GFDL. !> @author ORIGINATOR - GFDL/CAMPANA !> @date ??-??-?? \n !> @author LUCCI !> @date 18-03-20 \n !> @version V1.1.0 !> @details MODERNIZATION OF THE CODE, INCLUDING: !! * F77 TO F90/F95 !! * INDENTATION & UNIFORMIZATION CODE !! * REPLACEMENT OF COMMONS BLOCK FOR MODULES !! * DOCUMENTATION WITH DOXYGEN !! * OPENMP FUNCTIONALITY !< !> @details Use Module: !! @arg @c F77KINDS !! @arg @c MPPSTAFF !! @arg @c O3DATA !! @arg @c SEASO3 !> @details Driver: !! @arg @c INIT !> @details Calls: !-------------------------------------------------------------------------------------------------- SUBROUTINE O3CLIM !-------------------------------------------------------------------------------------------------- ! SUBROUTINE O3CLIM ! ! SUBROUTINE: O3CLIM - GENERATE SEASONAL OZONE DISTRIBUTION ! PROGRAMMER: GFDL/CAMPANA ! ORG: W/NP22 ! DATE: ??-??-?? ! ! ABSTRACT: ! O3CLIM COMPUTES THE SEASONAL CLIMATOLOGY OF OZONE USING 81-LAYER DATA FROM GFDL. ! ! PROGRAM HISTORY LOG: ! ??-??-?? GFDL/KC - ORIGINATOR ! 96-07-26 BLACK - MODIFIED FOR ETA MODEL ! 18-01-15 LUCCI - MODERNIZATION OF THE CODE, INCLUDING: ! * F77 TO F90/F95 ! * INDENTATION & UNIFORMIZATION CODE ! * REPLACEMENT OF COMMONS BLOCK FOR MODULES ! * DOCUMENTATION WITH DOXYGEN ! * OPENMP FUNCTIONALITY ! ! INPUT ARGUMENT LIST: ! NONE ! ! OUTPUT ARGUMENT LIST: ! NONE ! ! INPUT/OUTPUT ARGUMENT LIST: ! NONE ! ! USE MODULES: F77KINDS ! MPPSTAFF ! O3DATA ! SEASO3 ! ! DRIVER : INIT ! ! CALLS : ----- !-------------------------------------------------------------------------------------------------- USE F77KINDS USE MPPSTAFF USE O3DATA USE SEASO3 ! IMPLICIT NONE ! REAL (KIND=R4) ::& & AVG , A1 , B1 , B2 ! INTEGER(KIND=I4) ::& & I , JJ , K , N , NCASE , KI , KK , KQ , NKM , & & IPLACE , NKMM , KEN INTEGER(KIND=I4) ,PARAMETER :: NLGTH = 37 * NL INTEGER(KIND=I4) ,PARAMETER :: NKK = 41 INTEGER(KIND=I4) ,PARAMETER :: NK = 81 INTEGER(KIND=I4) ,PARAMETER :: NKP = NK + 1 ! REAL (KIND=R4) , DIMENSION(NL) ::& & PSTD , RBAR , PHALF ! REAL (KIND=R4) , DIMENSION(19) ::& & TEMPN ! REAL (KIND=R4) , DIMENSION(37, NL) ::& & O35DEG ! REAL (KIND=R4) , DIMENSION(19, NL) ::& & DDUO3N ! REAL (KIND=R4) , DIMENSION(19, 41) ::& & DUO3N ! REAL (KIND=R4) , DIMENSION(10, 41) ::& & RO3 , RO31 , RO32 ! REAL (KIND=R4) , DIMENSION(10, 40) ::& & RO3M ! REAL (KIND=R4) , DIMENSION(81) ::& & RSTD , RDATA !-------------------------------------------------------------------------------------------------- ! COMPUTE DETAILED O3 PROFILE FROM THE ORIGINAL GFDL PRESSURES WHERE OUTPUT FROM O3INT (PSTD) IS ! TOP DOWN IN MB*1.E3 AND PSFC=1013.25 MB ! K.A.C. DEC94 !-------------------------------------------------------------------------------------------------- DO K=1,NK PH(K) = PH(K) * 1013250. P(K) = P(K) * 1013250. END DO ! PH(NKP) = PH(NKP) * 1013250. ! DO K=1,NL PSTD(K) = P(K) END DO ! DO K=1,25 DO N=1,10 RO31(N,K) = O3HI(N,K) RO32(N,K) = O3HI(N,K) END DO END DO ! DO 100 NCASE=1,4 !-------------------------- ! NCASE=1: SPRING (IN N.H.) ! NCASE=2: FALL (IN N.H.) ! NCASE=3: WINTER (IN N.H.) ! NCASE=4: SUMMER (IN N.H.) !-------------------------- IPLACE = 2 ! IF (NCASE == 2) IPLACE = 4 IF (NCASE == 3) IPLACE = 1 IF (NCASE == 4) IPLACE = 3 ! IF (NCASE == 1 .OR. NCASE == 2) THEN DO K=26,41 DO N=1,10 RO31(N,K) = O3LO1(N,K-25) RO32(N,K) = O3LO2(N,K-25) END DO END DO END IF ! IF (NCASE == 3 .OR. NCASE == 4) THEN DO K=26,41 DO N=1,10 RO31(N,K) = O3LO3(N,K-25) RO32(N,K) = O3LO4(N,K-25) END DO END DO END IF ! DO 25 KK=1,NKK DO N=1,10 DUO3N(N ,KK) = RO31(11-N,KK) DUO3N(N+9,KK) = RO32( N,KK) END DO DUO3N(10,KK) = 0.5 * (RO31(1,KK) + RO32(1,KK)) 25 END DO !-------------------------------------------------------------------- ! FOR NCASE=2 OR NCASE=4,REVERSE LATITUDE ARRANGEMENT OF CORR. SEASON !-------------------------------------------------------------------- IF (NCASE == 2 .OR. NCASE == 4) THEN DO 50 KK=1,NKK DO N=1,19 TEMPN(N) = DUO3N(20-N,KK) END DO DO N=1,19 DUO3N(N,KK) = TEMPN(N) END DO 50 END DO END IF !---------------------------------------------------------------------- ! DUO3N NOW IS O3 PROFILE FOR APPROPRIATE SEASON AT STD PRESSURE LEVELS !---------------------------------------------------------------------- ! !----------------------------- ! BEGIN LATITUDE (10 DEG) LOOP !----------------------------- DO 75 N=1,19 ! DO KK=1,NKK RSTD(KK) = DUO3N(N,KK) END DO ! NKM = NK-1 NKMM = NK-3 !----------------------------------------- ! BESSELS HALF-POINT INTERPOLATION FORMULA !----------------------------------------- DO K=4,NKMM,2 KI = K / 2 RDATA(K) = 0.5 * (RSTD(KI ) + RSTD(KI+1)) - (RSTD(KI+2) & & - RSTD(KI+1) - RSTD(KI ) + RSTD(KI-1)) / 16. END DO ! RDATA(2) = 0.5 * (RSTD(2) + RSTD(1)) RDATA(NKM) = 0.5 * (RSTD(NKK) + RSTD(NKK-1)) !---------------------------------- ! PUT UNCHANGED DATA INTO NEW ARRAY !---------------------------------- DO K=1,NK,2 KQ = (K+1) / 2 RDATA(K) = RSTD(KQ) END DO ! DO KK=1,NL DDUO3N(N,KK) = RDATA(KK) * .01 END DO ! 75 END DO !--------------------- ! END OF LATITUDE LOOP !--------------------- ! !----------------------------------------------------------------------- ! CREATE 5 DEG OZONE QUANTITIES BY LINEAR INTERPOLATION OF 10 DEG VALUES !----------------------------------------------------------------------- DO 90 KK=1,NL ! DO N=1,19 O35DEG(2*N-1,KK) = DDUO3N(N,KK) END DO ! DO N=1,18 O35DEG(2*N ,KK) = 0.5 * (DDUO3N(N,KK) + DDUO3N(N+1,KK)) END DO ! 90 CONTINUE ! DO JJ=1,37 DO KEN=1,NL O3O3(JJ,KEN,IPLACE) = O35DEG(JJ,KEN) END DO END DO ! 100 CONTINUE !----------------------- ! END OF LOOP OVER CASES !----------------------- ! !-------------------------------------------------------------------------------------------------- ! AVERAGE CLIMATOLOGICAL VALUS OF O3 FROM 5 DEG LAT MEANS, SO THAT TIME AND SPACE INTERPOLATION ! WILL WORK (SEE SUBR OZON2D). !-------------------------------------------------------------------------------------------------- DO I=1,NLGTH AVG = 0.25 * ( XRAD1(I) + XRAD2(I) + XRAD3(I) + XRAD4(I)) A1 = 0.50 * ( XRAD2(I) - XRAD4(I)) B1 = 0.50 * ( XRAD1(I) - XRAD3(I)) B2 = 0.25 * ((XRAD1(I) + XRAD3(I)) - (XRAD2(I) + XRAD4(I))) ! XRAD1(I) = AVG XRAD2(I) = A1 XRAD3(I) = B1 XRAD4(I) = B2 ! END DO !---------------------------------------- ! CONVERT GFDL PRESSURE (MICROBARS) TO PA !---------------------------------------- DO N=1,NL PRGFDL(N) = PSTD(N) * 1.E-1 END DO ! RETURN ! END SUBROUTINE O3CLIM