Changeset 12


Ignore:
Timestamp:
Jan 24, 2011, 1:57:13 PM (14 years ago)
Author:
aslmd
Message:

spiga@svn-planeto:les copains j ai fait ce changement

Location:
trunk/mesoscale/LMD_MM_MARS
Files:
4 added
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/LMD_MM_MARS/SRC/WPS/geogrid/GEOGRID.TBL.ARW

    r11 r12  
    77        df_dx=SLPX
    88        df_dy=SLPY
    9         smooth_option = smth-desmth_special; smooth_passes=1
    10         fill_missing=0.
    11         interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
    12         interp_option =      2m:four_pt
    13         interp_option =      5m:four_pt
    14         interp_option =     10m:four_pt
    15         interp_option = default:four_pt
    16         rel_path=     30s:topo_30s/
    17         rel_path=      2m:topo_2m/
    18         rel_path=      5m:topo_5m/
    19         rel_path=     10m:topo_10m/
    20         rel_path= default:topo_2m/
     9        fill_missing=9000.
     10        smooth_option = 1-2-1; smooth_passes=5
     11#        smooth_option = smth-desmth_special; smooth_passes=1
     12        interp_option =         default:four_pt
     13        interp_option =         32ppd:sixteen_pt+four_pt
     14        interp_option =         32ppd_HRalb:four_pt
     15        interp_option =         gcm:four_pt
     16        interp_option =         64ppd:four_pt
     17        interp_option =         64ppd_noHRti:four_pt
     18        interp_option =         128ppd:average_gcell(4.0)+four_pt+average_4pt
     19        interp_option =         custom:four_pt
     20        interp_option =         omega:four_pt
     21        rel_path=               default:mola_topo32/
     22        rel_path=               32ppd:mola_topo32/
     23        rel_path=               32ppd_HRalb:mola_topo32/
     24        rel_path=               gcm:mola_GCM/
     25        rel_path=               64ppd:mola_topo64/
     26        rel_path=               64ppd_noHRti:mola_topo64/
     27        rel_path=               custom:crater2/
     28        rel_path=               omega:mola_topo64/
    2129===============================
    2230name = HGT_U
    … …  
    2432        priority = 1
    2533        dest_type = continuous
    26         interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
    27         interp_option =      2m:four_pt
    28         interp_option =      5m:four_pt
    29         interp_option =     10m:four_pt
    30         interp_option = default:four_pt
    31         smooth_option = smth-desmth_special; smooth_passes=1
    32         fill_missing=0.
    33         rel_path=     30s:topo_30s/
    34         rel_path=      2m:topo_2m/
    35         rel_path=      5m:topo_5m/
    36         rel_path=     10m:topo_10m/
    37         rel_path= default:topo_2m/
     34        fill_missing=9000.
     35        smooth_option = 1-2-1; smooth_passes=5
     36        interp_option =         default:four_pt
     37        interp_option =         32ppd:four_pt
     38        interp_option =         32ppd_HRalb:four_pt
     39        interp_option =         64ppd:four_pt
     40        interp_option =         64ppd_noHRti:four_pt
     41        interp_option =         128ppd:average_gcell(4.0)+four_pt+average_4pt
     42        interp_option =         custom:four_pt
     43        interp_option =         omega:four_pt
     44        interp_option =         gcm:four_pt
     45        rel_path=               default:mola_topo32/
     46        rel_path=               32ppd:mola_topo32/
     47        rel_path=               32ppd_HRalb:mola_topo32/
     48        rel_path=               gcm:mola_GCM/
     49        rel_path=               64ppd:mola_topo64/
     50        rel_path=               64ppd_noHRti:mola_topo64/
     51        rel_path=               custom:crater2/
     52        rel_path=               omega:mola_topo64/
    3853===============================
    3954name = HGT_V
    … …  
    4156        priority = 1
    4257        dest_type = continuous
    43         interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
    44         interp_option =      2m:four_pt
    45         interp_option =      5m:four_pt
    46         interp_option =     10m:four_pt
    47         interp_option = default:four_pt
    48         smooth_option = smth-desmth_special; smooth_passes=1
    49         fill_missing=0.
    50         rel_path=     30s:topo_30s/
    51         rel_path=      2m:topo_2m/
    52         rel_path=      5m:topo_5m/
    53         rel_path=     10m:topo_10m/
    54         rel_path= default:topo_2m/
    55 ===============================
    56 name=LANDUSEF
    57         priority=1
    58         dest_type=categorical
    59         z_dim_name=land_cat
    60         landmask_water=16               # Calculate a landmask from this field
    61         dominant=LU_INDEX
    62         interp_option =     30s:nearest_neighbor
    63         interp_option =      2m:four_pt
    64         interp_option =      5m:four_pt
    65         interp_option =     10m:four_pt
    66         interp_option = default:four_pt
    67         rel_path=     30s:landuse_30s/
    68         rel_path=      2m:landuse_2m/
    69         rel_path=      5m:landuse_5m/
    70         rel_path=     10m:landuse_10m/
    71         rel_path= default:landuse_2m/
    72 ===============================
    73 name=SOILTEMP
    74         priority=1
    75         dest_type=continuous
    76         interp_option=default:sixteen_pt+four_pt+average_4pt+average_16pt+search
    77         masked=water
    78         fill_missing=0.
    79         rel_path=default:soiltemp_1deg/
    80 ===============================
    81 name=SOILCTOP
    82         dominant=SOILCAT
    83         priority=1
    84         dest_type=categorical
    85         z_dim_name=soil_cat
    86         interp_option =     30s:nearest_neighbor
    87         interp_option =      2m:four_pt
    88         interp_option =      5m:four_pt
    89         interp_option =     10m:four_pt
    90         interp_option = default:four_pt
    91         rel_path=     30s:soiltype_top_30s/
    92         rel_path=      2m:soiltype_top_2m/
    93         rel_path=      5m:soiltype_top_5m/
    94         rel_path=     10m:soiltype_top_10m/
    95         rel_path= default:soiltype_top_2m/
    96 ===============================
    97 name=SOILCBOT
    98         priority=1
    99         dest_type=categorical
    100         z_dim_name=soil_cat
    101         interp_option =     30s:nearest_neighbor
    102         interp_option =      2m:four_pt
    103         interp_option =      5m:four_pt
    104         interp_option =     10m:four_pt
    105         interp_option = default:four_pt
    106         rel_path=     30s:soiltype_bot_30s/
    107         rel_path=      2m:soiltype_bot_2m/
    108         rel_path=      5m:soiltype_bot_5m/
    109         rel_path=     10m:soiltype_bot_10m/
    110         rel_path= default:soiltype_bot_2m/
    111 ===============================
    112 name=ALBEDO12M
    113         priority=1
    114         dest_type=continuous
    115         z_dim_name=month
    116         masked = water
    117         fill_missing = 8.
    118         interp_option=default:four_pt+average_4pt+average_16pt+search
    119         rel_path=default:albedo_ncep/
    120 ===============================
    121 name=GREENFRAC
    122         priority=1
    123         dest_type=continuous
    124         interp_option=default:four_pt+average_4pt+average_16pt+search
    125         z_dim_name=month
    126         masked = water
    127         fill_missing = 0.
    128         rel_path=default:greenfrac/
    129 ===============================
    130 name=SNOALB
    131         priority=1
    132         dest_type=continuous
    133         interp_option=default:four_pt+average_4pt+average_16pt+search
    134         masked = water
    135         fill_missing = 0.
    136         rel_path=default:maxsnowalb/
    137 ===============================
    138 name=SLOPECAT
    139         priority=1
    140         dominant_only=SLOPECAT
    141         dest_type=categorical
    142         z_dim_name=slope_cat
    143         masked = water
    144         fill_missing = 0.
    145         interp_option=default:nearest_neighbor+average_16pt+search
    146         rel_path=default:islope/
    147 ===============================
     58        fill_missing=9000.
     59        smooth_option = 1-2-1; smooth_passes=5
     60        interp_option =         default:four_pt
     61        interp_option =         32ppd:four_pt
     62        interp_option =         32ppd_HRalb:four_pt
     63        interp_option =         64ppd:four_pt
     64        interp_option =         64ppd_noHRti:four_pt
     65        interp_option =         128ppd:average_gcell(4.0)+four_pt+average_4pt
     66        interp_option =         custom:four_pt
     67        interp_option =         omega:four_pt
     68        interp_option =         gcm:four_pt
     69        rel_path=               default:mola_topo32/
     70        rel_path=               32ppd:mola_topo32/
     71        rel_path=               32ppd_HRalb:mola_topo32/
     72        rel_path=               gcm:mola_GCM/
     73        rel_path=               64ppd:mola_topo64/
     74        rel_path=               64ppd_noHRti:mola_topo64/
     75        rel_path=               custom:crater2/
     76        rel_path=               omega:mola_topo64/
     77===============================
     78name = THERMAL_INERTIA
     79        priority = 1
     80        dest_type = continuous
     81        smooth_option = 1-2-1; smooth_passes=5
     82        interp_option =         32ppd:sixteen_pt+four_pt+average_4pt+average_16pt+search
     83        rel_path=               32ppd:thermal_GCM/
     84        interp_option =         32ppd_HRalb:sixteen_pt+four_pt+average_4pt+average_16pt+search
     85        rel_path=               32ppd_HRalb:thermal_GCM/
     86        interp_option =         64ppd:sixteen_pt+four_pt
     87        rel_path=               64ppd:thermal_TES/
     88        interp_option =         64ppd_noHRti:sixteen_pt+four_pt+average_4pt+average_16pt+search
     89        rel_path=               64ppd_noHRti:thermal_GCM/
     90        interp_option =         custom:sixteen_pt+four_pt
     91#       rel_path=               custom:thermal_TES/
     92        rel_path=               custom:thermal_GCM/
     93        interp_option =         omega:sixteen_pt+four_pt
     94        rel_path=               omega:thermal_GCM/
     95        interp_option =         gcm:four_pt
     96        rel_path=               gcm:thermal_GCM/
     97#       fill_missing=400. #north polar TI (needed with thermal_TES)
     98#        fill_missing=180. #polar projection artefacts
     99===============================
     100name = ALBEDO_GCM
     101        priority = 1
     102        dest_type = continuous
     103        smooth_option = 1-2-1; smooth_passes=5
     104#        fill_missing=0.25
     105#        interp_option = default:sixteen_pt+four_pt+average_4pt+average_16pt+search
     106        interp_option = 32ppd:four_pt+average_4pt+average_16pt+search
     107        rel_path =      32ppd:albedo_GCM/
     108        interp_option = 32ppd_HRalb:sixteen_pt+four_pt
     109        rel_path=       32ppd_HRalb:albedo_TES/
     110        interp_option = 64ppd:sixteen_pt+four_pt
     111        rel_path=       64ppd:albedo_TES/
     112        interp_option = 64ppd_noHRti:sixteen_pt+four_pt
     113        rel_path=       64ppd_noHRti:albedo_TES/
     114#        fill_missing=0.35 #north polar alb (needed with albedo_TES - albedo_GCM/1.3 to be consistent)
     115        interp_option = custom:sixteen_pt+four_pt
     116        rel_path=       custom:albedo_TES/
     117        interp_option = omega:sixteen_pt+four_pt
     118        rel_path=       omega:albedo_OMG/
     119        interp_option = gcm:four_pt
     120        rel_path =      gcm:albedo_GCM/
     121#===============================
     122#name=LANDUSEF
     123#        priority=1
     124#        dest_type=categorical
     125#        z_dim_name=land_cat
     126#        landmask_water=16               # Calculate a landmask from this field
     127#        dominant=LU_INDEX
     128#        interp_option =     30s:nearest_neighbor
     129#        interp_option =      2m:four_pt
     130#        interp_option =      5m:four_pt
     131#        interp_option =     10m:four_pt
     132#        interp_option = default:four_pt
     133#        rel_path=     30s:landuse_30s/
     134#        rel_path=      2m:landuse_2m/
     135#        rel_path=      5m:landuse_5m/
     136#        rel_path=     10m:landuse_10m/
     137#        rel_path= default:landuse_2m/
     138#===============================
     139#name=SOILTEMP
     140#        priority=1
     141#        dest_type=continuous
     142#        interp_option=default:sixteen_pt+four_pt+average_4pt+average_16pt+search
     143#        masked=water
     144#        fill_missing=0.
     145#        rel_path=default:soiltemp_1deg/
     146#===============================
     147#name=SOILCTOP
     148#        dominant=SOILCAT
     149#        priority=1
     150#        dest_type=categorical
     151#        z_dim_name=soil_cat
     152#        interp_option =     30s:nearest_neighbor
     153#        interp_option =      2m:four_pt
     154#        interp_option =      5m:four_pt
     155#        interp_option =     10m:four_pt
     156#        interp_option = default:four_pt
     157#        rel_path=     30s:soiltype_top_30s/
     158#        rel_path=      2m:soiltype_top_2m/
     159#        rel_path=      5m:soiltype_top_5m/
     160#        rel_path=     10m:soiltype_top_10m/
     161#        rel_path= default:soiltype_top_2m/
     162#===============================
     163#name=SOILCBOT
     164#        priority=1
     165#        dest_type=categorical
     166#        z_dim_name=soil_cat
     167#        interp_option =     30s:nearest_neighbor
     168#        interp_option =      2m:four_pt
     169#        interp_option =      5m:four_pt
     170#        interp_option =     10m:four_pt
     171#        interp_option = default:four_pt
     172#        rel_path=     30s:soiltype_bot_30s/
     173#        rel_path=      2m:soiltype_bot_2m/
     174#        rel_path=      5m:soiltype_bot_5m/
     175#        rel_path=     10m:soiltype_bot_10m/
     176#        rel_path= default:soiltype_bot_2m/
     177#===============================
     178#name=ALBEDO12M
     179#        priority=1
     180#        dest_type=continuous
     181#        z_dim_name=month
     182#        masked = water
     183#        fill_missing = 8.
     184#        interp_option=default:four_pt+average_4pt+average_16pt+search
     185#        rel_path=default:albedo_ncep/
     186#===============================
     187#name=GREENFRAC
     188#        priority=1
     189#        dest_type=continuous
     190#        interp_option=default:four_pt+average_4pt+average_16pt+search
     191#        z_dim_name=month
     192#        masked = water
     193#        fill_missing = 0.
     194#        rel_path=default:greenfrac/
     195#===============================
     196#name=SNOALB
     197#        priority=1
     198#        dest_type=continuous
     199#        interp_option=default:four_pt+average_4pt+average_16pt+search
     200#        masked = water
     201#        fill_missing = 0.
     202#        rel_path=default:maxsnowalb/
     203#===============================
     204#name=SLOPECAT
     205#        priority=1
     206#        dominant_only=SLOPECAT
     207#        dest_type=categorical
     208#        z_dim_name=slope_cat
     209#        masked = water
     210#        fill_missing = 0.
     211#        interp_option=default:nearest_neighbor+average_16pt+search
     212#        rel_path=default:islope/
     213#===============================
     214

  • trunk/mesoscale/LMD_MM_MARS/SRC/WPS/geogrid/src/constants_module.f90

    r11 r12  
    77
    88   real, parameter :: PI = 3.141592653589793
    9    real, parameter :: OMEGA_E = 7.292e-5 ! Angular rotation rate of the earth
     9!   real, parameter :: OMEGA_E = 7.292e-5 ! Angular rotation rate of the earth
     10!*** Mars
     11   real, parameter :: OMEGA_E = 7.0721E-5 ! Angular rotation rate of Mars
     12!*** Mars
    1013
    1114   real, parameter :: DEG_PER_RAD = 180./PI
    … …  
    1417   ! Mean Earth Radius in m.  The value below is consistent
    1518   ! with NCEP's routines and grids.
     19!************************************************************************
     20! Earth ellipsoid: only used on specific polar stereographic projections
     21! - equatorial axis
    1622   real, parameter :: A_WGS84 = 6378137.
     23! - polar axis
    1724   real, parameter :: B_WGS84 = 6356752.314
    1825   real, parameter :: RE_WGS84 = A_WGS84
    1926   real, parameter :: E_WGS84 = 0.081819192
    20    real, parameter :: EARTH_RADIUS_M = 6370000.   ! same as MM5 system
     27!************************************************************************
     28!   real, parameter :: EARTH_RADIUS_M = 6370000.   ! same as MM5 system
     29!*** Mars
     30!!necessary if WRFSI intermediate file is used
     31!!instead of WPS intermediate file format
     32   real, parameter :: EARTH_RADIUS_M = 3397200.   ! same as MM5 system
     33!*** Mars
    2134   real, parameter :: EARTH_CIRC_M = 2.*PI*EARTH_RADIUS_M
    2235

  • trunk/mesoscale/LMD_MM_MARS/SRC/WPS/metgrid/METGRID.TBL.ARW

    r11 r12  
    5454        fill_lev = 4 : SOILM040(200100)     
    5555        fill_lev = 5 : SOILM160(200100)     
    56         fill_lev = 6 : SOILM300(200100)     
     56        fill_lev = 6 : SOILM300(200100)     
    5757========================================
    5858name=SOILT
    … …  
    9797========================================
    9898name=SKINTEMP
    99         interp_option=sixteen_pt+four_pt+wt_average_4pt+wt_average_16pt+search
     99#        interp_option=sixteen_pt+four_pt+wt_average_4pt+wt_average_16pt+search
     100#        interp_option=four_pt
     101#        mandatory=yes    # MUST HAVE THIS FIELD
     102## commenter les trois prochaines lignes ?
     103#        masked=both
     104#        interp_land_mask  = LANDSEA(1)
     105#        interp_water_mask = LANDSEA(0)
     106#        fill_missing=0.
    100107        mandatory=yes    # MUST HAVE THIS FIELD
    101         masked=both
    102         interp_land_mask  = LANDSEA(1)
    103         interp_water_mask = LANDSEA(0)
    104         fill_missing=0.
     108        interp_option=four_pt+average_4pt
     109        #careful: inputs from the GCMs are in a "polygonal" world !
     110        #        interp_option=sixteen_pt+four_pt+average_4pt
     111        interp_option=sixteen_pt+four_pt
     112        fill_missing=-1.0e+30
     113        missing_value=-9.99e+33
    105114========================================
    106115name=PSFC
    107         interp_option=four_pt+average_4pt
     116#        interp_option=four_pt+average_4pt
     117#        interp_option=sixteen_pt+average_16pt
     118#careful: inputs from the GCMs are in a "polygonal" world !
     119        interp_option=sixteen_pt+four_pt+average_4pt
    108120        fill_lev=200100:const(200100.)
    109121        flag_in_output=FLAG_PSFC
    … …  
    122134========================================
    123135name=SW000010
    124         interp_option=sixteen_pt+four_pt+average_4pt+search
     136#        interp_option=sixteen_pt+four_pt+average_4pt+search
    125137        masked=water
    126138        interp_mask=LANDSEA(0)
    … …  
    129141========================================
    130142name=SW010040
    131         interp_option=sixteen_pt+four_pt+average_4pt+search
     143#        interp_option=sixteen_pt+four_pt+average_4pt+search
    132144        masked=water
    133145        interp_mask=LANDSEA(0)
    … …  
    136148========================================
    137149name=SW040100
    138         interp_option=sixteen_pt+four_pt+average_4pt+search
     150#        interp_option=sixteen_pt+four_pt+average_4pt+search
    139151        masked=water
    140152        interp_mask=LANDSEA(0)
    … …  
    143155========================================
    144156name=SW100200
    145         interp_option=sixteen_pt+four_pt+average_4pt+search
     157#        interp_option=sixteen_pt+four_pt+average_4pt+search
    146158        masked=water
    147159        interp_mask=LANDSEA(0)
    … …  
    150162========================================
    151163name=SW010200
    152         interp_option=sixteen_pt+four_pt+average_4pt+search
     164#        interp_option=sixteen_pt+four_pt+average_4pt+search
    153165        masked=water
    154166        interp_mask=LANDSEA(0)
    … …  
    157169========================================
    158170name=SM000010
    159         interp_option=sixteen_pt+four_pt+average_4pt+search
    160         masked=water
    161         interp_mask=LANDSEA(0)
     171#        interp_option=sixteen_pt+four_pt+average_4pt+search
     172        interp_option=four_pt
     173#        masked=water
     174#        interp_mask=LANDSEA(0)
    162175        fill_missing=1.
    163176        flag_in_output=FLAG_SM000010
    164177========================================
    165178name=SM010040
    166         interp_option=sixteen_pt+four_pt+average_4pt+search
    167         masked=water
    168         interp_mask=LANDSEA(0)
     179#        interp_option=sixteen_pt+four_pt+average_4pt+search
     180        interp_option=four_pt
     181#        masked=water
     182#        interp_mask=LANDSEA(0)
    169183        fill_missing=1.
    170184        flag_in_output=FLAG_SM010040
    171185========================================
    172186name=SM040100
    173         interp_option=sixteen_pt+four_pt+average_4pt+search
    174         masked=water
    175         interp_mask=LANDSEA(0)
     187#        interp_option=sixteen_pt+four_pt+average_4pt+search
     188        interp_option=four_pt
     189#        masked=water
     190#        interp_mask=LANDSEA(0)
    176191        fill_missing=1.
    177192        flag_in_output=FLAG_SM040100
    178193========================================
    179194name=SM100200
    180         interp_option=sixteen_pt+four_pt+average_4pt+search
    181         masked=water
    182         interp_mask=LANDSEA(0)
     195#        interp_option=sixteen_pt+four_pt+average_4pt+search
     196        interp_option=four_pt
     197#        masked=water
     198#        interp_mask=LANDSEA(0)
    183199        fill_missing=1.
    184200        flag_in_output=FLAG_SM100200
    185201========================================
    186202name=SM010200
    187         interp_option=sixteen_pt+four_pt+average_4pt+search
    188         masked=water
    189         interp_mask=LANDSEA(0)
     203#        interp_option=sixteen_pt+four_pt+average_4pt+search
     204        interp_option=four_pt
     205#        masked=water
     206#        interp_mask=LANDSEA(0)
    190207        fill_missing=1.
    191208        flag_in_output=FLAG_SM010200
    192209========================================
     210# EMISSIVITY
    193211name=ST000010
    194         interp_option=sixteen_pt+four_pt+average_4pt+search
    195         masked=water
    196         interp_mask=LANDSEA(0)
    197         fill_missing=285.
     212#        interp_option=sixteen_pt+four_pt+average_4pt+search
     213        interp_option=sixteen_pt+four_pt
     214#        masked=water
     215#        interp_mask=LANDSEA(0)
     216        fill_missing=0.
    198217        flag_in_output=FLAG_ST000010
    199218========================================
     219# CO2 ICE
    200220name=ST010040
    201         interp_option=sixteen_pt+four_pt+average_4pt+search
    202         masked=water
    203         interp_mask=LANDSEA(0)
    204         fill_missing=285.
     221#        interp_option=sixteen_pt+four_pt+average_4pt+search
     222#        interp_option=sixteen_pt+four_pt
     223        interp_option=four_pt
     224###ATTENTION sixteen_pt donne des valeurs negatives
     225#        masked=water
     226#        interp_mask=LANDSEA(0)
     227        fill_missing=0.
    205228        flag_in_output=FLAG_ST010040
    206229========================================
    207230name=ST040100
    208         interp_option=sixteen_pt+four_pt+average_4pt+search
    209         masked=water
    210         interp_mask=LANDSEA(0)
     231#        interp_option=sixteen_pt+four_pt+average_4pt+search
     232        interp_option=four_pt
     233#        masked=water
     234#        interp_mask=LANDSEA(0)
    211235        fill_missing=285.
    212236        flag_in_output=FLAG_ST040100
    213237========================================
    214238name=ST100200
    215         interp_option=sixteen_pt+four_pt+average_4pt+search
    216         masked=water
    217         interp_mask=LANDSEA(0)
     239#        interp_option=sixteen_pt+four_pt+average_4pt+search
     240        interp_option=four_pt
     241#        masked=water
     242#        interp_mask=LANDSEA(0)
    218243        fill_missing=285.
    219244        flag_in_output=FLAG_ST100200
    220245========================================
    221246name=ST010200
    222         interp_option=sixteen_pt+four_pt+average_4pt+search
    223         masked=water
    224         interp_mask=LANDSEA(0)
     247#        interp_option=sixteen_pt+four_pt+average_4pt+search
     248        interp_option=four_pt
     249#        masked=water
     250#        interp_mask=LANDSEA(0)
    225251        fill_missing=285.
    226252        flag_in_output=FLAG_ST010200
    227253========================================
    228254name=SM000007
    229         interp_option=sixteen_pt+four_pt+average_4pt+search
     255#        interp_option=sixteen_pt+four_pt+average_4pt+search
     256        interp_option=four_pt
    230257        masked=water
    231258        interp_mask=LANDSEA(0)
    … …  
    371398name=RH
    372399        interp_option=sixteen_pt+four_pt+average_4pt
    373         fill_missing=0.
    374         fill_lev=200100:const(0.)
     400        fill_missing=100.
     401        fill_lev=200100:const(0.)       
    375402========================================
    376403name=TAVGSFC
    … …  
    383410                         #   write the field out as TT
    384411========================================
     412#average_4pt, average_16pt, wt_average_4pt, wt_average_16pt,
     413#nearest_neighbor, four_pt, sixteen_pt, search, average_gcell(r)
     414#
     415#wt_average: bof
     416#average_gcell: bof aussi
    385417name=TT
    386418        mandatory=yes    # MUST HAVE THIS FIELD
    387         interp_option=four_pt+average_4pt
    388         fill_missing=0.
     419#        interp_option=four_pt+average_4pt
     420#careful: inputs from the GCMs are in a "polygonal" world !
     421        interp_option=sixteen_pt+four_pt+average_4pt
     422#        interp_option=sixteen_pt+four_pt
     423        fill_missing=-1.0e+30
     424        missing_value=-9.99e+33
    389425========================================
    390426name=U ; output_name=UU  # If we get U, use entry from UU and
    … …  
    394430        mandatory=yes    # MUST HAVE THIS FIELD
    395431        interp_option=sixteen_pt+four_pt+average_4pt
     432#        interp_option=sixteen_pt+four_pt
    396433        is_u_field=yes
    397434        output_stagger=U
    398         fill_missing=0.
     435        fill_missing=-1.0e+30
     436        missing_value=-9.99e+33
    399437========================================
    400438name=V ; output_name=VV  # If we get V, use entry from VV and
    … …  
    404442        mandatory=yes    # MUST HAVE THIS FIELD
    405443        interp_option=sixteen_pt+four_pt+average_4pt
     444#        interp_option=sixteen_pt+four_pt
    406445        is_v_field=yes
    407446        output_stagger=V
    408         fill_missing=0.
     447        fill_missing=-1.0e+30
     448        missing_value=-9.99e+33
    409449========================================
    410450name=SST
    … …  
    474514========================================
    475515name=SOILHGT
    476         interp_option=four_pt+average_4pt
     516#        interp_option=four_pt+average_4pt
     517#careful: inputs from the GCMs are in a "polygonal" world !
     518        interp_option=sixteen_pt+four_pt+average_4pt
    477519        flag_in_output=FLAG_SOILHGT
    478520========================================
     521name=HV
     522#        mandatory=yes    # MUST HAVE THIS FIELD
     523#        interp_option=sixteen_pt+four_pt+average_4pt
     524#        interp_option=sixteen_pt+four_pt
     525        interp_option=four_pt
     526#        fill_lev=200100:const(0.)
     527        fill_missing=-1.0e+30
     528        missing_value=-9.99e+33
     529========================================
     530name=HI
     531#        mandatory=yes    # MUST HAVE THIS FIELD
     532#        interp_option=sixteen_pt+four_pt+average_4pt
     533#        interp_option=sixteen_pt+four_pt
     534        interp_option=four_pt       
     535#        fill_lev=200100:const(0.)
     536        fill_missing=-1.0e+30
     537        missing_value=-9.99e+33
     538========================================
     539name=TSOIL
     540#        mandatory=yes    # MUST HAVE THIS FIELD
     541#         interp_option=sixteen_pt+four_pt+average_4pt
     542         interp_option=sixteen_pt+four_pt
     543#        fill_lev=200100:const(0.)
     544         fill_missing=-1.0e+30
     545         missing_value=-9.99e+33
     546========================================

  • trunk/mesoscale/LMD_MM_MARS/SRC/WPS/metgrid/src/module_date_pack.f90

    r11 r12  
    1111!  3. split_date_char ( date , century_year , month , day , hour , minute , second )
    1212!  Given the date, return the integer components.
     13
     14!****************
     15! subroutines are similar to the ones in module_date_time in the WRF share directory
     16!****************
    1317
    1418use module_debug
    … …  
    6872     
    6973      !  Assign the number of days in a months
    70      
    71       mday( 1) = 31
    72       mday( 2) = 28
    73       mday( 3) = 31
    74       mday( 4) = 30
    75       mday( 5) = 31
    76       mday( 6) = 30
    77       mday( 7) = 31
    78       mday( 8) = 31
    79       mday( 9) = 30
    80       mday(10) = 31
    81       mday(11) = 30
    82       mday(12) = 31
     74
     75! months as defined in Martian ESMF routines
     76 
     77      mday( 1) = 61
     78      mday( 2) = 66
     79      mday( 3) = 66
     80      mday( 4) = 65
     81      mday( 5) = 60
     82      mday( 6) = 54
     83      mday( 7) = 50
     84      mday( 8) = 46
     85      mday( 9) = 47
     86      mday(10) = 47
     87      mday(11) = 51
     88      mday(12) = 56
    8389     
    8490      !  Break down old hdate into parts
    … …  
    143149      !  Check that the day of NDATE makes sense.
    144150     
    145       IF (monew.ne.2) THEN
     151!      IF (monew.ne.2) THEN
    146152      ! ...... For all months but February
    147153         IF ((dynew.GT.mday(monew)).or.(dynew.LT.1)) THEN
    … …  
    149155            npass = .false.
    150156         END IF
    151       ELSE IF (monew.eq.2) THEN
    152       ! ...... For February
    153          IF ((dynew.GT.nfeb(yrnew)).OR.(dynew.LT.1)) THEN
    154             PRINT*, 'GETH_IDTS:  Day of NDATE = ', dynew
    155             npass = .false.
    156          END IF
    157       END IF
     157!      ELSE IF (monew.eq.2) THEN
     158!      ! ...... For February
     159!         IF ((dynew.GT.nfeb(yrnew)).OR.(dynew.LT.1)) THEN
     160!            PRINT*, 'GETH_IDTS:  Day of NDATE = ', dynew
     161!            npass = .false.
     162!         END IF
     163!      END IF
    158164     
    159165      !  Check that the day of ODATE makes sense.
    160166     
    161       IF (moold.ne.2) THEN
     167!      IF (moold.ne.2) THEN
    162168      ! ...... For all months but February
    163169         IF ((dyold.GT.mday(moold)).or.(dyold.LT.1)) THEN
    … …  
    165171            opass = .false.
    166172         END IF
    167       ELSE IF (moold.eq.2) THEN
    168       ! ....... For February
    169          IF ((dyold.GT.nfeb(yrold)).or.(dyold.LT.1)) THEN
    170             PRINT*, 'GETH_IDTS:  Day of ODATE = ', dyold
    171             opass = .false.
    172          END IF
    173       END IF
     173!      ELSE IF (moold.eq.2) THEN
     174!      ! ....... For February
     175!         IF ((dyold.GT.nfeb(yrold)).or.(dyold.LT.1)) THEN
     176!            PRINT*, 'GETH_IDTS:  Day of ODATE = ', dyold
     177!            opass = .false.
     178!         END IF
     179!      END IF
    174180     
    175181      !  Check that the hour of NDATE makes sense.
    … …  
    189195      !  Check that the minute of NDATE makes sense.
    190196     
    191       IF ((minew.GT.59).or.(minew.LT.0)) THEN
     197      IF ((minew.GT.36).or.(minew.LT.0)) THEN
    192198         PRINT*, 'GETH_IDTS:  Minute of NDATE = ', minew
    193199         npass = .false.
    … …  
    196202      !  Check that the minute of ODATE makes sense.
    197203     
    198       IF ((miold.GT.59).or.(miold.LT.0)) THEN
     204      IF ((miold.GT.36).or.(miold.LT.0)) THEN
    199205         PRINT*, 'GETH_IDTS:  Minute of ODATE = ', miold
    200206         opass = .false.
    … …  
    203209      !  Check that the second of NDATE makes sense.
    204210     
    205       IF ((scnew.GT.59).or.(scnew.LT.0)) THEN
     211      IF ((scnew.GT.99).or.(scnew.LT.0)) THEN
    206212         PRINT*, 'GETH_IDTS:  SECOND of NDATE = ', scnew
    207213         npass = .false.
    … …  
    210216      !  Check that the second of ODATE makes sense.
    211217     
    212       IF ((scold.GT.59).or.(scold.LT.0)) THEN
     218      IF ((scold.GT.99).or.(scold.LT.0)) THEN
    213219         PRINT*, 'GETH_IDTS:  Second of ODATE = ', scold
    214220         opass = .false.
    … …  
    216222     
    217223      IF (.not. npass) THEN
     224         !WRITE( wrf_err_message , * ) 'module_date_time: geth_idts: Bad NDATE: ', ndate(1:nlen)
     225         !CALL wrf_error_fatal ( TRIM ( wrf_err_message ) )
     226         !****difference with module_date_time
    218227         call mprintf(.true.,ERROR,'Screwy NDATE: %s',s1=ndate(1:nlen))
    219228      END IF
    220229     
    221230      IF (.not. opass) THEN
     231         !WRITE( wrf_err_message , * ) 'module_date_time: geth_idts: Bad ODATE: ', odate(1:olen)
     232         !CALL wrf_error_fatal ( TRIM ( wrf_err_message ) )
     233         !****difference with module_date_time
    222234         call mprintf(.true.,ERROR,'Screwy ODATE: %s',s1=odate(1:olen))
    223235      END IF
    … …  
    231243      newdys = 0
    232244      DO i = yrold, yrnew - 1
    233          newdys = newdys + (365 + (nfeb(i)-28))
     245!         newdys = newdys + (365 + (nfeb(i)-28))
     246! so cool to deal with a planet without any february month :)
     247         newdys = newdys + 669
    234248      END DO
    235249     
    236250      IF (monew .GT. 1) THEN
    237          mday(2) = nfeb(yrnew)
     251!         mday(2) = nfeb(yrnew)
    238252         DO i = 1, monew - 1
    239253            newdys = newdys + mday(i)
    240254         END DO
    241          mday(2) = 28
     255!         mday(2) = 28
    242256      END IF
    243257     
    … …  
    250264     
    251265      IF (moold .GT. 1) THEN
    252          mday(2) = nfeb(yrold)
     266!         mday(2) = nfeb(yrold)
    253267         DO i = 1, moold - 1
    254268            olddys = olddys + mday(i)
    255269         END DO
    256          mday(2) = 28
     270!         mday(2) = 28
    257271      END IF
    258272     
    … …  
    261275      !  Determine the time difference in seconds
    262276     
    263       idts = (newdys - olddys) * 86400
    264       idts = idts + (hrnew - hrold) * 3600
    265       idts = idts + (minew - miold) * 60
     277      idts = (newdys - olddys) * 88800    !martian sols
     278      idts = idts + (hrnew - hrold) * 3700  !hours
     279      idts = idts + (minew - miold) * 100  !minutes
    266280      idts = idts + (scnew - scold)
    267281     
    … …  
    328342      CHARACTER (LEN=10) :: hfrc
    329343      CHARACTER (LEN=1) :: sp
    330       ! INTEGER, EXTERNAL :: nfeb  ! in the same module now
    331 
     344!      INTEGER, EXTERNAL :: nfeb  ! in the same module now
     345     
    332346      !  Assign the number of days in a months
    333      
    334       mday( 1) = 31
    335       mday( 2) = 28
    336       mday( 3) = 31
    337       mday( 4) = 30
    338       mday( 5) = 31
    339       mday( 6) = 30
    340       mday( 7) = 31
    341       mday( 8) = 31
    342       mday( 9) = 30
    343       mday(10) = 31
    344       mday(11) = 30
    345       mday(12) = 31
     347
     348      mday( 1) = 61
     349      mday( 2) = 66
     350      mday( 3) = 66
     351      mday( 4) = 65
     352      mday( 5) = 60
     353      mday( 6) = 54
     354      mday( 7) = 50
     355      mday( 8) = 46
     356      mday( 9) = 47
     357      mday(10) = 47
     358      mday(11) = 51
     359      mday(12) = 56
    346360     
    347361      !  Break down old hdate into parts
    … …  
    377391      END IF
    378392     
    379       !  Set the number of days in February for that year.
    380      
    381       mday(2) = nfeb(yrold)
     393!      !  Set the number of days in February for that year.
     394     
     395!      mday(2) = nfeb(yrold)
    382396     
    383397      !  Check that ODATE makes sense.
    … …  
    408422      !  Check that the minute of ODATE makes sense.
    409423     
    410       IF ((miold.GT.59).or.(miold.LT.0)) THEN
     424      IF ((miold.GT.36).or.(miold.LT.0)) THEN
    411425         WRITE(*,*) 'GETH_NEWDATE:  Minute of ODATE = ', miold
    412426         opass = .FALSE.
    … …  
    415429      !  Check that the second of ODATE makes sense.
    416430     
    417       IF ((scold.GT.59).or.(scold.LT.0)) THEN
     431      IF ((scold.GT.99).or.(scold.LT.0)) THEN
    418432         WRITE(*,*) 'GETH_NEWDATE:  Second of ODATE = ', scold
    419433         opass = .FALSE.
    … …  
    424438     
    425439      IF (.not.opass) THEN
     440         !WRITE( wrf_err_message , * ) 'module_date_time: GETH_NEWDATE: Bad ODATE: ', odate(1:olen), olen
     441         !CALL wrf_error_fatal ( TRIM ( wrf_err_message ) )
     442         !**** difference avec module_date_time
    426443         call mprintf(.true.,ERROR,'GETH_NEWDATE: Crazy ODATE: %s %i',s1=odate(1:olen),i1=olen)
    427444      END IF
    … …  
    435452         ifrc = olen-20
    436453         ifrc = 10**ifrc
    437          nday   = ABS(idt)/(86400*ifrc)
    438          nhour  = MOD(ABS(idt),86400*ifrc)/(3600*ifrc)
    439          nmin   = MOD(ABS(idt),3600*ifrc)/(60*ifrc)
    440          nsec   = MOD(ABS(idt),60*ifrc)/(ifrc)
     454         nday   = ABS(idt)/(88800*ifrc)
     455         nhour  = MOD(ABS(idt),88800*ifrc)/(3700*ifrc)
     456         nmin   = MOD(ABS(idt),3700*ifrc)/(100*ifrc)
     457         nsec   = MOD(ABS(idt),100*ifrc)/(ifrc)
    441458         nfrac = MOD(ABS(idt), ifrc)
    442459      ELSE IF (olen.eq.19) THEN  !idt should be in seconds
    443460         ifrc = 1
    444          nday   = ABS(idt)/86400 ! Integer number of days in delta-time
    445          nhour  = MOD(ABS(idt),86400)/3600
    446          nmin   = MOD(ABS(idt),3600)/60
    447          nsec   = MOD(ABS(idt),60)
     461         nday   = ABS(idt)/88800 ! Integer number of days in delta-time
     462         nhour  = MOD(ABS(idt),88800)/3700
     463         nmin   = MOD(ABS(idt),3700)/100
     464         nsec   = MOD(ABS(idt),100)
    448465         nfrac  = 0
    449466      ELSE IF (olen.eq.16) THEN !idt should be in minutes
    450467         ifrc = 1
    451          nday   = ABS(idt)/1440 ! Integer number of days in delta-time
    452          nhour  = MOD(ABS(idt),1440)/60
    453          nmin   = MOD(ABS(idt),60)
     468         nday   = ABS(idt)/2400 ! Integer number of days in delta-time
     469         nhour  = MOD(ABS(idt),2400)/100
     470         nmin   = MOD(ABS(idt),100)
    454471         nsec   = 0
    455472         nfrac  = 0
    … …  
    469486         nfrac  = 0
    470487      ELSE
     488         !WRITE( wrf_err_message , * ) 'module_date_time: GETH_NEWDATE: Strange length for ODATE: ',olen
     489         !CALL wrf_error_fatal ( TRIM ( wrf_err_message ) )
     490         !**** difference with module_date_time
    471491         call mprintf(.true.,ERROR,'GETH_NEWDATE: Strange length for ODATE: %i',i1=olen)
    472492      END IF
    … …  
    481501     
    482502         scnew = scold + nsec
    483          IF (scnew .GE. 60) THEN
    484             scnew = scnew - 60
     503         IF (scnew .GE. 100) THEN
     504            scnew = scnew - 100
    485505            nmin  = nmin + 1
    486506         END IF
    487507     
    488508         minew = miold + nmin
    489          IF (minew .GE. 60) THEN
    490             minew = minew - 60
     509         IF (minew .GE. 37) THEN
     510            minew = minew - 37
    491511            nhour  = nhour + 1
    492512         END IF
    … …  
    509529                  monew = 1
    510530                  yrnew = yrnew + 1
    511                   ! If the year changes, recompute the number of days in February
    512                   mday(2) = nfeb(yrnew)
     531                  !! If the year changes, recompute the number of days in February
     532                  !mday(2) = nfeb(yrnew)
    513533               END IF
    514534            END IF
    … …  
    525545         scnew = scold - nsec
    526546         IF (scnew .LT. 00) THEN
    527             scnew = scnew + 60
     547            scnew = scnew + 100
    528548            nmin  = nmin + 1
    529549         END IF
    … …  
    531551         minew = miold - nmin
    532552         IF (minew .LT. 00) THEN
    533             minew = minew + 60
     553            minew = minew + 37
    534554            nhour  = nhour + 1
    535555         END IF
    … …  
    551571                  monew = 12
    552572                  yrnew = yrnew - 1
    553                   ! If the year changes, recompute the number of days in February
    554                   mday(2) = nfeb(yrnew)
     573                  !! If the year changes, recompute the number of days in February
     574                  !mday(2) = nfeb(yrnew)
    555575               END IF
    556576               dynew = mday(monew)
    … …  
    593613!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    594614
    595    FUNCTION nfeb ( year ) RESULT (num_days)
    596    
    597       ! Compute the number of days in February for the given year
    598    
    599       IMPLICIT NONE
    600    
    601       INTEGER :: year
    602       INTEGER :: num_days
    603    
    604       num_days = 28 ! By default, February has 28 days ...
    605       IF (MOD(year,4).eq.0) THEN 
    606          num_days = 29  ! But every four years, it has 29 days ...
    607          IF (MOD(year,100).eq.0) THEN
    608             num_days = 28  ! Except every 100 years, when it has 28 days ...
    609             IF (MOD(year,400).eq.0) THEN
    610                num_days = 29  ! Except every 400 years, when it has 29 days.
    611             END IF
    612          END IF
    613       END IF
    614    
    615    END FUNCTION nfeb
     615!   FUNCTION nfeb ( year ) RESULT (num_days)
     616!   
     617!      ! Compute the number of days in February for the given year
     618!   
     619!      IMPLICIT NONE
     620!   
     621!      INTEGER :: year
     622!      INTEGER :: num_days
     623!
     624!        num_days = 99999
     625!        PRINT *, 'WARNING !' 
     626!! normally never called ...       
     627! 
     628!!      num_days = 28 ! By default, February has 28 days ...
     629!!      IF (MOD(year,4).eq.0) THEN 
     630!!         num_days = 29  ! But every four years, it has 29 days ...
     631!!         IF (MOD(year,100).eq.0) THEN
     632!!            num_days = 28  ! Except every 100 years, when it has 28 days ...
     633!!            IF (MOD(year,400).eq.0) THEN
     634!!               num_days = 29  ! Except every 400 years, when it has 29 days.
     635!!            END IF
     636!!         END IF
     637!!      END IF
     638!   
     639!   END FUNCTION nfeb
    616640
    617641!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    623647   
    624648      CHARACTER(LEN=19) , INTENT(IN) :: date
    625   
     649 
    626650      !  Output data.
    627651   
    … …  
    636660   
    637661   END SUBROUTINE split_date_char
    638    
     662
    639663END MODULE date_pack

  • trunk/mesoscale/LMD_MM_MARS/SRC/WPS/metgrid/src/process_domain_module.F90

    r11 r12  
    106106   
    107107         call geth_newdate(valid_date, trim(start_date(n)), t*interval_seconds)
     108
     109         print *,'now working on ...'       
     110         print *,valid_date
     111
    108112         temp_date = ' '
    109113
    110          if (mod(interval_seconds,3600) == 0) then
    111             write(temp_date,'(a13)') valid_date(1:10)//'_'//valid_date(12:13)
    112          else if (mod(interval_seconds,60) == 0) then
    113             write(temp_date,'(a16)') valid_date(1:10)//'_'//valid_date(12:16)
    114          else
     114!*****MARS !!!
     115!         if (mod(interval_seconds,3700) == 0) then
     116!            write(temp_date,'(a13)') valid_date(1:10)//'_'//valid_date(12:13)
     117!         else if (mod(interval_seconds,100) == 0) then
     118!            write(temp_date,'(a16)') valid_date(1:10)//'_'//valid_date(12:16)
     119!         else
    115120            write(temp_date,'(a19)') valid_date(1:10)//'_'//valid_date(12:19)
    116          end if
     121!         end if
    117122   
    118123         call mprintf(.true.,STDOUT, ' Processing %s', s1=trim(temp_date))
    … …  
    726731                           met_dy = abs(deltalat)
    727732                        else
     733!
     734! MARS MARS: BEWARE WITH THAT PART !
     735!
    728736! BUG: Need to more correctly handle dx/dy in meters.
    729737                           met_dx = met_dx / 111000.  ! Convert meters to approximate degrees
Note: See TracChangeset for help on using the changeset viewer.