|
Last change
on this file since 2831 was
2779,
checked in by romain.vande, 3 years ago
|
|
LMDZ_COMMON:
First commit of PEM. Small writting improvements will follow.
Only works with GCM slopes parametrisation of Lucas L.
Example to compile : ./makelmdz_fcm -arch X64_OCCIGEN -arch_path ../ARCH -p mars -d 32x24x26 pem
RV
|
|
File size:
1.1 KB
|
| Line | |
|---|
| 1 | ! |
|---|
| 2 | ! $Id: conf_evol.F $ |
|---|
| 3 | |
|---|
| 4 | SUBROUTINE conf_evol( tapedef, etatinit ) |
|---|
| 5 | |
|---|
| 6 | #ifdef CPP_IOIPSL |
|---|
| 7 | use IOIPSL |
|---|
| 8 | #else |
|---|
| 9 | ! if not using IOIPSL, we still need to use (a local version of) getin |
|---|
| 10 | use ioipsl_getincom |
|---|
| 11 | #endif |
|---|
| 12 | |
|---|
| 13 | USE temps_mod_evol, ONLY: nyear, dt_pem, alpha_criterion |
|---|
| 14 | |
|---|
| 15 | IMPLICIT NONE |
|---|
| 16 | !----------------------------------------------------------------------- |
|---|
| 17 | ! Read the run_pem.def file |
|---|
| 18 | ! |
|---|
| 19 | ! Arguments : |
|---|
| 20 | ! |
|---|
| 21 | ! nyear : Maximun number of year over which the PEM can interpolate |
|---|
| 22 | ! dt_pem : in years, the time step used by the PEM |
|---|
| 23 | ! alpha_criterion : percentage of change before stopping the PEM |
|---|
| 24 | ! |
|---|
| 25 | LOGICAL,INTENT(IN) :: etatinit |
|---|
| 26 | INTEGER,INTENT(IN) :: tapedef |
|---|
| 27 | |
|---|
| 28 | ! Declarations : |
|---|
| 29 | ! -------------- |
|---|
| 30 | include "dimensions.h" |
|---|
| 31 | include "paramet.h" |
|---|
| 32 | include "comdissnew.h" |
|---|
| 33 | include "iniprint.h" |
|---|
| 34 | |
|---|
| 35 | ! local: |
|---|
| 36 | ! ------ |
|---|
| 37 | |
|---|
| 38 | REAL clonn,clatt,grossismxx,grossismyy |
|---|
| 39 | REAL dzoomxx,dzoomyy, tauxx,tauyy |
|---|
| 40 | LOGICAL fxyhypbb, ysinuss |
|---|
| 41 | LOGICAL use_filtre_fft |
|---|
| 42 | |
|---|
| 43 | |
|---|
| 44 | nyear=100 |
|---|
| 45 | CALL getin('nyear', nyear) |
|---|
| 46 | |
|---|
| 47 | dt_pem=1 |
|---|
| 48 | CALL getin('dt_pem', dt_pem) |
|---|
| 49 | |
|---|
| 50 | alpha_criterion=0.2 |
|---|
| 51 | CALL getin('alpha_criterion', alpha_criterion) |
|---|
| 52 | |
|---|
| 53 | |
|---|
| 54 | END SUBROUTINE conf_evol |
|---|
Note: See
TracBrowser
for help on using the repository browser.
