Changeset 194 for trunk/MESOSCALE


Ignore:
Timestamp:
Jul 4, 2011, 6:13:23 PM (14 years ago)
Author:
aslmd
Message:

MESOSCALE: corrected graphics script to display Ls. and added a figure in the first page of the user manual.

Location:
trunk/MESOSCALE
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/DOC/SRC/these_spiga.sty

    r168 r194  
    44%\documentclass[12pt,a4paper,twoside]{report}
    55%
    6 \documentclass[10pt,a4paper,twoside,openright]{book}
     6\documentclass[10pt,a4paper,twoside,openright,titlepage]{book}
    77%\documentclass[10pt,a4paper,twoside]{report}
    88%\documentclass[10pt,a4paper,twoside]{book}

  • trunk/MESOSCALE/DOC/SRC/user_manual.tex

    r168 r194  
    1212%
    1313
    14 %
    15 \author{
    16 \huge{A. Spiga}\\
    17 \\
    18 Laboratoire de M\'et\'eorologie Dynamique\\
    19 Institut Pierre Simon Laplace\\
    20 Universit\'e Pierre et Marie Curie\\
    21 Paris, France\\
    22 \\
    23 Contact: \email{aymeric.spiga}{upmc.fr}\\
    24 }
    25 \title{
    26 %\huge{ESA/CNES contracts ``Mars Environment Models"\\
    27 %WP11.1 Deliverable Report}\\
    28 %\mbox{}\\
    29 %\Large{Ref:~ESA 11369/95/NL/JG(SC) \\
    30 %CNES ``Base de donn\'ees atmosph\'eriques martiennes"}\\
    31 \mbox{}\\
    32 \mbox{}\\
    33 \Ovalbox{
    34 \begin{minipage}{1\textwidth}
    35 \begin{center}
    36 \Huge{User Manual for the LMD Martian Mesoscale Model}
    37 \end{center}
    38 \end{minipage}
    39 }
    40 }
    41 %
    4214\begin{document}
    4315
    … …  
    4820%\pagestyle{fancy}
    4921\selectlanguage{english}
    50 \maketitle
     22%\maketitle
     23
     24\begin{titlepage}
     25\mbox{}\\
     26\Ovalbox{
     27\begin{minipage}{1\textwidth} \begin{center}
     28\Huge{LMD Martian Mesoscale Model}\\\huge{User Manual}
     29\end{center} \end{minipage}
     30}
     31\mbox{}\\
     32\begin{center}
     33\includegraphics[width=0.85\textwidth]{domain_100.png}
     34\end{center}
     35\mbox{}\\
     36\begin{center}
     37\huge{A. Spiga}\\
     38\email{aymeric.spiga}{upmc.fr}\\
     39Laboratoire de M\'et\'eorologie Dynamique\\
     40%Institut Pierre Simon Laplace\\
     41Universit\'e Pierre et Marie Curie, Paris, France\\
     42\end{center}
     43\end{titlepage}
     44
     45
    5146\clearemptydoublepage
    5247

  • trunk/MESOSCALE/PLOT/PYTHON/mylib/api_g95.sh

    r191 r194  
    99\rm loge
    1010f2py -c -m api ../../../LMD_MM_MARS/SRC/POSTPROC/api.F90 --fcompiler=g95 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include --f90flags="-Wall -Wno=112,141,137,155 -fno-second-underscore -ffree-form" > logc 2> loge
     11f2py -c -m timestuff time.F --fcompiler=g95

  • trunk/MESOSCALE/PLOT/PYTHON/mylib/api_pgf90.sh

    r191 r194  
    1313\rm loge
    1414f2py -c -m api ../../../LMD_MM_MARS/SRC/POSTPROC/api.F90 --fcompiler=pg -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include --f90flags="-Mfree" > logc 2> loge
     15f2py -c -m timestuff time.F --fcompiler=pg

  • trunk/MESOSCALE/PLOT/PYTHON/mylib/myplot.py

    r191 r194  
    180180    else:                  kcolor=color
    181181    if key: p = plt.quiverkey(q,posx,posy,scale,\
    182                    str(int(scale)),coordinates='data',color=kcolor,labelpos='S')
     182                   str(int(scale)),coordinates='data',color=kcolor,labelpos='S',labelsep = 0.03)
    183183    return
    184184

  • trunk/MESOSCALE/PLOT/PYTHON/scripts/winds.py

    r193 r194  
    2626    from netCDF4 import Dataset
    2727    from myplot import getcoord2d,define_proj,makeplotpng,simplinterv,vectorfield,ptitle,latinterv,getproj,wrfinterv,dumpbdy
     28    from mymath import deg
    2829    from matplotlib.pyplot import contourf, subplot, figure, rcParams, savefig, colorbar
    2930    from matplotlib.cm import get_cmap
     31    from timestuff import sol2ls
    3032    import numpy as np
    3133
    … …  
    7375    else:                    [wlon,wlat] = simplinterv(lon2d,lat2d)
    7476
    75     ##################
    76     ### Get local time
    77     if typefile in ['mesoapi','meso']:  ltst = int(getattr(nc, 'GMT') + 0.5*(wlon[0]+wlon[1])/15.)
    78     elif typefile in ['gcm']:           ltst = 0
     77    #########################
     78    ### Get local time and Ls
     79    if typefile in ['mesoapi','meso']: 
     80        ltst = int(getattr(nc, 'GMT') + 0.5*(wlon[0]+wlon[1])/15.)
     81        zetime = nc.variables['Times'][0]
     82        zetimestart = getattr(nc, 'START_DATE')
     83        zeday = int(zetime[8]+zetime[9]) - int(zetimestart[8]+zetimestart[9])
     84        if zeday < 0:    lschar=""  ## might have crossed a month... fix soon
     85        else:            lschar="_Ls"+str( int( sol2ls ( getattr( nc, 'JULDAY' ) + zeday ) ) )
     86    elif typefile in ['gcm']:           
     87        ltst = 0
     88        lschar = ""
    7989
    8090    ##############################################################################
    … …  
    217227       
    218228       ### Next subplot
    219        ptitle( basename + "_LT"+str((ltst+i)%24) )
     229       ptitle( basename + "_LT"+str((ltst+i)%24) + lschar )
    220230       sub += 1
    221231
Note: See TracChangeset for help on using the changeset viewer.