slope_radiative_effect.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
SUBROUTINE slope_radiative_effect(PTSTEP,PZENITH,PAZIM,PPS,PTA,PRAIN,PDIR_SW,PLW, &
                                  PZS,PZS_XY,PSLOPANG,PSLOPAZI,PSURF_TRIANGLE     )
!##############################################################
!
!!**** *SLOPE_RADIATIVE_EFFECT* compute direct short-wave radiation modified by slopes and shadows, 
!                               BUT renormalized on the horizontal surface of the grid mesh  
!                               to serve as input for ISBA
!
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!     Based on original package to correct shortwave radiation included in
!!     MESO-NH (developed by V. Masson)
!!
!!    AUTHOR
!!    ------
!!
!!    V. Vionnet        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    04/11
!!             29/04/11  : VV implementation of adjustment of incoming longwave 
!!                         radiation as a function of the slope (routine
!!                         originally implemented in operational chain SCM) 
!!                         (routine meteo.f90 of Crocus) 
!!                03/14  : M Lafaysse, modifs for optimization and parallelization
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modi_surf_rad_modif
!
USE modd_csts,       ONLY : xpi,xstefan
USE modd_slope_effect
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
REAL,                 INTENT(IN) :: PTSTEP   ! surface time step (s)
REAL, DIMENSION(:),   INTENT(IN) :: PZENITH  ! zenithal angle       (radian from the vertical)
REAL, DIMENSION(:),   INTENT(IN) :: PAZIM    ! azimuthal angle      (radian from North, clockwise)
REAL, DIMENSION(:),   INTENT(IN) :: PPS      ! pressure at atmospheric model surface (Pa)
REAL, DIMENSION(:),   INTENT(IN) :: PTA      ! air temperature forcing               (K)
REAL, DIMENSION(:),   INTENT(IN) :: PRAIN    ! liquid precipitation                  (kg/m2/s)

! OROGRAPHY OF THE SUBDOMAIN (THIS PROCESSOR)
REAL, DIMENSION(:,:),     INTENT(IN) :: PZS        ! (resolved) model orography
REAL, DIMENSION(:,:),     INTENT(IN) :: PZS_XY     ! orography at vort. points
REAL, DIMENSION(:,:,:),     INTENT(IN) :: PSLOPAZI   ! azimuthal slope angle of triangles
REAL, DIMENSION(:,:,:),     INTENT(IN) :: PSLOPANG   ! vertical slope angle of triangles
REAL, DIMENSION(:,:,:),     INTENT(IN) :: PSURF_TRIANGLE ! surface of triangles

REAL, DIMENSION(:,:),   INTENT(INOUT) :: PDIR_SW  !   IN : input down (direct) short-wave radiation
!                                                   OUT : down (direct) short-wave radiation modified by slopes and shadows, 
!                                                        BUT renormalized on the horizontal surface of the grid mesh
REAL, DIMENSION(:),   INTENT(INOUT) :: PLW      !   IN : longwave radiation (on horizontal surf.)
!                                                   OUT : longwave radiation modified by slopee and valley effects, 
!                                                        BUT renormalized on the horizontal surface of the grid mesh
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER :: JX       ! loop counter
INTEGER :: JY       ! loop counter
INTEGER :: JB       ! loop counter
INTEGER :: INBANDS   ! number of radiative bands
INTEGER :: IINDY

REAL, DIMENSION(:,:), ALLOCATABLE :: ZMAP   ! map factor
REAL, DIMENSION(:,:), ALLOCATABLE :: ZCOSZEN ! cosine of solar zenithal angle (=1 at zenith)
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSINZEN ! sinus of solar zenithal angle (=0 at zenith)
REAL, DIMENSION(:,:), ALLOCATABLE :: ZAZIMSOL ! solar azimuthal angle
!                                 ! 0     ==> Sun from the South
!                                 ! Pi/2  ==> Sun from the East
!                                 ! Pi    ==> Sun from the North
!                                 ! 3Pi/2 ==> Sun from the West
! REAL, DIMENSION(:,:), ALLOCATABLE :: ZSSO_SURF    ! ratio between sloping orography and horizontal surface
! REAL, DIMENSION(:), ALLOCATABLE :: ZSSO_SURF_1D   ! ratio between sloping orography and horizontal surface (1D vector)
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZDIRFLASWD ! input down (direct) short-wave radiation
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZDIRSRFSWD ! down (direct) short-wave radiation modified by slopes and shadows, 
!                                      ! BUT renormalized on the horizontal surface of the grid mesh
!                                      ! For both variables, 3rd dimension is wavelength
REAL, DIMENSION(:), ALLOCATABLE :: ZALPHA ! intermediate variable 
!
REAL, PARAMETER :: VPRES1 = 87000.
!

inbands = SIZE(pdir_sw,2)

!-------------------------------------------------------------------------------
!
!*    1.    Downwards solar radiation and solar angles
!     ------------------------------------------
!
! PRINT*,"SIZE(PTA)=",SIZE(PTA)

!
! solar zenithal angle    (cosine=1 : zenith)
ALLOCATE(zcoszen(nnxloc,nnyloc))
ALLOCATE(zsinzen(nnxloc,nnyloc))
ALLOCATE(zazimsol(nnxloc,nnyloc))
ALLOCATE(zmap(nnxloc,nnyloc))
! ALLOCATE(ZSSO_SURF(NNX,NNY))
ALLOCATE(zdirflaswd(nnxloc,nnyloc,inbands))
ALLOCATE(zdirsrfswd(nnxloc,nnyloc,inbands))
! ALLOCATE(ZSSO_SURF_1D(IX*IY))
ALLOCATE(zalpha(nix*niy))

! PRINT*,SIZE(PZENITH)
! PRINT*,NNX,NNY,IX,IY

