Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/ALT_POLAR_top0.1_levels =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/ALT_POLAR_top0.1_levels (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/ALT_POLAR_top0.1_levels (revision 212) @@ -0,0 +1,101 @@ + 1.00000 + 0.999500 + 0.998000 + 0.995000 + 0.985000 + 0.970000 + 0.948382 + 0.899410 + 0.852949 + 0.808870 + 0.767051 + 0.727376 + 0.689735 + 0.654024 + 0.620144 + 0.588001 + 0.557507 + 0.528575 + 0.501127 + 0.475086 + 0.450380 + 0.426941 + 0.404704 + 0.383607 + 0.363591 + 0.344602 + 0.326586 + 0.309494 + 0.293279 + 0.277894 + 0.263299 + 0.249451 + 0.236314 + 0.223850 + 0.212026 + 0.200807 + 0.190164 + 0.180066 + 0.170486 + 0.161398 + 0.152775 + 0.144594 + 0.136833 + 0.129470 + 0.122484 + 0.115856 + 0.109568 + 0.103603 + 0.0979435 + 0.0925741 + 0.0874799 + 0.0826470 + 0.0780618 + 0.0737118 + 0.0695847 + 0.0656693 + 0.0619546 + 0.0584303 + 0.0550868 + 0.0519146 + 0.0489052 + 0.0460500 + 0.0433411 + 0.0407712 + 0.0383331 + 0.0360199 + 0.0338253 + 0.0317433 + 0.0297680 + 0.0278940 + 0.0261160 + 0.0244293 + 0.0228290 + 0.0213107 + 0.0198703 + 0.0185037 + 0.0172072 + 0.0159772 + 0.0148102 + 0.0137031 + 0.0126527 + 0.0116562 + 0.0107108 + 0.00981387 + 0.00896290 + 0.00815556 + 0.00738962 + 0.00666295 + 0.00597353 + 0.00531947 + 0.00469894 + 0.00411022 + 0.00355168 + 0.00302178 + 0.00251905 + 0.00204210 + 0.00158960 + 0.00116030 + 0.000753006 + 0.000366596 + 0.00000 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_GEOGRID.TBL =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_GEOGRID.TBL (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_GEOGRID.TBL (revision 212) @@ -0,0 +1,240 @@ +# See options.txt for a (somewhat up to date) list of the +# options that may be specified here. +=============================== +name = HGT_M + priority = 1 + dest_type = continuous + df_dx=SLPX + df_dy=SLPY + fill_missing=9000. +# smooth_option = 1-2-1; smooth_passes=5 #10000 +# smooth_option = 1-2-1; smooth_passes=2 +# smooth_option = smth-desmth_special; smooth_passes=5 + smooth_option = smth-desmth_special; smooth_passes=2 + interp_option = default:average_gcell(4.0)+four_pt+average_4pt +# smooth_option = 1-2-1; smooth_passes=1 +# smooth_option = smth-desmth; smooth_passes=5 +## smooth_option = 1-2-1; smooth_passes=5 +#### tous a four_pt normalement +# interp_option = default:sixteen_pt+four_pt +# #interp_option = 32ppd:sixteen_pt+four_pt +# interp_option = 32ppd_HRalb:sixteen_pt+four_pt +# interp_option = gcm:four_pt +# interp_option = 64ppd:sixteen_pt+four_pt +# interp_option = 64ppd_noHRti:sixteen_pt+four_pt +# interp_option = 128ppd:average_gcell(4.0)+four_pt+average_4pt +# interp_option = custom:four_pt +# interp_option = omega:four_pt +# interp_option = 32ppd:average_gcell(4.0)+four_pt+average_4pt +# #interp_option = 32ppd:average_gcell(16.0)+sixteen_pt+average_16pt + rel_path= default:mola_topo32/ + rel_path= 32ppd:mola_topo32/ + rel_path= 32ppd_HRalb:mola_topo32/ + rel_path= gcm:mola_GCM/ + rel_path= 64ppd:mola_topo64/ + rel_path= 64ppd_noHRti:mola_topo64/ + rel_path= custom:crater2/ + rel_path= omega:mola_topo64/ +=============================== +name = HGT_U + output_stagger = U + priority = 1 + dest_type = continuous + fill_missing=9000. +# smooth_option = 1-2-1; smooth_passes=5 +# smooth_option = 1-2-1; smooth_passes=2 +# smooth_option = smth-desmth_special; smooth_passes=5 + smooth_option = smth-desmth_special; smooth_passes=2 + interp_option = default:average_gcell(4.0)+four_pt+average_4pt +##### tous a four_pt normalement +# interp_option = default:sixteen_pt+four_pt +# interp_option = 32ppd:sixteen_pt+four_pt +# interp_option = 32ppd_HRalb:sixteen_pt+four_pt +# interp_option = 64ppd:sixteen_pt+four_pt +# interp_option = 64ppd_noHRti:sixteen_pt+four_pt +# interp_option = 128ppd:average_gcell(4.0)+four_pt+average_4pt +# interp_option = custom:four_pt +# interp_option = omega:four_pt +# interp_option = gcm:four_pt + rel_path= default:mola_topo32/ + rel_path= 32ppd:mola_topo32/ + rel_path= 32ppd_HRalb:mola_topo32/ + rel_path= gcm:mola_GCM/ + rel_path= 64ppd:mola_topo64/ + rel_path= 64ppd_noHRti:mola_topo64/ + rel_path= custom:crater2/ + rel_path= omega:mola_topo64/ +=============================== +name = HGT_V + output_stagger = V + priority = 1 + dest_type = continuous + fill_missing=9000. +# smooth_option = 1-2-1; smooth_passes=5 +# smooth_option = 1-2-1; smooth_passes=2 +# smooth_option = smth-desmth_special; smooth_passes=5 + smooth_option = smth-desmth_special; smooth_passes=2 + interp_option = default:average_gcell(4.0)+four_pt+average_4pt +##### tous a four_pt normalement +# interp_option = default:sixteen_pt+four_pt +# interp_option = 32ppd:sixteen_pt+four_pt +# interp_option = 32ppd_HRalb:sixteen_pt+four_pt +# interp_option = 64ppd:sixteen_pt+four_pt +# interp_option = 64ppd_noHRti:sixteen_pt+four_pt +# interp_option = 128ppd:average_gcell(4.0)+four_pt+average_4pt +# interp_option = custom:four_pt +# interp_option = omega:four_pt +# interp_option = gcm:four_pt + rel_path= default:mola_topo32/ + rel_path= 32ppd:mola_topo32/ + rel_path= 32ppd_HRalb:mola_topo32/ + rel_path= gcm:mola_GCM/ + rel_path= 64ppd:mola_topo64/ + rel_path= 64ppd_noHRti:mola_topo64/ + rel_path= custom:crater2/ + rel_path= omega:mola_topo64/ +=============================== +name = THERMAL_INERTIA + priority = 1 + dest_type = continuous + smooth_option = 1-2-1; smooth_passes=5 +# smooth_option = smth-desmth_special; smooth_passes=5 +# smooth_option = smth-desmth_special; smooth_passes=1 + interp_option = 32ppd:sixteen_pt+four_pt+average_4pt+average_16pt+search + rel_path= 32ppd:thermal_GCM/ + interp_option = 32ppd_HRalb:sixteen_pt+four_pt+average_4pt+average_16pt+search + rel_path= 32ppd_HRalb:thermal_GCM/ + interp_option = 64ppd:sixteen_pt+four_pt + rel_path= 64ppd:thermal_TES/ + interp_option = 64ppd_noHRti:sixteen_pt+four_pt+average_4pt+average_16pt+search + rel_path= 64ppd_noHRti:thermal_GCM/ + interp_option = custom:sixteen_pt+four_pt +# rel_path= custom:thermal_TES/ + rel_path= custom:thermal_GCM/ + interp_option = omega:sixteen_pt+four_pt + rel_path= omega:thermal_GCM/ + interp_option = gcm:four_pt + rel_path= gcm:thermal_GCM/ + fill_missing=400. #north polar TI (needed with thermal_TES) + fill_missing=4000. #north polar TI (needed with thermal_TES 64ppd) + fill_missing=4100. #north polar TI (needed with thermal_TES 64ppd) +# fill_missing=180. #polar projection artefacts +=============================== +name = ALBEDO_GCM + priority = 1 + dest_type = continuous + smooth_option = 1-2-1; smooth_passes=5 +# smooth_option = smth-desmth_special; smooth_passes=5 #3 +# fill_missing=0.25 +# interp_option = default:sixteen_pt+four_pt+average_4pt+average_16pt+search + interp_option = 32ppd:four_pt+average_4pt+average_16pt+search + rel_path = 32ppd:albedo_GCM/ + interp_option = 32ppd_HRalb:sixteen_pt+four_pt + rel_path= 32ppd_HRalb:albedo_TES/ + interp_option = 64ppd:sixteen_pt+four_pt + rel_path= 64ppd:albedo_TES/ + interp_option = 64ppd_noHRti:sixteen_pt+four_pt + rel_path= 64ppd_noHRti:albedo_TES/ + fill_missing=0.35 #north polar alb (needed with albedo_TES - albedo_GCM/1.3 to be consistent) + fill_missing=0.33 #north polar alb (needed with albedo_TES 64 ppd) + fill_missing=0.325 #north polar alb (needed with albedo_TES 64 ppd) + interp_option = custom:sixteen_pt+four_pt + rel_path= custom:albedo_TES/ + interp_option = omega:sixteen_pt+four_pt + rel_path= omega:albedo_OMG/ + interp_option = gcm:four_pt + rel_path = gcm:albedo_GCM/ +#=============================== +#name=LANDUSEF +# priority=1 +# dest_type=categorical +# z_dim_name=land_cat +# landmask_water=16 # Calculate a landmask from this field +# dominant=LU_INDEX +# interp_option = 30s:nearest_neighbor +# interp_option = 2m:four_pt +# interp_option = 5m:four_pt +# interp_option = 10m:four_pt +# interp_option = default:four_pt +# rel_path= 30s:landuse_30s/ +# rel_path= 2m:landuse_2m/ +# rel_path= 5m:landuse_5m/ +# rel_path= 10m:landuse_10m/ +# rel_path= default:landuse_2m/ +#=============================== +#name=SOILTEMP +# priority=1 +# dest_type=continuous +# interp_option=default:sixteen_pt+four_pt+average_4pt+average_16pt+search +# masked=water +# fill_missing=0. +# rel_path=default:soiltemp_1deg/ +#=============================== +#name=SOILCTOP +# dominant=SOILCAT +# priority=1 +# dest_type=categorical +# z_dim_name=soil_cat +# interp_option = 30s:nearest_neighbor +# interp_option = 2m:four_pt +# interp_option = 5m:four_pt +# interp_option = 10m:four_pt +# interp_option = default:four_pt +# rel_path= 30s:soiltype_top_30s/ +# rel_path= 2m:soiltype_top_2m/ +# rel_path= 5m:soiltype_top_5m/ +# rel_path= 10m:soiltype_top_10m/ +# rel_path= default:soiltype_top_2m/ +#=============================== +#name=SOILCBOT +# priority=1 +# dest_type=categorical +# z_dim_name=soil_cat +# interp_option = 30s:nearest_neighbor +# interp_option = 2m:four_pt +# interp_option = 5m:four_pt +# interp_option = 10m:four_pt +# interp_option = default:four_pt +# rel_path= 30s:soiltype_bot_30s/ +# rel_path= 2m:soiltype_bot_2m/ +# rel_path= 5m:soiltype_bot_5m/ +# rel_path= 10m:soiltype_bot_10m/ +# rel_path= default:soiltype_bot_2m/ +#=============================== +#name=ALBEDO12M +# priority=1 +# dest_type=continuous +# z_dim_name=month +# masked = water +# fill_missing = 8. +# interp_option=default:four_pt+average_4pt+average_16pt+search +# rel_path=default:albedo_ncep/ +#=============================== +#name=GREENFRAC +# priority=1 +# dest_type=continuous +# interp_option=default:four_pt+average_4pt+average_16pt+search +# z_dim_name=month +# masked = water +# fill_missing = 0. +# rel_path=default:greenfrac/ +#=============================== +#name=SNOALB +# priority=1 +# dest_type=continuous +# interp_option=default:four_pt+average_4pt+average_16pt+search +# masked = water +# fill_missing = 0. +# rel_path=default:maxsnowalb/ +#=============================== +#name=SLOPECAT +# priority=1 +# dominant_only=SLOPECAT +# dest_type=categorical +# z_dim_name=slope_cat +# masked = water +# fill_missing = 0. +# interp_option=default:nearest_neighbor+average_16pt+search +# rel_path=default:islope/ +#=============================== + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_levels =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_levels (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_levels (revision 212) @@ -0,0 +1,61 @@ + 1.00000 + 0.999500 + 0.998000 + 0.995000 + 0.985000 + 0.970000 + 0.928546 + 0.862104 + 0.800322 + 0.742874 + 0.689456 + 0.639785 + 0.593598 + 0.550651 + 0.510717 + 0.473583 + 0.439055 + 0.406948 + 0.377094 + 0.349334 + 0.323521 + 0.299519 + 0.277200 + 0.256447 + 0.237149 + 0.219206 + 0.202521 + 0.187006 + 0.172580 + 0.159165 + 0.146692 + 0.135093 + 0.124308 + 0.114280 + 0.104955 + 0.0962840 + 0.0882214 + 0.0807244 + 0.0737531 + 0.0672710 + 0.0612435 + 0.0556388 + 0.0504272 + 0.0455812 + 0.0410752 + 0.0368852 + 0.0329891 + 0.0293664 + 0.0259977 + 0.0228654 + 0.0199527 + 0.0172444 + 0.0147260 + 0.0123843 + 0.0102069 + 0.00818220 + 0.00629953 + 0.00454891 + 0.00292109 + 0.00140746 + 0.00000 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/callphys.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +10 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/makemeso_input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/makemeso_input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/makemeso_input (revision 212) @@ -0,0 +1,7 @@ +3 +4 +201 +201 +101 +1 +1 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.input (revision 212) @@ -0,0 +1,60 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 02, !! (p1) Start Martian Month + start_day = 48, !! (p1) Start Martian Day + start_hour = 06, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 02, !! (p1) End Martian Month + end_day = 53, !! (p1) End Martian Day + end_hour = 06, !! (p1) End Martian Hour (at longitude 0) + history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) + frames_per_outfile = 24, !! Size of time dimension in files + restart = .false. !! (*) Output restart files ? + restart_interval = 8880 !! (*) Frequency of output restart files ? + io_form_history = 2 !! (*) Choice of NETCDF for ouputs + io_form_restart = 2 !! (*) Choice of NETCDF for ouputs + io_form_input = 2 !! (*) Choice of NETCDF for ouputs + io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs + debug_level = 0 !! (*) Verbose level + !! + !! OPTIONAL + !! + interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) + / + + + &domains + time_step = 37 !! Dynamical timestep + dx = 10000, !! (p2) Horizontal resolution + dy = 10000, !! (p2) Horizontal resolution (should be equal to dx) + e_we = 201, !! (r)(p2) Number of longitude grid points + e_sn = 201, !! (r)(p2) Number of latitude grid points + e_vert = 101, !! (r)(p2) Number of vertical levels + p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model + !! + !! OPTIONAL + !! + force_sfc_in_vinterp = 6 !! (p3) Number of levels hardwired in the PBL + !! NB: decrease this parameter when low model top + / + + &physics + / + + &dynamics + tiso = 100. + / + + &bdy_control + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI_vincent/LS53/namelist.wps (revision 212) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 1 !! (r)(n)(p2) Total number of domains + start_date = '2024-02-48_06:00:00','2024-02-48_00:00:00','2024-02-48_00:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '2024-02-53_06:00:00','2024-02-53_00:00:00','2024-02-53_00:00:00' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs +/ + + +&geogrid + parent_id = 0, !! (n)(p2) Associated parent domain + parent_grid_ratio = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest + i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner + j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner + e_we = 201, !! (r)(p2) Number of longitude grid points + e_sn = 201, !! (r)(p2) Number of latitude grid points + dx = 10000 + dy = 10000 + geog_data_res = '64ppd', !! choice of static data sources + !! NB: possible: '64ppd', '32ppd', ... + !! NB: please glance at geogrid/GEOGRID.TBL + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = -1.88, !! north latitude of the center of the domain + ref_lon = 353.87, !! east longitude of the center of the domain + truelat1 = -1.88, !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = 353.87, !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] +/ + +&metgrid + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0 (revision 212) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.false. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0.ideal =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0.ideal (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars0.ideal (revision 212) @@ -0,0 +1,119 @@ +##General options +##~~~~~~~~~~~~~~~ +#Directory where external input files are +datadir = /san/work/colaitis/datafile/ + +#Run with or without tracer transport ? +tracer=.false. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 1 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis = 0.5 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .false. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .false. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .false. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars1 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars1 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars1 (revision 212) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .true. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11 (revision 212) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 2 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +# DUST: lifted by GCM surface winds ? +lifting = .true. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .true. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. #.true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11_radact =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11_radact (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars11_radact (revision 212) @@ -0,0 +1,119 @@ +##General options +##~~~~~~~~~~~~~~~ +#Directory where external input files are +datadir = /san/home/spiga/DATAFILE/ + +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 2 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. +# DUST: lifted by GCM surface winds ? +lifting = .true. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .true. +# WATER: Compute water cycle +water = .true. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars2 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars2 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars2 (revision 212) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 1 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .true. #.false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. #.true. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. #.true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars3 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars3 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/callphys.def.mars3 (revision 212) @@ -0,0 +1,119 @@ +##General options +##~~~~~~~~~~~~~~~ +#Directory where external input files are +datadir = /san/home/spiga/DATAFILE/ + +#Run with or without tracer transport ? +tracer=.true. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +callstats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 26 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .false. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. +# Impose polar cap surface albedos as observed by TES? +TESicealbedo = .true. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 6 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 2 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .true. #.false. +# DUST: use mass and number mixing ratios to predict dust size ? +doubleq = .true. #.false. +# DUST: lifted by GCM surface winds ? +lifting = .true. #.false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .true. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. #.true. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. #.true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input (revision 212) @@ -0,0 +1,66 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 05, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 01, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 05, !! (p1) End Martian Month + end_day = 09, !! (p1) End Martian Day + end_hour = 01, !! (p1) End Martian Hour (at longitude 0) + history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) + frames_per_outfile = 24, !! Size of time dimension in files + restart = .false. !! (*) Output restart files ? + restart_interval = 8880 !! (*) Frequency of output restart files ? + io_form_history = 2 !! (*) Choice of NETCDF for ouputs + io_form_restart = 2 !! (*) Choice of NETCDF for ouputs + io_form_input = 2 !! (*) Choice of NETCDF for ouputs + io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs + debug_level = 0 !! (*) Verbose level + / + + &domains + time_step = 20 !! Dynamical timestep + dx = 30000, !! (p2) Horizontal resolution + dy = 30000, !! (p2) Horizontal resolution (should be equal to dx) + e_we = 121, !! (r)(p2) Number of longitude grid points + e_sn = 121, !! (r)(p2) Number of latitude grid points + e_vert = 61, !! (r)(p2) Number of vertical levels + p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model + !! + !! OPTIONAL + !! + num_metgrid_levels = 26 !! (p1) number of vertical levels in GCM inputs (+1) + / + + &physics + mars = 11 + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + !! + !! OPTIONAL + !! + / + + &dynamics + tiso = 100. + pd_scalar = T !!! tres IMPORTANT, sinon valeurs negatives + !!! + !!! OPTIONAL + !!! + !non_hydrostatic = F, !! Integrate in non-hydrostatic/hydrostatic mode + ! !! --- necessaire pour garder un pas de tps long + ! !! --- et parce qu'un bout d'Arsia est au bord + / + + &bdy_control + / + + &grib2 + / + + &fdda + / + + &namelist_quilt !! (*) + nio_tasks_per_group = 0, !! (*) + nio_groups = 1, !! (*) + / !! (*) Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.wps (revision 212) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 1 !! (r)(n)(p2) Total number of domains + start_date = '2024-05-01_01:00:00','2024-05-01_01:00:00','2024-05-01_01:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '2024-05-09_01:00:00','2024-05-01_01:00:00','2024-05-01_01:00:00' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs +/ + + +&geogrid + parent_id = 0, !! (n)(p2) Associated parent domain + parent_grid_ratio = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest + i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner + j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner + e_we = 121, !! (r)(p2) Number of longitude grid points + e_sn = 121, !! (r)(p2) Number of latitude grid points + dx = 30000 + dy = 30000 + geog_data_res = '64ppd_noHRti', !! choice of static data sources + !! NB: possible: '64ppd', '32ppd', ... + !! NB: please glance at geogrid/GEOGRID.TBL + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = 10., !! north latitude of the center of the domain + ref_lon = -125., !! east longitude of the center of the domain + truelat1 = 10., !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = -125., !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] +/ + +&metgrid + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/run.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/run.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/run.def (revision 212) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.input (revision 212) @@ -0,0 +1,54 @@ + &time_control + start_year = 2024, + start_month = 09, + start_day = 23, + start_hour = 00, + end_year = 2024, + end_month = 09, + end_day = 26, + end_hour = 00, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 10, !5 + dx = 2500, + dy = 2500, + e_we = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 31, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 50 + / + + &physics + radt = 2, !5 ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUb_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00' + end_date = '2024-09-26_00:00:00', '2026-01-47_00:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 12 + j_parent_start = 1, 12 + e_we = 150 !275, 25 + e_sn = 150 !275, 25 + geog_data_res = '64ppd_noHRti' ! 'custom', 'custom' !'32ppd','32ppd' + dx = 2500, !1500, + dy = 2500, !1500, + map_proj = 'mercator', + ref_lat = 11 !14.76, + ref_lon = 159.34, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './LMD/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='/tmp15/aslmd/wrf/MET_EM' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.input (revision 212) @@ -0,0 +1,57 @@ + &time_control + start_year = 2024, + start_month = 09, + start_day = 23, + start_hour = 00, + end_year = 2024, + end_month = 09, + end_day = 26, + end_hour = 00, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains +!!! time_step = 1, ! 2 crashe, meme en baissant le top, 1 aussi au bout d'un moment + time_step = 0 + time_step_fract_num = 1, + time_step_fract_den = 2, + dx = 500, !1000, + dy = 500, !1000, + e_we = 181, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 181, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 31, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 50 + / + + &physics + radt = 20, !5 ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_finer_namelist.wps (revision 212) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00' + end_date = '2024-09-26_00:00:00', '2026-01-47_00:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 12 + j_parent_start = 1, 12 + e_we = 181 !150 !150 !275, 25 + e_sn = 181 !150 !150 !275, 25 + geog_data_res = '64ppd_noHRti' ! 'custom', 'custom' !'32ppd','32ppd' + dx = 500, !1500, + dy = 500, !1500, + map_proj = 'mercator', +! ref_lat = 13.4 !13.04, +! ref_lon = 158.6 !158.8, + ref_lat = 13.04, + ref_lon = 158.95, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.input (revision 212) @@ -0,0 +1,54 @@ + &time_control + start_year = 2024, + start_month = 09, + start_day = 23, + start_hour = 00, + end_year = 2024, + end_month = 09, + end_day = 26, + end_hour = 00, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 5 + dx = 1000, + dy = 1000, + e_we = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 31, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 50 + / + + &physics + radt = 4, !5 ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/DUs_namelist.wps (revision 212) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-09-23_00:00:00', '2026-01-47_00:00:00' + end_date = '2024-09-26_00:00:00', '2026-01-47_00:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 12 + j_parent_start = 1, 12 + e_we = 150 !150 !150 !275, 25 + e_sn = 150 !150 !150 !275, 25 + geog_data_res = '64ppd_noHRti' ! 'custom', 'custom' !'32ppd','32ppd' + dx = 1000, !1500, + dy = 1000, !1500, + map_proj = 'mercator', + ref_lat = 13.4 !13.04, + ref_lon = 158.6 !158.8, +! ref_lat = -4.87, +! ref_lon = 170.14, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphys.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +20 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd2_.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd2_.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd2_.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +40 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd3_.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd3_.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_callphysd3_.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +80 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_launch_simu =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_launch_simu (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_launch_simu (revision 212) @@ -0,0 +1,141 @@ +#! /bin/bash + + +# +# USE: launch_simu & +# + +# +# must set env. var. WRF_current +# + +# PRELIMINARY STEP +# - model must be compiled and linked (wrf.exe, real.exe) +# - domain must be correctly set with geogrid.exe +# - put namelist.input, namelist.wps, callphys*.def + +# 'ls' must yield: +# callphys*.def launch_simu namelist.input namelist.wps real.exe wrf.exe +# + + +############################################################### +gcmstep=0 #1 +wpsstep=0 #1 +realstep=1 +wrfstep=1 +############################################################### +# GCM STEP + daygcm=371 +# WPS STEP + name='em2' # un nom de dossier non pris dans WRFFEED +# WRF STEP + proc=4 +############################################################### +# FOLDERS + wheregcm=$GCM + whereprep=$WRF_current/PREP_MARS + wherewps=$WRF_current/WPS +############################################################### + + +# +# LOG FILE HEADER +# + +where=$PWD +echo "check what's happening in " ${where}/log_launch + +rm -rf ${where}/log_launch 2> /dev/null +echo SIMU folder is ${where} > ${where}/log_launch +echo SIMU WILL BE LAUNCHED on the machine >> ${where}/log_launch +hostname >> ${where}/log_launch + +# +# 0. GCM step +# +if [ ${gcmstep} -eq 1 ] +then + + echo "0. GCM run and output preparation for WPS" >> ${where}/log_launch + + cd ${wheregcm} + ln -sf ../STARTBASE_64_48_25_t2/start/start${daygcm}.nc start.nc + ln -sf ../STARTBASE_64_48_25_t2/startfi/startfi${daygcm}.nc startfi.nc + + rm -rf restart.nc + rm -rf restartfi.nc + rm -rf log_gcm + + echo running the GCM - check in ${wheregcm}/log_gcm >> ${where}/log_launch + nohup gcm.e > log_gcm + + echo preparing outputs for WPS >> ${where}/log_launch + cd ${whereprep} + compile_and_exec >> ${where}/log_launch + +fi + +# +# 1. WPS step +# +if [ ${wpsstep} -eq 1 ] +then + + echo "1. WPS: horizontal interpolations of static and meteorological data" >> ${where}/log_launch + +# rm -rf ${where}/met_em + mkdir ${where}/met_em + cd ${where}/met_em + ln -sf ../namelist.wps . + + cd ${wherewps} + rm -rf WHERE_SIMU_IS + ln -sf ${where}/met_em WHERE_SIMU_IS + mkdir ./WRFFEED/current + nohup geogrid.exe >> ${where}/log_launch + nohup metgrid.exe >> ${where}/log_launch + mv ./WRFFEED/current ./WRFFEED/${name} + + ln -sf $FEED/WRFFEED/${name}/* ${where}/met_em/ + +fi + +# +# 2. REAL.EXE +# +if [ ${realstep} -eq 1 ] +then + + echo "2. REAL.EXE: vertical interpolation and preparation of boundary conditions" >> ${where}/log_launch + + cd ${where}/met_em + rm -rf rsl.* 2> /dev/null + rm -rf wrfbdy* 2> /dev/null + rm -rf wrfinput* 2> /dev/null + + ln -sf ../namelist.input . + ln -sf ../real.exe . + real.exe + cd .. + ln -sf ./met_em/wrf*d0* . + +fi + + +# +# 3. WRF.EXE +# +if [ ${wrfstep} -eq 1 ] +then + + echo "3. WRF.EXE: launch the mesoscale simulation" >> ${where}/log_launch + + rm -rf rsl.* + rm -rf wrfout* + + $WHERE_MPI/mpirun -np ${proc} wrf.exe < /dev/null & + echo OK, simulation launched ! + +fi + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.input (revision 212) @@ -0,0 +1,138 @@ + &time_control + start_year = 2024, 2024, 2024, + start_month = 07, 07, 07, + start_day = 01, 01, 01, + start_hour = 01, 01, 01, + end_year = 2024, 2024, 2024, + end_month = 07, 07, 07, + end_day = 04, 04, 04, + end_hour = 01, 01, 01, + history_interval = 37, 37, 37, + frames_per_outfile = 24, 24, 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, T, T, ! F en second pour que d1 init. d2 (smoother) + / + + &domains + time_step = 30 + dx = 30000, 10000, 3333.333, + dy = 30000, 10000, 3333.333, + e_we = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, 61, 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 +!!!!! NESTING STUFF + max_dom = 3 !!!! PLEASE DO NOT FORGET callphysd*_.def + grid_id = 1, 2, 3, + parent_id = 0, 1, 2, + i_parent_start = 0, 51, 60, ! attention au raccord avec WPS ! + j_parent_start = 0, 55, 45, ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 3, 3, + parent_time_step_ratio = 1, 3, 3, + feedback = 0 + smooth_option = 2 + blend_width = 5 +!!!!! NESTING STUFF +! p_top_requested = 0.5 +! eta_levels = 1.00000, +! 0.999500, +! 0.998000, +! 0.995000, +! 0.985000, +! 0.970000, +! 0.940000, +! 0.900000, +! 0.800000, +! 0.700000, +! 0.592538, +! 0.501572, +! 0.424571, +! 0.359392, +! 0.304219, +! 0.257516, +! 0.217982, +! 0.184518, +! 0.156191, +! 0.132213, +! 0.111916, +! 0.0947347, +! 0.0801912, +! 0.0678804, +! 0.0574595, +! 0.0486384, +! 0.0411715, +! 0.0348509, +! 0.0295007, +! 0.0249718, +! 0.0211382, +! 0.0178931, +! 0.0151462, +! 0.0128209, +! 0.0108527, +! 0.00918661, +! 0.00777630, +! 0.00658249, +! 0.00557196, +! 0.00471656, +! 0.00399248, +! 0.00337956, +! 0.00286074, +! 0.00242156, +! 0.00204981, +! 0.00173513, +! 0.00146875, +! 0.00124327, +! 0.00105241, +! 0.000890844, +! 0.000754083, +! 0.000638317, +! 0.000540324, +! 0.000457374, +! 0.000387159, +! 0.000327723, +! 0.000277412, +! 0.000234824, +! 0.000198774, +! 0.000168259, +! 0.00000 + / + + &physics + radt = 1, 1, 1 ! 1 is the best (change iradia in callphys) + mars = 0, 0, 0 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, 2, 2, !! NECESSAIRE + diff_6th_factor = 0.2, 0.2, 0.2 !! NECESSAIRE + diff_opt = 1, + km_opt = 4, + pd_scalar = F !! no water cycle (mars = 0) + time_step_sound = 10, 10, 10 + !!damp_opt = 1 + !!dampcoef = 0.1, 0.1, 0.1 + / + + &grib2 + / + + &bdy_control + specified = T, F, F + nested = F, T, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_3NEST_namelist.wps (revision 212) @@ -0,0 +1,35 @@ +&share + wrf_core = 'ARW', + max_dom = 3, + start_date = '2024-07-01_01:00:00', '2024-07-01_01:00:00', '2024-07-01_01:00:00' + end_date = '2024-07-04_01:00:00', '2024-07-01_01:00:00', '2024-07-01_01:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + +&geogrid + parent_id = 0, 1, 2 + parent_grid_ratio = 1, 3, 3 + i_parent_start = 0, 51, 60 + j_parent_start = 0, 55, 45 + e_we = 153, 157, 157 !!! ATTENTION + e_sn = 153, 157, 157 !!! ATTENTION + geog_data_res = '32ppd_HRalb', '32ppd_HRalb', '32ppd_HRalb' + dx = 30000 + dy = 30000 + map_proj = 'mercator', + ref_lat = 24.66, + ref_lon = 146.76, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.input (revision 212) @@ -0,0 +1,59 @@ + &time_control + start_year = 2024, + start_month = 07, + start_day = 01, + start_hour = 01, + end_year = 2024, + end_month = 07, + end_day = 05, + end_hour = 01, + history_interval = 37 ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 10, + dx = 7000, + dy = 7000, + e_we = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 7 + / + + &physics + radt = 1, ! 1 is the best (change iradia in callphys) + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + non_hydrostatic = F + diff_6th_opt = 2, !! NECESSAIRE + diff_6th_factor = 0.2, !! NECESSAIRE + diff_opt = 1, + km_opt = 4, + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/EM_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-07-01_01:00:00' + end_date = '2024-07-05_01:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 201, + e_sn = 201, + geog_data_res = '64ppd_noHRti' + dx = 7000, + dy = 7000, + map_proj = 'mercator', + ref_lat = 24.8, + ref_lon = 146.9, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.input (revision 212) @@ -0,0 +1,78 @@ + &time_control + start_year = 2024, 2024 + start_month = 07, 04 + start_day = 01, 01 + start_hour = 01, 06 + end_year = 2024, 2024 + end_month = 07, 04 + end_day = 11, 05 + end_hour = 00, 06 + history_interval = 37, 74 ! *100s : 74 is each 2 MH, 37 is each MH !! attention, compatibilite avec time_step + frames_per_outfile = 24, 12 + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, F ! F en second pour que d1 init. d2 (smoother) + / + + &domains + time_step = 185 + max_dz = 4000 + dx = 296461.63, 20000 + dy = 296461.63, 20000 + e_we = 74, 211 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 50, 211 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 26, 61 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 1. +!!!!!! NESTING STUFF +! max_dom = 2 +! grid_id = 1, 2 +! parent_id = 0, 1 +! i_parent_start = 0, 69 ! attention au raccord avec WPS ! +! j_parent_start = 0, 69 ! attention au raccord avec WPS ! +! parent_grid_ratio = 1, 3 +! parent_time_step_ratio = 1, 5 !3 !4 !5 +! feedback = 1 +! smooth_option = 2 +! blend_width = 5 +!!!!!! NESTING STUFF + / + + &physics + radt = 1, !10 ! call each half hour with radt = 1850/dt + mars = 1, 1 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + non_hydrostatic = F, T + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES !! OPTION 2 SINON PB QH2O + diff_6th_factor = 0.5, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES !! FAIBLE SINON PB QH2O_ICE (si pas BDY) + !!diff_opt = 1, !! non sinon crash + !!km_opt = 4, + smdiv = 0., + emdiv = 0., + epssm = 0., + / + + &grib2 + / + + &bdy_control + specified = T, F + periodic_x = T, F + nested = F, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/GCM_namelist.wps (revision 212) @@ -0,0 +1,40 @@ +&share + wrf_core = 'ARW', + max_dom = 1 !!2, + start_date = '2024-10-01_01:00:00', '2024-10-01_01:00:00' + end_date = '2024-10-11_00:00:00', '2024-10-01_01:00:00' +! start_date = '2024-01-01_01:00:00', '2024-10-01_01:00:00' +! end_date = '2024-01-11_00:00:00', '2024-10-01_01:00:00' +! start_date = '2024-07-01_01:00:00', '2024-10-01_01:00:00' +! end_date = '2024-07-11_00:00:00', '2024-10-01_01:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 69 + j_parent_start = 1, 69 + e_we = 74, 211 !!! ATTENTION + e_sn = 50, 211 !!! ATTENTION + geog_data_res = '32ppd','64ppd_noHRti' + dx = 296461.63 !!1deglat 59292.325 + dy = 296461.63 !!1deglat 59292.325 + map_proj = 'mercator', + ref_lat = 0.00, + ref_lon = 175., + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.input (revision 212) @@ -0,0 +1,59 @@ + &time_control + start_year = 2025, + start_month = 07, + start_day = 02, + start_hour = 00 + end_year = 2025, + end_month = 07, + end_day = 03 !!15, + end_hour = 00, + history_interval = 1 !!74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + interval_seconds = 7400 ! assim assim + / + + &domains + time_step = 18 + dx = 50000, + dy = 50000, + e_we = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 31, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 2 + max_dz = 2000 + / + + &physics + radt = 2, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + h_mom_adv_order = 3 !! for better stability + h_sca_adv_order = 3 !! for better stability + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/HE_namelist.wps (revision 212) @@ -0,0 +1,37 @@ +&share + wrf_core = 'ARW', + max_dom = 2, + start_date = '2025-07-02_00:00:00' '2026-01-47_00:00:00' +! end_date = '2025-07-15_00:00:00', '2026-01-47_00:00:00' + end_date = '2025-07-03_00:00:00' + interval_seconds = 7400 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 40 + j_parent_start = 1, 40 + e_we = 150, 181 + e_sn = 150, 181 + geog_data_res = '32ppd','32ppd' + dx = 50000, + dy = 50000, + map_proj = 'mercator', + ref_lat = -10., + ref_lon = 100., + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/OXF' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.input (revision 212) @@ -0,0 +1,63 @@ + &time_control + start_year = 2024, + start_month = 01, + start_day = 05, + start_hour = 21, + end_year = 2024, + end_month = 01, + end_day = 06, + end_hour = 06, + history_interval = 1, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 37, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 0 + time_step_fract_num = 1, + time_step_fract_den = 2, + dx = 100, + dy = 100, + e_we = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 150, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 71, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 250 + force_sfc_in_vinterp = 5, + / + + &physics + radt = 60, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + init_MU=0.01 ! ... quasi no-wind + init_MV=0.01 ! ... quasi no-wind + init_LES=T ! LES MODE + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + specified = F + periodic_x = T + periodic_y = T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/LES_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-01-05_18:00:00' + end_date = '2024-01-06_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 150, + e_sn = 150, + geog_data_res = '32ppd' + dx = 100, + dy = 100, + map_proj = 'mercator', + ref_lat = -14.5 !!-14.57, + ref_lon = 175.4 !!175.47, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.input (revision 212) @@ -0,0 +1,75 @@ + &time_control + start_year = 2024, 2024, 2024, + start_month = 01, 01, 01, + start_day = 49, 49, 49, + start_hour = 06, 06, 06, + end_year = 2024, 2024, 2024, + end_month = 01, 01, 01, + end_day = 54, 54, 54, + end_hour = 06, 06, 06, + history_interval = 74, 74, 74, + frames_per_outfile = 12, 12, 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, T, T, ! F en second pour que d1 init. d2 (smoother) + / + + &domains + time_step = 40 + dx = 36000, 9000, 2250, !!alt 48000 12000 3000 + dy = 36000, 9000, 2250, + e_we = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, 61, 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 +!!!!! NESTING STUFF + max_dom = 3 !!!! PLEASE DO NOT FORGET callphysd*_.def + grid_id = 1, 2, 3, + parent_id = 0, 1, 2, + i_parent_start = 0, 57, 57, ! attention au raccord avec WPS ! + j_parent_start = 0, 57, 57, ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 4, 4, + parent_time_step_ratio = 1, 2, 2, + feedback = 0 + smooth_option = 2 + blend_width = 5 +!!!!! NESTING STUFF + / + + &physics + radt = 1, 1, 1 ! 1 is the best (change iradia in callphys) + mars = 0, 0, 0 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) +!!! init_MU=3. ! ... quasi no-wind +!!! init_MV=3. ! ... quasi no-wind + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, 2, 2, !! NECESSAIRE + diff_6th_factor = 0.2, 0.2, 0.2 !! NECESSAIRE + diff_opt = 1, + km_opt = 4, + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + specified = T, F, F + nested = F, T, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OMG_namelist.wps (revision 212) @@ -0,0 +1,35 @@ +&share + wrf_core = 'ARW', + max_dom = 3, + start_date = '2024-01-49_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00' + end_date = '2024-01-54_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + +&geogrid + parent_id = 0, 1, 2 + parent_grid_ratio = 1, 4, 4 + i_parent_start = 0, 57, 57 + j_parent_start = 0, 57, 57 + e_we = 153, 157, 157 !!! ATTENTION + e_sn = 153, 157, 157 !!! ATTENTION + geog_data_res = '64ppd', '64ppd', '64ppd' + dx = 36000 + dy = 36000 + map_proj = 'mercator', + ref_lat = 11.5, + ref_lon = 15.5, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input (revision 212) @@ -0,0 +1,59 @@ + &time_control + start_year = 2024, + start_month = 06, + start_day = 37, + start_hour = 18, + end_year = 2024, + end_month = 06, + end_day = 44, + end_hour = 18, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 10, + dx = 10000, + dy = 10000, + e_we = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) +! init_TI = 85., + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + pd_scalar = F, + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input_before_improv =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input_before_improv (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input_before_improv (revision 212) @@ -0,0 +1,55 @@ + &time_control + start_year = 2024, + start_month = 06, + start_day = 37, + start_hour = 18, + end_year = 2024, + end_month = 06, + end_day = 44, + end_hour = 18, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 10, + dx = 10000, + dy = 10000, + e_we = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 + / + + &physics + radt = 5, ! call each half hour with radt = 1850/dt + mars = 0 !!1, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, !!0 !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-06-37_18:00:00' + end_date = '2024-06-44_18:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 201, + e_sn = 201, + geog_data_res = '64ppd_noHRti' + dx = 10000, + dy = 10000, + map_proj = 'mercator', + ref_lat = 18.00, + ref_lon = -134.00, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.input (revision 212) @@ -0,0 +1,58 @@ + &time_control + start_year = 2024, + start_month = 06, + start_day = 02, + start_hour = 00, + end_year = 2024, + end_month = 06, + end_day = 07, + end_hour = 00, + history_interval = 37, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 24, + time_step_fract_num = 2, + time_step_fract_den = 3, + dx = 20000, + dy = 20000, + e_we = 181, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 181, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 5 + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + / + + &fdda + / + + &dynamics + pd_scalar = F !! no water cycle (mars = 0) + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PATH_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-06-02_00:00:00', '2026-01-47_00:00:00' + end_date = '2024-06-07_00:00:00', '2026-01-47_00:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 12 + j_parent_start = 1, 12 + e_we = 181, 25 + e_sn = 181, 25 + geog_data_res = '64ppd_noHRti' !'32ppd','32ppd' + dx = 20000, + dy = 20000, + map_proj = 'lambert' !!'mercator', + ref_lat = 32. !!19.0996, + ref_lon = -40. !!-33.2526, + truelat1 = 30., + truelat2 = 0.0, + stand_lon = -33., + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/callphys.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +50 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.input (revision 212) @@ -0,0 +1,60 @@ + &time_control + start_year = 2024, + start_month = 04, + start_day = 01, + start_hour = 02, + end_year = 2024, + end_month = 04, + end_day = 12, + end_hour = 00, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + interval_seconds = 7400 + / + + &domains + time_step = 37, + dx = 40000, + dy = 40000, + e_we = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 10 + / + + &physics + radt = 1, ! 1 is the best (change iradia in callphys) + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + non_hydrostatic = F + diff_6th_opt = 2, !! NECESSAIRE + diff_6th_factor = 0.2, !! NECESSAIRE + diff_opt = 1, + km_opt = 4, + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/PO/namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-04-01_02:00:00' + end_date = '2024-04-12_00:00:00' + interval_seconds = 7400 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 201, + e_sn = 201, + geog_data_res = '32ppd' + dx = 40000 + dy = 40000 + map_proj = 'polar', + ref_lat = -90, + ref_lon = 0, + truelat1 = -90.0, + truelat2 = 0.0, + stand_lon = 180.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.input (revision 212) @@ -0,0 +1,66 @@ + &time_control + start_year = 2024, + start_month = 07, + start_day = 01, + start_hour = 06, + end_year = 2024, + end_month = 07, + end_day = 02, + end_hour = 06, + history_interval = 37, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + interval_seconds = 7400 !! OPT OPT + / + + &domains + time_step = 25, + dx = 20000, + dy = 20000, + e_we = 51, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 51, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 5 + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + init_TI = 100., !! OPT OPT + init_AL = 0.27, !! OPT OPT + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + pd_scalar = F, !! OPT OPT + non_hydrostatic = F !! OPT OPT + !!time_step_sound = 10 + / + + &grib2 + / + + &bdy_control + spec_bdy_width = 3, !! OPT OPT + spec_zone = 1, !! OPT OPT 1+1 est instable + relax_zone = 2, !! OPT OPT + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TESTARSIA_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-07-01_06:00:00' + end_date = '2024-07-02_06:00:00' + interval_seconds = 7400 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 51, + e_sn = 51, + geog_data_res = '64ppd_noHRti' + dx = 20000, + dy = 20000, + map_proj = 'mercator', + ref_lat = -9.5, + ref_lon = 239.5, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.input (revision 212) @@ -0,0 +1,73 @@ + &time_control + start_year = 2024, 2024 + start_month = 05, 05, !!06, 06 + start_day = 01, 01, !!37, 37 + start_hour = 06, 06 + end_year = 2024, 2024 + end_month = 05, 05, !!06, 06 + end_day = 05, 05, !!43, 43 + end_hour = 06, 06 + history_interval = 74, 74 ! *100s : 74 is each 2 MH, 37 is each MH !! attention, compatibilite avec time_step + frames_per_outfile = 12, 12 + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, T ! F en second pour que d1 init. d2 (smoother) + / + fine_input_stream = 2, 2 + + &domains + time_step = 20 + dx = 20000, 5000 + dy = 20000, 5000 + e_we = 205, 209 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 205, 209 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, 61 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 +!!!!! NESTING STUFF + max_dom = 2 !! PLEASE DO NOT FORGET callphysd2_.def + grid_id = 1, 2 + parent_id = 0, 1 + i_parent_start = 0, 75 ! attention au raccord avec WPS ! + j_parent_start = 0, 33 ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 4 + parent_time_step_ratio = 1, 2 !!4 + feedback = 0 !!1 + smooth_option = 2 + blend_width = 5 +!!!!! NESTING STUFF + / + + &physics + radt = 1, 1 ! call each half hour with radt = 1850/dt + mars = 1, 1 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + / + + &grib2 + / + + &bdy_control + specified = T, F + nested = F, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_ARSIA_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 2, + start_date = '2024-05-01_06:00:00','2024-05-01_06:00:00' + end_date = '2024-05-05_06:00:00','2024-05-01_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 0, 4 + i_parent_start = 0, 75 + j_parent_start = 0, 33 + e_we = 205, 209 !!! ATTENTION + e_sn = 205, 209 !!! ATTENTION + geog_data_res = '32ppd_HRalb','32ppd_HRalb' + dx = 20000, + dy = 20000, + map_proj = 'mercator', + ref_lat = 5.00, + ref_lon = -120.00, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.input (revision 212) @@ -0,0 +1,72 @@ + &time_control + start_year = 2024, 2024 + start_month = 04, 04 + start_day = 01, 01 + start_hour = 06, 06 + end_year = 2024, 2024 + end_month = 04, 04 + end_day = 01, 01 !!05, 05 + end_hour = 12, 12 !!06, 06 + history_interval = 74, 74 ! *100s : 74 is each 2 MH, 37 is each MH !! attention, compatibilite avec time_step + frames_per_outfile = 12, 12 + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, T !!F ! F en second pour que d1 init. d2 (smoother) + / + fine_input_stream = 2, 2 + + &domains + time_step = 100 !!60 + dx = 60000, 20000 + dy = 60000, 20000 + e_we = 207, 211 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 207, 211 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, 61 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 +!!!!! NESTING STUFF + max_dom = 2 !! PLEASE DO NOT FORGET callphysd2_.def + grid_id = 1, 2 + parent_id = 0, 1 + i_parent_start = 0, 69 ! attention au raccord avec WPS ! + j_parent_start = 0, 69 ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 3 + parent_time_step_ratio = 1, 5 !3 !4 !5 + feedback = 1 + smooth_option = 2 + blend_width = 5 +!!!!! NESTING STUFF + / + + &physics + radt = 5, 5 !10, 10 ! call each half hour with radt = 1850/dt + mars = 1, 1 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + non_hydrostatic = F, T + diff_6th_opt = 1, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.8, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + / + + &grib2 + / + + &bdy_control + specified = T, F + nested = F, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_NEST_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 2, + start_date = '2024-04-01_06:00:00','2024-04-01_06:00:00' + end_date = '2024-04-05_06:00:00','2024-04-01_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, 1 + parent_grid_ratio = 1, 3 + i_parent_start = 1, 69 + j_parent_start = 1, 69 + e_we = 207,211 !!! ATTENTION + e_sn = 207,211 !!! ATTENTION + geog_data_res = '32ppd','64ppd_noHRti' + dx = 60000, !40000, + dy = 60000, !40000, + map_proj = 'mercator', + ref_lat = 5.00, + ref_lon = -120.00, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.input (revision 212) @@ -0,0 +1,57 @@ + &time_control + start_year = 2024, + start_month = 04, + start_day = 01, + start_hour = 06, + end_year = 2024, + end_month = 04, + end_day = 05, + end_hour = 06, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 20 + dx = 20000, + dy = 20000, + e_we = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 201, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 1, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + diff_6th_opt = 2 !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/TH_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-04-01_06:00:00' + end_date = '2024-04-05_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 201, + e_sn = 201, + geog_data_res = '64ppd_noHRti' + dx = 20000, + dy = 20000, + map_proj = 'mercator', + ref_lat = 5.00, + ref_lon = -120.00, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_callphys.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +100 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_launch_simu =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_launch_simu (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_launch_simu (revision 212) @@ -0,0 +1,141 @@ +#! /bin/bash + + +# +# USE: launch_simu & +# + +# +# must set env. var. WRF_current +# + +# PRELIMINARY STEP +# - model must be compiled and linked (wrf.exe, real.exe) +# - domain must be correctly set with geogrid.exe +# - put namelist.input, namelist.wps, callphys*.def + +# 'ls' must yield: +# callphys*.def launch_simu namelist.input namelist.wps real.exe wrf.exe +# + + +############################################################### +gcmstep=0 +wpsstep=1 +realstep=1 +wrfstep=1 +############################################################### +# GCM STEP + daygcm=593 +# WPS STEP + name='vm' # un nom de dossier non pris dans WRFFEED +# WRF STEP + proc=4 +############################################################### +# FOLDERS + wheregcm=$GCM + whereprep=$WRF_current/PREP_MARS + wherewps=$WRF_current/WPS +############################################################### + + +# +# LOG FILE HEADER +# + +where=$PWD +echo "check what's happening in " ${where}/log_launch + +rm -rf ${where}/log_launch 2> /dev/null +echo SIMU folder is ${where} > ${where}/log_launch +echo SIMU WILL BE LAUNCHED on the machine >> ${where}/log_launch +hostname >> ${where}/log_launch + +# +# 0. GCM step +# +if [ ${gcmstep} -eq 1 ] +then + + echo "0. GCM run and output preparation for WPS" >> ${where}/log_launch + + cd ${wheregcm} + ln -sf ../STARTBASE_64_48_25_t2/start/start${daygcm}.nc start.nc + ln -sf ../STARTBASE_64_48_25_t2/startfi/startfi${daygcm}.nc startfi.nc + + rm -rf restart.nc + rm -rf restartfi.nc + rm -rf log_gcm + + echo running the GCM - check in ${wheregcm}/log_gcm >> ${where}/log_launch + nohup gcm.e > log_gcm + + echo preparing outputs for WPS >> ${where}/log_launch + cd ${whereprep} + compile_and_exec >> ${where}/log_launch + +fi + +# +# 1. WPS step +# +if [ ${wpsstep} -eq 1 ] +then + + echo "1. WPS: horizontal interpolations of static and meteorological data" >> ${where}/log_launch + +# rm -rf ${where}/met_em + mkdir ${where}/met_em + cd ${where}/met_em + ln -sf ../namelist.wps . + + cd ${wherewps} + rm -rf WHERE_SIMU_IS + ln -sf ${where}/met_em WHERE_SIMU_IS + mkdir ./WRFFEED/current + nohup geogrid.exe >> ${where}/log_launch + nohup metgrid.exe >> ${where}/log_launch + mv ./WRFFEED/current ./WRFFEED/${name} + + ln -sf $FEED/WRFFEED/${name}/* ${where}/met_em/ + +fi + +# +# 2. REAL.EXE +# +if [ ${realstep} -eq 1 ] +then + + echo "2. REAL.EXE: vertical interpolation and preparation of boundary conditions" >> ${where}/log_launch + + cd ${where}/met_em + rm -rf rsl.* 2> /dev/null + rm -rf wrfbdy* 2> /dev/null + rm -rf wrfinput* 2> /dev/null + + ln -sf ../namelist.input . + ln -sf ../real.exe . + real.exe + cd .. + ln -sf ./met_em/wrf*d0* . + +fi + + +# +# 3. WRF.EXE +# +if [ ${wrfstep} -eq 1 ] +then + + echo "3. WRF.EXE: launch the mesoscale simulation" >> ${where}/log_launch + + rm -rf rsl.* + rm -rf wrfout* + + $WHERE_MPI/mpirun -np ${proc} wrf.exe < /dev/null & + echo OK, simulation launched ! + +fi + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.input (revision 212) @@ -0,0 +1,58 @@ + &time_control + start_year = 2024, + start_month = 11, + start_day = 34, + start_hour = 00, + end_year = 2024, + end_month = 11, + end_day = 38, + end_hour = 00, + history_interval = 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 10, + dx = 12000, + dy = 12000, + e_we = 401, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 121, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 !! 1 et 2 ne passent pas + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 0, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + pd_scalar = F, + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/VM_namelist.wps (revision 212) @@ -0,0 +1,36 @@ +&share + wrf_core = 'ARW', + max_dom = 1, + start_date = '2024-11-34_00:00:00' + end_date = '2024-11-38_00:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + + +&geogrid + parent_id = 1, + parent_grid_ratio = 1, + i_parent_start = 1, + j_parent_start = 1, + e_we = 401, + e_sn = 121, + geog_data_res = '32ppd_HRalb', + dx = 12000, + dy = 12000, + map_proj = 'mercator', + ref_lat = -8.00, + ref_lon = -68.00, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 0.0, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.assim =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.assim (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.assim (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +1352 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.gcm.mars1 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.gcm.mars1 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.gcm.mars1 (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +T +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +T +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +10 1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +T +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +T +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +T +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0 (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_PATH =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_PATH (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_PATH (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +15 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_noconvadj =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_noconvadj (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars0_noconvadj (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1 =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1 (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1 (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +T +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +T +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +T +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +T +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1.MGSdust =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1.MGSdust (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1.MGSdust (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +T +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +3 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +T +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +T +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +T +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1_noconvadj =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1_noconvadj (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.mars1_noconvadj (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +T +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +T +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +T +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +T +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.microscale =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.microscale (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/callphys.def.microscale (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +1 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +F +callcond (call CO2 condensation ?) +F +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +1 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/dustopacity.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/dustopacity.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/dustopacity.def (revision 212) @@ -0,0 +1,1 @@ +0.3 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/old_namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/old_namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/old_namelist.input (revision 212) @@ -0,0 +1,133 @@ + &time_control + run_days = 0, + run_hours = 0, + run_minutes = 0, + run_seconds = 0, + start_year = 2025, 2026, + start_month = 07 !06, 01, + start_day = 02 !49, 47, + start_hour = 00, 00, + start_minute = 00, 00, + start_second = 00, 00, + end_year = 2025, 2026, + end_month = 07, 07, + end_day = 15, 09, + end_hour = 00, 22, + end_minute = 00, 00, + end_second = 00, 00, + interval_seconds = 7400 + input_from_file = .true., .false. ! .false. en second pour que d1 init. d2 + fine_input_stream = 0 + history_interval = 74, 74, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 12, 12, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 !200 + / + + + &domains + time_step = 18 !37, + time_step_fract_num = 0, + time_step_fract_den = 1, +max_dom = 1 +s_we = 1,1, +e_we = 150,49, +s_sn = 1,1, +e_sn = 150,49, +s_vert = 1,1, +e_vert = 31,50, + num_metgrid_levels = 26, + max_dz = 2000. + sfcp_to_sfcp = .true. + lowest_lev_from_sfc = .false. + p_top_requested = 2, + dx = 50000, 5000, + dy = 50000, 5000, + grid_id = 1, 2, + parent_id = 0, 1, + i_parent_start = 0, 12, ! attention au raccord avec WPS ! + j_parent_start = 0, 12, ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 3, + parent_time_step_ratio = 1, 2, + feedback = 1, + smooth_option = 2 + blend_width = 5 + interp_type = 2, ! pas d'influence sur la stabilite + zap_close_levels = 0, ! set it to >0 when mountaneous domain + force_sfc_in_vinterp = 9 ! default: 8 (number of first layers prescribed for PBL) + / + eta_levels = 1.0,0.95,0.9,0.85,0.8,0.75,0.7,0.65,0.6,0.55,0.5,0.45,0.4,0.35,0.3,0.25,0.2,0.15,0.1,0.05,0 + lagrange_order = 1 + interp_type = 2 + adjust_heights = .false. + ztop = 60000. + + + + + &physics + mp_physics = 0, 0, + ra_lw_physics = 0, 0, + ra_sw_physics = 0, 0, + radt = 2, 10, + sf_sfclay_physics = 0, 0, + sf_surface_physics = 0, 0, + bl_pbl_physics = 0, 0, + bldt = 0, 0, + cu_physics = 0, 0, + cudt = 400 !200, 1, + isfflx = 1, + ifsnow = 0, + icloud = 1, + num_soil_layers = 5, + / + + &fdda + / + + &dynamics + dyn_opt = 2, + rk_ord = 3, + diff_opt = 0, + km_opt = 4, + diff_6th_opt = 0, + diff_6th_factor = 0.12, + base_temp = 220. + base_pres = 610. + base_lapse = 30. + damp_opt = 0, + zdamp = 3000, + dampcoef = 0.05, + w_damping = 0 + smdiv = 0.1, 0.1, + emdiv = 0.01, 0.01, + epssm = 0.1, 0.1, + non_hydrostatic = .true., .true., + time_step_sound = 4, 4, + h_mom_adv_order = 3, !5, 5, + v_mom_adv_order = 3, 3, + h_sca_adv_order = 3, !5, 5, + v_sca_adv_order = 3, 3, + + + &grib2 + / + + &bdy_control + spec_bdy_width = 5, + spec_zone = 1, + relax_zone = 4, + specified = .true., .false., + nested = .false., .true., + / + + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/callphys.def_dust =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/callphys.def_dust (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/callphys.def_dust (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +T +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +75 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +1 0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +T F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +T F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/namelist.input_dust =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/namelist.input_dust (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/namelist.input_dust (revision 212) @@ -0,0 +1,66 @@ + &time_control + start_year = 2024, + start_month = 07, + start_day = 01, + start_hour = 06, + end_year = 2024, + end_month = 07, + end_day = 02, + end_hour = 06, + history_interval = 37, ! *100s : 74 is each 2 MH, 37 is each MH + frames_per_outfile = 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + interval_seconds = 7400 !! OPT OPT + / + + &domains + time_step = 25, + dx = 20000, + dy = 20000, + e_we = 51, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 51, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 5 + / + + &physics + radt = 1, ! call each half hour with radt = 1850/dt + mars = 2, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) + init_TI = 100., !! OPT OPT + init_AL = 0.27, !! OPT OPT + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_6th_factor = 0.2, !! NECESSAIRE DANS LES ZONES DE VENT FAIBLES + diff_opt = 1, + km_opt = 4, + pd_scalar = F, !! OPT OPT + non_hydrostatic = F !! OPT OPT + !!time_step_sound = 10 + / + + &grib2 + / + + &bdy_control + spec_bdy_width = 3, !! OPT OPT + spec_zone = 1, !! OPT OPT 1+1 est instable + relax_zone = 2, !! OPT OPT + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/stress.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/stress.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/dust/stress.def (revision 212) @@ -0,0 +1,2 @@ +0.01 +1. Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/callphys.def (revision 212) @@ -0,0 +1,116 @@ +##General options +##~~~~~~~~~~~~~~~ +#Run with or without tracer transport ? +tracer=.false. + +#Diurnal cycle ? if diurnal=False, diurnal averaged solar heating +diurnal=.true. + +#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" +season = .true. + +#write some more output on the screen ? +lwrite = .false. + +#Save statistics in file "stats.nc" ? +stats =.false. + +#Save EOF profiles in file "profiles" for Climate Database? +calleofdump = .false. + +## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +## ~~~~~~~~~~~~~ +# =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, +# =4 Mars Year 24 from TES assimilation (same as =24 for now) +# =24 Mars Year 24 from TES assimilation (ie: MCD reference case) +# =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) +# =26 Mars Year 26 from TES assimilation +iaervar = 24 +# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case) +tauvis=0.2 +# Dust vertical distribution: +# (=0: old distrib. (Pollack90), =1: top set by "topdustref", +# =2: Viking scenario; =3 MGS scenario) +iddist = 3 +# Dust top altitude (km). (Matters only if iddist=1) +topdustref = 55. + +## Physical Parameterizations : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call radiative transfer ? +callrad = .true. +# call NLTE radiative schemes ? matters only if callrad=T +callnlte = .true. +# call CO2 NIR absorption ? matters only if callrad=T +callnirco2 = .true. +# call turbulent vertical diffusion ? +calldifv = .true. +# call convective adjustment ? +calladj = .true. +# call CO2 condensation ? +callcond =.true. +# call thermal conduction in the soil ? +callsoil = .true. +# call Lott's gravity wave/subgrid topography scheme ? +calllott = .false. + +## Radiative transfer options : +## ~~~~~~~~~~~~~~~~~~~~~~~~~~ +# the rad.transfer is computed every "iradia" physical timestep +iradia = 1 +# Output of the exchange coefficient mattrix ? for diagnostic only +callg2d = .false. +# Rayleigh scattering : (should be .false. for now) +rayleigh = .false. + +## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): +## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# DUST: Transported dust ? (if >0, use 'dustbin' dust bins) +dustbin = 0 +# DUST: Radiatively active dust ? (matters if dustbin>0) +active = .false. +# DUST: use mass and number mixing ratios to predict dust size ? +# (must also have dustbin=1) +doubleq = .false. +# DUST: lifted by GCM surface winds ? +lifting = .false. +# DUST: lifted by dust devils ? +callddevil = .false. +# DUST: Scavenging by CO2 snowfall ? +scavenging = .false. +# DUST/WATERICE: Gravitationnal sedimentation ? +sedimentation = .false. +# WATERICE: Radiatively active transported atmospheric water ice ? +activice = .false. +# WATER: Compute water cycle +water = .false. +# WATER: current permanent caps at both poles. True IS RECOMMENDED +# (with .true., North cap is a source of water and South pole +# is a cold trap) +caps = .false. +#caps = .true. +# PHOTOCHEMISTRY: include chemical species +photochem = .false. + +## Thermospheric options (relevant if tracer=T) : +##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +# call thermosphere ? +callthermos = .false. +# WATER: included without cycle (only if water=.false.) +thermoswater = .false. +# call thermal conduction ? (only if callthermos=.true.) +callconduct = .false. +# call EUV heating ? (only if callthermos=.true.) +calleuv=.false. +# call molecular viscosity ? (only if callthermos=.true.) +callmolvis = .false. +# call molecular diffusion ? (only if callthermos=.true.) +callmoldiff = .false. +# call thermospheric photochemistry ? (only if callthermos=.true.) +thermochem = .false. +# date for solar flux calculation: (1985 < date < 2002) +## (Solar min=1996.4 ave=1993.4 max=1990.6) +solarcondate = 1993.4 + + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.input (revision 212) @@ -0,0 +1,55 @@ + &time_control + start_year = 2024, !! (p1) Start Martian Year (20XX for MY XX) + start_month = 01, !! (p1) Start Martian Month + start_day = 01, !! (p1) Start Martian Day + start_hour = 06, !! (p1) Start Martian Hour (at longitude 0) + end_year = 2024, !! (p1) End Martian Year (20XX for MY XX) + end_month = 01, !! (p1) End Martian Month + end_day = 04, !! (p1) End Martian Day + end_hour = 06, !! (p1) End Martian Hour (at longitude 0) + history_interval = 37, + frames_per_outfile = 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + / + + &domains + time_step = 25, + dx = 20000, + dy = 20000, + e_we = 61, + e_sn = 61, + e_vert = 61, !! (r)(p2) Number of vertical levels + p_top_requested = 5. !! (p3) Chosen value of pressure at the top of the model + num_metgrid_levels = 19 + / + + &physics + !init_TI = 85., !! Define constant thermal inertia value + num_soil_layers = 18 !!! NOUVELLE PHYSIQUE + / + + &fdda + / + + &dynamics + tiso = 100. + / + + &grib2 + / + + &bdy_control + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/namelist.wps (revision 212) @@ -0,0 +1,38 @@ +&share + wrf_core = 'ARW', !! [do not modify: choice of dynamical core] + max_dom = 1 !! (r)(n)(p2) Total number of domains + start_date = '2024-01-01_06:00:00','2024-05-01_06:00:00','2024-05-01_06:00:00' !! YYYY-MM-DD_HH:mm:ss start date + end_date = '2024-01-04_06:00:00','2024-05-01_06:00:00','2024-05-01_06:00:00' !! YYYY-MM-DD_HH:mm:ss end date + interval_seconds = 3700 !! (p2) Frequency of large-scale fields update (s) + io_form_geogrid = 2, !! [do not modify: choice of NETCDF outputs] + debug_level = 0, !! verbose level of the programs + opt_output_from_geogrid_path='./' !! location of the geogrid outputs +/ + + +&geogrid + parent_id = 0, !! (n)(p2) Associated parent domain + parent_grid_ratio = 1, !! (n)(p2) Ratio of horizontal resolution parent/nest + i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner + j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner + e_we = 61, + e_sn = 61, + dx = 20000 + dy = 20000 + geog_data_res = '64ppd', !! choice of static data sources + !! NB: possible: '64ppd', '32ppd', ... + !! NB: please glance at geogrid/GEOGRID.TBL + map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' + ref_lat = 18, !! north latitude of the center of the domain + ref_lon = -136, !! east longitude of the center of the domain + truelat1 = 18, !! (lambert or polar) lat position of projection cone + truelat2 = 0.0, !! [do not modify] + stand_lon = -136, !! (lambert or polar) lon position of projection cone + geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] +/ + +&metgrid + fg_name = './WPSFEED/LMD' !! [do not modify: symbolic link in the WPS folder] + io_form_metgrid = 2, !! [do not modify: choice of NETCDF outputs] + opt_output_from_metgrid_path='./WRFFEED/current' !! [do not modify: symbolic link] +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/run.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/run.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_test/run.def (revision 212) @@ -0,0 +1,2 @@ +# some definitions for the physics, in file 'callphys.def' +INCLUDEDEF=callphys.def Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphys.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphys.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphys.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +20 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd2_.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd2_.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd2_.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +40 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd3_.def =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd3_.def (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/callphysd3_.def (revision 212) @@ -0,0 +1,94 @@ +General options +~~~~~~~~~~~~~~~ +tracer (Run with or without tracer transport ?) +F +diurnal (Diurnal cycle ? if diurnal=F, diurnal averaged solar heating) +T +season (Seasonal cycle ? if season=F, Ls stays constant like in "start") +T +lwrite (want some more output on the screen ?) +F +stats (Saving statistics in file "cumul" ?) -- noWRF +F +calleofdump (Saving EOF profiles in file "profiles" for Climate Database ?) -- noWRF +F +Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) +~~~~~~~~~~~~~ +iaervar (=1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario +4 (=4 Mars Year 24 from TES assimilation) +iddist (Dust vertical distribution: =0: old distrib. (Pollack90) +3 (=1: top set by "topdustref"; =2: Viking scenario; =3 MGS scenario ) +topdustref (Dust top altitude (km). Matter only if iddist=1) +55. +Physical Parameterizations : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +callrad (call radiative transfer ?) +T +callnlte (call NLTE radiative schemes ? matter only if callrad=T) +F +callnirco2 (call CO2 NIR absorption ? matter only if callrad=T) +T +calldifv (call turbulent vertical diffusion ?) +T +calladj (call convective adjustment ?) +T F +callcond (call CO2 condensation ?) +T +callsoil (call thermal conduction in the soil ?) +T +calllott (call Lott's gravity wave/subgrid topography scheme ?) -- noWRF +F +Radiative transfer options : +~~~~~~~~~~~~~~~~~~~~~~~~~~ +iradia (the rad.transfer is computed every "iradia" physical timestep) +80 +callg2d (Output of the exchange coefficient mattrix ? for diagnostic only) +F +rayleigh (Rayleigh scattering : should be =F for now) +F +Tracer (dust water, ice and/or chemical species) options (use if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +dustbin (DUST: Transported dust ? (if >0, uses q(1) to q(dustbin)) +0 +active (DUST: Radiatively active dust ? (uses q(1) to q(dustbin)) +F +doubleq (DUST: needs dustbin=1, use mass (q(1)) and number (q(2)) mixing ratio to predict dust size ?) +F +lifting (DUST: lifted by GCM surface winds ?) +F +dustdevil (DUST: lifted by dust devils ?) +F +scavenging (DUST: Scavenging by CO2 snowfall ?) +F +sedimentation (DUST/WATERICE: Gravitationnal sedimentation ?) +F +iceparty (WATERICE: Water cycle includes water ice mixing ratio q(nqmx-1)) +F +activice (WATERICE: Radiatively active transported atmospheric water ice ?) +F +water (WATER: Compute water cycle using q(nqmx) ) +F +caps (WATER: put the current permanent caps at both poles) +F +photochem (PHOTOCHEMISTRY: chemical species included) +F +Thermospheric options (relevant if tracer=T) : +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +callthermos (call thermosphere ?) +F +thermoswater (WATER: included without cycle only if water=F) +F +callconduct (call thermal conduction ? matter only if callthermos=T) +F +calleuv (call EUV heating ? matter only if callthermos=T) +F +callmolvis (call molecular viscosity ? matter only if callthermos=T) +F +callmoldiff (call molecular diffusion ? matter only if callthermos=T) +F +thermochem (call thermospheric photochemistry ? matter only if callthermos=T) +F +solarcondate (date for solar flux calculation: 1985 < date < 2002)) +1993.4 (Solar min=1996.4 ave=1993.4 max=1990.6) + + Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/makemeso_input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/makemeso_input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/makemeso_input (revision 212) @@ -0,0 +1,6 @@ +1 +153 +153 +61 +3 +1 Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.input =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.input (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.input (revision 212) @@ -0,0 +1,76 @@ + &time_control + start_year = 2024, 2024, 2024, + start_month = 01, 01, 01, + start_day = 49, 49, 49, + start_hour = 06, 06, 06, + end_year = 2024, 2024, 2024, + end_month = 01, 01, 01, + end_day = 54, 54, 54, + end_hour = 06, 06, 06, + history_interval = 37, 37, 37, +! history_interval_s = 1,1,1 + frames_per_outfile = 24, 24, 24, + restart = .false. + restart_interval = 8880, + io_form_history = 2 + io_form_restart = 2 + io_form_input = 2 + io_form_boundary = 2 + debug_level = 0 + input_from_file = T, T, T, ! F en second pour que d1 init. d2 (smoother) + / + + &domains + time_step = 40 + dx = 36000, 9000, 2250, !!alt 48000 12000 3000 + dy = 36000, 9000, 2250, + e_we = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_sn = 153, 157, 157, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + e_vert = 61, 61, 61, ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS + p_top_requested = 3 +!!!!! NESTING STUFF + max_dom = 3 !!!! PLEASE DO NOT FORGET callphysd*_.def + grid_id = 1, 2, 3, + parent_id = 0, 1, 2, + i_parent_start = 0, 57, 57, ! attention au raccord avec WPS ! + j_parent_start = 0, 57, 57, ! attention au raccord avec WPS ! + parent_grid_ratio = 1, 4, 4, + parent_time_step_ratio = 1, 2, 2, + feedback = 0 + smooth_option = 2 + blend_width = 5 +!!!!! NESTING STUFF + / + + &physics + radt = 1, 1, 1 ! 1 is the best (change iradia in callphys) + mars = 0, 0, 0 ! IF YOU CHANGE THIS YOU MUST RECOMPILE THE PHYSICS (AND CHANGE CALLPHYS.DEF) +!!! init_MU=3. ! ... quasi no-wind +!!! init_MV=3. ! ... quasi no-wind + / + + &fdda + / + + &dynamics + diff_6th_opt = 2, 2, 2, !! NECESSAIRE + diff_6th_factor = 0.2, 0.2, 0.2 !! NECESSAIRE + diff_opt = 1, + km_opt = 4, + pd_scalar = F !! no water cycle (mars = 0) + / + + &grib2 + / + + &bdy_control + specified = T, F, F + nested = F, T, T + / + +!! DO NOT MODIFY + &namelist_quilt + nio_tasks_per_group = 0, + nio_groups = 1, + / +!! DO NOT MODIFY Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.wps =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.wps (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/yeyeye/namelist.wps (revision 212) @@ -0,0 +1,35 @@ +&share + wrf_core = 'ARW', + max_dom = 3, + start_date = '2024-01-49_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00' + end_date = '2024-01-54_06:00:00', '2024-01-49_06:00:00', '2024-01-49_06:00:00' + interval_seconds = 3700 + io_form_geogrid = 2, + debug_level = 0, + opt_output_from_geogrid_path='./' +/ + +&geogrid + parent_id = 0, 1, 2 + parent_grid_ratio = 1, 4, 4 + i_parent_start = 0, 57, 57 + j_parent_start = 0, 57, 57 + e_we = 153, 157, 157 !!! ATTENTION + e_sn = 153, 157, 157 !!! ATTENTION + geog_data_res = '64ppd', '64ppd', '64ppd' + dx = 36000 + dy = 36000 + map_proj = 'mercator', + ref_lat = 11.5, + ref_lon = 15.5, + truelat1 = 0.0, + truelat2 = 0.0, + stand_lon = 00.00, + geog_data_path = './WPS_GEOG', +/ + +&metgrid + fg_name = './WPSFEED/LMD' + io_form_metgrid = 2, + opt_output_from_metgrid_path='./WRFFEED/current' +/ Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/elapsed =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/elapsed (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/elapsed (revision 212) @@ -0,0 +1,8 @@ +#! /bin/bash + +if [[ -f rsl.error.0000 ]] +then + more rsl.error.0000 | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}' +else + more log_wrf | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}' +fi Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome (revision 212) @@ -0,0 +1,61 @@ +####################################################################### +####################################################################### nom +# @ job_name = LMD_LES_MARS +####################################################################### MPICH / serial +# @ job_type = MPICH +####################################################################### 8 / 16 / 20 / 32 / 64 +# @ total_tasks = 8 +####################################################################### A / B / C / AP / BP / CP +# @ class = AP +####################################################################### 24gb (cf. llps pendant fonctionnement) +# @ resources = ConsumableCpus(1)ConsumableMemory(1gb) +####################################################################### real ou ideal ??? +# @ environment = real=0 ; tmpdir=$(job_name).$(jobid) ; para=$(job_type) +####################################################################### + +####################################################################### +####################################################################### +# @ output = $(job_name).$(jobid)/LOG +# @ error = $(output) +# @ blocking = unlimited +# @ notification = never +# @ queue +grep "^# @" $0 +set -vx +####################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR +\rm zeLAST +ln -sf ${tmpdir} zeLAST +cd zeLAST +####################################################################### +cp $LOADL_STEP_INITDIR/launch.gnome . +cp $LOADL_STEP_INITDIR/*.def . +cp $LOADL_STEP_INITDIR/*.exe . +cp $LOADL_STEP_INITDIR/namelist* . +cp $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale + ### au prealable !!! +####################################################################### +if [[ "${real}" == "0" ]] +then + cp $LOADL_STEP_INITDIR/input_* . + cp $LOADL_STEP_INITDIR/ze_hill . ## topo + cp $LOADL_STEP_INITDIR/levels . ## custom vert levels + ##mpirun ideal.exe ##marche pas + ##mkdir ideal + ##mv rsl.* ideal/ +else + cp $LOADL_STEP_INITDIR/wrfb* . +fi +####################################################################### +#### avant wrf.exe, necessaire ??? +#declare -x MP_STACK_SIZE=640000000 +#declare -x MP_STACK_SIZE=64000000 +####################################################################### +if [[ "${para}" == "MPICH" ]] +then + mpirun wrf.exe +else + ./wrf.exe +fi +####################################################################### Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome.ideal =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome.ideal (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/launch.gnome.ideal (revision 212) @@ -0,0 +1,46 @@ +####################################################################### +####################################################################### nom +# @ job_name = PREP_LES +####################################################################### real ou ideal ??? +# @ environment = real=0 ; tmpdir=$(job_name).$(jobid) +####################################################################### MPICH / serial +# @ job_type = serial +####################################################################### A / B / C / AP / BP / CP +# @ class = A +####################################################################### 16gb OK jusque 281 (non teste plus loin) +# @ resources = ConsumableCpus(1)ConsumableMemory(16gb) +####################################################################### + + +####################################################################### +####################################################################### +# @ output = $(job_name).$(jobid)/LOG +# @ error = $(output) +# @ notification = never +# @ queue +grep "^# @" $0 +set -vx +####################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR +\rm zeLAST +ln -sf ${tmpdir} zeLAST +cd zeLAST +####################################################################### +cp $LOADL_STEP_INITDIR/*.def . +cp $LOADL_STEP_INITDIR/*.exe . +cp $LOADL_STEP_INITDIR/namelist* . +cp $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale + ### au prealable !!! +####################################################################### +if [[ "${real}" == "0" ]] +then + cp $LOADL_STEP_INITDIR/input_* . + cp $LOADL_STEP_INITDIR/ze_hill . ## topo + cp $LOADL_STEP_INITDIR/levels . ## custom vert levels +else + ln -sf $LOADL_STEP_INITDIR/wrfb* . +fi +####################################################################### +ideal.exe ### marche pas en parallele, pb avec execution ensuite +####################################################################### Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/mpd.hosts =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/mpd.hosts (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/mpd.hosts (revision 212) @@ -0,0 +1,7 @@ +auric.lmd.jussieu.fr +smith.lmd.jussieu.fr + +auric.lmd.jussieu.fr +smith.lmd.jussieu.fr +widor.lmd.jussieu.fr +parker.lmd.jussieu.fr Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/save_this_simu =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/save_this_simu (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/save_this_simu (revision 212) @@ -0,0 +1,48 @@ +#! /bin/bash + +\rm simtime 2> /dev/null + +mkdir results + +if [[ -f rsl.error.0000 ]] +then + + timesec=$(more rsl.error.0000 | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}') + touch simtime + echo ${timesec} seconds + echo ${timesec} seconds >> simtime + + head -10000 rsl.out.0000 > red_rsl.out.0000 + head -4000 rsl.error.0000 > red_rsl.error.0000 + cp met_em/rsl.out.0000 real_rsl.out.0000 + mv red_rsl.* results/ + mv real_rsl.out.0000 results/ + mv simtime results/ + +else + + timesec=$(more log_wrf | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}') + touch simtime + echo ${timesec} seconds + echo ${timesec} seconds >> simtime + + head -10000 log_wrf > red_log_wrf + mv red_log_wrf results/ + mv simtime results/ + mv met_em/log_real results/ + +fi + +tar czvf def.tar.gz namelist.input namelist.wps callphys*.def makemeso* #launch_simu makemeso* +cp def.tar.gz results/ + +mv *.def results/ +mv results/callphys*def ./ + +mv wrfout* results/ +ls -lh results +du -bh results +echo '*****************************************************' +echo '*** BEWARE ! *** output files were moved in the folder named results !!!' +echo 'what you may keep: results folder and def.tar.gz file' +echo '*****************************************************' Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/levspe.pro (revision 212) @@ -0,0 +1,94 @@ +pro levspe + + +; +; TWEAK PARAM +; +@param.idl + + +; +; HARD PARAM +; +!p.multi=[0,2,2] +set_plot, 'ps' +device, file='plot.ps', /landscape +tinv=5 ;; intervalle plot +hache=10. ;; scale height +psurf=610. ;; pression surface -- ne change pas les eta levels + +; +; PRELIM +; +param = 'e_vert='+string(nlev,'(I0)')+' ztop='+string(altmax,'(I0)')+' e='+string(epsilon,'(I0)')+' c='+string(elong_cos,'(F4.2)') +; +epsilon = epsilon / (nlev-1) / 100. +; +exposant = !pi/nlev/elong_cos +; +; alpha (plus grand ecart en km) est determine pour que max(altitudes)=altmax +; +alpha = altmax / ( (nlev-1)/2. - sin(2.*exposant*(nlev-1))/4./exposant + epsilon*(nlev-1)^2/2. ) +print, 'alpha (km)', alpha +; +x=[findgen(nlev)] + + +; +; CALC +; +ecart = alpha - alpha*cos(exposant*x)^2 + epsilon*alpha*x +altitudes = alpha*x/2. - alpha*sin(2.*exposant*x)/4./exposant + epsilon*alpha*x^2/2. +logpressions = alog(psurf) - altitudes/hache +pressions=exp(logpressions) +ptop=psurf*exp(-altmax/hache) & print, ptop +etas=(pressions-ptop)/(psurf-ptop) +etas(nlev-1)=0. + + +; +; PRINT +; +print, 'eta levels' +print, etas +; +print, 'pressure' +pi=etas*(psurf-ptop)+ptop +print, pi +print, 'ptop', ptop +; +print, 'pseudo-altitude' +pseudo=10.*alog(610./pi) +print, 1000.*10.*alog(610./pi) + + +; +; PLOT +; +plot, x, etas, title='ETA '+param, xtickinterval=tinv +!psym=7 & oplot, x, etas +!psym=0 +plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv +!psym=7 & oplot, x, pseudo +!psym=0 +plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv +!psym=7 & oplot, x, pi +!psym=0 +; +diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2) +; +plot, x, ecart, title='ECART (km, L=theor x=calc)', xtickinterval=tinv +!psym=7 & oplot, x, diff +!psym=0 + + +; +; PRINT +; +openw, 1, 'levels' +for i=0, nlev-1 do begin + printf, 1, etas(i);,',' +endfor +close, 1 + +end Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/param.idl =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/param.idl (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/param.idl (revision 212) @@ -0,0 +1,44 @@ +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=120. ;; altitude maximale en km +;; +epsilon=30. ;; augmentation totale en pourcentage de l'ecart max (modifie inflexion) + ;; -- epsilon = 0 : point d'inflexion parfaitement plat + ;; -- epsilon = 100 : pas de point d'inflexion +elong_cos=2./3. ;; plus petit rapproche l'inflexion du sol ;; 2/3 parfait pour un nlev divisible par 3 +elong_cos=0.64 +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;nlev=61 +;altmax=150. +;epsilon=30. +;elong_cos=2./3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=150. +epsilon=100. +elong_cos=2. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=201 +altmax=150. +epsilon=0. +elong_cos=2. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +;nlev=101 +;altmax=80. ;;100. +;epsilon=100. +;elong_cos=60. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=60. ;;100. +epsilon=100. +elong_cos=2.*3.14/3. +;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; +nlev=61 +altmax=20. ;;100. +epsilon=100. +elong_cos=2.*3.14/3. +elong_cos=20. Index: /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/simple.pro =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/simple.pro (revision 212) +++ /trunk/MESOSCALE/LMD_MM_MARS/SIMU/RUN/vert_level/simple.pro (revision 212) @@ -0,0 +1,94 @@ +pro simple, z_scale + + +nlev=101 +; +; z-scale doit etre suffisamment bas pour que l avant-dernier niveau soit assez haut +; mais pas trop sinon instable (il faut avoir des niveaux espaces en haut semble-t-il) +; +z_scale = .25 +;z_scale = .15 marche pas (niveaux pas assez espaces en haut ?) +ptop = 0.2 +npbl = 5 + ;z_scale = .15 + ;ptop = 0.05 ;marche pas + +z_scale = .20 & ptop = 0.1 & npbl = 6 +z_scale = .25 & ptop = 0.5 & npbl = 6 & nlev=61 ;; HOLDEN NESTS + +psurf=610. + +pbl = [1.0000,0.9995,0.9980,0.9950,0.9850,0.9700,0.9400,0.9000] +x=[findgen(nlev)] +nlev = nlev - npbl + 1 + +etas=[findgen(nlev)] +FOR k=0, nlev-1 DO BEGIN + etas(k) = (exp(-float(k)/float(nlev-1)/z_scale) - exp(-1./z_scale)) / (1.-exp(-1./z_scale)) +ENDFOR +etas(nlev-1)=0. + +etas = [pbl(0:npbl-1),etas(1:nlev-1)] +nlev = nlev + npbl - 1 +print, nlev + +;pressions = ptop + etas * (psurf-ptop) +;zzz = 10. * alog ( psurf / pressions ) +;print, zzz + +; +; HARD PARAM +; +!p.multi=[0,2,2] +;set_plot, 'ps' +;device, file='plot.ps', /landscape + + + +; +; PRINT +; +print, 'eta levels' +print, etas +; +print, 'pressure' +pi=etas*(psurf-ptop)+ptop +print, pi +print, 'ptop', ptop +; +print, 'pseudo-altitude' +pseudo=10.*alog(610./pi) +print, 1000.*10.*alog(610./pi) + + +; +; PLOT +; +plot, x, etas, title='ETA', xtickinterval=tinv +!psym=7 & oplot, x, etas +!psym=0 +plot, x, pseudo, title='ALTITUDE (km)', xtickinterval=tinv +!psym=7 & oplot, x, pseudo +!psym=0 +plot, x, pi, title='PRESSURE (Pa)', ylog=1, xtickinterval=tinv +!psym=7 & oplot, x, pi +!psym=0 +; +diff=pseudo - shift(pseudo,-1) & diff=-diff(0:nlev-2) +; +plot, x, diff, title='ECART (km, L=theor x=calc)', xtickinterval=tinv +!psym=7 & oplot, x, diff +!psym=0 + +;stop + +; +; PRINT +; +openw, 1, 'levels' +for i=0, nlev-1 do begin + printf, 1, etas(i);,',' +endfor +close, 1 + +end Index: unk/MESOSCALE/LMD_MM_MARS/SIMU/elapsed =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/elapsed (revision 211) +++ (revision ) @@ -1,8 +1,0 @@ -#! /bin/bash - -if [[ -f rsl.error.0000 ]] -then - more rsl.error.0000 | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}' -else - more log_wrf | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}' -fi Index: unk/MESOSCALE/LMD_MM_MARS/SIMU/launch.gnome =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/launch.gnome (revision 211) +++ (revision ) @@ -1,61 +1,0 @@ -####################################################################### -####################################################################### nom -# @ job_name = LMD_LES_MARS -####################################################################### MPICH / serial -# @ job_type = MPICH -####################################################################### 8 / 16 / 20 / 32 / 64 -# @ total_tasks = 8 -####################################################################### A / B / C / AP / BP / CP -# @ class = AP -####################################################################### 24gb (cf. llps pendant fonctionnement) -# @ resources = ConsumableCpus(1)ConsumableMemory(1gb) -####################################################################### real ou ideal ??? -# @ environment = real=0 ; tmpdir=$(job_name).$(jobid) ; para=$(job_type) -####################################################################### - -####################################################################### -####################################################################### -# @ output = $(job_name).$(jobid)/LOG -# @ error = $(output) -# @ blocking = unlimited -# @ notification = never -# @ queue -grep "^# @" $0 -set -vx -####################################################################### -mkdir $LOADL_STEP_INITDIR/${tmpdir} -cd $LOADL_STEP_INITDIR -\rm zeLAST -ln -sf ${tmpdir} zeLAST -cd zeLAST -####################################################################### -cp $LOADL_STEP_INITDIR/launch.gnome . -cp $LOADL_STEP_INITDIR/*.def . -cp $LOADL_STEP_INITDIR/*.exe . -cp $LOADL_STEP_INITDIR/namelist* . -cp $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale - ### au prealable !!! -####################################################################### -if [[ "${real}" == "0" ]] -then - cp $LOADL_STEP_INITDIR/input_* . - cp $LOADL_STEP_INITDIR/ze_hill . ## topo - cp $LOADL_STEP_INITDIR/levels . ## custom vert levels - ##mpirun ideal.exe ##marche pas - ##mkdir ideal - ##mv rsl.* ideal/ -else - cp $LOADL_STEP_INITDIR/wrfb* . -fi -####################################################################### -#### avant wrf.exe, necessaire ??? -#declare -x MP_STACK_SIZE=640000000 -#declare -x MP_STACK_SIZE=64000000 -####################################################################### -if [[ "${para}" == "MPICH" ]] -then - mpirun wrf.exe -else - ./wrf.exe -fi -####################################################################### Index: unk/MESOSCALE/LMD_MM_MARS/SIMU/launch.gnome.ideal =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/launch.gnome.ideal (revision 211) +++ (revision ) @@ -1,46 +1,0 @@ -####################################################################### -####################################################################### nom -# @ job_name = PREP_LES -####################################################################### real ou ideal ??? -# @ environment = real=0 ; tmpdir=$(job_name).$(jobid) -####################################################################### MPICH / serial -# @ job_type = serial -####################################################################### A / B / C / AP / BP / CP -# @ class = A -####################################################################### 16gb OK jusque 281 (non teste plus loin) -# @ resources = ConsumableCpus(1)ConsumableMemory(16gb) -####################################################################### - - -####################################################################### -####################################################################### -# @ output = $(job_name).$(jobid)/LOG -# @ error = $(output) -# @ notification = never -# @ queue -grep "^# @" $0 -set -vx -####################################################################### -mkdir $LOADL_STEP_INITDIR/${tmpdir} -cd $LOADL_STEP_INITDIR -\rm zeLAST -ln -sf ${tmpdir} zeLAST -cd zeLAST -####################################################################### -cp $LOADL_STEP_INITDIR/*.def . -cp $LOADL_STEP_INITDIR/*.exe . -cp $LOADL_STEP_INITDIR/namelist* . -cp $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale - ### au prealable !!! -####################################################################### -if [[ "${real}" == "0" ]] -then - cp $LOADL_STEP_INITDIR/input_* . - cp $LOADL_STEP_INITDIR/ze_hill . ## topo - cp $LOADL_STEP_INITDIR/levels . ## custom vert levels -else - ln -sf $LOADL_STEP_INITDIR/wrfb* . -fi -####################################################################### -ideal.exe ### marche pas en parallele, pb avec execution ensuite -####################################################################### Index: unk/MESOSCALE/LMD_MM_MARS/SIMU/mpd.hosts =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/mpd.hosts (revision 211) +++ (revision ) @@ -1,7 +1,0 @@ -auric.lmd.jussieu.fr -smith.lmd.jussieu.fr - -auric.lmd.jussieu.fr -smith.lmd.jussieu.fr -widor.lmd.jussieu.fr -parker.lmd.jussieu.fr Index: unk/MESOSCALE/LMD_MM_MARS/SIMU/save_this_simu =================================================================== --- /trunk/MESOSCALE/LMD_MM_MARS/SIMU/save_this_simu (revision 211) +++ (revision ) @@ -1,48 +1,0 @@ -#! /bin/bash - -\rm simtime 2> /dev/null - -mkdir results - -if [[ -f rsl.error.0000 ]] -then - - timesec=$(more rsl.error.0000 | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}') - touch simtime - echo ${timesec} seconds - echo ${timesec} seconds >> simtime - - head -10000 rsl.out.0000 > red_rsl.out.0000 - head -4000 rsl.error.0000 > red_rsl.error.0000 - cp met_em/rsl.out.0000 real_rsl.out.0000 - mv red_rsl.* results/ - mv real_rsl.out.0000 results/ - mv simtime results/ - -else - - timesec=$(more log_wrf | grep elapsed | grep 'Timing for main' | awk '{s+=$9} END{print s}') - touch simtime - echo ${timesec} seconds - echo ${timesec} seconds >> simtime - - head -10000 log_wrf > red_log_wrf - mv red_log_wrf results/ - mv simtime results/ - mv met_em/log_real results/ - -fi - -tar czvf def.tar.gz namelist.input namelist.wps callphys*.def makemeso* #launch_simu makemeso* -cp def.tar.gz results/ - -mv *.def results/ -mv results/callphys*def ./ - -mv wrfout* results/ -ls -lh results -du -bh results -echo '*****************************************************' -echo '*** BEWARE ! *** output files were moved in the folder named results !!!' -echo 'what you may keep: results folder and def.tar.gz file' -echo '*****************************************************'