Ignore:
Timestamp:
Mar 8, 2011, 5:59:28 PM (14 years ago)
Author:
aslmd
Message:

mars: test outliers [dans initracer.F, commente] LMD_MM_MARS: modifications mineures [retrocompatibilite ancienne physique, callphys.def test pour nouvelle physique] PLOT: generalisation de la routine map_uvt pour pouvoir tracer des figures en projection polaire complete

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/PLOT/SPEC/MAP/defs/map_uvt_inc.pro

    r85 r91  
    1 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; FOLDER
    2 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; FOLDER
    3         folder  = '/d5/aslmd/LMD_MM_MARS_SIMUS/OM/'
    4         coord2d = 'false'       ;; for non-regular projections 
    5 filename        = folder + 'OM6_TI85/wrfout_d01_2024-06-43_06:00:00_zabg'
    6 save_ps         = 'OM6_TI85_tsurf'
    7 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TIME
    8 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TIME
    9 ini_utc         = 6     ;; cf. name of file
    10 freq            = 1     ;; cf. 1 output per ... hour
    11 utc_to_lt       = -9    ;; cf. longitude -- LT = UTC + utc_to_lt
    12         use_lt  = 2     ;; cf. what user wants
    13         use_utc = use_lt - utc_to_lt
    14         ntime   = (use_utc - ini_utc)/freq      ;; TRUE IDL SUBSCRIPT ... ajouter un modulo...
    15         print, 'CHECK ... ', filename, ntime
    16 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; LVL & FLD
    17 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; LVL & FLD
    18 nlevel          = 0
    19 field1          = 'TSURF'       ;; comment to trace horizontal velocity
    20 no3d            = 'true'
    21 overvector_x    = 0             ;; comment out to get rid of vectors
    22 overvector_y    = 0             ;; comment out to get rid of vectors
    23 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TWEAK VAR
    24 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TWEAK VAR
    25         ;a0 = 65.052165 & a1 = 3.1228993 & a2 = 0.0053787417 ;; Ls ~ 173
    26         ;;a0 = 64.039300 & a1 = 3.1378104 & a2 = 0.0055369148 ;; Ls ~ 120
    27         ;what_I_plot = - ( alog ( a1 - what_I_plot / a0 ) ) /  a2
    28         ;print, max(what_I_plot), min(what_I_plot)
    29         ;overvector_x=0
    30         ;overvector_y=0
    31 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    32 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     1;;;;;*************************************** FOLDER
     2        folder  = '/tmp7/aslmd/'
     3        coord2d = 'false'       ;; for non-regular projections  ;;CHANGER windowx and windowy
     4        coord2d = 'polar'       ;; for polar projection with map_set ;; utiliser use_lt=99
     5        filename        = folder + 'wrfout_d01_2024-05-02_01:00:00'
     6        ini_utc         = 1           ;; cf. name of file
     7        freq            = 1           ;; cf. 1 output per ... hour
     8        save_ps         = 'POLAR_newphys_RAC_TAUTES'
    339
     10;;;;;*************************************** LVL & FLD & TIME
     11        nlevel          = 0
     12        use_lt          = 99               ;; cf. what user wants (99 pour tous)
     13        field1          = 'TAU_ICE'        ;; comment to trace horizontal velocity
     14        ;field2         = 'XLAT'           ;; contour
     15        no3d            = 'true'
     16        ;overvector_x    = 0               ;; comment out to get rid of vectors
     17        ;overvector_y    = 0               ;; comment out to get rid of vectors
     18
     19;;;;;*************************************** TWEAK VAR
     20        ;what_I_plot = what_I_plot * 0.
     21        ;print, max(what_I_plot), min(what_I_plot)
    3422
    3523;;;;;*************************************** PLOT TITLES
    36         ;title_user     =       'Surface temperature (K) and Winds 10m ABG (m s!U-1!N)'
    37         ;title_user     =       'Wind-induced apparent thermal inertia (J m!U-2!N s!U-0.5!N K!U-1!N)'
    38         ;title_user     =       'Winds 10m ABG (m s!U-1!N)'
    39         title_user      =       'Surface temperature (K)'
    40         ;subtitle_user   =       'LT = 03:00am  /  Ls = 173!Uo!N  /  dx = 10km [single]  /  Uniform TI = 85 J m!U-2!N s!U-0.5!N K!U-1!N'
    41         subtitle_user   =       'LT = 02:00am  /  Ls = 173!Uo!N  /  dx = 6km [single]  /  Uniform TI = 85 J m!U-2!N s!U-0.5!N K!U-1!N'
     24        title_user      =       'Water ice cloud optical depth at 825 cm!U-1!N'
    4225        title_axis      =       ['East longitude','North latitude']
     26        ;subtitle_user   = 'LMD_MM'
     27        ;;subtitle_user   = subtitle_user + ' / LT = '+string(use_lt,'(I0)')+'h'
     28        ;;subtitle_user   = subtitle_user + ' / UTC!D0!N = 06:00am' ;; pour les series
     29        ;subtitle_user   = subtitle_user + ' / Ls = 120!Uo!N'
     30        ;subtitle_user   = subtitle_user + ' / dx = 10km [S]'
    4331
    4432;;;;;*************************************** COLOR TABLES
    4533        flag_cb         =       'true'
    46         poscb           =       0.80;0.96
    47         format          =       '(I0)'
    48         ;format          =       '(F4.1)'
    49         ;colors          =       128
     34        poscb           =       0.85
     35        format          =       '(F5.2)'
    5036        colors          =       32
    5137        pal             =       22              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
    5238
    5339;;;;;*************************************** FILL LIMITS
    54         minfield_init   =       150.
    55         maxfield_init   =       180.
    56         ;minfield_init  =       50.
    57         ;maxfield_init  =       250.
     40        minfield_init   =       0.
     41        maxfield_init   =       0.55
     42        ndiv            =       11
    5843
    5944        ;;;;;*************************************** LIMIT TRICKS
    … …  
    6550
    6651;;;;;*************************************** WINDS
    67         ;windex          =       10.            ;; DEF: 20.
    6852        windex          =       30.             ;; DEF: 20.
    69         ;stride          =       1.             ;; DEF: 5.
    70         ;stride          =       2.              ;; DEF: 5.
    7153        stride          =       3.              ;; DEF: 5.
    7254
    7355;;;;;*************************************** CONTOUR
    7456        overcontour     =       overcontour/1000.
    75         ;lev             =       50. + 50.*findgen(20)
    76         ;lev             =       -10. + 0.2*findgen(20)
    77         lev             =       -10. + 2.*findgen(20)
    78         ;lev             =       -10. + 1.*findgen(40)
     57        lev             =       -10. + 0.5*findgen(80)
    7958
    8059;;;;;*************************************** AXIS
    8160        isotropic       =       'true'
    82         ;intervalx       =       0.5
    83         intervalx       =       2.0
    84         ;intervaly       =      intervalx
    85         intervaly       =       1.0
     61        intervalx       =       30.0
     62        intervaly       =       05.0
    8663
    8764;;;;;*************************************** MAP LIMITS
    88         ;windowx         =       [-144.,-126.]
    89         ;windowy         =       [10.,26.]
    90         windowx         =       [-146.,-126.]       
    91         ;windowy         =      [10.,28.]
    92         ;windowx         =      [-146.,-128.]
    93         ;windowy         =      [12.,26.]
    94         windowy         =       [11.,27.]
     65        windowx         =       [-180.,180.]       
     66        windowy         =       [65.,90.]
     67
     68;;;;;***************************************;;;;;
     69;;;;;***************************************;;;;;
     70;;;;;***************************************;;;;;
     71;;;;;***************************************;;;;;
     72
     73;;;;;*************************************** SETTING TIME (do not modify)
     74        if (n_elements(windowx) eq 0) then windowx = 0
     75        utc_to_lt       = mean(windowx) / 15.    ;; cf. longitude -- LT = UTC + utc_to_lt
     76        use_utc         = use_lt - utc_to_lt
     77        zentime         = floor(((24 + use_utc - ini_utc) MOD 24)/freq)        ;; TRUE IDL SUBSCRIPT...
     78        if (use_lt ne 99) then ntime = zentime else ntime = 99    ;; ou commenter pour avoir tous les pas de temps
     79        print, zentime, ntime
    9580
    9681;;;;;***************************************
Note: See TracChangeset for help on using the changeset viewer.