!--------------------------------------- LICENCE BEGIN -----------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
! version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
module ozoneclim 1,1
!************************************************************************
!*
!* MODULE OZONECLIM
!* -----------------
!*
!* PURPOSE: CLIMATOLOGICAL OZONE (1998)
!*
!* AUTHOR: A. BEAULNE (CMDA/SMC) May 2006
!*
!* REVISION:
!*
!************************************************************************
use obsSpaceData_mod
implicit none
save
private
! public procedures
public :: OZO_GET_PROFILE, OZO_READ_CLIMATOLOGY
! Number of latitudes and vertical levels in climatology file
INTEGER, PARAMETER :: NLATO3=19, NLEVO3=28
! Climatological ozone field (ppmv) and total ozone
REAL :: FOZO_r4(NLATO3,NLEVO3), TOTOZO_r4(NLATO3,12)
! Pressure height of climatology file vertical levels (mb)
REAL(8) :: PO3(NLEVO3)
DATA PO3 / 0.010D0, 0.015D0, 0.022D0, 0.032D0, 0.046D0, 0.068D0, 0.100D0, &
0.150D0, 0.200D0, 0.300D0, 0.500D0, 1.000D0, 2.000D0, 3.000D0, &
5.000D0, 7.000D0, 10.00D0, 20.00D0, 30.00D0, 50.00D0, 70.00D0, &
100.0D0, 150.0D0, 200.0D0, 300.0D0, 500.0D0, 700.0D0, 1000.D0 /
contains
SUBROUTINE OZO_GET_PROFILE(O3P,toto3obs,zlat,plev,nlev,nprf,datestamp) 1
!***********************************************************************
!
!**ID OZO_GET_PROFILE -- GET OZONE PROFILE FROM CLIMATOLOGY INTERPOLATED TO DESIRED P LEVELS
!
! AUTHOR: A. BEAULNE (CMDA/SMC) March 2006
!
! REVISION:
!
! OBJECT: GET OZONE PROFILE FROM CLIMATOLOGY INTERPOLATED TO DESIRED P LEVELS
!
! ARGUMENTS:
! INPUT:
! -TOTO3OBS(NPRF) : TOTAL OZONE IF KNOWN FROM OUTSIDE SOURCE SUCH AS TOMS
! -ZLAT(NPRF) : ARRAY OF LATITUDE (-90S TO 90N)
! -PLEV(NLEV) : PRESSURE LEVELS (HPA)
! -NLEV : NUMBER OF VERTICAL LEVELS
! -NPRF : NUMBER OF PROFILES
!
! OUTPUT:
! -O3P(NLEV,NPRF) : OZONE PROFILES (PPMV)
!
!
!***********************************************************************
IMPLICIT NONE
integer ,intent(in) :: nlev,nprf
REAL(8),intent(in) :: ZLAT(NPRF), TOTO3OBS(NPRF),PLEV(NLEV)
REAL(8),intent(out) :: O3P(NLEV,NPRF)
integer, intent(in) :: datestamp
!***************************************************************
INTEGER :: JN, K, NUMLAT(NPRF)
INTEGER :: IJOUR,ITIME,IJ,IMONTH,IER
REAL(8) :: QO3B(NLEVO3,NPRF)
REAL(8) :: PRO3(NLEVO3,NPRF)
INTEGER :: NEWDATE
!* assign default qgas values if need be
DO JN = 1, NPRF
NUMLAT(JN) = NINT( (ZLAT(JN)+90.D0) / (180.D0/(REAL(NLATO3-1,8))) ) + 1
DO K = 1, NLEVO3
QO3B(K,JN) = FOZO_r4(NUMLAT(JN),K)
END DO
END DO
!* interpolation of field QO3B at NLEVO3 levels of height PO3mbb
!* into field O3P at NLEV levels of height PLEV
FORALL(K=1:NLEVO3) PRO3(K,:) = PO3(K)
CALL LINTV2(pro3,qo3b,nlevo3,nlevo3,nprf,nlev,plev,O3P)
!* if total climatological ozone is known from outside source
!* then set the climatological profile so that it matches that total
DO JN = 1, NPRF
IF ( NINT(TOTO3OBS(JN)) /= 0 ) THEN
! IJOUR = obs_headElem_i(lobsSpaceData,OBS_DAT,1)
ier = newdate(datestamp,ijour,itime,-3)
IJ = IJOUR/100
IMONTH = IJ - (IJ/100)*100
O3P(:,JN) = O3P(:,JN) * TOTO3OBS(JN) / TOTOZO_r4(NUMLAT(JN),IMONTH)
END IF
ENDDO
END SUBROUTINE OZO_GET_PROFILE
SUBROUTINE OZO_READ_CLIMATOLOGY(datestamp) 1,1
!***********************************************************************
!
!**ID OZONECLIM -- READ OZONE CLIMATOLOGICAL FIELDS
!
! AUTHOR: A. BEAULNE (CMDA/SMC) March 2006
!
! REVISION:
!
! OBJECT: READ OZONE CLIMATOLOGICAL FIELDS AND PUT IN COMMON DECK
!
! ARGUMENTS:
! INPUT: NONE
! OUTPUT: NONE
!
!
!***********************************************************************
IMPLICIT NONE
integer :: datestamp
INTEGER :: IJOUR,ITIME,IMONTH,IJ,IV,K,IER
CHARACTER(len=100) :: CFILE
INTEGER :: NIOZO,NJOZO,NKOZO
INTEGER, EXTERNAL :: FNOM,FSTOUV,FSTLIR,FSTFRM,FCLOS,NEWDATE
integer :: IOZTEST
integer :: iv1,iv2,iv3,iv4,iv5,iv6
! IJOUR = obs_headElem_i(lobsSpaceData,OBS_DAT,1)
ier = newdate(datestamp,ijour,itime,-3)
IJ= IJOUR/100
IMONTH = IJ - (IJ/100)*100
ioztest=0
CFILE='ozoneclim98'
IV1=FNOM(IOZTEST,CFILE,'RND+R/O',0)
IV2=FSTOUV(IOZTEST,'RND')
IV3=FSTLIR(FOZO_r4,IOZTEST,NIOZO,NJOZO,NKOZO,-1,' ',-1,-1,IMONTH,' ','O3')
IV4=FSTLIR(TOTOZO_r4,IOZTEST,NIOZO,NJOZO,NKOZO,-1,' ',-1,-1,-1,' ','TO')
IV5=FSTFRM(IOZTEST)
IV6=FCLOS(IOZTEST)
if(iv1.lt.0.or.iv2.lt.0.or.iv3.lt.0.or.iv4.lt.0.or.iv5.lt.0.or.iv6.lt.0) then
write(*,*) 'LES IV DE OZO_READ_CLIMATOLOGY ',iv1,iv2,iv3,iv4,iv5,iv6
write(*,*) 'THESE NUMBERS SHOULD NOT BE NEGATIVE'
write(*,*) 'datestamp,ijour,itime,imonth = ',datestamp,ijour,itime,imonth
call abort3d('Problem with file in ozo_read_climatology (ozoneclim_mod)')
endif
END SUBROUTINE OZO_READ_CLIMATOLOGY
End module ozoneclim