Changeset 3570 for palm/trunk/SOURCE/chemistry_model_mod.f90
- Timestamp:
- Nov 27, 2018 5:44:21 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r3543 r3570 21 21 ! 22 22 ! Current revisions: 23 ! ----------------- 23 ! ------------------ 24 24 ! 25 25 ! … … 27 27 ! ----------------- 28 28 ! $Id$ 29 ! resler: 30 ! Break lines at 132 characters 31 ! 32 ! 3543 2018-11-20 17:06:15Z suehring 29 33 ! Remove tabs 30 34 ! … … 323 327 ! 324 328 ! Set alpha for f_light (4.57 is conversion factor from 1./(mumol m-2 s-1) naar W m-2 325 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: alpha =(/0.009,0.009, 0.009,0.006,0.006, -999., -999.,0.009,-999.,-999.,0.009,0.006,-999.,0.009,0.008/)*4.57 329 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: alpha = & 330 (/0.009,0.009, 0.009,0.006,0.006, -999., -999.,0.009,-999.,-999.,0.009,0.006,-999.,0.009,0.008/)*4.57 326 331 ! 327 332 ! Set temperatures per land use for F_temp 328 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmin = (/12.0, 12.0, 12.0, 0.0, 0.0, -999., -999., 12.0, -999., -999., 12.0, 0.0, -999., 12.0, 8.0/) 329 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Topt = (/26.0, 26.0, 26.0, 18.0, 20.0, -999., -999., 26.0, -999., -999., 26.0, 20.0, -999., 26.0, 24.0 /) 330 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmax = (/40.0, 40.0, 40.0, 36.0, 35.0, -999., -999., 40.0, -999., -999., 40.0, 35.0, -999., 40.0, 39.0 /) 333 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmin = & 334 (/12.0, 12.0, 12.0, 0.0, 0.0, -999., -999., 12.0, -999., -999., 12.0, 0.0, -999., 12.0, 8.0/) 335 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Topt = & 336 (/26.0, 26.0, 26.0, 18.0, 20.0, -999., -999., 26.0, -999., -999., 26.0, 20.0, -999., 26.0, 24.0 /) 337 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmax = & 338 (/40.0, 40.0, 40.0, 36.0, 35.0, -999., -999., 40.0, -999., -999., 40.0, 35.0, -999., 40.0, 39.0 /) 331 339 ! 332 340 ! Set F_min: 333 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: F_min =(/0.01, 0.01, 0.01, 0.1, 0.1, -999., -999.,0.01, -999.,-999.,0.01,0.1,-999.,0.01, 0.04/) 341 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: F_min = & 342 (/0.01, 0.01, 0.01, 0.1, 0.1, -999., -999.,0.01, -999.,-999.,0.01,0.1,-999.,0.01, 0.04/) 334 343 335 344 ! Set maximal conductance (m/s) 336 345 ! (R T/P) = 1/41000 mmol/m3 is given for 20 deg C to go from mmol O3/m2/s to m/s 337 346 ! Could be refined to a function of T and P. in Jones 338 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: g_max =(/270., 300., 300., 140., 150., -999., -999.,270., -999.,-999.,270., 150.,-999.,300., 422./)/41000 347 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: g_max = & 348 (/270., 300., 300., 140., 150., -999., -999.,270., -999.,-999.,270., 150.,-999.,300., 422./)/41000 339 349 ! 340 350 ! Set max, min for vapour pressure deficit vpd; 341 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_max =(/1.3, 0.9, 0.9, 0.5, 1.0, -999., -999.,1.3, -999.,-999.,1.3,1.0, -999.,0.9, 2.8/) 342 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_min =(/3.0, 2.8, 2.8, 3.0, 3.25, -999., -999.,3.0, -999.,-999.,3.0,3.25, -999.,2.8, 4.5/) 351 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_max = & 352 (/1.3, 0.9, 0.9, 0.5, 1.0, -999., -999.,1.3, -999.,-999.,1.3,1.0, -999.,0.9, 2.8/) 353 REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_min = & 354 (/3.0, 2.8, 2.8, 3.0, 3.25, -999., -999.,3.0, -999.,-999.,3.0,3.25, -999.,2.8, 4.5/) 343 355 ! 344 356 ! … … 4958 4970 4959 4971 4960 REAL(kind=wp), PARAMETER :: rsoil_wet(ncmp) = (/2000., 10., 2000., -999., 10., -999., -999., -999., -999., -999./) 4961 REAL(kind=wp), PARAMETER :: rsoil_frozen(ncmp) = (/2000., 500., 2000., -999., 1000., -999., -999., -999., -999., -999./) 4972 REAL(kind=wp), PARAMETER :: rsoil_wet(ncmp) = & 4973 (/2000., 10., 2000., -999., 10., -999., -999., -999., -999., -999./) 4974 REAL(kind=wp), PARAMETER :: rsoil_frozen(ncmp) = & 4975 (/2000., 500., 2000., -999., 1000., -999., -999., -999., -999., -999./) 4962 4976 4963 4977 … … 5027 5041 ! b = empirical constant for computation of rinc (in canopy resistance) (= 14 m-1 or -999 IF not applicable) 5028 5042 ! h = vegetation height (m) grass arabl crops conIF decid water urba othr desr ice sav trf wai med semi 5029 REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: b = (/ -999, 14, 14, 14, 14, -999, -999, -999, -999, -999, -999, 14, -999, 14, 14 /) 5030 REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: h = (/ -999, 1, 1, 20, 20, -999, -999, -999, -999, -999, -999, 20, -999, 1 , 1 /) 5043 REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: b = & 5044 (/ -999, 14, 14, 14, 14, -999, -999, -999, -999, -999, -999, 14, -999, 14, 14 /) 5045 REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: h = & 5046 (/ -999, 1, 1, 20, 20, -999, -999, -999, -999, -999, -999, 20, -999, 1 , 1 /) 5031 5047 5032 5048 ! Compute Rinc only for arable land, perm. crops, forest; otherwise Rinc = 0:
Note: See TracChangeset
for help on using the changeset viewer.