Changeset 207


Ignore:
Timestamp:
Jul 12, 2011, 11:51:13 AM (13 years ago)
Author:
aslmd
Message:

MESOSCALE: A GENERAL CLEAN-UP FOLLOWING UPDATING THE USER MANUAL. EVERYTHING ESSENTIAL IS IN MESOSCALE (much lighter than before). EVERYTHING FOR DEVELOPPERS OR EXPERTS IS IN MESOSCALE_DEV.

Location:
trunk
Files:
117 added
19 deleted
2 edited
7 copied
21 moved

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/copy_model

    r200 r207  
    1717
    1818mkdir WRFV2
     19ln -sf $WRFSRC/prepare_ini .
     20
    1921cd WRFV2
    2022

  • trunk/MESOSCALE/LMD_MM_MARS/makemeso

    r167 r207  
    174174              tra=1     
    175175            fi
     176            if [[ "${phys}" == "newphys_" ]]
     177            then
     178              echo "How many types of scatterer are you using ? 1 [dust] or 2 [dust + water ice]" ; read scat
     179            else
     180              scat=1 ## dummy stuff
     181            fi
    176182         fi
    177183  # folder
    … …  
    185191        if [[ "${scenario}" == "" ]]
    186192        then
    187           sed s+"mars_lmd"+"mars_lmd_new"+g             copy_model > copy_model_tmp
     193          sed s+"mars_lmd"+"mars_lmd_new"+g             SRC/copy_model > copy_model_tmp
    188194        else
    189195          echo SCENARIO ${scenario} DID YOU INCLUDE THE SPECIFIC folder ? ; read dummy
    190           sed s+"mars_lmd"+"mars_lmd_new_${scenario}"+g copy_model > copy_model_tmp
     196          sed s+"mars_lmd"+"mars_lmd_new_${scenario}"+g SRC/copy_model > copy_model_tmp
    191197             if [[ "${config}" == "les" ]]
    192198             then
    … …  
    217223     cd Registry ; Registry.bash ; cd ..
    218224    else
    219      #./copy_model
    220225         if [[ "${config}" == "les" ]]
    221226         then
    222227           ### NB: With physics: LES folder // Without physics : LESnophys_ folder
    223228           ### ---- because differences in 'modif' folder ----
    224            sed s+"PWD/SRC/"+"PWD/SRC/LES$phys/"+g copy_model > copy_model_tmp
     229           sed s+"PWD/SRC/"+"PWD/SRC/LES$phys/"+g SRC/copy_model > copy_model_tmp
    225230         else
    226            cp copy_model copy_model_tmp
     231           cp SRC/copy_model copy_model_tmp
    227232         fi   
    228233     chmod 755 copy_model_tmp
    … …  
    265270            fi
    266271    fi
    267   fi
    268   if [[ "${phys}" == "newphys_" ]]
    269   then
    270      echo "How many scatterers are you using : 1,2,..." ; read scat
    271   else
    272      scat=1 ## dummy stuff
    273272  fi
    274273

  • trunk/MESOSCALE_DEV/PLOT/PYTHON/mylib/api_wrapper.py

    r186 r207  
    1414                    extrapolate     = 0, \
    1515                    unstagger_grid  = False, \
    16                     onelevel        = 0.020  ):
     16                    onelevel        = 0.020, \
     17                    nocall          = False ):
    1718    import api
    1819    import numpy as np
    … …  
    2728    print input_name, output_name
    2829
    29     api.api_main ( path_to_input, input_name, path_to_output, output_name, \
     30    if nocall:     pass
     31    else:          api.api_main ( path_to_input, input_name, path_to_output, output_name, \
    3032                   process, fields, debug, bit64, oldvar, np.arange(299), \
    3133                   interp_method, extrapolate, unstagger_grid, onelevel )

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

    r204 r207  
    6868    return prefix
    6969
    70 def api_onelevel (  path_to_input   = None, \
    71                     input_name      = 'wrfout_d0?_????-??-??_??:00:00', \
    72                     path_to_output  = None, \
    73                     output_name     = 'output.nc', \
    74                     process         = 'list', \
    75                     fields          = 'tk,W,uvmet,HGT', \
    76                     debug           = False, \
    77                     bit64           = False, \
    78                     oldvar          = True, \
    79                     interp_method   = 4, \
    80                     extrapolate     = 0, \
    81                     unstagger_grid  = False, \
    82                     onelevel        = 0.020  ):
    83     import api
    84     import numpy as np
    85     if not path_to_input:   path_to_input  = './'
    86     if not path_to_output:  path_to_output = path_to_input
    87     api.api_main ( path_to_input, input_name, path_to_output, output_name, \
    88                    process, fields, debug, bit64, oldvar, np.arange (299), \
    89                    interp_method, extrapolate, unstagger_grid, onelevel )
    90     return
    91 
    9270def getproj (nc):
    9371    map_proj = getattr(nc, 'MAP_PROJ')
    … …  
    240218    elif wlat[0] <= -80.:  blat = -40.
    241219    else:                  blat = wlat[0]
    242     h = 2000.
     220    h = 50.  ## en km
    243221    radius = 3397200.
    244222    if   char == "cyl":     m = Basemap(rsphere=radius,projection='cyl',\
    … …  
    250228    elif char == "npstere": m = Basemap(rsphere=radius,projection='npstere', boundinglat=blat, lon_0=0.)
    251229    elif char == "spstere": m = Basemap(rsphere=radius,projection='spstere', boundinglat=blat, lon_0=0.)
     230    elif char == "nplaea":  m = Basemap(rsphere=radius,projection='nplaea', boundinglat=wlat[0], lon_0=meanlon)
     231    elif char == "laea":    m = Basemap(rsphere=radius,projection='laea',lon_0=meanlon,lat_0=meanlat,lat_ts=meanlat,\
     232                              llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])
    252233    elif char == "nsper":   m = Basemap(rsphere=radius,projection='nsper',lon_0=meanlon,lat_0=meanlat,satellite_height=h*1000.)
    253234    elif char == "merc":    m = Basemap(rsphere=radius,projection='merc',lat_ts=0.,\
    254235                              llcrnrlat=wlat[0],urcrnrlat=wlat[1],llcrnrlon=wlon[0],urcrnrlon=wlon[1])
    255236    fontsizemer = int(mpl.rcParams['font.size']*3./4.)
    256     if char in ["cyl","lcc","merc"]:   step = findstep(wlon)
    257     else:                              step = 10.
     237    if char in ["cyl","lcc","merc","nsper","laea"]:   step = findstep(wlon)
     238    else:                                             step = 10.
    258239    m.drawmeridians(np.r_[-180.:180.:step*2.], labels=[0,0,0,1], color='grey', fontsize=fontsizemer)
    259240    m.drawparallels(np.r_[-90.:90.:step], labels=[1,0,0,0], color='grey', fontsize=fontsizemer)

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

    r204 r207  
    7171        if not var:                                        back = "mola"    ## if no var:         draw mola
    7272        elif typefile in ['mesoapi','meso'] \
    73              and proj not in ['merc','lcc']:               back = "molabw"  ## if var but meso:   draw molabw
     73             and proj not in ['merc','lcc','nsper','laea']:       back = "molabw"  ## if var but meso:   draw molabw
    7474        else:                                              pass             ## else:              draw None
    7575
    … …  
    8282    elif typefile in ['gcm']:
    8383        [lon2d,lat2d] = getcoord2d(nc,nlat="latitude",nlon="longitude",is1d=True)
    84     if proj == "npstere":    [wlon,wlat] = latinterv("North_Pole")
    85     elif proj == "lcc":      [wlon,wlat] = wrfinterv(lon2d,lat2d)
    86     else:                    [wlon,wlat] = simplinterv(lon2d,lat2d)
     84    if proj == "npstere":             [wlon,wlat] = latinterv("North_Pole")
     85    elif proj in ["lcc","laea"]:      [wlon,wlat] = wrfinterv(lon2d,lat2d)
     86    else:                             [wlon,wlat] = simplinterv(lon2d,lat2d)
    8787    if zoom: 
    8888        dlon = abs(wlon[1]-wlon[0])/2.
    … …  
    125125        ### some already defined colormaps
    126126        if colorb is True:    colorb = defcolorb(var)
     127    else:
     128        dimension = 0
    127129 
    128130    ###########################
    … …  
    139141    if var:     basename = basename + var
    140142    ###
    141     if typefile is 'meso':                      stralt = "_lvl" + str(nvert)
    142     elif typefile is 'mesoapi': 
    143         zelevel = int(nc.variables['vert'][nvert])
    144         if 'altitude'       in nc.dimensions:   stralt = "_"+str(zelevel)+"m-AMR"
    145         elif 'altitude_abg' in nc.dimensions:   stralt = "_"+str(zelevel)+"m-ALS"
    146         elif 'bottom_top'   in nc.dimensions:   stralt = "_"+str(zelevel)+"m"
    147         elif 'pressure'     in nc.dimensions:   stralt = "_"+str(zelevel)+"Pa"
    148     else:                                       stralt = ""         
     143    if dimension == 4:
     144        if typefile is 'meso':                      stralt = "_lvl" + str(nvert)
     145        elif typefile is 'mesoapi': 
     146            zelevel = int(nc.variables['vert'][nvert])
     147            if 'altitude'       in nc.dimensions:   stralt = "_"+str(zelevel)+"m-AMR"
     148            elif 'altitude_abg' in nc.dimensions:   stralt = "_"+str(zelevel)+"m-ALS"
     149            elif 'bottom_top'   in nc.dimensions:   stralt = "_"+str(zelevel)+"m"
     150            elif 'pressure'     in nc.dimensions:   stralt = "_"+str(zelevel)+"Pa"
     151        else:                                       stralt = ""         
     152    else:
     153        stralt = ""
    149154    ###
    150155    basename = basename + stralt
    … …  
    175180            rcParams['font.size'] = int( rcParams['font.size'] * 1. / 2. )
    176181        elif numplot == 9:
    177             sub = 331 
     182            sub = 331
    178183            fig.subplots_adjust(wspace = 0.3, hspace = 0.3)
    179184            rcParams['font.size'] = int( rcParams['font.size'] * 1. / 2. )
    … …  
    181186            sub = 99999
    182187        else:
    183             print "supported: 1,2,3,4,6,8"
     188            print "supported: 1,2,3,4,6,8,9"
    184189            exit()
    185190        ### Prepare time loop
    … …  
    306311    parser.add_option('-i', action='store', dest='interp',      type="int",     default=None,  help='interpolation method (2: press, 3: z-amr, 4:z-als)')
    307312    parser.add_option('-c', action='store', dest='colorb',      type="string",  default=True,  help='change colormap')
    308     parser.add_option('-x', action='store_false', dest='winds',                 default=True,  help='flag: no wind vectors')
     313    parser.add_option('-x', action='store_false', dest='winds',                 default=True,  help='no wind vectors')
    309314    parser.add_option('-m', action='store', dest='vmin',        type="float",   default=None,  help='bounding minimum value for color plot')   
    310315    parser.add_option('-M', action='store', dest='vmax',        type="float",   default=None,  help='bounding maximum value for color plot')
    311316    parser.add_option('-T', action='store_true', dest='tile',                   default=False, help='draw a tiled plot (no blank zone)')
    312317    parser.add_option('-z', action='store', dest='zoom',        type="float",   default=None,  help='zoom factor in %')
     318    parser.add_option('-N', action='store_true', dest='nocall',                 default=False, help='do not recreate api file')
    313319    #parser.add_option('-V', action='store', dest='comb',        type="float",   default=None,  help='a defined combination of variables')
    314320    (opt,args) = parser.parse_args()
    … …  
    339345                                 fields          = zefields, \
    340346                                 interp_method   = opt.interp, \
    341                                  onelevel        = zelevel )
     347                                 onelevel        = zelevel, \
     348                                 nocall          = opt.nocall )
     349        print zefile
    342350        zelevel = 0 ## so that zelevel could play again the role of nvert
    343351
Note: See TracChangeset for help on using the changeset viewer.