radozc.F90

Go to the documentation of this file.

SUBROUTINE radozc ( KIDIA , KFDIA , KLON , KTDIA , KLEV &
                  &, krint , kdlon , kshift &
                  &, paprs , pgemu &
                  &, pozon                 )

!***********************************************************************
! CAUTION: THIS ROUTINE WORKS ONLY ON A NON-ROTATED, UNSTRETCHED GRID
!***********************************************************************

!**** *RADOZC* - COMPUTES DISTRIBUTION OF OZONE FROM CLIMATOLOGY

!     PURPOSE.
!     --------

!**   INTERFACE.
!     ----------
!        CALL *RADOZC* FROM *RADINT*

!        EXPLICIT ARGUMENTS :
!        --------------------
!     ==== INPUTS ===
!     ==== OUTPUTS ===

!        IMPLICIT ARGUMENTS :   NONE
!        --------------------

!     METHOD.
!     -------


!     EXTERNALS.
!     ----------

!          NONE

!     REFERENCE.
!     ----------

!        SEE RADIATION'S PART OF THE MODEL'S DOCUMENTATION AND
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE "I.F.S"

!     AUTHOR.
!     -------
!     J.-J. MORCRETTE  E.C.M.W.F.    95/01/25

!     MODIFICATIONS.
!     --------------
!          D.SALMOND   ECMWF    99-06-14   Optimisation

!-----------------------------------------------------------------------

#include "tsmbkind.h"

USE yomcst   , ONLY : r
USE yoeozoc  , ONLY : rsinc    ,rozt     ,rproc


IMPLICIT NONE


!     DUMMY INTEGER SCALARS
integer_m :: kdlon
integer_m :: kfdia
integer_m :: kidia
integer_m :: klev
integer_m :: klon
integer_m :: krint
integer_m :: kshift
integer_m :: ktdia




!     -----------------------------------------------------------------

!*       0.1   ARGUMENTS.
!              ----------

real_b :: paprs(klon,klev+1),pgemu(klon)

real_b :: pozon(kdlon,klev)

!     ----------------------------------------------------------------- 

!*       0.2   LOCAL ARRAYS.
!              -------------

real_b :: zozlt(0:35) , zozon(kdlon,klev+1)
real_b :: zrrr(0:34)

!     LOCAL INTEGER SCALARS
integer_m :: il, inla, jc, jk, jl

!     LOCAL REAL SCALARS
real_b :: zpmr, zsilat, zsin

!     ------------------------------------------------------------------
!     ------------------------------------------------------------------
!     ------------------------------------------------------------------

!*         1.     LATITUDE INDEX WITHIN OZONE CLIMATOLOGY
!                 ---------------------------------------



zsin=pgemu(kidia)
inla=0
zsilat=-9999._jprb
DO jl=18,1,-1
  IF (zsin <= rsinc(jl+1).AND.zsin >= rsinc(jl)) THEN
    inla=jl
  ENDIF
ENDDO
IF (inla == 0) THEN
!  CALL ABOR1(' Problem with lat. interpolation in radozc!')
  print *,' Problem with lat. interpolation in radozc!'
ENDIF
zsilat=(zsin-rsinc(inla))/(rsinc(inla+1)-rsinc(inla))

!     ------------------------------------------------------------------


!*         2.     LATITUDE INTERPOLATED FIELD
!                 ---------------------------


IF(inla == 18.OR.inla == 1) THEN
  DO jc=0,35
    zozlt(jc)=rozt(inla,jc)
  ENDDO
ELSE
  DO jc=0,35
    zozlt(jc)=rozt(inla,jc)+zsilat*(rozt(inla+1,jc)-rozt(inla,jc))
  ENDDO
ENDIF


!     ------------------------------------------------------------------

!*         3.     VERTICAL INTERPOLATION 
!                 ----------------------

DO jc=0,34
  zrrr(jc)=(_one_/(rproc(jc)-rproc(jc+1)))*(zozlt(jc)-zozlt(jc+1))
ENDDO


il=kshift
DO jl=kidia,kfdia,krint
  il=il+1
  DO jc=0,34
    DO jk=1,klev+1
    zpmr=paprs(jl,jk)
      IF(zpmr >= rproc(jc).AND.zpmr < rproc(jc+1)) &
&        zozon(il,jk)=zozlt(jc+1)+(zpmr-rproc(jc+1))*zrrr(jc)
    ENDDO
  ENDDO
ENDDO

il=kshift
DO jl=kidia,kfdia,krint
  il=il+1
  DO jk=1,klev+1
    zpmr=paprs(jl,jk)
    zpmr=paprs(jl,jk)
    IF(zpmr >= rproc(35)) zozon(il,jk)=zozlt(35)
  ENDDO
ENDDO

! INTEGRATION IN THE VERTICAL:
il=kshift
DO jl=kidia,kfdia,krint
  il=il+1
  DO jk=1,klev
    pozon(il,jk)=(paprs(jl,jk+1)-paprs(jl,jk))&
     &*(zozon(il,jk)+zozon(il,jk+1))*_half_
  ENDDO
ENDDO

!     -----------------------------------------------------------

RETURN
END SUBROUTINE radozc