|
Last change
on this file since 3026 was
2963,
checked in by romain.vande, 20 months ago
|
|
Mars PEM :
Adapt PEM to the subslope PCM configuration, it is now fully compatible.
Create a PEM folder in deftank that contains:
run_pem1: a bash file that runs chained simulation of PEM as well as running a parameterizable number of PCM simulation in between.
It also takes care of reshaping XIOS output as well as renaming outputs etc… in the spirit of run_month1.
It is written for Irene machine and the header needs to be adapted for other machines.
run_PEM.def: A text file that shows the possible parameters to choose before a PEM simulation.
It should be included at the end of run.def just like callphys.def
ob_ex_lsp.asc: An ascii file containing the obliquity, eccentricity, ls_peri data from Laskar in Martian year.
README: A txt file explaining the content of the folder
Adapt field_def_physics_mars.xml to consider the case with 7 subslopes in the PCM.
Change context_lmdz_physics.xml to be able to output the file needed by the PEM.
Correct a few other minor bugs.
RV
|
|
File size:
1.1 KB
|
| Rev | Line | |
|---|
| [2448] | 1 | c ******************************************************* |
|---|
| [2228] | 2 | subroutine call_dayperi(Lsperi,e_elips,dayperi,year_day) |
|---|
| 3 | IMPLICIT NONE |
|---|
| 4 | |
|---|
| 5 | c======================================================================= |
|---|
| 6 | c |
|---|
| 7 | c Objet: |
|---|
| 8 | c ------ |
|---|
| 9 | c Computing the Martian date (number of sols since Ls=0 = spring |
|---|
| 10 | c equinox) at perihelion |
|---|
| 11 | c |
|---|
| 12 | c Arguments: |
|---|
| 13 | c ---------- |
|---|
| 14 | c |
|---|
| 15 | c Input: |
|---|
| 16 | c ------ |
|---|
| 17 | c Lsperi solar longitude (Ls) of perohelion (rad) |
|---|
| 18 | c e_elips Excentricity |
|---|
| [2963] | 19 | c real year_day ! number of sols per Mars yar |
|---|
| [2228] | 20 | c |
|---|
| 21 | c output |
|---|
| 22 | c ------ |
|---|
| 23 | c dayperi Martian date at perihelion (sol) |
|---|
| [2963] | 24 | |
|---|
| [2228] | 25 | c----------------------------------------------------------------------- |
|---|
| 26 | |
|---|
| 27 | c arguments: |
|---|
| 28 | c ---------- |
|---|
| 29 | real Lsperi,e_elips,dayperi |
|---|
| 30 | real year_day ! number of sols per Mars yar |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | c Local |
|---|
| 34 | real x1, x2,pi |
|---|
| 35 | |
|---|
| 36 | |
|---|
| 37 | pi=2.*asin(1.) |
|---|
| 38 | |
|---|
| 39 | x1 = sqrt( (1-e_elips)/(1+e_elips) ) |
|---|
| 40 | x2 = e_elips*sqrt(1-e_elips**2) |
|---|
| 41 | |
|---|
| 42 | dayperi = 0.5*(year_day/pi)* |
|---|
| 43 | & ( 2*atan(x1*tan(0.5*Lsperi)) |
|---|
| 44 | & -x2*sin(Lsperi)/(1+e_elips*cos(Lsperi)) ) |
|---|
| 45 | if(dayperi.lt.0) dayperi=dayperi+year_day |
|---|
| 46 | return |
|---|
| 47 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
