- Timestamp:
- 11/28/11 16:25:13 (13 years ago)
- Location:
- CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES
- Files:
-
- 2 added
- 2 deleted
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/COMP/inca.card
r1554 r1601 36 36 37 37 [OutputFiles] 38 List= (inca_avgr.nc, ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc, NONE ) 39 38 List= (inca_avgr.nc, ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc, Post_1M_inca_avgr ) 40 39 41 40 [Post_1M_inca_avgr] -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/COMP/lmdz.card
r1554 r1601 5 5 [UserChoices] 6 6 LMDZ_NbPeriod_adjust=3 7 # LMDZ_NbPeriod_adjust=0 8 # To force usage of this Bands file, put LMDZ_NbPeriod_adjust=0 and replace XXXXXXX by Restart Job Name 9 # LMDZ_Bands_file_name=${ARCHIVE}/IGCM_OUT/IPSLCM5/XXXXXXX/ATM/Debug/XXXXXXX_Bands_96x95x19_3prc.dat_3 10 #LMDZ_Bands_file_name=/dmnfs11/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_NA_AR5WERF_co2/ATM/Debug/LOI_NA_AR5WERF_co2_20060501_20060531_Bands_96x95x19_4prc.dat 7 #LMDZ_NbPeriod_adjust=0 8 # To force usage of this Bands file, put LMDZ_NbPeriod_adjust=0 and replace XXXXXXX by Name/Resolution/NB Prco Job Restart 9 # LMDZ_Bands_file_name=${ARCHIVE}/IGCM_OUT/LMDZORINCA/XXXXXXX/ATM/Debug/XXXXXXX_Bands_XXXXXXX_XXXprc.dat_3 10 11 # Other specific parameters to this experience 12 aer_type=scenario 13 pmagic=0.0 14 15 # Nugding : ok_guide=y/n 16 ok_guide=n 11 17 12 18 … … 18 24 19 25 [BoundaryFiles] 20 List= (/dmnfs/cont003/p86cozic/INPUT/39LEVELS/limit_365.nc , limit.nc) 26 List= (${R_BC}/CHM/LMDZORINCA/LMD9695_leap/EL-LMD9695_${year}1231_limit.nc , limit.nc ) 27 21 28 ListNonDel= () 22 29 23 30 [ParametersFiles] 24 List= (${SUBMIT_DIR}/PARAM/physiq.def , .), \25 (${SUBMIT_DIR}/PARAM/gcm.def , gcm.def), \26 (${SUBMIT_DIR}/PARAM/run.def , .),\27 (${SUBMIT_DIR}/PARAM/ traceur.def, .)31 List= (${SUBMIT_DIR}/PARAM/physiq.def , physiq.def ), \ 32 (${SUBMIT_DIR}/PARAM/gcm.def , gcm.def ), \ 33 (${SUBMIT_DIR}/PARAM/run.def , . ),\ 34 (${SUBMIT_DIR}/PARAM/guide.def , . ) 28 35 29 36 [RestartFiles] … … 32 39 33 40 [OutputText] 34 List= (physiq.def, gcm.def, run.def, traceur.def, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)41 List= (physiq.def, gcm.def, run.def, guide.def, traceur.def, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat) 35 42 36 43 37 44 [OutputFiles] 38 List= (histmth.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc, NONE)45 List= (histmth.nc, ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc, Post_1M_histmth ) 39 46 40 47 [Post_1M_histmth] -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/COMP/lmdz.driver
r1555 r1601 3 3 4 4 #----------------------------------------------------------------- 5 function LMDZ_sed 6 { 7 IGCM_debug_PushStack "LMDZ_sed" 8 9 sed -e "s/^${2}\ *=.*/${2}= ${3}/" ${1} > ${1}.tmp 10 RET=$? 11 echo "LMDZ_sed : ${1} ${2} ${3}" 12 \mv ${1}.tmp ${1} 13 14 IGCM_debug_PopStack "LMDZ_sed" 15 return $RET 16 } 17 18 #----------------------------------------------------------------- 19 5 20 function ATM_Initialize 6 21 { … … 31 46 fi 32 47 33 ##--Frequency purpose .... 48 49 ##-- Output frequency purpose .... 34 50 ##-- Initialisation .... 35 OK_instan=n 36 ##-- 37 OK_journe=n 38 OK_mensuel=y 51 ok_mensuel=y 52 ok_journe=n 53 ok_hf=n 54 ok_hf3h=n 55 ok_hf3hm=n 56 ok_stn=n 39 57 40 case ${config_UserChoices_PeriodLength} in41 1Y|1y|1M|1m) OK_mensuel=y ;;42 5D|5d|1D|1d) OK_journe=n ;;43 esac44 45 for frequency in ${config_ATM_WriteFrequency} ; do46 case ${frequency} in47 5D|5d|1D|1d) OK_journe=n ;;48 esac49 done50 51 ##--Variables used by LMDZ --52 PAT_INST=$( grep 'OK_instan' ${SUBMIT_DIR}/PARAM/physiq.def )53 PAT_JOUR=$( grep 'OK_journe' ${SUBMIT_DIR}/PARAM/physiq.def )54 PAT_MOIS=$( grep 'OK_mensuel' ${SUBMIT_DIR}/PARAM/physiq.def )55 56 PAT_iphysiq=$( grep 'iphysiq' ${SUBMIT_DIR}/PARAM/gcm.def )57 PAT_iperiod=$( grep 'iperiod' ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)58 PAT_day_step=$( grep 'day_step' ${SUBMIT_DIR}/PARAM/gcm.def )59 PAT_ecritphy=$( grep 'ecritphy' ${SUBMIT_DIR}/PARAM/gcm.def )60 58 61 59 … … 70 68 71 69 72 ##-- This could be define in lmdz.card, inside section [UserChoices]73 ##-- Otherwise we get the value in *.def74 iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)75 iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')76 # day_step : number of steps per day (multiple of iperiod)77 (( day_step = 48 * iphysiq ))78 79 70 IGCM_debug_PopStack "ATM_Initialize" 80 71 } … … 86 77 87 78 ecritphy=${PeriodLengthInDays} 79 LMDZ_ecrit_ISCCP=${PeriodLengthInDays}. 80 LMDZ_periodav=${PeriodLengthInDays}. 88 81 89 82 ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def … … 93 86 RAZ_DATE=0 94 87 fi 88 89 90 ## algo : For CumulPeriod=1 ; 91 ## IF NbPeriod_Adjust = 0, Get Bands file directly from server, store it with _0 suffix 92 ## ELSE LMDZ_adjust=y ; save Bands file on file server ; 93 ## For CumulPeriod=2 to LMDZ_NbPeriod_adjust ; Get Bands file from server (CumulPeriod-1) ; LMDZ_adjust=y ; save Bands file on file server ; 94 ## For CumulPeriod > LMDZ_NbPeriod_adjust ; LMDZ_adjust=n ; Get Bands file from server (CumulPeriod=LMDZ_NbPeriod_adjust 0 to force it) ; Use Bands ; 95 95 96 96 if [ ${CumulPeriod} -eq 1 ] ; then … … 111 111 [ ${CumulPeriod} -gt 1 ] && ( IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat; IGCM_sys_Chmod u+w Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ) 112 112 113 ## Mise en forme du fichier physiq.def114 sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/" \115 -e "s/OK_journe=.*/OK_journe=${OK_journe}/" \116 -e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" \117 physiq.def > physiq.def.tmp118 IGCM_sys_Mv physiq.def.tmp physiq.def119 113 120 ## Mise en forme du fichier gcm.def 121 sed -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" \ 122 -e "s/${PAT_day_step}/day_step=${day_step}/" \ 123 -e "s/${PAT_iperiod}/iperiod=${iperiod}/" \ 124 -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" \ 125 gcm.def > gcm.def.tmp 126 IGCM_sys_Mv gcm.def.tmp gcm.def 114 ## physiq.def parameters modified from initial physiq.def 115 # columns in phys_output_filekeys refer to the filenames given in phys_out_filenames in output.def : 116 # histmth histday histhf histfh3h histhf3hm histstn 117 LMDZ_sed physiq.def phys_out_filekeys "${ok_mensuel} ${ok_journe} ${ok_hf} ${ok_hf3h} ${ok_hf3hm} ${ok_stn}" 118 LMDZ_sed physiq.def ecrit_ISCCP ${LMDZ_ecrit_ISCCP} 127 119 128 greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 129 jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat ) 130 nbjour=$( expr \( $jul_dat \% 1000 \) ) 131 yractu=$( expr $greg_dat / 10000 ) 132 echo "dayref = " $nbjour "year actu = " $yractu 120 if [ ! X${lmdz_UserChoices_ok_guide} = X ] ; then 121 LMDZ_sed guide.def ok_guide ${lmdz_UserChoices_ok_guide} 122 fi 133 123 134 ## Mise en forme du fichier run.def 135 sed -e "s/_dayref_/${nbjour}/" \ 136 -e "s/_anneeref_/${yractu}/" \ 137 -e "s/_calend_/${CalendarTypeForLmdz}/" \ 138 -e "s/_nday_/${PeriodLengthInDays}/" \ 139 -e "s/_raz_date_/${RAZ_DATE}/" \ 140 -e "s/_adjust_/${LMDZ_adjust}/" \ 141 run.def > run.def.tmp 142 IGCM_sys_Mv run.def.tmp run.def 124 ## run.def parameters 125 LMDZ_sed run.def dayref ${InitDay} 126 LMDZ_sed run.def anneeref ${InitYear} 127 LMDZ_sed run.def calend ${CalendarTypeForLmdz} 128 LMDZ_sed run.def nday ${PeriodLengthInDays} 129 LMDZ_sed run.def raz_date ${RAZ_DATE} 130 LMDZ_sed run.def periodav ${LMDZ_periodav} 131 LMDZ_sed run.def adjust ${LMDZ_adjust} 132 133 ## Other specific parameters to this experience 134 if [ ! X${lmdz_UserChoices_aer_type} = X ] ; then 135 LMDZ_sed physiq.def aer_type ${lmdz_UserChoices_aer_type} 136 fi 137 if [ ! X${lmdz_UserChoices_pmagic} = X ] ; then 138 LMDZ_sed physiq.def pmagic ${lmdz_UserChoices_pmagic} 139 fi 143 140 144 141 IGCM_debug_PopStack "ATM_Update" -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/COMP/orchidee.card
r1073 r1601 4 4 [UserChoices] 5 5 NEWHYDROL=n 6 sechiba_LEVEL=10 6 sechiba_LEVEL=11 7 VEGET_UPDATE=1Y 8 LAND_COVER_CHANGE=y 7 9 8 10 [InitialStateFiles] 9 List= (${R_INIT}/SRF/LMDZORINCA/soils_param.nc, .), \10 (${R_INIT}/SRF/LMDZORINCA/routing.nc, .)11 List= (${R_INIT}/SRF/LMDZORINCA/soils_param.nc, . ), \ 12 (${R_INIT}/SRF/LMDZORINCA/routing.nc, . ) 11 13 12 14 [BoundaryFiles] 13 List= () 14 ListNonDel= (${R_BC}/SRF/LMDZORINCA/PFTmap.20C3M.nc, PFTmap.nc) 15 List= (${BC_SRF_PFTmap}, PFTmap.nc) 16 ListNonDel= () 17 #Tous les 5 ans : 18 #ListNonDel= (${BC_SRF_PFTmap}, PFTmap.nc) 19 20 [SmoothFiles] 21 List=() 15 22 16 23 [ParametersFiles] … … 22 29 23 30 [OutputText] 24 List= (out_orchidee )31 List= (out_orchidee, orchidee.def) 25 32 26 33 [OutputFiles] 27 List= (sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history) \28 (sechiba_out_2.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_out2.nc, NONE)\29 ( watchout.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_watchout.nc,NONE)34 List= (sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history), \ 35 (sechiba_out_2.nc, ${R_OUT_SRF_O_H}/${PREFIX}_HF_sechiba_out_2.nc, Post_HF_sechiba_out_2), \ 36 (orchidee_watchout.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_watchout.nc, NONE) 30 37 31 38 [Post_1M_sechiba_history] 32 Patches = () 33 GatherWithInternal = (lon, lat, veget, time_counter, time_counter_bnds, Areas) 34 TimeSeriesVars2D = (nobiofrac, alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, snownobio, snowf, subli, tair, temp_sol, tsol_max, tsol_min, drainage) 35 ChunckJob2D = NONE 36 TimeSeriesVars3D = (lai, maxvegetfrac, vegetfrac, CO2FLUX) 37 ChunckJob3D = NONE 39 Patches= () 40 GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac) 41 TimeSeriesVars2D= (nobiofrac, alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, snownobio, snowf, subli, tair, temp_sol, tsol_max, tsol_min, drainage,mrsos, mrso, mrros, mrro, prveg, evspsblveg, evspsblsoi, tran, treeFrac, grassFrac, cropFrac, baresoilFrac, residualFrac) 42 ChunckJob2D= NONE 43 TimeSeriesVars3D= (lai, maxvegetfrac, vegetfrac, CO2FLUX, ptn, nee) 44 ChunckJob3D= NONE 45 Seasonal= ON 46 47 [Post_HF_sechiba_out_2] 48 Patches= () 49 GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac) 50 TimeSeriesVars2D= (mrsos, mrro) 51 ChunckJob2D= 50Y 52 TimeSeriesVars3D= (ptn) 53 ChunckJob3D= 10Y 54 Seasonal= OFF -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/COMP/orchidee.driver
r858 r1601 24 24 RESOL_SRF=ALL 25 25 26 for frequency in ${config_SRF_WriteFrequency} ; do 27 case ${frequency} in 28 HF|hf) SRF_ok_hf=y ;; 29 esac 30 done 31 26 32 IGCM_debug_PopStack "SRF_Initialize" 33 } 34 35 #----------------------------------------------------------------- 36 function SRF_PeriodStart 37 { 38 39 IGCM_debug_PushStack "SRF_PeriodStart" 40 if [ ${year} -gt 2004 ] ; then 41 42 PFTmap="${R_BC}/SRF/LMDZORINCA/ORCHIDEE_PFTmap_2006to2100_AR5_RCP8.5_MESSAGE.nc" 43 if [ ${year} -eq 2005 ] ; then 44 IGCM_sys_Rm -f PFTmap.nc 45 fi 46 47 else 48 PFTmap="${R_BC}/SRF/LMDZORINCA/PFTmap_1850to2005_AR5_LUHa.rc2.nc" 49 fi 50 51 BC_SRF_PFTmap=${PFTmap} 52 53 IGCM_debug_PopStack "SRF_PeriodStart" 54 27 55 } 28 56 … … 34 62 typeset SECHIBA_WRITE_STEP 35 63 36 case ${config_SRF_WriteFrequency} in 64 SRF_WriteFrequency=$( echo ${config_SRF_WriteFrequency} | sed -e 's/\([0-9]*[yYmMdDs]\).*/\1/' ) 65 case ${SRF_WriteFrequency} in 37 66 *Y|*y) 38 WriteInYears=$( echo ${ 1} | awk -F '[yY]' "{print ${config_SRF_WriteFrequency}}")39 PeriodLengthInYears=$( echo ${ 1} | awk -F '[yY]' "{print ${config_UserChoices_PeriodLength}}")67 WriteInYears=$( echo ${SRF_WriteFrequency} | awk -F '[yY]' '{print $1}' ) 68 PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' ) 40 69 (( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInYears / PeriodLengthInYears * 86400 )) ;; 41 *M|*m) 42 WriteInMonths=$( echo ${1} | awk -F '[mM]' "{print ${config_SRF_WriteFrequency}}" ) 70 1M) 43 71 case ${config_UserChoices_PeriodLength} in 44 72 *Y|*y) 45 PeriodLengthInYears=$( echo ${1} | awk -F '[yY]' "{print ${config_UserChoices_PeriodLength}}" ) 73 SECHIBA_WRITE_STEP=-1. 74 ;; 75 *M|*m) 76 SECHIBA_WRITE_STEP=-1. 77 ;; 78 *) 79 (( SECHIBA_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) * 86400 )) 80 ;; 81 esac 82 ;; 83 *M|*m) 84 WriteInMonths=$( echo ${SRF_WriteFrequency} | awk -F '[mM]' '{print $1}' ) 85 case ${config_UserChoices_PeriodLength} in 86 *Y|*y) 87 PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' ) 46 88 (( SECHIBA_WRITE_STEP = PeriodLengthInDays * 86400 / PeriodLengthInYears / 12 )) 47 89 ;; 48 90 *M|*m) 49 PeriodLengthInMonths=$( echo ${ 1} | awk -F '[mM]' "{print ${config_UserChoices_PeriodLength}}")91 PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' ) 50 92 (( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInMonths * 86400 / PeriodLengthInMonths )) 51 93 ;; … … 60 102 (( SECHIBA_WRITE_STEP = 86400 )) ;; 61 103 *s) 62 WriteInSeconds=$( echo ${ 1} | awk -F '[s]' "{print ${config_SRF_WriteFrequency}}")104 WriteInSeconds=$( echo ${SRF_WriteFrequency} | awk -F '[s]' '{print $1}' ) 63 105 (( SECHIBA_WRITE_STEP = WriteInSeconds )) ;; 64 106 *) 65 IGCM_debug_Exit "SRF_Update " ${ config_SRF_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D."107 IGCM_debug_Exit "SRF_Update " ${SRF_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D." 66 108 IGCM_debug_Verif_Exit ;; 67 109 esac 68 110 111 # Dans le cas de test : 112 # if [ ${year} -eq 1979 ] ; then 113 # ORCHIDEE_sed VEGET_REINIT y 114 # ORCHIDEE_sed VEGET_YEAR 131 115 # fi 116 # 117 # # Dans le cas historique : 118 # if [ ${year} -eq 2005 ] ; then 119 # ORCHIDEE_sed VEGET_REINIT y 120 # ORCHIDEE_sed VEGET_YEAR 1 121 ## ORCHIDEE_sed VEGET_UPDATE '0Y' 122 # fi 123 124 # # Dans le cas sans changement de carte aprÚs 2005 125 # if [ ${year} -gt 2005 ] ; then 126 # ORCHIDEE_sed VEGET_UPDATE '0Y' 127 # fi 128 69 129 ORCHIDEE_sed HYDROL_CWRR ${orchidee_UserChoices_NEWHYDROL} 70 130 … … 72 132 ORCHIDEE_sed SECHIBA_HISTLEVEL ${orchidee_UserChoices_sechiba_LEVEL} 73 133 74 ## if [ ${year} -eq 1950 ] ; then 75 # ORCHIDEE_sed SECHIBA_HISTFILE2 y 76 # ORCHIDEE_sed SECHIBA_HISTLEVEL2 10 77 # ORCHIDEE_sed WRITE_STEP2 86400.0 78 ## fi 134 # Outputs HF in HISTFILE2 if required 135 [ X${SRF_ok_hf} = Xy ] && ORCHIDEE_sed SECHIBA_HISTFILE2 y 136 ORCHIDEE_sed SECHIBA_HISTLEVEL2 1 137 ORCHIDEE_sed WRITE_STEP2 10800.0 79 138 80 139 if ( [ ${CumulPeriod} -ne 1 ] || [ "${config_SRF_Restart}" != "n" ] ) ; then 81 140 ORCHIDEE_sed SECHIBA_restart_in sechiba_rest_in.nc 141 fi 142 143 # Desactaivate river routing for high resolution or very low resolution 144 if ( [ ${RESOL_ATM_3D} = "280x280x39-ZAf" ] || [ ${RESOL_ATM_3D} = "360x180x19-ZIn" ] || [ ${RESOL_ATM_3D} = "32x24x11" ] ) ; then 145 ORCHIDEE_sed RIVER_ROUTING n 146 fi 147 148 # Parameters set in orchidee.card 149 if [ ! X${orchidee_UserChoices_VEGET_UPDATE} = X ] ; then 150 ORCHIDEE_sed VEGET_UPDATE ${orchidee_UserChoices_VEGET_UPDATE} 151 fi 152 if [ ! X${orchidee_UserChoices_LAND_COVER_CHANGE} = X ] ; then 153 ORCHIDEE_sed LAND_COVER_CHANGE ${orchidee_UserChoices_LAND_COVER_CHANGE} 82 154 fi 83 155 -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/PARAM/gcm.def
r1553 r1601 4 4 ## periode pour le pas Matsuno (en pas) 5 5 iperiod=5 6 ## periode de la dissipation (en pas)7 idissip=108 6 ## choix de l'operateur de dissipation (star ou non star ) 9 7 lstardis=y … … 24 22 ## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 25 23 purmats=n 26 ## guide ou pas 27 ok_guide=.false. 28 ## avec ou sans physique 29 ##physic=n 30 iflag_phys=1 31 ## periode de la physique (en pas) 24 ## avec ou sans physique 25 physic=y 26 ## periode de la physique (en pas) 32 27 iphysiq=10 33 28 ## Avec ou sans strato 34 ok_strato=y 29 ok_strato=y 35 30 # Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche 36 31 iflag_top_bound=2 37 32 # Coefficient pour la couche eponge (valeur derniere couche) 38 33 tau_top_bound=5.e-5 39 ## frequence (en jours ) de l'ecriture du fichier histphy 40 ecritphy=1 41 ## frequence (en jours ) de l'ecriture du fichier histmth 42 ecrit_mth=1 43 ## Cycle diurne ou non 44 cycle_diurne=y 45 ## Soil Model ou non 46 soil_model=y 47 ## Choix ou non de New oliq 48 new_oliq=y 49 ## Orodr ou non pour l orographie 50 ok_orodr=y 51 ## Orolf ou non pour l orographie 52 ok_orolf=y 53 ## Si = .T. , lecture du fichier limit avec la bonne annee 54 ok_limitvrai=n 55 ## longitude en degres du centre du zoom 34 ## longitude en degres du centre du zoom 56 35 clon=0. 57 36 ## latitude en degres du centre du zoom -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/PARAM/orchidee.def
r858 r1601 248 248 STOMATE_HISTLEVEL = 10 249 249 # default = 10 250 251 #-------------------------------------------------------------------- 252 # STOMATE_IPCC_OUTPUT_FILE 253 # This file is going to be created by the model 254 # and will contain the output from the model. 255 # This file is a truly COADS compliant netCDF file. 256 # It will be generated by the hist software from 257 # the IOIPSL package. 258 # Name of file in which STOMATE's output is going 259 # to be written 260 STOMATE_IPCC_OUTPUT_FILE = stomate_ipcc_history.nc 261 # default = stomate_ipcc_history.nc 262 263 # STOMATE_IPCC_HIST_DT 264 # Time step of the STOMATE IPCC history file 265 # STOMATE IPCC history time step (d) 266 STOMATE_IPCC_HIST_DT = -1 267 # default = 0. 250 268 251 269 # Write frequency for output files (SECHIBA in seconds et … … 572 590 573 591 # Total depth of soil reservoir 574 HYDROL_SOIL_DEPTH = 2.592 HYDROL_SOIL_DEPTH = 4. 575 593 # default = 2. 594 595 # Root profile 596 # Default values were defined for 2 meters soil depth. 597 # For 4 meters soil depth, you may use those ones : 598 # 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1. 599 HYDROL_HUMCSTE = 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1. 600 # default = 5., .8, .8, 1., .8, .8, 1., 1., .8, 4., 4., 4., 4. 576 601 577 602 # Initial restart deep soil moisture if not found in restart … … 682 707 # Add from Nathalie - the 28 of March 2006 - advice from Fred Hourdin 683 708 # By PFT. 684 RVEG_PFT = 1., 0.5, 0.5, 1., 1., 1., 1., 1., 1., 1., 0.5, 1., 0.5709 RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 685 710 # default = 1. 686 711 … … 703 728 # Year of the land_use vegetation map readed 704 729 # year off the pft map 705 # If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc )706 VEGET_YEAR = 1 51730 # If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc, 1 for PFTmap_IPCC_2000.nc) 731 VEGET_YEAR = 1 707 732 # default = 282 733 734 # booleen to indicate that a new LAND USE file will be used (since 1.9.5 version). 735 # The parameter is used to bypass veget_year count 736 # and reinitialize it with VEGET_YEAR parameter. 737 # Then it is possible to change LAND USE file. 738 # If LAND_USE 739 VEGET_REINIT = n 740 # default = n 708 741 709 742 # Update vegetation frequency (since 2.0 version) 710 743 # The veget datas will be update each this time step. 711 744 # If LAND_USE 712 VEGET_UPDATE = 0Y745 VEGET_UPDATE = XY 713 746 # default = 1Y 714 747 … … 718 751 # deforestation. 719 752 # If LAND_USE 720 LAND_COVER_CHANGE = n753 LAND_COVER_CHANGE = x 721 754 # default = y 722 755 -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/PARAM/physiq.def
r1553 r1601 1 ##- $Id: physiq.def 71 2007-06-21 14:26:03Z acosce $ 2 type_ocean=force 3 VEGET=y 4 OK_journe=n 5 OK_mensuel=y 6 OK_instan=n 7 if_ebil=0 8 # 9 # parametres KE 10 # 11 epmax = .999 12 ok_adj_ema = n 13 iflag_clw = 1 1 ## $Id$ 2 # 3 # PARAMETRES ANCIENNEMENT DANS gcm.def 4 ## Cycle diurne ou non 5 cycle_diurne=y 6 ## Soil Model ou non 7 soil_model=y 8 ## Choix ou non de New oliq 9 new_oliq=y 14 10 ## Activation ou non de la parametrisation de Hines pour la strato 15 11 ok_hines=y 12 ## Orodr ou non pour l orographie 13 ok_orodr=y 14 ## Orolf ou non pour l orographie 15 ok_orolf=y 16 ## Si = .T. , lecture du fichier limit avec la bonne annee 17 ok_limitvrai=n 16 18 ## Nombre d'appels des routines de rayonnements ( par jour) 17 19 nbapp_rad=12 18 ## Flag pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)20 ## Flag pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect 19 21 iflag_con=30 22 ## Facteur multiplication des precip convectives dans KE 23 cvl_corr=1.0 24 ## Facteur additif pour l'albedo 25 # variable modife avec lmdz.card 26 pmagic=999.9 27 # 28 # 29 # Parametres fichiers de sortie 30 # 31 # Noms des fichiers 32 phys_out_filenames= histmth histday histhf histhf3h histhf3hm histstn 33 # Sortir ou non les fichiers 34 phys_out_filekeys= y y y y y n 35 # Niveaux de sorties 36 phys_out_filelevels= 5 5 5 5 5 5 37 ### Type de fichier : global (n) ou stations (y) 38 phys_out_filestations = n n n n n y 39 ### Frequences des sorties 40 phys_out_filetimesteps = 1.day, 1.day, 0.25day, 0.125day, 0.125day, 1800.s 41 42 # more variables for level 5 : vitw and pres 43 # Vertical wind 44 flag_vitw = 2, 3, 5, 5, 1 45 name_vitw = vitw 46 # Air pressure 47 flag_pres = 2, 3, 5, 5, 1 48 name_pres = pres 49 # 50 ### parametres pour hist*NMC.nc 51 #lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier 52 # histdayNMC: 8 ou 17 53 lev_histdayNMC=8 54 ### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence 55 ok_histNMC = y, n, n 56 57 ### freqin_isccp = frequence input en secondes du simulateur ISCCP 58 freq_ISCCP=10800. 59 ### freqout_isccp = frequence output en jours du simulateur ISCCP 60 ecrit_ISCCP=_ecrit-ISCCP_ 61 ### niveau du diagnostique de conservation d energie 62 if_ebil=0 63 # 64 # parametres KE 65 # 66 ### epmax = Efficacite precipitation maximale 67 epmax = .999 68 ### ok_adj_ema = ?? pas utilise 69 ok_adj_ema = n 70 ### iflag_clw Flag calcul eau liquide 71 # - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l) 72 # - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l) 73 # - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip 74 iflag_clw = 1 20 75 # 21 76 # parametres nuages 22 77 # 78 ### cld_lc_lsc contenu en eau liquide des nuages large-scale (fisrtilp) 23 79 cld_lc_lsc = 4.16e-4 80 ### cld_lc_con contenu en eau liquide des nuages convectifs (fisrtilp) 24 81 cld_lc_con = 4.16e-4 82 ### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 25 83 cld_tau_lsc = 1800. 84 ### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 26 85 cld_tau_con = 1800. 86 ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 27 87 ffallv_lsc = 0.5 88 ### ffallv_lsc cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 28 89 ffallv_con = 0.5 29 coef_eva = 0.00002 90 ### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?) 91 coef_eva = 2.e-5 92 ### reevap_ice reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 30 93 reevap_ice = y 94 ### iflag_cldcon flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 95 # - iflag_cldcon<=-1 diag. rain_Tiedtke 96 # - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k) 97 # - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc) 98 # - iflag_cldcon=3, ratqs=ratqss 31 99 iflag_cldcon = 3 100 ### iflag_pdf : flag calcul distribution sous-maille de l eau et des nuages 101 # - iflag_pdf=0, version ratqs, 102 # - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir 103 # - des PDFS de Sandrine Bony 32 104 iflag_pdf = 1 105 ### fact_cldcon constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 33 106 fact_cldcon = 1. 34 # facttemps = 1.e-4107 ### facttemps= facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3) 35 108 facttemps = 0.0001 109 ## ok_newmicro =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 36 110 ok_newmicro = y 111 ### iflag_ratqs=0 correspond a la version IPCC AR4 112 iflag_ratqs=0 113 ### ratqsbas ratqs en bas si iflag_cldcon=1 37 114 ratqsbas = 0.005 115 ### ratqshaut ratqs en haut pour ratqss "stables" 38 116 ratqshaut = 0.33 117 ### rad_froid rayon cristaux des nuages de glace "froids" 39 118 rad_froid = 35 119 ### rad_chau1 rayon goutelettes d eau chauds", en haut: k=4-klev 40 120 rad_chau1=12 121 ### rad_chau2 rayon goutelettes d eau chauds", en bas: k=1-3 41 122 rad_chau2=11 42 #ksta_ter=1.e-7 43 ksta_ter=0.0000001 44 # 45 # parametres outputs 46 # 47 #niveau de sortie "hf" lev_histhf avec 48 # - lev_histhf=3 => defaut 49 # - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele 50 lev_histhf=3 51 #niveau de sortie "day" lev_histday 52 # - lev_histday=2 => defaut 53 # - lev_histday=3 => + champs 3D => F. Lott 54 # - lev_histday=4 => + champs sous-surfaces 55 lev_histday=2 56 #niveau de sortie "mth" lev_histmth avec 57 # - lev_histmth=2 => defaut 58 # - lev_histmth=3 => albedo, rugosite sous-surfaces 59 # - lev_histmth=4 => champs tendances 3d 60 lev_histmth=4 61 # 62 # parametres climatique 63 # 123 # 124 # Coefficient et parametres sur les drags 125 # 126 f_cdrag_ter=0.8 127 f_cdrag_stable=1. 128 f_cdrag_oce=0.8 129 f_rugoro=0. 130 ### ksta_ter coef.diffusion minimale sur terre/sic/lic 131 ksta_ter=1.e-7 132 ### cdmmax = cdrag maximum pour le moment 133 cdmmax = 2.5E-3 134 ### cdhmax = cdrag maximum pour l energie 135 cdhmax = 2.0E-3 136 # 137 # Parametres "orbitaux/ ere geologique" 138 # 139 ### R_ecc = Excentricite 64 140 R_ecc = 0.016715 141 ### R_peri = Equinoxe 65 142 R_peri = 102.7 143 ### R_incl = Inclinaison 66 144 R_incl = 23.441 145 ### solaire = Constante solaire 67 146 solaire = 1365. 147 # 148 # Taux gaz a effet de serre 149 # 150 ### co2_ppm = taux CO2 en ppm 68 151 co2_ppm = 348. 152 ### CH4_ppb = taux CH4 en ppb 69 153 CH4_ppb = 1650. 154 ### N2O_ppb = taux N2O en ppb 70 155 N2O_ppb = 306. 156 ### CFC11_ppt = taux CFC11 en ppt 71 157 CFC11_ppt = 280. 72 CFC12_ppt = 484. 73 # 74 # parametres simulateur ISCCP 75 # 158 ### CFC12_ppt = taux CFC12 en ppt 159 CFC12_ppt = 484 160 # 161 # Parametres effets directs/indirects des "aerosols" 162 # 163 ### ok_ade=y/n flag Aerosol direct effect 164 ok_ade=n 165 ### ok_aie=y/n flag Aerosol indirect effect 166 ok_aie=n 167 ### aer_type = Aerosol variation type : actuel / preind / scenario / annuel 168 # variable modife avec lmdz.card 169 aer_type=xxx 170 ### type of coupled aerosol =1 (default) =2 => bc only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol 171 flag_aerosol=1 172 ### bl95_b0 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) 173 bl95_b0=2 174 ### bl95_b1 = Parameter in CDNC-maer link (Boucher&Lohmann 1995) 175 bl95_b1=0.2 176 # 177 # Parametre de lecture de l'ozone 178 # 179 # Allowed values are 0, 1 and 2 180 # 0: do not read an ozone climatology 181 # 1: read a single ozone climatology that will be used day and night 182 # 2: read two ozone climatologies, the average day and night climatology and the daylight climatology 183 read_climoz=0 184 # 185 # Parametres simulateur COSP (CFMIP Observational Simulator Package) 186 # 187 ### ok_cosp=y/n flag simulateur COSP 188 ok_cosp=n 189 ### freq_COSP = frequence d'appel de COSP en secondes 190 freq_COSP=10800. 191 ### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc 192 ok_mensuelCOSP=y 193 ### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc 194 ok_journeCOSP=y 195 ### ok_hfCOSP=y/n, ecrire sorties haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc 196 ok_hfCOSP=n 197 # 198 # Parametres simulateur ISCCP 199 # 200 ### ok_isccp=y/n flag simulateur ISCCP 201 ok_isccp=n 202 ### top_height = flag choix calcul nuages par le simulateur en utilisant 203 # - les donnees IR et/ou VIS et l algorithme ISCCP-D1 204 # - top_height = 1 -> algo IR-VIS 205 # - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)" 206 # - top_height = 3 -> algo IR 76 207 top_height = 3 77 #overlap = 1, 2 ou 3 208 ### overlap = Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP 209 # - overlap=1 Max overlap 210 # - overlap=2 Random overlap 211 # - overlap=3 Max/Random overlap 78 212 overlap = 3 79 #cdmmax80 #cdmmax = 2.5E-381 cdmmax = 0.002582 #cdhmax83 #cdhmax = 2.0E-384 cdhmax = 0.00285 #autres flags86 iflag_pbl = 187 iflag_thermals = 088 iflag_ratqs=089 iflag_clos=090 #91 # parametre chimie/aerosols92 #93 213 config_inca=chem 94 214 aerosol_couple=n 95 flag_aerosol=1 96 new_aod=.true. 97 ok_ade=n 98 ok_aie=n 99 215 ## frequence (en jours ) de l'ecriture du fichier histphy 216 ecritphy=1 217 ## frequence (en jours ) de l'ecriture du fichier histmth 218 ecrit_mth=1 -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/PARAM/run.def
r1553 r1601 4 4 INCLUDEDEF=orchidee.def 5 5 INCLUDEDEF=guide.def 6 ## activation du calcul d equilibrage de charge7 adjust=_adjust_8 6 ## Type de calendrier utilise 9 7 ## valeur possible: earth_360d (defaut), earth_365d, earth_366d … … 13 11 ## Annee de l'etat initial ( avec 4 chiffres ) 14 12 anneeref=_anneeref_ 15 ## raz_date16 raz_date=_raz_date_17 13 ## Nombre de jours d'integration 18 14 nday=_nday_ 15 ## Remise a zero de la date initiale 16 raz_date=_raz_date_ 19 17 ## periode de sortie des variables de controle (en pas) 20 iconser= 18 iconser=480 21 19 ## periode d'ecriture du fichier histoire (en jour) 22 iecri=10 23 ## periode de stockage fichier histmoy (en jour) 24 periodav=1. 25 ## unite de sortie des impressions 26 lunout=6 20 iecri=1 21 ## flag de sortie dynzon 22 ok_dynzon=n 23 ## periode de stockage fichier dynzon (en jour) 24 periodav=_periodav_ 25 ## activation du calcul d equilibrage de charge 26 adjust=_adjust_ 27 ## activation du filtre fft 28 use_filtre_fft=y 29 ## niveau d impression de controle 30 prt_level=0 31 ### 32 ### Informations sur la configuration utilisee 33 ### 34 ### type_ocean = force / slab /couple 35 type_ocean=force 36 ### VEGET= y si ORCHIDEE, =n si bucket 37 VEGET=y 38 ### Parameters read only by ce0l : 39 ### create file grilles_gcm.nc 40 grilles_gcm_netcdf=y 41 -
CONFIG/LMDZORINCA/branches/LMDZORINCA3_39/EXP_GES/config.card
r1553 r1601 1 1 #- $Id: config.card 69 2007-06-21 14:08:22Z acosce $ 2 2 #======================================================================== 3 #D-- Carte pour un run GES en resolution 96x 723 #D-- Carte pour un run GES en resolution 96x95x39 4 4 #======================================================================== 5 5 #D-- Compatibility - … … 12 12 #D-- (<8 chars MAX for JobName) 13 13 JobName=LOIG3_39_00 14 LongName=' LMDz-ORCHIDEE-INCA with New scripts'14 LongName='Ajout de Stomate (NEW) - Fichiers sflx / aircraft / start / startphy créés par S.Szopa ' 15 15 TagName=LMDZORINCA/GES 16 16 #============================ … … 40 40 #D- For each component, Name of component, Tag of component 41 41 ATM= (lmdz, LMDZ4) 42 SRF= (orchidee, orchidee_1_9_2) 42 SRF= (orchidee, orchidee_1_9_5) 43 SBG= (stomate, ORCHIDEE_1_9_5) 43 44 CHM= (inca, INCA3) 44 45 … … 50 51 ATM= (gcm.e, gcm.e) 51 52 SRF= ("", "") 53 SBG= ("", "") 52 54 CHM= (inca.dat, inca.dat) 53 55 … … 58 60 OverRule=n 59 61 #D- Last day of the experience used as restart 60 RestartDate= 1869-12-3062 RestartDate=2005-12-31 61 63 #D- Define restart simulation name (=> JOB_OS) 62 RestartJobName= JobName64 RestartJobName=LOI_AR5WERF_E0_2006 63 65 #D- Path Server Group Login (=> PSGL) 64 RestartPath= ${ARCHIVE}/IGCM_OUT/LMDZORINCA/GES66 RestartPath=/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER 65 67 # Attention login depend de la machine 66 68 … … 72 74 RebuildFrequency=1D 73 75 #D- Do we rebuild parallel output from archive 74 RebuildFromArchive= none76 RebuildFromArchive=NONE 75 77 #D- If you want to produce time series, this flag determines 76 78 #D- frequency of post-processing submission … … 108 110 OldName= 109 111 112 #======================================================================== 113 #D-- SBG - STOMATE 114 [SBG] 115 WriteFrequency="1M" 116 # If config_Restarts_OverRule == 'n' next 4 params are read 117 Restart=n 118 #-- Last day of the experience used as restart for this component if Restart=y 119 RestartDate=2000-01-31 120 # Define restart simulation name for this component 121 RestartJobName=EXP00 122 RestartPath=${ARCHIVE}/IGCM_OUT/LMDZOR/DEVT/amip 123 #-- Old component name for restart (if empty, use new name) 124 OldName= 110 125 111 126 #========================================================================
Note: See TracChangeset
for help on using the changeset viewer.