[219] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
| 4 | |
---|
[373] | 5 | subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, & |
---|
| 6 | & fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, & |
---|
| 7 | & ratqsbas,ratqshaut) |
---|
[219] | 8 | |
---|
[235] | 9 | use IOIPSL |
---|
| 10 | implicit none |
---|
| 11 | |
---|
[373] | 12 | #include "conema3.h" |
---|
| 13 | #include "fisrtilp.inc" |
---|
[386] | 14 | #include "nuage.h" |
---|
[219] | 15 | ! |
---|
[241] | 16 | ! Configuration de la "physique" de LMDZ a l'aide de la fonction |
---|
[219] | 17 | ! GETIN de IOIPSL |
---|
| 18 | ! |
---|
| 19 | ! LF 05/2001 |
---|
| 20 | ! |
---|
| 21 | |
---|
| 22 | ! |
---|
| 23 | ! ocean: type d'ocean (force, slab, couple) |
---|
| 24 | ! ok_veget: type de modele de vegetation |
---|
| 25 | ! ok_journe: sorties journalieres |
---|
| 26 | ! ok_mensuel: sorties mensuelles |
---|
| 27 | ! ok_instan: sorties instantanees |
---|
| 28 | |
---|
| 29 | |
---|
| 30 | |
---|
| 31 | ! Sortie: |
---|
| 32 | character (len = 6) :: ocean |
---|
[373] | 33 | logical :: ok_veget, ok_newmicro |
---|
[219] | 34 | logical :: ok_journe, ok_mensuel, ok_instan |
---|
[373] | 35 | real :: fact_cldcon, facttemps,ratqsbas,ratqshaut |
---|
| 36 | integer :: iflag_cldcon |
---|
[223] | 37 | ! Local |
---|
| 38 | integer :: numout = 6 |
---|
[219] | 39 | |
---|
| 40 | ! |
---|
| 41 | ! |
---|
| 42 | ! |
---|
| 43 | |
---|
| 44 | |
---|
| 45 | !Config Key = OCEAN |
---|
| 46 | !Config Desc = Type d'ocean |
---|
| 47 | !Config Def = force |
---|
| 48 | !Config Help = Type d'ocean utilise: force, slab,couple |
---|
| 49 | ! |
---|
| 50 | ocean = 'force ' |
---|
| 51 | call getin('OCEAN', ocean) |
---|
| 52 | ! |
---|
| 53 | !Config Key = VEGET |
---|
| 54 | !Config Desc = Type de modele de vegetation |
---|
| 55 | !Config Def = .false. |
---|
| 56 | !Config Help = Type de modele de vegetation utilise |
---|
| 57 | ! |
---|
| 58 | ok_veget = .false. |
---|
| 59 | call getin('VEGET', ok_veget) |
---|
| 60 | ! |
---|
| 61 | !Config Key = OK_journe |
---|
| 62 | !Config Desc = Pour des sorties journalieres |
---|
| 63 | !Config Def = .false. |
---|
| 64 | !Config Help = Pour creer le fichier histday contenant les sorties |
---|
| 65 | ! journalieres |
---|
| 66 | ! |
---|
| 67 | ok_journe = .false. |
---|
| 68 | call getin('OK_journe', ok_journe) |
---|
| 69 | ! |
---|
| 70 | !Config Key = OK_mensuel |
---|
| 71 | !Config Desc = Pour des sorties mensuelles |
---|
| 72 | !Config Def = .true. |
---|
| 73 | !Config Help = Pour creer le fichier histmth contenant les sorties |
---|
| 74 | ! mensuelles |
---|
| 75 | ! |
---|
| 76 | ok_mensuel = .true. |
---|
| 77 | call getin('OK_mensuel', ok_mensuel) |
---|
| 78 | ! |
---|
| 79 | !Config Key = OK_instan |
---|
| 80 | !Config Desc = Pour des sorties instantanees |
---|
| 81 | !Config Def = .false. |
---|
| 82 | !Config Help = Pour creer le fichier histins contenant les sorties |
---|
| 83 | ! instantanees |
---|
| 84 | ! |
---|
| 85 | ok_instan = .false. |
---|
| 86 | call getin('OK_instan', ok_instan) |
---|
[373] | 87 | !! |
---|
| 88 | !! KE |
---|
| 89 | ! |
---|
| 90 | !Config Key = epmax |
---|
| 91 | !Config Desc = Efficacite precip |
---|
| 92 | !Config Def = 0.993 |
---|
| 93 | !Config Help = |
---|
| 94 | ! |
---|
| 95 | epmax = .993 |
---|
| 96 | call getin('epmax', epmax) |
---|
| 97 | ! |
---|
| 98 | !Config Key = ok_adj_ema |
---|
| 99 | !Config Desc = |
---|
| 100 | !Config Def = false |
---|
| 101 | !Config Help = |
---|
| 102 | ! |
---|
| 103 | ok_adj_ema = .false. |
---|
| 104 | call getin('ok_adj_ema',ok_adj_ema) |
---|
| 105 | ! |
---|
| 106 | !Config Key = iflag_clw |
---|
| 107 | !Config Desc = |
---|
| 108 | !Config Def = 0 |
---|
| 109 | !Config Help = |
---|
| 110 | ! |
---|
| 111 | iflag_clw = 0 |
---|
| 112 | call getin('iflag_clw',iflag_clw) |
---|
| 113 | ! |
---|
| 114 | !Config Key = cld_lc_lsc |
---|
| 115 | !Config Desc = |
---|
| 116 | !Config Def = 2.6e-4 |
---|
| 117 | !Config Help = |
---|
| 118 | ! |
---|
| 119 | cld_lc_lsc = 2.6e-4 |
---|
| 120 | call getin('cld_lc_lsc',cld_lc_lsc) |
---|
| 121 | ! |
---|
| 122 | !Config Key = cld_lc_con |
---|
| 123 | !Config Desc = |
---|
| 124 | !Config Def = 2.6e-4 |
---|
| 125 | !Config Help = |
---|
| 126 | ! |
---|
| 127 | cld_lc_con = 2.6e-4 |
---|
| 128 | call getin('cld_lc_con',cld_lc_con) |
---|
| 129 | ! |
---|
| 130 | !Config Key = cld_tau_lsc |
---|
| 131 | !Config Desc = |
---|
| 132 | !Config Def = 3600. |
---|
| 133 | !Config Help = |
---|
| 134 | ! |
---|
| 135 | cld_tau_lsc = 3600. |
---|
| 136 | call getin('cld_tau_lsc',cld_tau_lsc) |
---|
| 137 | ! |
---|
| 138 | !Config Key = cld_tau_con |
---|
| 139 | !Config Desc = |
---|
| 140 | !Config Def = 3600. |
---|
| 141 | !Config Help = |
---|
| 142 | ! |
---|
| 143 | cld_tau_con = 3600. |
---|
| 144 | call getin('cld_tau_con',cld_tau_con) |
---|
| 145 | ! |
---|
| 146 | !Config Key = ffallv_lsc |
---|
| 147 | !Config Desc = |
---|
| 148 | !Config Def = 1. |
---|
| 149 | !Config Help = |
---|
| 150 | ! |
---|
| 151 | ffallv_lsc = 1. |
---|
| 152 | call getin('ffallv_lsc',ffallv_lsc) |
---|
| 153 | ! |
---|
| 154 | !Config Key = ffallv_con |
---|
| 155 | !Config Desc = |
---|
| 156 | !Config Def = 1. |
---|
| 157 | !Config Help = |
---|
| 158 | ! |
---|
| 159 | ffallv_con = 1. |
---|
| 160 | call getin('ffallv_con',ffallv_con) |
---|
| 161 | ! |
---|
| 162 | !Config Key = coef_eva |
---|
| 163 | !Config Desc = |
---|
| 164 | !Config Def = 2.e-5 |
---|
| 165 | !Config Help = |
---|
| 166 | ! |
---|
| 167 | coef_eva = 2.e-5 |
---|
| 168 | call getin('coef_eva',coef_eva) |
---|
| 169 | ! |
---|
| 170 | !Config Key = reevap_ice |
---|
| 171 | !Config Desc = |
---|
| 172 | !Config Def = .false. |
---|
| 173 | !Config Help = |
---|
| 174 | ! |
---|
| 175 | reevap_ice = .false. |
---|
| 176 | call getin('reevap_ice',reevap_ice) |
---|
| 177 | ! |
---|
| 178 | !Config Key = iflag_cldcon |
---|
| 179 | !Config Desc = |
---|
| 180 | !Config Def = 1 |
---|
| 181 | !Config Help = |
---|
| 182 | ! |
---|
| 183 | iflag_cldcon = 1 |
---|
| 184 | call getin('iflag_cldcon',iflag_cldcon) |
---|
[219] | 185 | |
---|
| 186 | ! |
---|
[373] | 187 | !Config Key = iflag_pdf |
---|
| 188 | !Config Desc = |
---|
| 189 | !Config Def = 0 |
---|
| 190 | !Config Help = |
---|
[219] | 191 | ! |
---|
[373] | 192 | iflag_pdf = 0 |
---|
| 193 | call getin('iflag_pdf',iflag_pdf) |
---|
[219] | 194 | ! |
---|
[373] | 195 | !Config Key = fact_cldcon |
---|
| 196 | !Config Desc = |
---|
| 197 | !Config Def = 0.375 |
---|
| 198 | !Config Help = |
---|
| 199 | ! |
---|
| 200 | fact_cldcon = 0.375 |
---|
| 201 | call getin('fact_cldcon',fact_cldcon) |
---|
[219] | 202 | |
---|
[373] | 203 | ! |
---|
| 204 | !Config Key = facttemps |
---|
| 205 | !Config Desc = |
---|
| 206 | !Config Def = 1.e-4 |
---|
| 207 | !Config Help = |
---|
| 208 | ! |
---|
| 209 | facttemps = 1.e-4 |
---|
| 210 | call getin('facttemps',facttemps) |
---|
| 211 | |
---|
| 212 | ! |
---|
| 213 | !Config Key = ok_newmicro |
---|
| 214 | !Config Desc = |
---|
| 215 | !Config Def = .true. |
---|
| 216 | !Config Help = |
---|
| 217 | ! |
---|
| 218 | ok_newmicro = .true. |
---|
| 219 | call getin('ok_newmicro',ok_newmicro) |
---|
| 220 | ! |
---|
| 221 | !Config Key = ratqsbas |
---|
| 222 | !Config Desc = |
---|
| 223 | !Config Def = 0.01 |
---|
| 224 | !Config Help = |
---|
| 225 | ! |
---|
| 226 | ratqsbas = 0.01 |
---|
| 227 | call getin('ratqsbas',ratqsbas) |
---|
| 228 | ! |
---|
| 229 | !Config Key = ratqshaut |
---|
| 230 | !Config Desc = |
---|
| 231 | !Config Def = 0.3 |
---|
| 232 | !Config Help = |
---|
| 233 | ! |
---|
| 234 | ratqshaut = 0.3 |
---|
| 235 | call getin('ratqshaut',ratqshaut) |
---|
| 236 | |
---|
| 237 | ! |
---|
[386] | 238 | !Config Key = rad_froid |
---|
| 239 | !Config Desc = |
---|
| 240 | !Config Def = 35.0 |
---|
| 241 | !Config Help = |
---|
[373] | 242 | ! |
---|
[386] | 243 | rad_froid = 35.0 |
---|
| 244 | call getin('rad_froid',rad_froid) |
---|
| 245 | |
---|
[373] | 246 | ! |
---|
[386] | 247 | !Config Key = rad_chau1 |
---|
| 248 | !Config Desc = |
---|
| 249 | !Config Def = 13.0 |
---|
| 250 | !Config Help = |
---|
| 251 | ! |
---|
| 252 | rad_chau1 = 13.0 |
---|
| 253 | call getin('rad_chau1',rad_chau1) |
---|
[373] | 254 | |
---|
[386] | 255 | ! |
---|
| 256 | !Config Key = rad_chau2 |
---|
| 257 | !Config Desc = |
---|
| 258 | !Config Def = 9.0 |
---|
| 259 | !Config Help = |
---|
| 260 | ! |
---|
| 261 | rad_chau2 = 9.0 |
---|
| 262 | call getin('rad_chau2',rad_chau2) |
---|
| 263 | |
---|
| 264 | ! |
---|
| 265 | !Config Key = |
---|
| 266 | !Config Desc = |
---|
| 267 | !Config Def = |
---|
| 268 | !Config Help = |
---|
| 269 | ! |
---|
| 270 | ! = |
---|
| 271 | ! call getin('',) |
---|
| 272 | ! |
---|
| 273 | ! |
---|
| 274 | ! |
---|
| 275 | ! |
---|
| 276 | |
---|
[223] | 277 | write(numout,*)' ##############################################' |
---|
| 278 | write(numout,*)' Configuration des parametres de la physique: ' |
---|
| 279 | write(numout,*)' Config ocean = ', ocean |
---|
| 280 | write(numout,*)' Config veget = ', ok_veget |
---|
| 281 | write(numout,*)' Sortie journaliere = ', ok_journe |
---|
| 282 | write(numout,*)' Sortie mensuelle = ', ok_mensuel |
---|
| 283 | write(numout,*)' Sortie instantanee = ', ok_instan |
---|
[373] | 284 | write(numout,*)' epmax = ', epmax |
---|
| 285 | write(numout,*)' ok_adj_ema = ', ok_adj_ema |
---|
| 286 | write(numout,*)' iflag_clw = ', iflag_clw |
---|
| 287 | write(numout,*)' cld_lc_lsc = ', cld_lc_lsc |
---|
| 288 | write(numout,*)' cld_lc_con = ', cld_lc_con |
---|
| 289 | write(numout,*)' cld_tau_lsc = ', cld_tau_lsc |
---|
| 290 | write(numout,*)' cld_tau_con = ', cld_tau_con |
---|
| 291 | write(numout,*)' ffallv_lsc = ', ffallv_lsc |
---|
| 292 | write(numout,*)' ffallv_con = ', ffallv_con |
---|
| 293 | write(numout,*)' coef_eva = ', coef_eva |
---|
| 294 | write(numout,*)' reevap_ice = ', reevap_ice |
---|
| 295 | write(numout,*)' iflag_pdf = ', iflag_pdf |
---|
| 296 | write(numout,*)' iflag_cldcon = ', iflag_cldcon |
---|
| 297 | write(numout,*)' fact_cldcon = ', fact_cldcon |
---|
| 298 | write(numout,*)' facttemps = ', facttemps |
---|
| 299 | write(numout,*)' ok_newmicro = ',ok_newmicro |
---|
| 300 | write(numout,*)' ratqsbas = ',ratqsbas |
---|
| 301 | write(numout,*)' ratqshaut = ',ratqshaut |
---|
[219] | 302 | |
---|
| 303 | return |
---|
| 304 | |
---|
| 305 | end subroutine conf_phys |
---|
| 306 | |
---|