Changeset 131 for trunk/mesoscale


Ignore:
Timestamp:
May 25, 2011, 4:19:36 PM (14 years ago)
Author:
aslmd
Message:

mars

M 130 mars/libf/phymars/callradite.F
M 130 mars/README
no more need to modify callradite.F prior to compilation
[but still dimradmars.h must be modified]

--> in callradite.F we now have

-- DEFAULT name_iaer(1) is "dust_conrath"
-- IF (doubleq.AND.active) name_iaer(1) = "dust_doubleq"
-- IF (water.AND.activice) name_iaer(2) = "h2o_ice"

M 130 000-USERS
small meeting about SVN for "planeto" team

M 130 mars/libf/phymars/meso_physiq.F
purely cosmetic changes

M 130 mesoscale/PLOT/MINIMAL/map_latlon.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwave2.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro
M 130 mesoscale/PLOT/SPEC/GW/gravitwave_inc.pro
A 0 mesoscale/PLOT/SPEC/MAP/defs/THARSIS_newphys_bugnoRAC_RACgcm_TAUTES.map_uvt_inc.pro
M 130 mesoscale/PLOT/SPEC/MAP/map_uvt.pro
graphic stuff for GW studies and H2O cycle

Location:
trunk/mesoscale/PLOT
Files:
1 added
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/PLOT/MINIMAL/map_latlon.pro

    r114 r131  
    375375;yeah=UNIQ(lev) & lev=lev[yeah]
    376376
     377colcont = 75
     378if (pal eq 16) then colcont=0 ;noir
     379if (pal eq 4) then colcont=255 ;blanc
    377380
    378381if (isotropic eq 'true') then begin
    … …  
    387390        max_value=max_contour, $
    388391        min_value=min_contour, $       
    389         c_labels=findgen(n_elements(lev))*0., $
     392        c_labels=findgen(n_elements(lev))*1., $
    390393;        c_charsize=!P.charsize/2, $
    391         color=75;0;50;75;100tropclair;150;0;, $
     394        color=colcont;75;0;50;75;100tropclair;150;0;, $
    392395;       /noerase, $
    393396;       xtitle=xtitleset, $
    … …  
    410413        min_value=min_contour, $
    411414        c_labels=findgen(n_elements(lev))*0., $
    412         color=0;, $
     415        color=colcont;, $
    413416;        /noerase, $
    414417;        xtitle=xtitleset, $

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro

    r117 r131  
    1313folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_2D.morepoints.160818/'
    1414folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544/'
    15 folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440/'
    16 folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468/'
     15folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints.161542/'
     16;folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440/'
     17;folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468/'
    1718file=folder+'wrfout_d01_9999-09-09_09:00:00_z'
    1819charvar='W' & charvarc='W'
    19 charvar='tk' & charvarc='W'
     20charvar='W' & charvarc='tk' & cond=0
     21;charvar='tk' & charvarc='W'
    2022;charvar='tk' & charvarc='W' & cond=1
    21 charvar='tk' & charvarc='tk' & cond=1
     23;charvar='tk' & charvarc='tk' & cond=1
     24;charvar='tk' & charvarc='tk' & cond=0
    2225;;charvar='tpot' & charvarc='tpot'
    23 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    24 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     26charvar='HR_NLTE' & charvarc='tk' & cond=0
     27;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     28        ;folder='/home/aslmd/GRAVITWAVE/'
     29        ;file=folder+'GW_MARS_highwind_3D_loweropacity_morepoints_diff.nc'
     30;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    2531;
    2632;
    … …  
    3945        charvar=charvarc, $
    4046        invar=invarc
     47;getcdf, $
     48;        file=file, $
     49;        charvar='vert', $
     50;        invar=vert
    4151getcdf, $
    42         file=file, $
    43         charvar='vert', $
    44         invar=vert
     52         file='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D.157077/wrfout_d01_9999-09-09_09:00:00_z', $
     53         charvar='vert', $
     54         invar=vert
    4555getcdf, $
    4656        file=file, $
    … …  
    5363;for nt=27,27 do begin
    5464for nt=1,40 do begin
     65;for nt=9,9 do begin
     66;for nt=0,20 do begin
    5567zefile='gravitwave_'+charvar+'_'+charvarc+'_'+string(nt+100,'(I0)')
    56 if (n_elements(cond) ne 0) then zefile = 'cond_' + zefile
     68if (n_elements(cond) ne 0) then if (cond ne 0) then zefile = 'cond_' + zefile
    5769;
    5870;
    5971;
    60   PS_Start, filename=zefile+'.ps'
    61   print, zefile+'.ps'
    62   !P.Charsize = 1.2
    63   !p.charthick = 2.0
    64   !p.thick = 2.0
    65   !x.thick = 2.0
    66   !y.thick = 2.0
     72  PS_Start, filename=zefile+'.ps' & print, zefile+'.ps'
    6773;
    6874;
    6975;
    7076s = size(invar) & middle = floor(s[2]/2) & print, 'PLOT at y subs ', middle
    71 what_I_plot = reform(invar(*,middle,*,nt))
     77what_I_plot = reform(invar(*,middle,*,nt)) & print, min(what_I_plot), max(what_I_plot)
    7278overcontour = reform(invarc(*,middle,*,nt))
    7379   if (n_elements(cond) ne 0) then begin
     80   if (cond ne 0) then begin
    7481     column = columnp
    7582     ;column = reform(columnp(*,1,*,1))
    … …  
    8087     if (charvarc eq 'tk') then overcontour = overcontour - yeye
    8188   endif
     89   endif
    8290xx = findgen(n_elements(what_I_plot(*,0)))
    8391zz = vert / 1000.
    84 
    85         ;w = where(what_I_plot lt 0.)
    86         w = where(what_I_plot lt -1.)
    87         ;w = where(what_I_plot lt -3.)
    88         if (w(0) ne -1) then begin
    89         for i=0,n_elements(w)-1 do begin
    90          ind = ARRAY_INDICES(what_I_plot, w(i))
    91          print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
    92         endfor
    93         endif
     92;
     93        ;;w = where(what_I_plot lt 2.)
     94        ;;w = where(what_I_plot lt 0.)
     95        ;w = where(what_I_plot lt -1.)
     96        ;;w = where(what_I_plot lt -3.)
     97        ;if (w(0) ne -1) then begin
     98        ;for i=0,n_elements(w)-1 do begin
     99        ; ind = ARRAY_INDICES(what_I_plot, w(i))
     100        ; print, 'x '+string(ind(0),'(I0)')+' z '+string(zz(ind(1)),'(I0)')+' t '+string(nt*925./3700.,'(F4.1)'), what_I_plot(ind[0],ind[1])
     101        ;endfor
     102        ;endif
     103;
     104;
     105;
     106if (charvar eq 'HR_NLTE') then begin
     107        what_I_plot = what_I_plot * 3700. ;; en K/hour
     108        ;w = where((zz le 60.) and (zz ge 50.)) ;; pour enlever la ligne d'altitude ou le NLTE est active
     109        ;if (w(0) ne -1) then what_I_plot(*,w) = 0.
     110endif
    94111;
    95112;
    … …  
    99116        xx, $                                   ; 1D latitude
    100117        zz, $                                   ; 1D longitude
    101 ;        minfield=minfield_init, $               ; minimum value of plotted field (=0: calculate)
    102 ;        maxfield=maxfield_init, $               ; maximum value of plotted field (=0: calculate)
     118;        minfield_init, $               ; minimum value of plotted field (=0: calculate)
     119;        maxfield_init, $               ; maximum value of plotted field (=0: calculate)
    103120        overcontour=overcontour, $              ; another 2D field to overplot with contour lines (=0: no)
    104121;        overvector_x=overvector_x, $            ; wind vector - x component (=0: no)

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwave_inc.pro

    r114 r131  
    22;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    33;;;;;;; W
    4 colors        = 32
    5 minfield_init = -6.
    6 maxfield_init =  6.
    7 format        = '(F4.1)'
    8 pal           = 4
     4colors        = 64
     5minfield_init = -10.
     6maxfield_init = 10.
     7ndiv          = 10
     8format        = '(I0)'
     9pal           = 17 ;4 ;16 ;18 ;4
    910title_user    = 'Vertical velocity (m s!U-1!N)'
     11lev=80.+10.*findgen(30)
    1012;;;;;;;; tk
    11 minfield_init =  110.
    12 maxfield_init =  210.
    13 format        = '(I0)'
    14 pal            = 33
    15 title_user    = 'Temperature (K)'
    16 ;;;;;;; tk - tsat
    17 minfield_init =  -4. ;-8.
    18 maxfield_init =  4. ;+8.
    19 ndiv          =  8     
    20 format        = '(F4.1)'
    21 pal           =  0
    22 title_user    = 'T!Datm!N - T!Dsat!N (K)'
    23 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     13;minfield_init = 110.
     14;maxfield_init = 230.
     15;ndiv         = 12
     16;format        = '(I0)'
     17;pal          = 33
     18;title_user    = 'Temperature (K)'
     19;lev = [-10.,-8.,-6.,-4.,-2.,2.,4.,6.,8.,10.]
     20;lev=80.+10.*findgen(30)
     21;;;;;;;; tk - tsat
     22;minfield_init =  -4. ;-8.
     23;maxfield_init =  4. ;+8.
     24;ndiv         =  8     
     25;format        = '(F4.1)'
     26;pal           =  0
     27;title_user    = 'T!Datm!N - T!Dsat!N (K)'
     28;lev = [-6.,-5.,-4.,-3.,-2.,-1.,0.]
     29;;;;;;;;; hr_nlte
     30;       minfield_init =  -6.
     31;       maxfield_init =  6.
     32;       ndiv          =  6
     33;       format        = '(F4.1)'
     34;       pal           =  33
     35;       title_user    = 'NLTE heating rate (K / hour)'
     36lev=[-20.,-18.,-16.,-14.,-12.,-10.,-8.,-6.,-4.,-2.,2.,4.,6.,8.,10.,12.,14.,16.,18.,20.]
     37title_user    = 'NLTE cooling rate (K hour!U-1!N)'
     38w=where(abs(what_I_plot) le 1e5) & if (w[0] ne -1) then what_I_plot[w] = - what_I_plot[w]
     39minfield_init = 0.
     40maxfield_init = 20.
     41ndiv          = 10
     42pal           = 1 ;4 ;22
     43lev=80.+10.*findgen(30)
     44format        = '(I0)'
     45;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    2446missing_value=1.e30
    2547lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max
    … …  
    3456windowy = [0.,130.]
    3557windowx = [10., 110.]
     58;windowy = [50.,130.]
    3659intervalx = 10.
    3760intervaly = 10.
    … …  
    3962poscb=0.55
    4063;isotropic='false'
    41 lev = 10.*findgen(100)
    42 lev = [-6.,-4.,-2.,2.,4.,6.]
    43 lev = [-8.,-5.,-1.,1.,5.,8.]
    44 lev = [-6.,-5.,-4.,-3.,-2.,-1.,0.]
     64;lev = 10.*findgen(100)
     65;lev = [-6.,-4.,-2.,2.,4.,6.]
     66;;lev = [-8.,-5.,-1.,1.,5.,8.]
     67;;lev = [-6.,-5.,-4.,-3.,-2.,-1.,0.]
    4568;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    4669;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    47 

  • trunk/mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro

    r117 r131  
    2828experiment='GW_MARS_highwind_3D_loweropacity_lotspoints.161468'
    2929experiment='GW_MARS_highwind_3D_loweropacity_morepoints.161542'
    30 experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544'
    31 experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440'
    32 experiment='GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468'
     30;experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544'
     31;experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440'
     32;experiment='GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468'
    3333        file=path+experiment+'_wrfout_d01_9999-09-09_09:00:00_z'
    3434file=path+experiment+'/wrfout_d01_9999-09-09_09:00:00_z'
    … …  
    5757;nts=1  & nte=1  & nxs=6   & nxe=6
    5858nts=1  & nte=100  & nxs=50   & nxe=50
     59nts=1  & nte=100  & nxs=60   & nxe=60
     60nts=1  & nte=100  & nxs=20   & nxe=20
    5961;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    6062;

  • trunk/mesoscale/PLOT/SPEC/MAP/map_uvt.pro

    r96 r131  
    148148endelse
    149149help, what_I_plot, lon, lat
    150 print, min(what_I_plot)
    151 print, max(what_I_plot)
    152 print, min(overcontour)
    153 print, max(overcontour)
    154 print, min(overvector_x)
    155 print, min(overvector_y)
     150print, 'COLOR PLOT min/max ', min(what_I_plot), max(what_I_plot)
     151print, 'CONTOUR PLOT min/max ', min(overcontour), max(overcontour)
     152print, 'VECTOR PLOT min/max ', min(overvector_x), min(overvector_y)
    156153;
    157154;
Note: See TracChangeset for help on using the changeset viewer.