Ignore:
Timestamp:
Jul 24, 2017, 5:06:29 PM (7 years ago)
Author:
aslmd
Message:

MESOSCALE. makemeso module_initialize_les module_initialize_real reverted back to r1723 because commit r1724 was based on corrupted versions.

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

    r1606 r1746  
    1414#psurf=92.e5 # Venus
    1515psurf=1.212862e6
     16
     17hache = 10.
     18psurf = 610.
     19
     20
    1621print hache, psurf
    1722#read paramlevspe
     23
     24
     25
     26
    1827param=np.loadtxt('paramlevspe')
    1928nlev=param[0]
    … …  
    5261#print ptop
    5362etas=(pressions-ptop)/(psurf-ptop)
    54 etas[nlev-1]=0.
     63etas[int(nlev)-1]=0.
    5564#print etas
    5665press=etas*(psurf-ptop)+ptop
    … …  
    6877np.savetxt('levels',etas)
    6978
    70 plt.figure(figsize=(15, 15))
     79plt.figure(figsize=(12, 6))
    7180plt.subplot(221)
    72 plt.plot(x, etas)
     81plt.plot(x, etas, 'b.')
    7382plt.xlabel('levels')
    7483plt.ylabel('etas')
    … …  
    7685#plt.title('a) NINO3 Sea Surface Temperature (seasonal)')
    7786#plt.hold(False)
     87axes = plt.gca()
     88axes.set_ylim([0,1])
    7889
    7990plt.subplot(222)
    80 plt.plot(x, pseudo)
     91plt.plot(x, pseudo, 'b.')
    8192plt.xlabel('levels')
    8293plt.grid()
    8394plt.ylabel('pseudo-altitude (km)')
     95plt.semilogy()
     96axes = plt.gca()
     97axes.set_ylim([0,60])
     98
    8499
    85100plt3 = plt.subplot(223)
    86 plt.semilogy(x, press)
     101plt.semilogy(x, press, 'b.')
    87102plt.xlabel('levels')
    88103plt.grid()
    89 plt.ylabel('pression (Pa)')
     104plt.ylabel('pressure (Pa)')
    90105
    91106plt4 = plt.subplot(224)
    92 plt.plot(x, res)
     107plt.plot(x, res, 'b.')
    93108plt.xlabel('levels')
    94109plt.grid()

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

    r1606 r1746  
    33#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
    44#elong_cos plus petit rapproche l'inflexion du sol ;; 2/3 parfait pour un nlev divisible par 3
    5 201
    6 280.
    7 0
    8 0.000000044
     561
     660.
     7100
     81.75
Note: See TracChangeset for help on using the changeset viewer.