|
Last change
on this file since 4191 was
4176,
checked in by dubos, 6 years ago
|
|
simple_physics : copy code from FH
|
|
File size:
419 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE lwtr(ngrid,coef,lstrong,dup,transm) |
|---|
| 2 | IMPLICIT NONE |
|---|
| 3 | INTEGER ngrid |
|---|
| 4 | REAL coef |
|---|
| 5 | LOGICAL lstrong |
|---|
| 6 | REAL dup(ngrid),transm(ngrid) |
|---|
| 7 | INTEGER ig |
|---|
| 8 | IF(lstrong) THEN |
|---|
| 9 | DO ig=1,ngrid |
|---|
| 10 | transm(ig)=exp(-coef*sqrt(dup(ig))) |
|---|
| 11 | ENDDO |
|---|
| 12 | ELSE |
|---|
| 13 | DO ig=1,ngrid |
|---|
| 14 | transm(ig)=exp(-coef*dup(ig)) |
|---|
| 15 | ENDDO |
|---|
| 16 | ENDIF |
|---|
| 17 | RETURN |
|---|
| 18 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
