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 345) +++ (revision ) @@ -1,101 +1,0 @@ - 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/GWPOLAR/ALT_POLAR_top0.1_levels =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/GWPOLAR/ALT_POLAR_top0.1_levels (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/GWPOLAR/ALT_POLAR_top0.1_levels (revision 346) @@ -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/HOLDEN_last_GEOGRID.TBL =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN/HOLDEN_last_GEOGRID.TBL (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN/HOLDEN_last_GEOGRID.TBL (revision 346) @@ -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/HOLDEN_last_levels =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN/HOLDEN_last_levels (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN/HOLDEN_last_levels (revision 346) @@ -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/HOLDEN_last_GEOGRID.TBL =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/HOLDEN_last_GEOGRID.TBL (revision 345) +++ (revision ) @@ -1,240 +1,0 @@ -# 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 345) +++ (revision ) @@ -1,61 +1,0 @@ - 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/MERIDIANI_vincent/LS53/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/callphys.def (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/callphys.def (revision 346) @@ -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/MERIDIANI_vincent/LS53/makemeso_input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/makemeso_input (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/makemeso_input (revision 346) @@ -0,0 +1,7 @@ +3 +4 +201 +201 +101 +1 +1 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/namelist.input (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/namelist.input (revision 346) @@ -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/MERIDIANI_vincent/LS53/namelist.wps =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/namelist.wps (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/MERIDIANI/MERIDIANI_vincent/LS53/namelist.wps (revision 346) @@ -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/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input (revision 345) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/NEWPHYS_TRACERS/namelist.input (revision 346) @@ -20,5 +20,5 @@ &domains - time_step = 20 !! Dynamical timestep + time_step = 15 !! Dynamical timestep dx = 30000, !! (p2) Horizontal resolution dy = 30000, !! (p2) Horizontal resolution (should be equal to dx) @@ -36,19 +36,9 @@ mars = 11 num_soil_layers = 18 !!! NOUVELLE PHYSIQUE - !! - !! OPTIONAL - !! / &dynamics tiso = 100. - pd_scalar = T !!! tres IMPORTANT, sinon valeurs negatives diff_6th_opt = 0 - !!! - !!! 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 / Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.input (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.input (revision 346) @@ -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/OLYMPUS/OM_namelist.input_before_improv =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.input_before_improv (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.input_before_improv (revision 346) @@ -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/OLYMPUS/OM_namelist.wps =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.wps (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/OLYMPUS/OM_namelist.wps (revision 346) @@ -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/OM_namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/OM_namelist.input (revision 345) +++ (revision ) @@ -1,59 +1,0 @@ - &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 345) +++ (revision ) @@ -1,55 +1,0 @@ - &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 345) +++ (revision ) @@ -1,36 +1,0 @@ -&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/yeyeye/callphys.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphys.def (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphys.def (revision 346) @@ -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/RES/yeyeye/callphysd2_.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphysd2_.def (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphysd2_.def (revision 346) @@ -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/RES/yeyeye/callphysd3_.def =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphysd3_.def (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/callphysd3_.def (revision 346) @@ -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/RES/yeyeye/makemeso_input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/makemeso_input (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/makemeso_input (revision 346) @@ -0,0 +1,6 @@ +1 +153 +153 +61 +3 +1 Index: trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/namelist.input =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/namelist.input (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/namelist.input (revision 346) @@ -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/RES/yeyeye/namelist.wps =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/namelist.wps (revision 346) +++ trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/RES/yeyeye/namelist.wps (revision 346) @@ -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' +/