val923.F90

Go to the documentation of this file.

SUBROUTINE val923(LDNEW)

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

USE yomcli   , ONLY : ntpmer   ,ntpgla   ,ntpdes   ,ntplac   ,&
 & smask    ,smanq    ,sther    ,salbn    ,salbx    ,&
 & salbm    ,salbg    ,salbb    ,salbd    ,semin    ,&
 & semix    ,semim    ,semig    ,semib    ,semid    ,&
 & sdepn    ,sdepx    ,sdepd    ,sargn    ,sargx    ,&
 & sargd    ,ssabn    ,ssabx    ,ssabd    ,srsmn    ,&
 & srsmx    ,srsmd    ,szz0n    ,szz0m    ,szz0b    ,&
 & szz0u    ,szz0d  
USE yomlun   , ONLY : nulout

#ifdef DOC
!**** *GEO923*

!     PURPOSE.
!     --------
!      Compute the constants (YOMCLI) which are used by configuration 923.

!     INTERFACE.
!     ----------
!      CALL VAL923(LDNEW)
!        LDNEW = .FALSE. if old fields required
!        Results in YOMCLI.

!     AUTHORS.
!     --------
!      D. Giard 97-05-06

!     MODIFICATIONS.
!     --------------
!      D. Giard   99-07-02 add values for desert and prints
!      S. Kertesz 99-12-16 new type for lakes
!        M.Hamrud      01-Oct-2003 CY28 Cleaning

#endif

IMPLICIT NONE

LOGICAL           ,INTENT(IN)    :: LDNEW 
REAL(KIND=JPRB) :: ZHOOK_HANDLE

!  Threshold defining the mask
IF (lhook) CALL dr_hook('VAL923',0,zhook_handle)
smask= 0.5_jprb
!  Value for missing data + 1
smanq=-9998._jprb
!  Land-use types for sea, ice-cap, desert, lakes
ntpmer= 1
ntpgla= 2
ntpdes= 3
ntplac= 5
!  Roughness length : minimum, sea, sea-ice, urban areas, desert
szz0n= 0.001_jprb
szz0m= 0.001_jprb
szz0b= 0.001_jprb
szz0u= 2.500_jprb
szz0d= 0.001_jprb
!  Ration of thermal to kinetic roughness length
sther= 0.10_jprb
!  Albedo : minimum, maximum, sea, ice-cap, sea-ice, desert
IF (ldnew) THEN
  salbn= 0.05_jprb
  salbx= 0.80_jprb
ELSE
  salbn= 0.07_jprb
  salbx= 0.70_jprb
ENDIF
salbm= 0.07_jprb
salbg= 0.75_jprb
salbb= 0.65_jprb
salbd= 0.10_jprb
!  Emissivity : minimum, maximum, sea, ice-cap, sea-ice, desert
semin= 0.90_jprb
semix= 1.00_jprb
semim= 0.96_jprb
semig= 0.98_jprb
semib= 0.97_jprb
semid= 0.943_jprb
!  Soil depth : minimum, maximum, desert
sdepn= 0.10_jprb
sdepx= 8.00_jprb
sdepd= 0.10_jprb
!  Percentage of clay : minimum, maximum, desert
sargn=  3._jprb
sargx= 58._jprb
sargd=  3._jprb
!  Percentage of sand : minimum, maximum, desert
ssabn=  6._jprb
ssabx= 92._jprb
ssabd= 92._jprb
!  Minimum surface resistance : minimum, maximum, desert
srsmx=5000._jprb
srsmn=   1.0_jprb
srsmd=5000._jprb

WRITE(unit=nulout,fmt=111) smask,smanq,sther,&
 & ntpmer,ntpgla,ntpdes,ntplac  
WRITE(unit=nulout,fmt=112) szz0n,szz0m,szz0b,szz0u,szz0d
WRITE(unit=nulout,fmt=113) salbn,salbx,salbm,salbg,salbb,salbd,&
 & semin,semix,semim,semig,semib,semid  
WRITE(unit=nulout,fmt=114) sdepn,sdepx,sdepd,sargn,sargx,sargd,&
 & ssabn,ssabx,ssabd,srsmn,srsmx,srsmd  
111 FORMAT(' COMMON YOMCLI',/,&
 & ' SMASK=',f4.2,' SMANQ=',f6.0,' STHER=',f4.2,/&
 & ' NTPMER=',i2,' NTPGLA=',i2,' NTPDES=',i2,' NTPLAC=',i2)  
112 FORMAT(' LONGUEUR DE RUGOSITE :',/,&
 & ' minimum    mer    banquise  villes   desert ',&
 & /,5f9.3)  
113 FORMAT(' ALBEDO ET EMISSIVITE :',/,&
 & ' minimum  maximum    mer    glacier  banquise  desert ',&
 & 2(/,6f9.3))  
114 FORMAT(' PROFONDEUR, % ARGILE, % SABLE, RESIS. MIN. :',/,&
 & ' minimum  maximum   desert ',&
 & 4(/,3f9.3))  

IF (lhook) CALL dr_hook('VAL923',1,zhook_handle)
END SUBROUTINE val923