Ignore:
Timestamp:
Feb 1, 2019, 3:18:29 PM (6 years ago)
Author:
mlefevre
Message:

MESOSCALE etas level generation. Surface pressure and scale height as argument

Location:
trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level_python
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level_python/levspe.py

    r1747 r2087  
    88import matplotlib.pyplot as plt
    99
     10param=np.loadtxt('paramlevspe')
     11hache=param[0]
     12psurf=param[1]
     13nlev=param[2]
     14altmax=param[3]
     15epsilon=param[4]
     16elong_cos=param[5]
    1017
    11 #hache=10.   # Venus scale height (km)
    12 hache=(3461.277*1305.70)/79.43000/1e3
    13 #psurf=610.  # Mars pression surface -- ne change pas les eta levels
    14 #psurf=92.e5 # Venus
    15 psurf=1.212862e6
    16 print hache, psurf
    17 #read paramlevspe
    18 param=np.loadtxt('paramlevspe')
    19 nlev=param[0]
    20 altmax=param[1]
    21 epsilon=param[2]
    22 elong_cos=param[3]
    23 
    24 print nlev, altmax, epsilon, elong_cos
     18print nlev, altmax, hache, psurf, epsilon, elong_cos
    2519
    2620epsilon = epsilon / (nlev-1) / 100.

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level_python/paramlevspe

    r1747 r2087  
     1#hache : scale height (km)
     2#psurf : surface pressure (Pa)
    13#nlev : nombre de niveaux.
    24#altmax : altitude maximale en km
    35#epsilon augmentation totale en pourcentage de l'ecart max (modifie inflexion). epsilon = 0   : point d'inflexion parfaitement plat. epsilon = 100 : pas de point d'inflexion
    46#elong_cos plus petit rapproche l'inflexion du sol ;; 2/3 parfait pour un nlev divisible par 3
     710.
     81.e5
    59201
    610280.
Note: See TracChangeset for help on using the changeset viewer.