Changeset 1933 for LMDZ5/trunk/libf/phylmd
- Timestamp:
- Jan 17, 2014, 7:06:50 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/radiation_ar4_param.F90
r1907 r1933 1 1 MODULE radiation_ar4_param 2 2 3 REAL *8, parameter :: ZPDH2O = 0.84 REAL *8, parameter :: ZPDUMG = 0.755 REAL *8, parameter :: ZPRH2O = 30000.06 REAL *8, parameter :: ZPRUMG = 30000.07 REAL *8, parameter :: RTDH2O = 0.408 REAL *8, parameter :: RTDUMG = 0.3759 REAL *8, parameter :: RTH2O = 240.010 REAL *8, parameter :: RTUMG = 240.011 12 REAL *8, dimension(2), parameter :: WG1 = (/1.0, 1.0/)13 REAL *8, dimension(11), parameter :: TINTP = (/ 187.5, 200., 212.5, 225., 237.5, 250., 262.5, 275., 287.5, 300., 312.5 /)14 15 real *8, dimension(11,16,3), parameter :: GA = reshape ( (/&3 REAL(KIND=8), parameter :: ZPDH2O = 0.8 4 REAL(KIND=8), parameter :: ZPDUMG = 0.75 5 REAL(KIND=8), parameter :: ZPRH2O = 30000.0 6 REAL(KIND=8), parameter :: ZPRUMG = 30000.0 7 REAL(KIND=8), parameter :: RTDH2O = 0.40 8 REAL(KIND=8), parameter :: RTDUMG = 0.375 9 REAL(KIND=8), parameter :: RTH2O = 240.0 10 REAL(KIND=8), parameter :: RTUMG = 240.0 11 12 REAL(KIND=8), dimension(2), parameter :: WG1 = (/1.0, 1.0/) 13 REAL(KIND=8), dimension(11), parameter :: TINTP = (/ 187.5, 200., 212.5, 225., 237.5, 250., 262.5, 275., 287.5, 300., 312.5 /) 14 15 real(kind=8), dimension(11,16,3), parameter :: GA = reshape ( (/ & 16 16 0.63499072E-02, 0.65566348E-02, 0.67849730E-02, 0.70481947E-02, 0.73585943E-02, 0.77242818E-02, 0.81472693E-02, 0.86227527E-02,& 17 17 0.91396814E-02, 0.96825438E-02, 0.10233955E-01, 0.77266491E-02, 0.81323287E-02, 0.86507620E-02, 0.92776391E-02, 0.99806312E-02,& … … 82 82 , (/ 11,16,3 /) ) 83 83 84 real *8, dimension(11,16,3), parameter :: GB = reshape ( (/&84 real(kind=8), dimension(11,16,3), parameter :: GB = reshape ( (/ & 85 85 0.63499072E-02, 0.65566348E-02, 0.67849730E-02, 0.70481947E-02, 0.73585943E-02, 0.77242818E-02, 0.81472693E-02, 0.86227527E-02,& 86 86 0.91396814E-02, 0.96825438E-02, 0.10233955E-01, 0.77266491E-02, 0.81323287E-02, 0.86507620E-02, 0.92776391E-02, 0.99806312E-02,& … … 151 151 , (/ 11,16,3 /) ) 152 152 153 real *8, dimension(6,6), parameter :: XP = reshape ( (/ &153 real(kind=8), dimension(6,6), parameter :: XP = reshape ( (/ & 154 154 0.46430621E+02, 0.12928299E+03, 0.20732648E+03, 0.31398411E+03, 0.18373177E+03, -0.11412303E+03, & 155 155 0.73604774E+02, 0.27887914E+03, 0.27076947E+03, -0.57322111E+02, -0.64742459E+02, 0.87238280E+02, & … … 160 160 161 161 ! If running with Reporbus type_trac=repr, values of RSUN might be overritten in radiation_AR4 162 REAL *8, dimension(2) :: RSUN = (/ 0.441676 , 0.558324 /)163 REAL *8, dimension(2,6), parameter :: RRAY = reshape ( &162 REAL(KIND=8), dimension(2) :: RSUN = (/ 0.441676 , 0.558324 /) 163 REAL(KIND=8), dimension(2,6), parameter :: RRAY = reshape ( & 164 164 (/ .428937E-01, .697200E-02,& 165 165 .890743E+00, .173297E-01,& … … 169 169 .161645E+01, .129662E+00 /) , (/2,6/) ) 170 170 171 REAL *8, dimension(2,5), parameter :: TAUA = reshape ( &171 REAL(KIND=8), dimension(2,5), parameter :: TAUA = reshape ( & 172 172 (/ 0.730719, 0.730719, 0.912819, 0.912819, 0.725059, & 173 173 0.725059, 0.745405, 0.745405, 0.682188, 0.682188 /),(/2,5/) ) 174 REAL *8, dimension(2,5), parameter :: RPIZA = reshape ( &174 REAL(KIND=8), dimension(2,5), parameter :: RPIZA = reshape ( & 175 175 (/ 0.872212, 0.872212, 0.982545, 0.982545, 0.623143, & 176 176 0.623143, 0.944887, 0.944887, 0.997975, 0.997975 /),(/2,5/) ) 177 REAL *8, dimension(2,5), parameter :: RCGA = reshape ( &177 REAL(KIND=8), dimension(2,5), parameter :: RCGA = reshape ( & 178 178 (/ 0.647596, 0.647596, 0.739002, 0.739002, 0.580845, & 179 179 0.580845, 0.662657, 0.662657, 0.624246, 0.624246 /),(/2,5/) ) 180 180 181 REAL *8, dimension(2,3,7), parameter :: APAD = reshape ( &181 REAL(KIND=8), dimension(2,3,7), parameter :: APAD = reshape ( & 182 182 (/ 0.912418292E+05, 0.376655383E-08, 0.000000000E-00,& 183 183 0.739646016E-08, 0.925887084E-04, 0.410177786E+03,& … … 194 194 0.000000000E-00, 0.000000000E+00, 0.000000000E-00,& 195 195 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /) , (/2,3,7/) ) 196 REAL *8, dimension(2,3,7), parameter :: BPAD = reshape ( &196 REAL(KIND=8), dimension(2,3,7), parameter :: BPAD = reshape ( & 197 197 (/ 0.912418292E+05, 0.376655383E-08, 0.000000000E-00,& 198 198 0.739646016E-08, 0.925887084E-04, 0.410177786E+03,& … … 209 209 0.000000000E-00, 0.100000000E+01, 0.000000000E-00,& 210 210 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /) , (/2,3,7/) ) 211 REAL *8, dimension(2,3), parameter :: D = reshape ( &211 REAL(KIND=8), dimension(2,3), parameter :: D = reshape ( & 212 212 (/ 0.0, 0.0, 0.0, 0.0, 0.0, 0.8 /) , (/2,3/) ) 213 213 214 REAL *8, parameter :: TREF = 250.0215 REAL *8, dimension(2), parameter :: RT1 = (/ -0.577350269, +0.577350269 /)216 REAL *8, dimension(5,5), parameter :: RAER= reshape ( &214 REAL(KIND=8), parameter :: TREF = 250.0 215 REAL(KIND=8), dimension(2), parameter :: RT1 = (/ -0.577350269, +0.577350269 /) 216 REAL(KIND=8), dimension(5,5), parameter :: RAER= reshape ( & 217 217 (/ .038520, .037196, .040532, .054934, .038520 & 218 218 , .12613 , .18313 , .10357 , .064106, .126130 & … … 221 221 , .013792, .026810, .052203, .066338, .013792 /) , (/5,5/) ) 222 222 223 REAL *8, dimension(8,3), parameter :: AT= reshape ( &223 REAL(KIND=8), dimension(8,3), parameter :: AT= reshape ( & 224 224 (/ 0.298199E-02,0.143676E-01,0.197861E-01,0.289560E-01,& 225 225 0.103800E-01,0.868859E-02,0.250073E-03,0.307423E-01,& … … 229 229 -.161431E-02,0.000000E-00,0.109242E-03,-.322172E-03 /) , (/8,3/) ) 230 230 231 REAL *8, dimension(8,3), parameter :: BT= reshape ( &231 REAL(KIND=8), dimension(8,3), parameter :: BT= reshape ( & 232 232 (/ -0.106432E-04,-0.553979E-04,-0.877012E-04,-0.165960E-03,& 233 233 -0.276744E-04,-0.278412E-04, 0.199846E-05,-0.108482E-03,& … … 237 237 0.127646E-04 ,0.117469E-05, 0.175991E-06,-0.814575E-06 /) , (/8,3/) ) 238 238 239 REAL *8, dimension(4), parameter :: OCT = (/ -.326E-03, -.102E-05, .137E-02, -.535E-05 /)239 REAL(KIND=8), dimension(4), parameter :: OCT = (/ -.326E-03, -.102E-05, .137E-02, -.535E-05 /) 240 240 241 241 end module radiation_AR4_param
Note: See TracChangeset
for help on using the changeset viewer.