!PRINT*,NIYLOC,NIXLOC

IF (lrevertgrid) THEN
  DO jy=1,niyloc
    iindy=niyloc-jy+1
    DO jx=1,nixloc
      zcoszen(jx+1,iindy+1) = cos(pzenith( jx + (jy-1)*nixloc )) 
      zsinzen(jx+1,iindy+1) = sin(pzenith( jx + (jy-1)*nixloc )) 
      zazimsol(jx+1,iindy+1) = pazim( jx + (jy-1)*nixloc )      
    END DO
  END DO
ELSE
  DO jy=1,niyloc
    DO jx=1,nixloc
      zcoszen(jx+1,jy+1) = cos(pzenith( jx + (jy-1)*nixloc )) 
      zsinzen(jx+1,jy+1) = sin(pzenith( jx + (jy-1)*nixloc ))    
      zazimsol(jx+1,jy+1) = pazim( jx + (jy-1)*nixloc )
    END DO
  END DO
ENDIF


zcoszen(1,:) = zcoszen(2,:)
zcoszen(nnxloc,:) = zcoszen(nnxloc-1,:)
zcoszen(:,1) = zcoszen(:,2)
zcoszen(:,nnyloc) = zcoszen(:,nnyloc-1)
zsinzen(1,:) = zsinzen(2,:)
zsinzen(nnxloc,:) = zsinzen(nnxloc-1,:)
zsinzen(:,1) = zsinzen(:,2)
zsinzen(:,nnyloc) = zsinzen(:,nnyloc-1)
zazimsol(1,:) = zazimsol(2,:)
zazimsol(nnxloc,:) = zazimsol(nnxloc-1,:)
zazimsol(:,1) = zazimsol(:,2)
zazimsol(:,nnyloc) = zazimsol(:,nnyloc-1)
!
! Downwards solar radiation
IF (lrevertgrid) THEN
  DO jb = 1,inbands
    DO jy=1,niyloc
      iindy=niyloc-jy+1
      DO jx=1,nixloc
        zdirflaswd(jx+1,iindy+1,jb) = pdir_sw( jx + (jy-1)*nixloc,jb)
      END DO
    END DO
  END DO
ELSE
  DO jb = 1,inbands
    DO jy=1,niyloc
      DO jx=1,nixloc
        zdirflaswd(jx+1,jy+1,jb) = pdir_sw( jx + (jy-1)*nixloc,jb)
      END DO
    END DO
  END DO
ENDIF


zdirflaswd(1,:,:) = zdirflaswd(2,:,:)
zdirflaswd(nnxloc,:,:) = zdirflaswd(nnxloc-1,:,:)
zdirflaswd(:,1,:) = zdirflaswd(:,2,:)
zdirflaswd(:,nnyloc,:) = zdirflaswd(:,nnyloc-1,:)

!
!  Map factor
!
zmap(:,:) = 1.

!------------------------------------------------------------------------------------------
!
!*        2. Calls radiative computations
!     ----------------------------
!
! PRINT*,XZSL
! PRINT*,"in slope radiative effect",SIZE(XZSL,1),SIZE(XZSL,2)

 CALL  surf_rad_modif ( zmap, xxhat_thread, xyhat_thread,                    &
                       zcoszen, zsinzen, zazimsol,pzs,pzs_xy,               &
                       pslopang,pslopazi,psurf_triangle,                    &
                       xxhat,xyhat,nindx1_x,nindx1_y,xzsl,xzs_xy,           &
                       zdirflaswd, zdirsrfswd                    )
!
!
!-------------------------------------------------------------------------------
!
!*    3.     Output field comes back into 1D vector
!            --------------------------------------
!

IF (lrevertgrid) THEN
  DO jb=1,inbands
    DO jy=1,niyloc
      iindy=niyloc-jy+1
      DO jx=1,nixloc
        pdir_sw( jx + (jy-1)*nixloc,jb ) = zdirsrfswd(jx+1,iindy+1,jb)
!     ZSSO_SURF_1D(JX + (JY-1)*IX)  = ZSSO_SURF(JX+1,JY+1)
      END DO
    END DO
  END DO
ELSE
  DO jb=1,inbands
    DO jy=1,niyloc
      DO jx=1,nixloc
        pdir_sw( jx + (jy-1)*nixloc,jb ) = zdirsrfswd(jx+1,jy+1,jb)
!     ZSSO_SURF_1D(JX + (JY-1)*IX)  = ZSSO_SURF(JX+1,JY+1)
      END DO
    END DO
  END DO
ENDIF
!
!
!
!-------------------------------------------------------------------------------
!
! !*    4.     Modify longwave incoming radiation due to account for the influence
! !            of opposite slope and correction in case of rain
! !            --------------------------------------------------------
! ZALPHA(:) = MAX(0.25*(ZSSO_SURF_1D(:)/SQRT(1+ZSSO_SURF_1D(:)**2.))**2., 0.1*PPS(:)/VPRES1)
! 
! PLW(:) = (1-ZALPHA(:))*PLW(:) +ZALPHA(:)*XSTEFAN*PTA(:)**4.
! 
! WHERE(PRAIN(:)*PTSTEP>0.001)
!         PLW(:) = MAX(PLW(:),0.95*XSTEFAN*PTA(:)**4.)
! END WHERE


DEALLOCATE(zcoszen)
DEALLOCATE(zsinzen)
DEALLOCATE(zazimsol)
DEALLOCATE(zmap)
! DEALLOCATE(ZSSO_SURF)
DEALLOCATE(zdirflaswd)
DEALLOCATE(zdirsrfswd)
! DEALLOCATE(ZSSO_SURF_1D)
DEALLOCATE(zalpha)
END SUBROUTINE slope_radiative_effect