[3990] | 1 | SUBROUTINE VAL923(LDNEW) |
---|
| 2 | |
---|
| 3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
| 4 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
| 5 | |
---|
| 6 | USE YOMCLI , ONLY : NTPMER ,NTPGLA ,NTPDES ,NTPLAC ,& |
---|
| 7 | & SMASK ,SMANQ ,STHER ,SALBN ,SALBX ,& |
---|
| 8 | & SALBM ,SALBG ,SALBB ,SALBD ,SEMIN ,& |
---|
| 9 | & SEMIX ,SEMIM ,SEMIG ,SEMIB ,SEMID ,& |
---|
| 10 | & SDEPN ,SDEPX ,SDEPD ,SARGN ,SARGX ,& |
---|
| 11 | & SARGD ,SSABN ,SSABX ,SSABD ,SRSMN ,& |
---|
| 12 | & SRSMX ,SRSMD ,SZZ0N ,SZZ0M ,SZZ0B ,& |
---|
| 13 | & SZZ0U ,SZZ0D |
---|
| 14 | USE YOMLUN , ONLY : NULOUT |
---|
| 15 | |
---|
| 16 | #ifdef DOC |
---|
| 17 | !**** *GEO923* |
---|
| 18 | |
---|
| 19 | ! PURPOSE. |
---|
| 20 | ! -------- |
---|
| 21 | ! Compute the constants (YOMCLI) which are used by configuration 923. |
---|
| 22 | |
---|
| 23 | ! INTERFACE. |
---|
| 24 | ! ---------- |
---|
| 25 | ! CALL VAL923(LDNEW) |
---|
| 26 | ! LDNEW = .FALSE. if old fields required |
---|
| 27 | ! Results in YOMCLI. |
---|
| 28 | |
---|
| 29 | ! AUTHORS. |
---|
| 30 | ! -------- |
---|
| 31 | ! D. Giard 97-05-06 |
---|
| 32 | |
---|
| 33 | ! MODIFICATIONS. |
---|
| 34 | ! -------------- |
---|
| 35 | ! D. Giard 99-07-02 add values for desert and prints |
---|
| 36 | ! S. Kertesz 99-12-16 new type for lakes |
---|
| 37 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning |
---|
| 38 | |
---|
| 39 | #endif |
---|
| 40 | |
---|
| 41 | IMPLICIT NONE |
---|
| 42 | |
---|
| 43 | LOGICAL ,INTENT(IN) :: LDNEW |
---|
| 44 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
| 45 | |
---|
| 46 | ! Threshold defining the mask |
---|
| 47 | IF (LHOOK) CALL DR_HOOK('VAL923',0,ZHOOK_HANDLE) |
---|
| 48 | SMASK= 0.5_JPRB |
---|
| 49 | ! Value for missing data + 1 |
---|
| 50 | SMANQ=-9998._JPRB |
---|
| 51 | ! Land-use types for sea, ice-cap, desert, lakes |
---|
| 52 | NTPMER= 1 |
---|
| 53 | NTPGLA= 2 |
---|
| 54 | NTPDES= 3 |
---|
| 55 | NTPLAC= 5 |
---|
| 56 | ! Roughness length : minimum, sea, sea-ice, urban areas, desert |
---|
| 57 | SZZ0N= 0.001_JPRB |
---|
| 58 | SZZ0M= 0.001_JPRB |
---|
| 59 | SZZ0B= 0.001_JPRB |
---|
| 60 | SZZ0U= 2.500_JPRB |
---|
| 61 | SZZ0D= 0.001_JPRB |
---|
| 62 | ! Ration of thermal to kinetic roughness length |
---|
| 63 | STHER= 0.10_JPRB |
---|
| 64 | ! Albedo : minimum, maximum, sea, ice-cap, sea-ice, desert |
---|
| 65 | IF (LDNEW) THEN |
---|
| 66 | SALBN= 0.05_JPRB |
---|
| 67 | SALBX= 0.80_JPRB |
---|
| 68 | ELSE |
---|
| 69 | SALBN= 0.07_JPRB |
---|
| 70 | SALBX= 0.70_JPRB |
---|
| 71 | ENDIF |
---|
| 72 | SALBM= 0.07_JPRB |
---|
| 73 | SALBG= 0.75_JPRB |
---|
| 74 | SALBB= 0.65_JPRB |
---|
| 75 | SALBD= 0.10_JPRB |
---|
| 76 | ! Emissivity : minimum, maximum, sea, ice-cap, sea-ice, desert |
---|
| 77 | SEMIN= 0.90_JPRB |
---|
| 78 | SEMIX= 1.00_JPRB |
---|
| 79 | SEMIM= 0.96_JPRB |
---|
| 80 | SEMIG= 0.98_JPRB |
---|
| 81 | SEMIB= 0.97_JPRB |
---|
| 82 | SEMID= 0.943_JPRB |
---|
| 83 | ! Soil depth : minimum, maximum, desert |
---|
| 84 | SDEPN= 0.10_JPRB |
---|
| 85 | SDEPX= 8.00_JPRB |
---|
| 86 | SDEPD= 0.10_JPRB |
---|
| 87 | ! Percentage of clay : minimum, maximum, desert |
---|
| 88 | SARGN= 3._JPRB |
---|
| 89 | SARGX= 58._JPRB |
---|
| 90 | SARGD= 3._JPRB |
---|
| 91 | ! Percentage of sand : minimum, maximum, desert |
---|
| 92 | SSABN= 6._JPRB |
---|
| 93 | SSABX= 92._JPRB |
---|
| 94 | SSABD= 92._JPRB |
---|
| 95 | ! Minimum surface resistance : minimum, maximum, desert |
---|
| 96 | SRSMX=5000._JPRB |
---|
| 97 | SRSMN= 1.0_JPRB |
---|
| 98 | SRSMD=5000._JPRB |
---|
| 99 | |
---|
| 100 | WRITE(UNIT=NULOUT,FMT=111) SMASK,SMANQ,STHER,& |
---|
| 101 | & NTPMER,NTPGLA,NTPDES,NTPLAC |
---|
| 102 | WRITE(UNIT=NULOUT,FMT=112) SZZ0N,SZZ0M,SZZ0B,SZZ0U,SZZ0D |
---|
| 103 | WRITE(UNIT=NULOUT,FMT=113) SALBN,SALBX,SALBM,SALBG,SALBB,SALBD,& |
---|
| 104 | & SEMIN,SEMIX,SEMIM,SEMIG,SEMIB,SEMID |
---|
| 105 | WRITE(UNIT=NULOUT,FMT=114) SDEPN,SDEPX,SDEPD,SARGN,SARGX,SARGD,& |
---|
| 106 | & SSABN,SSABX,SSABD,SRSMN,SRSMX,SRSMD |
---|
| 107 | 111 FORMAT(' COMMON YOMCLI',/,& |
---|
| 108 | & ' SMASK=',F4.2,' SMANQ=',F6.0,' STHER=',F4.2,/& |
---|
| 109 | & ' NTPMER=',I2,' NTPGLA=',I2,' NTPDES=',I2,' NTPLAC=',I2) |
---|
| 110 | 112 FORMAT(' LONGUEUR DE RUGOSITE :',/,& |
---|
| 111 | & ' minimum mer banquise villes desert ',& |
---|
| 112 | & /,5F9.3) |
---|
| 113 | 113 FORMAT(' ALBEDO ET EMISSIVITE :',/,& |
---|
| 114 | & ' minimum maximum mer glacier banquise desert ',& |
---|
| 115 | & 2(/,6F9.3)) |
---|
| 116 | 114 FORMAT(' PROFONDEUR, % ARGILE, % SABLE, RESIS. MIN. :',/,& |
---|
| 117 | & ' minimum maximum desert ',& |
---|
| 118 | & 4(/,3F9.3)) |
---|
| 119 | |
---|
| 120 | IF (LHOOK) CALL DR_HOOK('VAL923',1,ZHOOK_HANDLE) |
---|
| 121 | END SUBROUTINE VAL923 |
---|