| 1 | #! /usr/bin/env python |
|---|
| 2 | |
|---|
| 3 | from fmcd import call_mcd |
|---|
| 4 | import numpy as np |
|---|
| 5 | |
|---|
| 6 | lon = float(raw_input('Longitude?')) |
|---|
| 7 | lat = float(raw_input('Latitude?')) |
|---|
| 8 | xdate = float(raw_input('Ls?')) |
|---|
| 9 | loct = float(raw_input('Local time?')) |
|---|
| 10 | xz = float(raw_input('Altitude?')) |
|---|
| 11 | |
|---|
| 12 | dset = '/home/aymeric/Science/MCD_v4.3/data/' |
|---|
| 13 | |
|---|
| 14 | zkey = 3 # specify that xz is the altitude above surface (m) |
|---|
| 15 | datekey = 1 # 1 = "Mars date": xdate is the value of Ls |
|---|
| 16 | dust = 2 #our best guess MY24 scenario, with solar average conditions |
|---|
| 17 | hrkey = 1 #set high resolution mode on (hrkey=0 to set high resolution off) |
|---|
| 18 | perturkey = 0 #integer perturkey ! perturbation type (0: none) |
|---|
| 19 | seedin = 1 #random number generator seed (unused if perturkey=0) |
|---|
| 20 | gwlength = 0. #gravity Wave wavelength (unused if perturkey=0) |
|---|
| 21 | extvarkey = 1 |
|---|
| 22 | #extvarkeys = np.ones(100) |
|---|
| 23 | |
|---|
| 24 | (pres, dens, temp, zonwind, merwind, \ |
|---|
| 25 | meanvar, extvar, seedout, ierr) \ |
|---|
| 26 | = \ |
|---|
| 27 | call_mcd(zkey,xz,lon,lat,hrkey, \ |
|---|
| 28 | datekey,xdate,loct,dset,dust, \ |
|---|
| 29 | perturkey,seedin,gwlength,extvarkey ) |
|---|
| 30 | |
|---|
| 31 | print "temperature is %.0f K, pressure is %.0f Pa, density is %5.3e kg/m3, zonal wind is %.1f m/s, meridional wind is %.1f m/s" % (temp,pres,dens,zonwind,merwind) |
|---|
| 32 | |
|---|
