Last change
on this file since 4728 was
4728,
checked in by idelkadi, 10 months ago
|
Update of ecrad in the LMDZ_ECRad branch of LMDZ:
- version 1.6.1 of ecrad
- files are no longer grouped in the same ecrad directory.
- the structure of ecrad offline is preserved to facilitate updating in LMDZ
- cfg.bld modified to take into account the new added subdirectories.
- the interface routines and those added in ecrad are moved to the phylmd directory
|
File size:
652 bytes
|
Line | |
---|
1 | program test_fast_expm |
---|
2 | ! Test the fast matrix exponential exchange routine used for |
---|
3 | ! SPARTACUS shortwave entrapment calculations, but which can fail in |
---|
4 | ! single precision for very specific inputs leading to the matrix |
---|
5 | ! having two repeated eigenvalues |
---|
6 | |
---|
7 | use parkind1, only : jprb |
---|
8 | |
---|
9 | use radiation_matrix, only : fast_expm_exchange_3 |
---|
10 | |
---|
11 | real(jprb) :: a(1), b(1), c(1), d(1), r(1,3,3) |
---|
12 | |
---|
13 | a = 9.0408579E-02_jprb |
---|
14 | b = 9.2716664E-07_jprb |
---|
15 | c = 2.2503915E-03_jprb |
---|
16 | d = 8.8152386E-02_jprb |
---|
17 | |
---|
18 | call fast_expm_exchange_3(1,1,a(1:1),b(1:1),c(1:1),d(1:1),r) |
---|
19 | |
---|
20 | print *, r |
---|
21 | print *, 0.0000001 |
---|
22 | print *, epsilon(1.0_jprb) |
---|
23 | |
---|
24 | end program test_fast_expm |
---|
Note: See
TracBrowser
for help on using the repository browser.