Ignore:
Timestamp:
Jul 24, 2017, 5:10:39 PM (8 years ago)
Author:
aslmd
Message:

Amended previous commit. Committed more files than I wanted to.

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

    r1746 r1747  
    1414#psurf=92.e5 # Venus
    1515psurf=1.212862e6
    16 
    17 hache = 10.
    18 psurf = 610.
    19 
    20 
    2116print hache, psurf
    2217#read paramlevspe
    23 
    24 
    25 
    26 
    2718param=np.loadtxt('paramlevspe')
    2819nlev=param[0]
    … …  
    6152#print ptop
    6253etas=(pressions-ptop)/(psurf-ptop)
    63 etas[int(nlev)-1]=0.
     54etas[nlev-1]=0.
    6455#print etas
    6556press=etas*(psurf-ptop)+ptop
    … …  
    7768np.savetxt('levels',etas)
    7869
    79 plt.figure(figsize=(12, 6))
     70plt.figure(figsize=(15, 15))
    8071plt.subplot(221)
    81 plt.plot(x, etas, 'b.')
     72plt.plot(x, etas)
    8273plt.xlabel('levels')
    8374plt.ylabel('etas')
    … …  
    8576#plt.title('a) NINO3 Sea Surface Temperature (seasonal)')
    8677#plt.hold(False)
    87 axes = plt.gca()
    88 axes.set_ylim([0,1])
    8978
    9079plt.subplot(222)
    91 plt.plot(x, pseudo, 'b.')
     80plt.plot(x, pseudo)
    9281plt.xlabel('levels')
    9382plt.grid()
    9483plt.ylabel('pseudo-altitude (km)')
    95 plt.semilogy()
    96 axes = plt.gca()
    97 axes.set_ylim([0,60])
    98 
    9984
    10085plt3 = plt.subplot(223)
    101 plt.semilogy(x, press, 'b.')
     86plt.semilogy(x, press)
    10287plt.xlabel('levels')
    10388plt.grid()
    104 plt.ylabel('pressure (Pa)')
     89plt.ylabel('pression (Pa)')
    10590
    10691plt4 = plt.subplot(224)
    107 plt.plot(x, res, 'b.')
     92plt.plot(x, res)
    10893plt.xlabel('levels')
    10994plt.grid()

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

    r1746 r1747  
    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 61
    6 60.
    7 100
    8 1.75
     5201
     6280.
     70
     80.000000044
Note: See TracChangeset for help on using the changeset viewer.