| 1 | 2a3 |
|---|
| 2 | > !PLANET MARS |
|---|
| 3 | 17c18 |
|---|
| 4 | < REAL , PARAMETER :: g = 9.81 ! acceleration due to gravity (m {s}^-2) |
|---|
| 5 | --- |
|---|
| 6 | > REAL , PARAMETER :: g = 3.72 ! acceleration due to gravity (m {s}^-2) |
|---|
| 7 | 20,21c21,22 |
|---|
| 8 | < REAL , PARAMETER :: r_d = 287.04 |
|---|
| 9 | < REAL , PARAMETER :: cp = 1004.6 |
|---|
| 10 | --- |
|---|
| 11 | > REAL , PARAMETER :: r_d = 192. ! gas constant m2 s-2 K-1 |
|---|
| 12 | > REAL , PARAMETER :: cp = 844.6 ! r= 8.314511E+0 *1000.E+0/mugaz |
|---|
| 13 | 23,24c24,25 |
|---|
| 14 | < REAL , PARAMETER :: r_d = 287. |
|---|
| 15 | < REAL , PARAMETER :: cp = 7.*r_d/2. |
|---|
| 16 | --- |
|---|
| 17 | > REAL , PARAMETER :: r_d = 192. |
|---|
| 18 | > REAL , PARAMETER :: cp = 844.6 |
|---|
| 19 | 27c28 |
|---|
| 20 | < REAL , PARAMETER :: r_v = 461.6 |
|---|
| 21 | --- |
|---|
| 22 | > REAL , PARAMETER :: r_v = 461.6 ! gas constant for water vapor |
|---|
| 23 | 29c30 |
|---|
| 24 | < REAL , PARAMETER :: cpv = 4.*r_v |
|---|
| 25 | --- |
|---|
| 26 | > REAL , PARAMETER :: cpv = 4.4*r_v ! cp / R = 4.4 ?? ! rcp = r/cp ~0.256793 |
|---|
| 27 | 40,42c41,43 |
|---|
| 28 | < REAL , PARAMETER :: p1000mb = 100000. |
|---|
| 29 | < REAL , PARAMETER :: t0 = 300. |
|---|
| 30 | < REAL , PARAMETER :: p0 = p1000mb |
|---|
| 31 | --- |
|---|
| 32 | > REAL , PARAMETER :: p1000mb = 610. |
|---|
| 33 | > REAL , PARAMETER :: t0 = 220. ! earth : 300 ... à remplacer par 220 ? |
|---|
| 34 | > REAL , PARAMETER :: p0 = 610. |
|---|
| 35 | 47c48 |
|---|
| 36 | < REAL , PARAMETER :: reradius = 1./6370.0e03 |
|---|
| 37 | --- |
|---|
| 38 | > REAL , PARAMETER :: reradius = 1./3397200 ! Rayon de mars (m) ~3397200 m |
|---|
| 39 | 49,50c50,78 |
|---|
| 40 | < REAL , PARAMETER :: asselin = .025 |
|---|
| 41 | < ! REAL , PARAMETER :: asselin = .0 |
|---|
| 42 | --- |
|---|
| 43 | > !------------------------------- |
|---|
| 44 | > !****MARS |
|---|
| 45 | > ! here some LMD stuff |
|---|
| 46 | > ! >> TODO: possible modification with startfi.nc values ? |
|---|
| 47 | > REAL , PARAMETER :: wdaysec = 88800. ! duree du sol (s) ~88775 s |
|---|
| 48 | > REAL , PARAMETER :: wmugaz = 43.49 ! Masse molaire de l'atm (g.mol-1) ~43.49 |
|---|
| 49 | > REAL , PARAMETER :: womeg = 7.0721E-5 ! omega (rad.s-1) |
|---|
| 50 | > REAL , PARAMETER :: wyear_day = 669 ! Duree de l'annee (sols) ~668.6 |
|---|
| 51 | > REAL , PARAMETER :: wperiheli = 206.66 ! Dist.min. soleil-mars (Mkm) ~206.66 |
|---|
| 52 | > REAL , PARAMETER :: waphelie = 249.22 ! Dist.max. soleil-mars (Mkm) ~249.22 |
|---|
| 53 | > REAL , PARAMETER :: wperi_day = 485. ! Date du perihelie (sols depuis printemps) |
|---|
| 54 | > REAL , PARAMETER :: wobliquit = 25.2 ! Obliquite de la planete (deg) ~25.2 |
|---|
| 55 | > REAL , PARAMETER :: wz0 = 1.e-2 ! surface roughness (m) ~0.01 |
|---|
| 56 | > REAL , PARAMETER :: wlmixmin = 30 ! longueur de melange ~100 |
|---|
| 57 | > REAL , PARAMETER :: wemin_turb = 1.e-6 ! energie minimale ~1.e-8 |
|---|
| 58 | > REAL , PARAMETER :: wemissiv = 0.95 ! Emissivite du sol martien ~.95 |
|---|
| 59 | > REAL , PARAMETER :: wemissiceN = 0.95 ! Emissivite calotte nord |
|---|
| 60 | > REAL , PARAMETER :: wemissiceS = 0.95 ! Emissivite calotte sud |
|---|
| 61 | > REAL , PARAMETER :: walbediceN = 0.65 ! Albedo calotte nord !0.5 |
|---|
| 62 | > REAL , PARAMETER :: walbediceS = 0.65 ! Albedo calotte sud !0.5 |
|---|
| 63 | > REAL , PARAMETER :: wiceradiusN = 100.e-6 ! mean scat radius of CO2 snow (north) |
|---|
| 64 | > REAL , PARAMETER :: wiceradiusS = 100.e-6 ! mean scat radius of CO2 snow (south) |
|---|
| 65 | > REAL , PARAMETER :: wdtemisiceN = 0.4 ! time scale for snow metamorphism (north) !2 |
|---|
| 66 | > REAL , PARAMETER :: wdtemisiceS = 0.4 ! time scale for snow metamorphism (south) !2 |
|---|
| 67 | > ! |
|---|
| 68 | > !------------------------------- |
|---|
| 69 | > |
|---|
| 70 | > ! REAL , PARAMETER :: asselin = .025 |
|---|
| 71 | > REAL , PARAMETER :: asselin = .0 |
|---|
| 72 | 64,65c92,93 |
|---|
| 73 | < REAL , PARAMETER :: rhoair0 = 1.28 |
|---|
| 74 | < |
|---|
| 75 | --- |
|---|
| 76 | > REAL , PARAMETER :: rhoair0 = 0.02 !Earth Surface density: 1.217 kg/m3 |
|---|
| 77 | > !Mars Surface density: 0.020 kg/m3 |
|---|
| 78 | 67c95 |
|---|
| 79 | < REAL , PARAMETER :: DPD = 360./365. |
|---|
| 80 | --- |
|---|
| 81 | > REAL , PARAMETER :: DPD = 360./669. ! longitude solaire d'un jour? |
|---|
| 82 | 75,77c103,105 |
|---|
| 83 | < REAL , PARAMETER :: KARMAN=0.4 |
|---|
| 84 | < REAL , PARAMETER :: EOMEG=7.2921E-5 |
|---|
| 85 | < REAL , PARAMETER :: STBOLT=5.67051E-8 |
|---|
| 86 | --- |
|---|
| 87 | > REAL , PARAMETER :: KARMAN=0.4 ! von karman constant |
|---|
| 88 | > REAL , PARAMETER :: EOMEG=7.0721E-5 ! angular rotation rate (rad.s-1) |
|---|
| 89 | > REAL , PARAMETER :: STBOLT=5.67051E-8 ! stefan-boltzmann constant |
|---|
| 90 | 89c117,118 |
|---|
| 91 | < ! REAL , PARAMETER :: epsq2=0.02 |
|---|
| 92 | --- |
|---|
| 93 | > !try both ? |
|---|
| 94 | > ! REAL , PARAMETER :: epsq2=0.02 |
|---|
| 95 | 106a136 |
|---|
| 96 | > ! REAL , PARAMETER :: fcm=0.003 |
|---|
| 97 | 107a138 |
|---|
| 98 | > ! REAL , PARAMETER :: fcm=0.0 |
|---|
| 99 | 127a159 |
|---|
| 100 | > ! REAL , PARAMETER :: z0max=0.01 |
|---|