- Timestamp:
- 2009-09-21T11:28:25+02:00 (15 years ago)
- Location:
- branches/libIGCM/ORCA2_LIM_PISCES/IGCM00
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/COMP/opa9.card
r1396 r1631 4 4 [UserChoices] 5 5 OPA_NDT_DAY=15 6 #============================ 7 #--v3, v3_1 6 #============================================================== 7 #-- Name of directory for COMMON ACCOUNT at idris or CCRT 8 #-- available directories are: v3, v3_1 9 #-- related to NEMO_tag_revision (nemo_v3, nemo_v3_1) 10 #-- NOTE: if you use your own directory you have to comment it 8 11 OPA_version=v3_1 9 #============================ 12 #============================================================== 13 14 #-- if you run interannual is "y", if not is "n" 15 Interannual_Run=n 16 17 [Interannual] 18 #===================================================================== 19 #-- Surface Boundary Condition original files name ( get by the job ) 20 #===================================================================== 21 #-- NOTE: MANDATORY STRUCTURE OF FILE NAMES GET BY THE JOB : 22 #-- "basename_yyyy.nc" ( for interannual data ) 23 #-- "basename.nc" ( for climatologycal data ) 24 #-- MANDATORY: NAMES USED IN List_jobsbc VARIABLE MUST BE ONLY THE BASENAME 25 26 List_jobsbc=(taux_1m, tauy_1m, flx) 27 28 #========================================================================================== 29 #-- Surface Boundary Condition files name expected by NEMO ( same as the ones in namelist ) 30 #========================================================================================== 31 #-- NOTE: 2 possible cases : 32 #-- 1) put in List_runsbc file names different from List_jobsbc but identical to the ones in namelist 33 #-- 2) leave List_runsbc empty so NEMO will read file names specified in List_jobsbc variable above 34 35 List_runsbc=() 10 36 11 37 [InitialStateFiles] … … 21 47 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_potential_temperature_nomask.nc, .), \ 22 48 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_salinity_nomask.nc, .), \ 23 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/flx.nc, .),\24 49 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/geothermal_heating.nc, .),\ 25 50 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/runoff_1m_nomask.nc, .),\ 26 51 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/sss_data.nc, .),\ 27 52 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/sst_data.nc, .),\ 28 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/taux_1m.nc, .),\29 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/tauy_1m.nc, .),\30 53 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/subbasins.nc, .) 31 54 … … 44 67 (${PREFIX_NWRITE}_${DATE_OPA}_grid_V.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_V.nc, Post_1M_grid_V),\ 45 68 (${PREFIX_NWRITE}_${DATE_OPA}_grid_W.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_W.nc, Post_1M_grid_W),\ 46 (${PREFIX_NWRITE}_${DATE_OPA}_diaptr.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_diaptr.nc, Post_1M_diaptr),\47 (${PREFIX_NWRITE}_${DATE_OPA}_trends.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_trends.nc, NONE),\48 (${PREFIX_NWRITE}_${DATE_OPA}_diagap.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_diagap.nc, NONE),\49 69 (damping.coeff.nc , ${R_OUT_OCE_NWRITE}/${PREFIX}_damping.coeff.nc, NONE),\ 50 70 (mesh_mask.nc , ${R_OUT_OCE_O}/${config_UserChoices_JobName}_mesh_mask.nc, NONE),\ 51 71 (output.abort.nc , ${R_OUT_OCE_D}/${PREFIX}_output.abort.nc, NONE),\ 52 (output.init.nc , ${R_OUT_OCE_O_I}/${config_UserChoices_JobName}_${PeriodDateBegin}_output.init.nc, NONE),\ 72 (output.init.nc , ${R_OUT_OCE_O_I}/${config_UserChoices_JobName}_${PeriodDateBegin}_output.init.nc, NONE) 73 53 74 54 75 [Post_1M_grid_T] … … 71 92 GatherWithInternal = (nav_lon, nav_lat, depthw, time_counter) 72 93 TimeSeriesVars = 73 74 [Post_1M_diaptr]75 Patches = ()76 GatherWithInternal = (lat, deptht, depthw, time_counter)77 TimeSeriesVars = (zotemglo, zosalglo, zomsfglo, zotematl, zosalatl, zomsfatl, zotempac, zosalpac, zomsfpac, zotemind, zosalind, zomsfind, zotemipc, zosalipc, zomsfipc, sohtatl, sostatl, sohtpac, sostpac, sohtind, sostind, sohtipc, sostipc, sophtadv, sophtove, sophtldf, sopstadv, sopstove, zomsfeiv, sophteiv, sopsteiv) -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/COMP/opa9.driver
r1379 r1631 7 7 JOB_NAME=${config_UserChoices_JobName} 8 8 9 if [ -z "${opa_UserChoices_OPA_NDT_DAY}" ] ; then 10 OPA_NDT_DAY=15 11 else 12 OPA_NDT_DAY=${opa_UserChoices_OPA_NDT_DAY} 13 fi 9 OPA_NDT_DAY=${opa_UserChoices_OPA_NDT_DAY:=15} 10 14 11 15 12 ##--Variables used by OPA -- 16 13 17 # c experexperience name for vairmer format18 # n it000 number of the first time step19 # n itend number of the last time step20 # n leapy leap year calendar (0/1) (30 for 360d)21 # n write frequency of OUTPUT file14 # cn_exp experience name for vairmer format 15 # nn_it000 number of the first time step 16 # nn_itend number of the last time step 17 # nn_leapy leap year calendar (0/1) (30 for 360d) 18 # nn_write frequency of OUTPUT file 22 19 # ln_rstart boolean term for restart (true or false) 23 # nstock frequency of restart file 24 # nrstdt control of the time step (0, 1 or 2) 25 # ndate0 initial calendar date aammjj 26 # nmsh =1 create a mesh file (coordinates, scale factors, masks) 20 # nn_stock frequency of restart file 21 # nn_rstctl control of the time step (0, 1 or 2) 22 # nn_date0 initial calendar date aammjj 23 # nn_msh =1 create a mesh file (coordinates, scale factors, masks) 24 # rn_rdt time step in seconds (coming from namelist) 25 # nf_ptr_wri frequency of zonal means and transport output 26 27 27 28 28 # Local function to find namelists parameters … … 32 32 33 33 34 PAT_CEXPER=$( supergrep c exper${SUBMIT_DIR}/PARAM/namelist )35 PAT_NIT000=$( supergrep n it000 ${SUBMIT_DIR}/PARAM/namelist )36 PAT_NITEND=$( supergrep n itend ${SUBMIT_DIR}/PARAM/namelist )37 PAT_NLEAPY=$( supergrep n leapy ${SUBMIT_DIR}/PARAM/namelist )38 PAT_NWRITE=$( supergrep n write ${SUBMIT_DIR}/PARAM/namelist )34 PAT_CEXPER=$( supergrep cn_exp ${SUBMIT_DIR}/PARAM/namelist ) 35 PAT_NIT000=$( supergrep nn_it000 ${SUBMIT_DIR}/PARAM/namelist ) 36 PAT_NITEND=$( supergrep nn_itend ${SUBMIT_DIR}/PARAM/namelist ) 37 PAT_NLEAPY=$( supergrep nn_leapy ${SUBMIT_DIR}/PARAM/namelist ) 38 PAT_NWRITE=$( supergrep nn_write ${SUBMIT_DIR}/PARAM/namelist ) 39 39 PAT_RESTAR=$( supergrep ln_rstart ${SUBMIT_DIR}/PARAM/namelist ) 40 PAT_NSTOCK=$( supergrep nstock ${SUBMIT_DIR}/PARAM/namelist ) 41 PAT_NRSTAR=$( supergrep nrstdt ${SUBMIT_DIR}/PARAM/namelist ) 42 PAT_NDATE0=$( supergrep ndate0 ${SUBMIT_DIR}/PARAM/namelist ) 43 PAT_NMSH=$( supergrep nmsh ${SUBMIT_DIR}/PARAM/namelist ) 44 45 OPA_RDT=$( supergrep rdt ${SUBMIT_DIR}/PARAM/namelist | sed 's/ *rdt *=//' | sed 's/\. *,//' ) 40 PAT_NSTOCK=$( supergrep nn_stock ${SUBMIT_DIR}/PARAM/namelist ) 41 PAT_NRSTAR=$( supergrep nn_rstctl ${SUBMIT_DIR}/PARAM/namelist ) 42 PAT_NDATE0=$( supergrep nn_date0 ${SUBMIT_DIR}/PARAM/namelist ) 43 PAT_NMSH=$( supergrep nn_msh ${SUBMIT_DIR}/PARAM/namelist ) 44 PAT_NF_PTR=$( supergrep nf_ptr_wri ${SUBMIT_DIR}/PARAM/namelist ) 45 46 47 OPA_RDT=$( supergrep rn_rdt ${SUBMIT_DIR}/PARAM/namelist | sed 's/ *rn_rdt *=//' | sed 's/\. *,//' ) 46 48 47 49 # Period Length In Days between DateBegin and first day of calendar 0001 01 01 … … 49 51 (( DaysSinceJC = $( IGCM_date_DaysSinceJC ${DateBegin} ) + 1 )) 50 52 53 # Definition from opa9.card of List_jobsbc and List_runsbc used to do to interannual and climatological runs 54 IGCM_card_DefineArrayFromOption ${SUBMIT_DIR}/COMP/opa9.card Interannual List_jobsbc 55 IGCM_card_DefineArrayFromOption ${SUBMIT_DIR}/COMP/opa9.card Interannual List_runsbc 56 set -A ListFormulationJobsbc -- \${opa9_Interannual_Listjobsc${opa9_Interannual_List_jobsbc}[*]} 57 set -A ListFormulationRunsbc -- \${opa9_Interannual_Listrunsbc${opa9_Interannual_List_runsbc}[*]} 58 51 59 IGCM_debug_PopStack "OCE_Initialize" 52 60 } … … 55 63 function OCE_Update 56 64 { 57 ### set -vx58 65 IGCM_debug_PushStack "OCE_Update" 66 67 # Interannual run 68 # In this case job need to calculate previous and following year and give, and if every run needs to add "_y" before every year 69 if [ X${opa9_UserChoices_Interannual_Run} = Xy ] ; then 70 71 typeset file fileo 72 73 eval NbFileInter=${#opa9_Interannual_List_jobsbc[*]} 74 75 (( i = 0 )) 76 while [ $i -lt ${NbFileInter} ] ; do 77 eval file=${opa9_Interannual_List_jobsbc[$i]} 78 eval jobsbc_file_ym1=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_$(( year - 1 )).nc 79 eval jobsbc_file_y=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_${year}.nc 80 eval jobsbc_file_yp1=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_$(( year + 1 )).nc 81 # Copy of List_jobsbc in List_runsbc if this last one is empty 82 if [ X${opa9_Interannual_List_runsbc[0]} = X${NULL_STR} ] ; then 83 eval fileo=${opa9_Interannual_List_jobsbc[$i]} 84 eval runsbc_file_ym1=${file}_y$(( year - 1 )).nc 85 eval runsbc_file_y=${file}_y${year}.nc 86 eval runsbc_file_yp1=${file}_y$(( year + 1 )).nc 87 else 88 eval fileo=${opa9_Interannual_List_runsbc[$i]} 89 eval runsbc_file_ym1=${fileo}_y$(( year - 1 )).nc 90 eval runsbc_file_y=${fileo}_y${year}.nc 91 eval runsbc_file_yp1=${fileo}_y$(( year + 1 )).nc 92 fi 93 94 if [ X${Period} = X1 ] ; then 95 if [ ${month} -eq 01 ] ; then 96 IGCM_sys_Get ${jobsbc_file_ym1} ${runsbc_file_ym1} 97 fi 98 if [ ! -f ${runsbc_file_y} ] ; then 99 IGCM_sys_Get ${jobsbc_file_y} ${runsbc_file_y} 100 fi 101 if [ ! -f ${runsbc_file_yp1} ] ; then 102 IGCM_sys_Get ${jobsbc_file_yp1} ${runsbc_file_yp1} 103 fi 104 else 105 if [ ! -f ${runsbc_file_yp1} ] ; then 106 IGCM_sys_Get ${jobsbc_file_yp1} ${runsbc_file_yp1} 107 fi 108 fi 109 110 (( i = i + 1 )) 111 done 112 # End interannual 113 elif [ X${opa9_UserChoices_Interannual_Run} = Xn ] ; then 114 # Climatological run 115 typeset file fileo 116 117 eval NbFileInter=${#opa9_Interannual_List_jobsbc[*]} 118 119 (( i = 0 )) 120 while [ $i -lt ${NbFileInter} ] ; do 121 eval file=${opa9_Interannual_List_jobsbc[$i]} 122 eval jobsbc_file=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}.nc 123 if [ X${opa9_Interannual_List_runsbc[0]} = X${NULL_STR} ] ; then 124 eval fileo=${opa9_Interannual_List_jobsbc[$i]} 125 eval runsbc_file=${fileo}.nc 126 else 127 eval fileo=${opa9_Interannual_List_runsbc[$i]} 128 eval runsbc_file=${fileo}.nc 129 fi 130 131 IGCM_sys_Get ${jobsbc_file} ${runsbc_file} 132 (( i = i + 1 )) 133 done 134 fi 135 # End climatological 59 136 60 137 NbFreq=$( echo ${config_OCE_WriteFrequency} | wc -w ) … … 98 175 OPA_NSTOCK="${OPA_NITEND}" 99 176 100 ## Verification of number of time step in aday177 ## Verification of number of time steps per day 101 178 (( NB_SEC_DAY_MODEL = OPA_NDT_DAY * OPA_RDT )) 102 179 (( NB_SEC_DAY = 60 * 60 * 24 )) … … 104 181 if [ ${NB_SEC_DAY_MODEL} -ne ${NB_SEC_DAY} ] 105 182 then 106 echo " VERIF IER OPA_NDT_DAY dans le job ${JOB}"183 echo " VERIFY OPA_NDT_DAY in opa9.card " 107 184 exit 108 185 fi … … 111 188 if ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_OCE_Restart}" = "n" ] ) ; then 112 189 113 #echo " PAS DE RESTART OPA"190 #echo "NO OPA RESTART" 114 191 OPA_LRSTAR=.FALSE. 115 192 OPA_NRSTDT=0 116 # echo pas de meshmask en parallele193 #Put OPA_NMSH=0 when OPA runnig in parallel mode 117 194 OPA_NMSH=1 118 195 ( [ X${BATCH_NUM_PROC_TOT} != X ] && [ "${BATCH_NUM_PROC_TOT}" -gt 1 ] ) && OPA_NMSH=0 … … 120 197 elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_OCE_Restart}" = "y" ] ) ; then 121 198 122 #echo " RESTART OPA"199 #echo "OPA RESTART" 123 200 OPA_LRSTAR=.TRUE. 124 201 OPA_NRSTDT=1 … … 127 204 else 128 205 129 #echo " RESTART OPA"206 #echo "OPA RESTART" 130 207 OPA_LRSTAR=.TRUE. 131 208 OPA_NRSTDT=2 … … 133 210 134 211 fi 135 212 213 ## nleapy configuration 136 214 case ${config_UserChoices_CalendarType} in 137 215 leap) … … 147 225 NEMO_END=$( echo $( awk "BEGIN { printf \"%0${PRECIS}d\",${OPA_NITEND} }" ) ) 148 226 149 sed -e "s%${PAT_CEXPER}% cexper=\"${config_UserChoices_JobName}\"%" \ 150 -e "s%${PAT_NIT000}% nit000=${OPA_NIT000}%" \ 151 -e "s%${PAT_NITEND}% nitend=${OPA_NITEND}%" \ 152 -e "s%${PAT_NLEAPY}% nleapy=${OPA_NLEAPY}%" \ 153 -e "s%${PAT_NWRITE}% nwrite=${OPA_NWRITE}%" \ 227 sed -e "s%${PAT_CEXPER}% cn_exp=\"${config_UserChoices_JobName}\"%" \ 228 -e "s%${PAT_NIT000}% nn_it000=${OPA_NIT000}%" \ 229 -e "s%${PAT_NITEND}% nn_itend=${OPA_NITEND}%" \ 230 -e "s%${PAT_NWRITE}% nn_write=${OPA_NWRITE}%" \ 154 231 -e "s%${PAT_RESTAR}% ln_rstart=${OPA_LRSTAR}%" \ 155 -e "s%${PAT_NSTOCK}% nstock=${OPA_NSTOCK}%" \ 156 -e "s%${PAT_NRSTAR}% nrstdt=${OPA_NRSTDT}%" \ 157 -e "s%${PAT_NDATE0}% ndate0=${PeriodDateBegin}%" \ 158 -e "s%${PAT_NMSH}% nmsh=${OPA_NMSH}%" \ 232 -e "s%${PAT_NSTOCK}% nn_stock=${OPA_NSTOCK}%" \ 233 -e "s%${PAT_NRSTAR}% nn_rstctl=${OPA_NRSTDT}%" \ 234 -e "s%${PAT_NDATE0}% nn_date0=${PeriodDateBegin}%" \ 235 -e "s%${PAT_NMSH}% nn_msh=${OPA_NMSH}%" \ 236 -e "s%${PAT_NLEAPY}% nn_leapy=${OPA_NLEAPY}%" \ 237 -e "s%${PAT_NF_PTR}% nf_ptr_wri=${OPA_NWRITE}%" \ 159 238 namelist > namelist.tmp 160 239 … … 173 252 174 253 if [ -f date.file ] ; then 175 # Prefix use in opa.card AND in lim .card :254 # Prefix use in opa.card AND in lim2.card : 176 255 DATE_OPA=$( cat date.file | \ 177 256 sed "s/\ ${config_UserChoices_JobName}_[0-9]*[a-z]_\([0-9]*_[0-9]*\)_\ */\1/g" ) 178 ###DATE_OPA=$( $DATE_OPABF | cut -c1-30 )257 ###DATE_OPA=$( $DATE_OPABF | cut -c1-30 ) 179 258 MainPrefix=${config_UserChoices_JobName}_1d_${DATE_OPA} 180 259 SecondPrefix=${config_UserChoices_JobName}_5d_${DATE_OPA} -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/COMP/pisces.card
r1393 r1631 9 9 [BoundaryFiles] 10 10 List= () 11 ListNonDel= (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/kRGB61.txt, .), \ 11 ListNonDel= (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_DIC_nomask.nc, LEVITUS_DIC.nc), \ 12 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_Alkalini_nomask.nc, LEVITUS_Alkalini.nc), \ 13 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_O2_nomask.nc, LEVITUS_O2.nc), \ 14 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_NO3_nomask.nc, LEVITUS_NO3.nc), \ 15 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_PO4_nomask.nc, LEVITUS_PO4.nc), \ 16 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_Si_nomask.nc, LEVITUS_Si.nc), \ 17 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_DOC_nomask.nc, LEVITUS_DOC.nc), \ 18 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/data_1m_Fer_nomask.nc, LEVITUS_Fer.nc), \ 19 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/kRGB61.txt, .), \ 12 20 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/dust.orca.nc, .), \ 13 21 (${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/bathy.orca.nc, .), \ … … 23 31 24 32 [OutputText] 25 List= ( )33 List= (namelist_top, namelist_pisces) 26 34 27 35 [OutputFiles] 28 List= (${PREFIX_NWRITETRC}_${DATE_OPA}_ptrc_T.nc, ${R_OUT_ TOP_NWRITE}/${PREFIX}_${WFT1}_ptrc_T.nc, NONE), \29 (${PREFIX_NWRITETRC}_${DATE_OPA}_diad_T.nc, ${R_OUT_ TOP_NWRITE}/${PREFIX}_${WFT1}_diad_T.nc, NONE)36 List= (${PREFIX_NWRITETRC}_${DATE_OPA}_ptrc_T.nc, ${R_OUT_MBG_NWRITE}/${PREFIX}_${WFT1}_ptrc_T.nc, Post_1M_ptrc_T), \ 37 (${PREFIX_NWRITETRC}_${DATE_OPA}_diad_T.nc, ${R_OUT_MBG_NWRITE}/${PREFIX}_${WFT1}_diad_T.nc, Post_1M_diad_T) 30 38 39 [Post_1M_ptrc_T] 40 Patches= () 41 GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter) 42 TimeSeriesVars = (DCHL, NCHL, NO3, PO4, O2, Si, Fer, DIC, Alkalini, DOC) 43 44 [Post_1M_diad_T] 45 Patches = () 46 GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter) 47 TimeSeriesVars = (PMO, PMO2, Delc, Cflx, PPPHY, PPPHY2) 48 -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/COMP/pisces.driver
r1392 r1631 2 2 3 3 #----------------------------------- 4 function TOP_Initialize4 function MBG_Initialize 5 5 { 6 IGCM_debug_PushStack " TOP_Initialize"6 IGCM_debug_PushStack "MBG_Initialize" 7 7 8 ##--Variables used by TOP--8 ##--Variables used by MBG -- 9 9 10 10 # nwritetrc frequency of write in the tracer output file … … 24 24 25 25 PAT_TOP_NWRITETRC=$( supergrep nwritetrc ${SUBMIT_DIR}/PARAM/namelist_top ) 26 PAT_TOP_LRSTTR=$( supergrep l rsttr${SUBMIT_DIR}/PARAM/namelist_top )26 PAT_TOP_LRSTTR=$( supergrep ln_rsttr ${SUBMIT_DIR}/PARAM/namelist_top ) 27 27 PAT_TOP_NRSTTR=$( supergrep nrsttr ${SUBMIT_DIR}/PARAM/namelist_top ) 28 28 PAT_TOP_NWRITETRD=$( supergrep ntrd_trc ${SUBMIT_DIR}/PARAM/namelist_top ) … … 31 31 PAT_PIS_NWRITEDIA=$( supergrep nwritedia ${SUBMIT_DIR}/PARAM/namelist_pisces ) 32 32 33 IGCM_debug_PopStack " TOP_Initialize"33 IGCM_debug_PopStack "MBG_Initialize" 34 34 } 35 35 36 36 #----------------------------------- 37 function TOP_Update37 function MBG_Update 38 38 { 39 IGCM_debug_PushStack " TOP_Update"39 IGCM_debug_PushStack "MBG_Update" 40 40 41 41 ##--Write Frequency Purpose .... 42 topfrequency=$( echo ${config_ TOP_WriteFrequency} | awk "-F " '{print $1}' )42 topfrequency=$( echo ${config_MBG_WriteFrequency} | awk "-F " '{print $1}' ) 43 43 topfactor=$( echo ${topfrequency} | sed -e "s/[yYmMdD]//" ) 44 44 case ${topfrequency} in … … 46 46 (( TOP_NWRITE = OPA_NDT_DAY * topfactor * $( IGCM_date_DaysInYear ${year} ) )) ; 47 47 PREFIX_NWRITETRC=${config_UserChoices_JobName}_${topfactor}y ; 48 R_OUT_ TOP_NWRITE=${R_OUT_TOP_O_Y} ;48 R_OUT_MBG_NWRITE=${R_OUT_MBG_O_Y} ; 49 49 WFT1=${topfactor}Y ;; 50 50 1M|1m) 51 51 (( TOP_NWRITE = OPA_NDT_DAY * topfactor * $( IGCM_date_DaysInMonth ${year} ${month} ) )) ; 52 52 PREFIX_NWRITETRC=${config_UserChoices_JobName}_${topfactor}m ; 53 R_OUT_ TOP_NWRITE=${R_OUT_TOP_O_M} ;53 R_OUT_MBG_NWRITE=${R_OUT_MBG_O_M} ; 54 54 WFT1=${topfactor}M ;; 55 55 *D|*d) 56 56 (( TOP_NWRITE = OPA_NDT_DAY * topfactor )) ; 57 57 PREFIX_NWRITETRC=${config_UserChoices_JobName}_${topfactor}d ; 58 R_OUT_ TOP_NWRITE=${R_OUT_TOP_O_D} ;58 R_OUT_MBG_NWRITE=${R_OUT_MBG_O_D} ; 59 59 WFT1=${topfactor}D ;; 60 60 *) … … 62 62 63 63 ##-- Restart configuration 64 if ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_ TOP_Restart}" = "n" ] ) ; then64 if ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_MBG_Restart}" = "n" ] ) ; then 65 65 66 66 #echo "NO RESTART FOR TOP" … … 68 68 TOP_NRSTTR=0 69 69 70 elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_ TOP_Restart}" = "y" ] ) ; then70 elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_MBG_Restart}" = "y" ] ) ; then 71 71 72 72 #echo "RESTART TOP" … … 83 83 sed -e "s%${PAT_TOP_NWRITETRC}% nwritetrc=${TOP_NWRITE}%" \ 84 84 -e "s%${PAT_TOP_NWRITETRD}% ntrd_trc=${TOP_NWRITE}%" \ 85 -e "s%${PAT_TOP_LRSTTR}% l rsttr=${TOP_LRSTTR}%" \85 -e "s%${PAT_TOP_LRSTTR}% ln_rsttr=${TOP_LRSTTR}%" \ 86 86 -e "s%${PAT_TOP_NRSTTR}% nrsttr=${TOP_NRSTTR}%" \ 87 87 namelist_top > namelist_top.tmp … … 99 99 grep AUTO namelist* 100 100 101 IGCM_debug_PopStack " TOP_Update"101 IGCM_debug_PopStack "MBG_Update" 102 102 } 103 103 104 104 #----------------------------------- 105 function TOP_Finalize105 function MBG_Finalize 106 106 { 107 IGCM_debug_PushStack " TOP_Finalize"107 IGCM_debug_PushStack "MBG_Finalize" 108 108 109 echo FINALIZE TOP!!!109 echo FINALIZE MBG !!! 110 110 111 IGCM_debug_PopStack " TOP_Finalize"111 IGCM_debug_PopStack "MBG_Finalize" 112 112 } 113 113 -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/PARAM/namelist
r1379 r1631 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (nam _zgr, nam_zgr_sco, namdom)3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, nam alb)6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nam tide)5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 7 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam _traadv, nam_traldf, namtdp)9 !! 7 - dynamics (nam _dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)10 !! 8 - Verical physics (namzdf, nam npc, namric, namtke, namkpp, namddm)8 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_dmp) 9 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 10 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 11 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol, nam _mpp, nam_mpp_dyndist, namctl)12 !! 9 - miscellaneous (namsol, nammpp, nammpp_dyndist, namctl) 13 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. … … 23 23 &namrun ! parameters of the run 24 24 !----------------------------------------------------------------------- 25 no = 0 ! job number 26 cexper = "ORCA2" ! experience name 25 nn_no = 0 ! job number 26 cn_exp = "ORCA2P" ! experience name 27 nn_it000 = 1 ! first time step 28 nn_itend = 315 ! last time step (std 5475) 29 nn_date0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) 30 nn_leapy = 0 ! Leap year calendar (1) or not (0) 31 nn_istate = 0 ! output the initial state (1) or not (0) 32 nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) 33 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nit000) 34 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 35 ln_mskland = .true. ! mask land points in NetCDF outputs (costly: + ~15%) 36 ln_clobber = .false. ! clobber (overwrite) an existing file 37 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (working only with iom_nf90 routines) 38 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 39 nn_rstctl = 0 ! restart control = 0 nit000 is not compared to the restart file value 40 ! = 1 use ndate0 in namelist (not the value in the restart file) 41 ! = 2 calendar parameters read in the restart file 27 42 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 28 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 29 ln_rstart = .false. ! start from rest (F) or from a restart file (T)30 nrstdt = 0 ! restart control = 0 nit000 is not compared to the restart file value31 ! = 1 use ndate0 in namelist (not the value in the restart file)32 ! = 2 calendar parameters read in the restart file33 nit000 = 1 ! first time step34 nitend = 5475 ! last time step35 ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1)36 nleapy = 0 ! Leap year calendar (1) or not (0)37 ninist = 1 ! output the initial state (1) or not (0)38 nstock = 5475 ! frequency of creation of a restart file (modulo referenced to 1)39 nwrite = 5475 ! frequency of write in the output file (modulo referenced to nit000)40 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)41 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%)42 44 / 43 45 !!====================================================================== 44 46 !! *** Domain namelists *** 45 47 !!====================================================================== 46 !! nam _zgr vertical coordinate47 !! nam _zgr_sco s-coordinate or hybrid z-s-coordinate48 !! namdom 49 !!====================================================================== 50 51 !----------------------------------------------------------------------- 52 &nam _zgr! vertical coordinate48 !! namzgr vertical coordinate 49 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 50 !! namdom space and time domain (bathymetry, mesh, timestep) 51 !!====================================================================== 52 53 !----------------------------------------------------------------------- 54 &namzgr ! vertical coordinate 53 55 !----------------------------------------------------------------------- 54 56 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) … … 57 59 / 58 60 !----------------------------------------------------------------------- 59 &nam_zgr_sco ! s-coordinate or hybrid z-s-coordinate 60 !----------------------------------------------------------------------- 61 sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 62 sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 63 theta = 6.0 ! surface control parameter (0<=theta<=20) 64 thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 65 r_max = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 61 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 62 !----------------------------------------------------------------------- 63 rn_sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 64 rn_sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 65 rn_theta = 6.0 ! surface control parameter (0<=theta<=20) 66 rn_thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 67 rn_rmax = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 68 ln_s_sigma = .false. ! hybrid s-sigma coordinates 69 rn_bb = 0.8 ! stretching with s-sigma 70 rn_hc = 150.0 ! critical depth with s-sigma 66 71 / 67 72 !----------------------------------------------------------------------- 68 73 &namdom ! space and time domain (bathymetry, mesh, timestep) 69 74 !----------------------------------------------------------------------- 70 ntopo = 1 ! compute (=0) or read(=1) the bathymetry file 71 e3zps_min = 25. ! the thickness of the partial step is set larger than the minimum 72 e3zps_rat = 0.2 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 73 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 74 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 75 ! =0, no acceleration, rdt = rdttra 76 atfp = 0.1 ! asselin time filter parameter 77 rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) 78 rdtmin = 5760. ! minimum time step on tracers (used if nacc=1) 79 rdtmax = 5760. ! maximum time step on tracers (used if nacc=1) 80 rdth = 800. ! depth variation of tracer time step (used if nacc=1) 75 nn_bathy = 1 ! compute (=0) or read(=1) the bathymetry file 76 nn_closea = 0 ! closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 77 nn_msh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 78 rn_e3zps_min= 20. ! the thickness of the partial step is set larger than the minimum 79 rn_e3zps_rat= 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 80 ! 81 rn_rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) ==> 5760 81 82 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 82 nclosea = 0 ! = 0 no closed sea in the model domain 83 ! = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 83 rn_atfp = 0.1 ! asselin time filter parameter 84 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) 85 ! =0, not used, rdt = rdttra 86 rn_rdtmin = 28800. ! minimum time step on tracers (used if nacc=1) 87 rn_rdtmax = 28800. ! maximum time step on tracers (used if nacc=1) 88 rn_rdth = 800. ! depth variation of tracer time step (used if nacc=1) 84 89 / 85 90 !!====================================================================== … … 92 97 !! namsbc_core CORE bulk formulea formulation 93 98 !! namsbc_cpl CouPLed formulation ("key_coupled") 94 !! nam qsrpenetrative solar radiation99 !! namtra_qsr penetrative solar radiation 95 100 !! namsbc_rnf river runoffs 96 101 !! namsbc_ssr sea surface restoring term (for T and/or S) 97 !! nam albalbedo parameters102 !! namsbc_alb albedo parameters 98 103 !!====================================================================== 99 104 … … 117 122 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 118 123 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 119 nn_fwb = 0 ! FreshWater Budget: =0 unchecked ,120 ! =1 global mean of e-p-r set to zero at each nn_fsbc time step ,124 nn_fwb = 3 ! FreshWater Budget: =0 unchecked 125 ! =1 global mean of e-p-r set to zero at each time step 121 126 ! =2 annual global mean of e-p-r set to zero 127 ! =3 global emp set to zero and spread out over erp area 122 128 / 123 129 !----------------------------------------------------------------------- … … 134 140 &namsbc_flx ! surface boundary condition : flux formulation 135 141 !----------------------------------------------------------------------- 136 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !137 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !138 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. ,'yearly' , '' , ''139 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. ,'yearly' , '' , ''140 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. ,'yearly' , '' , ''141 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. ,'yearly' , '' , ''142 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. ,'yearly' , '' , ''142 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 143 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 144 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' , '' , '' 145 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' , '' , '' 146 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' , '' , '' 147 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' , '' , '' 148 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' , '' , '' 143 149 ! 144 150 cn_dir = './' ! root directory for the location of the flux files … … 147 153 &namsbc_clio ! namsbc_clio CLIO bulk formulea 148 154 !----------------------------------------------------------------------- 149 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !150 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !151 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. ,'yearly' , '' , ''152 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. ,'yearly' , '' , ''153 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. ,'yearly' , '' , ''154 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. ,'yearly' , '' , ''155 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. ,'yearly' , '' , ''156 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. ,'yearly' , '' , ''157 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. ,'yearly' , '' , ''155 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 156 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 157 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. , 'yearly' , '' , '' 158 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. , 'yearly' , '' , '' 159 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. , 'yearly' , '' , '' 160 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. , 'yearly' , '' , '' 161 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. , 'yearly' , '' , '' 162 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' , '' , '' 163 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' , '' , '' 158 164 ! 159 165 cn_dir = './' ! root directory for the location of the bulk files are … … 162 168 &namsbc_core ! namsbc_core CORE bulk formulea 163 169 !----------------------------------------------------------------------- 164 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! weights! rotation !165 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename! pairing !166 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'U1'167 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'V1'168 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''169 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''170 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''171 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''172 sn_prec = 'precip_core' , -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''173 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''170 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 171 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 172 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'U1' 173 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'V1' 174 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 175 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 176 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 177 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 178 sn_prec = 'precip_core', -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 179 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 174 180 ! 175 181 cn_dir = './' ! root directory for the location of the bulk files 176 182 ln_2m = .true. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 177 alpha_precip= 1. ! multiplicative factor for precipitation (total & snow)183 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 178 184 / 179 185 !----------------------------------------------------------------------- 180 186 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 181 187 !----------------------------------------------------------------------- 182 ! SEND 188 ! send 183 189 cn_snd_temperature= 'weighted oce and ice' ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 184 190 cn_snd_albedo = 'weighted ice' ! 'none' 'weighted ice' 'mixed oce-ice' … … 188 194 cn_snd_crt_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 189 195 cn_snd_crt_grid = 'T' ! 'T' 190 ! RECEIVE 196 ! receive 191 197 cn_rcv_w10m = 'coupled' ! 'none' 'coupled' 192 198 cn_rcv_tau_nature = 'oce only' ! 'oce only' 'oce and ice' 'mixed oce-ice' … … 202 208 / 203 209 !----------------------------------------------------------------------- 204 &namqsr ! penetrative solar radiation 205 !----------------------------------------------------------------------- 206 ln_traqsr = .true. ! penetrative solar radiation (T) or not (F) 207 rabs = 0.58 ! fraction of qsr associated with xsi1 208 xsi1 = 0.35 ! first depth of extinction 209 xsi2 = 23.0 ! second depth of extinction 210 &namsbc_cpl_co2 ! coupled ocean/biogeo/atmosphere model ("key_cpl_carbon_cycle") 211 !----------------------------------------------------------------------- 212 cn_snd_co2 = 'coupled' ! send : 'none' 'coupled' 213 cn_rcv_co2 = 'coupled' ! receive : 'none' 'coupled' 214 / 215 !----------------------------------------------------------------------- 216 &namtra_qsr ! penetrative solar radiation 217 !----------------------------------------------------------------------- 218 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 219 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 220 sn_chl = 'chlorophyll', -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 221 222 cn_dir = './' ! root directory for the location of the runoff files 223 ln_traqsr = .true. ! Light penetration (T) or not (F) 224 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration 225 ln_qsr_2bd = .false. ! 2 bands light penetration 226 ln_qsr_bio = .false. ! bio-model light penetration 227 nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) 228 rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1) 229 rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction 230 rn_si1 = 23.0 ! 2 bands: longest depth of extinction 231 rn_si2 = 62.0 ! 3 bands: longest depth of extinction (for blue waveband & 0.01 mg/m2 Chl) 210 232 / 211 233 !----------------------------------------------------------------------- 212 234 &namsbc_rnf ! runoffs namelist surface boundary condition 213 235 !----------------------------------------------------------------------- 214 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !215 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !216 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. ,'yearly' , '' , ''217 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. ,'yearly' , '' , ''218 ! 236 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 237 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 238 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' , '' , '' 239 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 240 219 241 cn_dir = './' ! root directory for the location of the runoff files 220 242 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) 221 243 ln_rnf_mouth = .false. ! specific treatment at rivers mouths 222 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used244 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used 223 245 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 246 rn_rfact = 1.e0 ! multiplicative factor for runoff 224 247 / 225 248 !----------------------------------------------------------------------- 226 249 &namsbc_ssr ! surface boundary condition : sea surface restoring 227 250 !----------------------------------------------------------------------- 228 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !229 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !230 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly', '' , ''231 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly', '' , ''232 !251 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 252 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 253 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' , '' , '' 254 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' , '' , '' 255 233 256 cn_dir = './' ! root directory for the location of the runoff files 234 257 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 235 nn_sssr = 1 ! add a damping term in the surface freshwater flux (=1) or not (=0) 236 dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 237 deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] 258 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 259 ! or to SSS only (=1) or no damping term (=0) 260 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 261 rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] 262 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 263 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 238 264 / 239 265 !----------------------------------------------------------------------- 240 &namalb ! albedo parameters 241 !----------------------------------------------------------------------- 242 cgren = 0.06 ! correction of the snow or ice albedo to take into account the 243 albice = 0.53 ! albedo of melting ice in the arctic and antarctic 244 alphd = 0.80 ! coefficients for linear interpolation used to 245 alphc = 0.65 ! compute albedo between two extremes values 246 alphdi = 0.72 ! (Pyane, 1972) 247 / 266 &namsbc_alb ! albedo parameters 267 !----------------------------------------------------------------------- 268 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 269 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 270 rn_alphd = 0.80 ! coefficients for linear interpolation used to 271 rn_alphc = 0.65 ! compute albedo between two extremes values 272 rn_alphdi = 0.72 ! (Pyane, 1972) 273 / 274 248 275 !!====================================================================== 249 276 !! *** Lateral boundary condition *** … … 260 287 &namlbc ! lateral momentum boundary condition 261 288 !----------------------------------------------------------------------- 262 shlat = 2. ! shlat = 0 : free slip 263 ! 0 < shlat < 2 : partial slip 264 ! shlat = 2 : no slip 265 ! 2 < shlat : strong slip 289 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 290 ! free slip ! partial slip ! no slip ! strong slip 266 291 / 267 292 !----------------------------------------------------------------------- 268 293 &namcla ! cross land advection 269 294 !----------------------------------------------------------------------- 270 n _cla= 0 ! advection between 2 ocean pts separates by land295 nn_cla = 0 ! advection between 2 ocean pts separates by land 271 296 / 272 297 !----------------------------------------------------------------------- 273 298 &namobc ! open boundaries parameters ("key_obc") 274 299 !----------------------------------------------------------------------- 275 nobc_dta = 1 ! = 0 the obc data are equal to the initial state276 ! = 1 the obc data are read in 'obc.dta' files277 cffile = 'annual' ! set to annual if obc datafile hold 1 year of data278 ! set to monthly if obc datafile hold 1 month of data279 rdpein = 1. ! ???280 rdpwin = 1. ! ???281 rdpnin = 1. ! ???282 rdpsin = 1. ! ???283 rdpeob = 3000. ! time relaxation (days) for the east open boundary284 rdpwob = 15. ! " " " west "285 rdpnob = 3000. ! " " " north "286 rdpsob = 15. ! " " " south "287 zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline288 zbsic2 = 1.e+6 ! " " second "289 zbsic3 = 0. ! " " thrid "290 300 ln_obc_clim= .false. ! climatological obc data files (T) or not (F) 291 301 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 292 / 293 !----------------------------------------------------------------------- 294 &namagrif ! ("key_agrif") 295 !----------------------------------------------------------------------- 296 nbclineupdate = 3 ! baroclinic update frequency 302 ln_obc_fla = .false. ! Flather open boundary condition 303 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 304 ! = 1 the obc data are read in 'obc.dta' files 305 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 306 ! set to monthly if obc datafile hold 1 month of data 307 rn_dpein = 1. ! damping time scale for inflow at east open boundary 308 rn_dpwin = 1. ! - - - west - - 309 rn_dpnin = 1. ! - - - north - - 310 rn_dpsin = 1. ! - - - south - - 311 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 312 rn_dpwob = 15. ! - - - west - - 313 rn_dpnob = 3000. ! - - - north - - 314 rn_dpsob = 15. ! - - - south - - 315 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 316 ! = 1 the total volume remains constant 317 / 318 !----------------------------------------------------------------------- 319 &namagrif ! AGRIF zoom ("key_agrif") 320 !----------------------------------------------------------------------- 321 nn_cln_update = 3 ! baroclinic update frequency 297 322 ln_spc_dyn = .true. ! use 0 as special value for dynamics 298 visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer299 visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer300 / 301 !----------------------------------------------------------------------- 302 &nambdy ! unstructured open boundaries parameters("key_bdy")323 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 324 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 325 / 326 !----------------------------------------------------------------------- 327 &nambdy ! unstructured open boundaries ("key_bdy") 303 328 !----------------------------------------------------------------------- 304 329 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) … … 320 345 / 321 346 !----------------------------------------------------------------------- 322 &nam tide! tidal forcing at unstructured boundaries347 &nambdy_tide ! tidal forcing at unstructured boundaries 323 348 !----------------------------------------------------------------------- 324 349 filtide = 'bdytide_' ! file name root of tidal forcing files … … 327 352 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 328 353 / 354 329 355 !!====================================================================== 330 356 !! *** Bottom boundary condition *** … … 338 364 &nambfr ! bottom friction 339 365 !----------------------------------------------------------------------- 340 n botfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction366 nn_bfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction 341 367 ! = 3 : free slip, = 1 : linear friction 342 bfri1= 4.e-4 ! bottom drag coefficient (linear case)343 bfri2= 1.e-3 ! bottom drag coefficient (non linear case)344 bfeb2= 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2)368 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 369 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 370 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 345 371 / 346 372 !----------------------------------------------------------------------- 347 373 &nambbc ! bottom temperature boundary condition 348 374 !----------------------------------------------------------------------- 349 n geo_flux = 2 ! geothermal heat flux = 0 no flux considered350 ! 351 ! 352 ngeo_flux_const = 86.4e-3! Constant value of geothermal heat flux [W/m2]375 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 376 ! = 1 constant flux 377 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 378 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 353 379 / 354 380 !----------------------------------------------------------------------- … … 357 383 ! ! diffusive bbl ("key_trabbl") 358 384 ! ! advective bbl ("key_trabbl_adv") 359 atrbbl= 10000. ! lateral mixing coefficient in the bbl [m2/s]385 rn_ahtbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 360 386 / 361 387 !!====================================================================== … … 363 389 !!====================================================================== 364 390 !! nameos equation of state 365 !! nam _traadv advection scheme366 !! nam _traldf lateral diffusion scheme367 !! namt dp tracer newtonian damping("key_tradmp")391 !! namtra_adv advection scheme 392 !! namtra_ldf lateral diffusion scheme 393 !! namtra_dmp T & S newtonian damping ("key_tradmp") 368 394 !!====================================================================== 369 395 … … 371 397 &nameos ! ocean physical parameters 372 398 !----------------------------------------------------------------------- 373 n eos= 0 ! type of equation of state and Brunt-Vaisala frequency399 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 374 400 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 375 401 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 376 402 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 377 r alpha= 2.e-4 ! thermal expension coefficient (neos= 1 or 2)378 r beta= 0.001 ! saline expension coefficient (neos= 2)379 / 380 !----------------------------------------------------------------------- 381 &nam _traadv ! advection scheme for tracer382 !----------------------------------------------------------------------- 383 ln_traadv_cen2 = . true.! 2nd order centered scheme384 ln_traadv_tvd = . false.! TVD scheme403 rn_alpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 404 rn_beta = 0.001 ! saline expension coefficient (neos= 2) 405 / 406 !----------------------------------------------------------------------- 407 &namtra_adv ! advection scheme for tracer 408 !----------------------------------------------------------------------- 409 ln_traadv_cen2 = .false. ! 2nd order centered scheme 410 ln_traadv_tvd = .true. ! TVD scheme 385 411 ln_traadv_muscl = .false. ! MUSCL scheme 386 412 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries … … 388 414 / 389 415 !----------------------------------------------------------------------- 390 &nam _traldf ! lateral diffusion scheme for tracer391 !----------------------------------------------------------------------- 392 !! Type of the operator :416 &namtra_ldf ! lateral diffusion scheme for tracer 417 !----------------------------------------------------------------------- 418 ! Type of the operator : 393 419 ln_traldf_lap = .true. ! laplacian operator 394 420 ln_traldf_bilap = .false. ! bilaplacian operator 395 421 ! Direction of action : 396 422 ln_traldf_level = .false. ! iso-level 397 423 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 398 424 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 399 ! ! Coefficient 400 aht0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 401 ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 402 aeiv0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 403 / 404 !----------------------------------------------------------------------- 405 &namtdp ! tracer newtonian damping ('key_tradmp') 406 !----------------------------------------------------------------------- 407 ndmp = -1 ! type of damping in temperature and salinity 408 ! ='latitude', damping poleward of 'ndmp' degrees and function 409 ! of the distance-to-coast. Red and Med Seas as ndmp=-1 410 ! =-1 damping only in Med and Red Seas 411 ndmpf = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 412 nmldmp = 1 ! type of damping: =0 damping throughout the water column 413 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 414 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 415 sdmp = 50. ! surface time scale for internal damping (days) 416 bdmp = 360. ! bottom time scale for internal damping (days) 417 hdmp = 800. ! depth of transition between sdmp and bdmp (meters) 425 ! Coefficient 426 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 427 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 428 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 429 / 430 !----------------------------------------------------------------------- 431 &namtra_dmp ! tracer: T & S newtonian damping ('key_tradmp') 432 !----------------------------------------------------------------------- 433 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only 434 ! =XX, damping poleward of XX degrees (XX>0) 435 ! + F(distance-to-coast) + Red and Med Seas 436 nn_zdmp = 1 ! vertical shape =0 damping throughout the water column 437 ! =1 no damping in the mixing layer (kz criteria) 438 ! =2 no damping in the mixed layer (rho crieria) 439 rn_surf = 50. ! surface time scale of damping [days] 440 rn_bot = 360. ! bottom time scale of damping [days] 441 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 442 nn_file = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 418 443 / 419 444 !!====================================================================== 420 445 !! *** Dynamics namelists *** 421 446 !!====================================================================== 422 !! nam_dynadv formulation of the momentum advection 423 !! nam_dynvor advection scheme 424 !! nam_dynhpg hydrostatic pressure gradient 425 !! namflg hydrostatic pressure gradient time stepping 426 !! nam_dynspg surface pressure gradient (CPP key only) 427 !! nam_dynldf lateral diffusion scheme 428 !!====================================================================== 429 430 !----------------------------------------------------------------------- 431 &nam_dynadv ! formulation of the momentum advection 447 !! namdyn_adv formulation of the momentum advection 448 !! namdyn_vor advection scheme 449 !! namdyn_hpg hydrostatic pressure gradient 450 !! namdyn_spg surface pressure gradient (CPP key only) 451 !! namdyn_ldf lateral diffusion scheme 452 !!====================================================================== 453 454 !----------------------------------------------------------------------- 455 &namdyn_adv ! formulation of the momentum advection 432 456 !----------------------------------------------------------------------- 433 457 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 436 460 / 437 461 !----------------------------------------------------------------------- 438 &nam _dynvor ! option of physics/algorithm (not control by CPP keys)462 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 439 463 !----------------------------------------------------------------------- 440 464 ln_dynvor_ene = .false. ! enstrophy conserving scheme 441 ln_dynvor_ens = . true.! energy conserving scheme465 ln_dynvor_ens = .false. ! energy conserving scheme 442 466 ln_dynvor_mix = .false. ! mixed scheme 443 ln_dynvor_een = . false.! energy & enstrophy scheme444 / 445 !----------------------------------------------------------------------- 446 &nam _dynhpg ! Hydrostatic pressure gradient option467 ln_dynvor_een = .true. ! energy & enstrophy scheme 468 / 469 !----------------------------------------------------------------------- 470 &namdyn_hpg ! Hydrostatic pressure gradient option 447 471 !----------------------------------------------------------------------- 448 472 ln_hpg_zco = .false. ! z-coordinate - full steps … … 453 477 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 454 478 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 455 gamm = 0.e0 ! weighting coefficient (wdj scheme) 456 / 457 !----------------------------------------------------------------------- 458 &namflg ! algorithm flags (algorithm not control by CPP keys) 459 !----------------------------------------------------------------------- 460 ln_dynhpg_imp = .false. ! hydrostatic pressure gradient: semi-implicit time scheme (T) 461 ! centered time scheme (F) 462 nn_dynhpg_rst = 0 ! add dynhpg implicit variables in restart ot not (1/0) 463 / 464 !----------------------------------------------------------------------- 465 !nam_dynspg ! surface pressure gradient (CPP key only) 479 rn_gamma = 0.e0 ! weighting coefficient (wdj scheme) 480 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 481 ! centered time scheme (F) 482 nn_dynhpg_rst = 0 ! =1 dynhpg restartable restart or not (=0) 483 / 484 !----------------------------------------------------------------------- 485 !namdyn_spg ! surface pressure gradient (CPP key only) 466 486 !----------------------------------------------------------------------- 467 487 ! ! explicit free surface ("key_dynspg_exp") 468 488 ! ! filtered free surface ("key_dynspg_flt") 469 489 ! ! split-explicit free surface ("key_dynspg_ts") 470 ! ! rigid-lid ("key_dynspg_rl") 471 472 !----------------------------------------------------------------------- 473 &nam_dynldf ! lateral diffusion on momentum 474 !----------------------------------------------------------------------- 475 ! ! Type of the operator : 476 ln_dynldf_lap = .true. ! laplacian operator 477 ln_dynldf_bilap = .false. ! bilaplacian operator 478 ! ! Direction of action : 479 ln_dynldf_level = .false. ! iso-level 480 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 481 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 482 ! Coefficient 483 ahm0 = 40000. ! horizontal eddy viscosity [m2/s] 484 ahmb0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 490 491 !----------------------------------------------------------------------- 492 &namdyn_ldf ! lateral diffusion on momentum 493 !----------------------------------------------------------------------- 494 ! Type of the operator : 495 ln_dynldf_lap = .true. ! laplacian operator 496 ln_dynldf_bilap = .false. ! bilaplacian operator 497 ! Direction of action : 498 ln_dynldf_level = .false. ! iso-level 499 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 500 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 501 ! Coefficient 502 rn_ahm_0 = 40000. ! horizontal eddy viscosity [m2/s] 503 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 485 504 / 486 505 !!====================================================================== … … 488 507 !!====================================================================== 489 508 !! namzdf vertical physics 490 !! nam npc non penetrative convection491 !! nam ric richardson number dependent vertical mixing ("key_zdfric" )492 !! nam tke TKE dependent vertical mixing ("key_zdftke" )493 !! nam kpp KPP dependent vertical mixing ("key_zdfkpp" )494 !! nam ddm double diffusive mixing parameterization ("key_zdfddm" )509 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric" ) 510 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke" ) 511 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp" ) 512 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm" ) 513 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx" ) 495 514 !!====================================================================== 496 515 … … 498 517 &namzdf ! vertical physics 499 518 !----------------------------------------------------------------------- 500 avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 501 avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 502 ln_zdfnpc = .false. ! convection: Non-Penetrative algorithm (T) or not (F) 503 ln_zdfevd = .true. ! convection: enhanced vertical diffusion (T) or not (F) 504 avevd = 100. ! vertical coefficient for enhanced diffusion scheme [m2/s] 505 n_evdm = 1 ! enhanced mixing apply on tracer (=0) or on tracer and momentum (=1) 506 ln_zdfexp = .false. ! split explicit (T) or implicit (F) time stepping 507 n_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 508 / 509 !----------------------------------------------------------------------- 510 &namnpc ! non penetrative convection 511 !----------------------------------------------------------------------- 512 nnpc1 = 1 ! non penetrative convective scheme computation frequency 513 nnpc2 = 365 ! non penetrative convective scheme print frequency 514 / 515 !----------------------------------------------------------------------- 516 &namric ! richardson number dependent vertical diffusion ("key_zdfric" ) 517 !----------------------------------------------------------------------- 518 avmri = 100.e-4 ! maximum value of the vertical viscosity 519 alp = 5. ! coefficient of the parameterization 520 nric = 2 ! coefficient of the parameterization 521 / 522 !----------------------------------------------------------------------- 523 &namtke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 524 !----------------------------------------------------------------------- 525 ln_rstke = .false. ! restart with tke from a run without tke (T) or not (F) 526 nn_itke = 50 ! number of iterative loops if ln_rstke=T 519 rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 520 rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 521 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 522 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 523 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) 524 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1) 525 rn_avevd = 100. ! evd mixing coefficient [m2/s] 526 ln_zdfnpc = .false. ! Non-Penetrative algorithm (T) or not (F) 527 nn_npc = 1 ! frequency of application of npc 528 nn_npcp = 365 ! npc control print frequency 529 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping 530 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 531 / 532 !----------------------------------------------------------------------- 533 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) 534 !----------------------------------------------------------------------- 535 rn_avmri = 100.e-4 ! maximum value of the vertical viscosity 536 rn_alp = 5. ! coefficient of the parameterization 537 nn_ric = 2 ! coefficient of the parameterization 538 / 539 !----------------------------------------------------------------------- 540 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 541 !----------------------------------------------------------------------- 527 542 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 528 543 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation 529 rn_ebb = 3.75 ! coef. of the surface input of tke 530 rn_efave = 1. ! boost of the tke diffusion ( avtke=rn_efave*avm ) 544 rn_ebb = 60. ! coef. of the surface input of tke 531 545 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 532 546 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 547 rn_bshear = 1.e-20 ! background shear (>0) 533 548 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 534 549 ! = 1 bounded by the local vertical scale factor … … 536 551 ! = 3 same criteria as case 2 but applied in a different way 537 552 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 538 nn_avb = 0 ! profile for constant background used on avt & avm (=1) or not (=0)539 nn_ave = 1 ! horizontal averaged on avt (=1) or not (=0)540 553 ln_mxl0 = .false. ! mixing length scale surface value as function of wind stress (T) or not (F) 541 rn_lmin = 0. 4! interior buoyancy lenght scale minimum value542 rn_lmin0 = 0. 4! surface buoyancy lenght scale minimum value554 rn_lmin = 0.001 ! interior buoyancy lenght scale minimum value 555 rn_lmin0 = 0.01 ! surface buoyancy lenght scale minimum value 543 556 nn_etau = 0 ! exponentially deceasing penetration of tke due to internal & intertial waves 544 557 ! = 0 no penetration ( O(2 km) resolution) 545 ! = 1 additional tke source 558 ! = 1 additional tke source 546 559 ! = 2 additional tke source applied only at the base of the mixed layer 547 nn_htau = 2! type of exponential decrease of tke penetration560 nn_htau = 1 ! type of exponential decrease of tke penetration 548 561 ! = 0 constant 10 m length scale 549 ! = 1 ??? 550 ! = 2 ??? 562 ! = 1 0.5m at the equator to 30m at high latitudes 551 563 rn_efr = 0.05 ! fraction of surface tke value which penetrates inside the ocean 552 ln_lc = .false. ! Langmuir cell effect564 ln_lc = .false. ! Langmuir cell parameterisation 553 565 rn_lc = 0.15 ! coef. associated to Langmuir cells 554 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)555 566 / 556 567 !------------------------------------------------------------------------ 557 &namkpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 558 ! "key_kppcustom" or "key_kpplktb") 559 !------------------------------------------------------------------------ 568 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 569 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 560 570 ln_kpprimix = .true. ! shear instability mixing 561 difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 562 difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 563 Riinfty = 0.8 ! local Richardson Number limit for shear instability 564 difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 565 bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 566 difcon = 1. ! maximum mixing in interior convection [m2/s] 567 navb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 568 nave = 1 ! constant (=0) or profile (=1) background on avt 569 / 570 !----------------------------------------------------------------------- 571 &namddm ! double diffusive mixing parameterization ("key_zdfddm") 572 !----------------------------------------------------------------------- 573 avts = 1.e-4 ! maximum avs (vertical mixing on salinity) 574 hsbfr = 1.6 ! heat/salt buoyancy flux ratio 571 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 572 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 573 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability 574 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 575 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 576 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 577 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 578 nn_ave = 1 ! constant (=0) or profile (=1) background on avt 579 / 580 !----------------------------------------------------------------------- 581 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") 582 !----------------------------------------------------------------------- 583 rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity) 584 rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio 585 / 586 !----------------------------------------------------------------------- 587 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") 588 !----------------------------------------------------------------------- 589 rn_htmx = 500. ! vertical decay scale for turbulence (meters) 590 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 591 rn_tfe = 0.333 ! tidal dissipation efficiency 592 rn_me = 0.2 ! mixing efficiency 593 ln_tmx_itf = .FALSE. ! ITF specific parameterisation 594 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 575 595 / 576 596 !!====================================================================== 577 597 !! *** Miscelaneous namelists *** 578 598 !!====================================================================== 579 !! nam _mppMassively Parallel Processing ("key_mpp_mpi)580 !! nam _mpp_dyndistMassively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist")599 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 600 !! nammpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist") 581 601 !! namctl Control prints & Benchmark 582 602 !! namsol elliptic solver / island / free surface … … 586 606 &namsol ! elliptic solver / island / free surface 587 607 !----------------------------------------------------------------------- 588 n solv= 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)608 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 589 609 ! =2 successive-over-relaxation (sor) 590 ! =3 FETI (fet) ("key_feti") 591 ! =4 sor with extra outer halo 592 nsol_arp = 0 ! absolute/relative (0/1) precision convergence test 593 nmin = 300 ! minimum of iterations for the SOR solver 594 nmax = 800 ! maximum of iterations for the SOR solver 595 nmod = 10 ! frequency of test for the SOR solver 596 eps = 1.e-6 ! absolute precision of the solver 597 resmax = 1.e-10 ! absolute precision for the SOR solver 598 sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 599 epsisl = 1.e-10 ! absolute precision on stream function solver 600 nmisl = 4000 ! maximum pcg iterations for island ("key_islands") 601 rnu = 1. ! strength of the additional force used in filtered free surface 602 / 603 !----------------------------------------------------------------------- 604 &nam_mpp ! Massively Parallel Processing ("key_mpp_mpi) 605 !----------------------------------------------------------------------- 606 c_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 610 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test 611 rn_eps = 1.e-6 ! absolute precision of the solver 612 nn_nmin = 300 ! minimum of iterations for the SOR solver 613 nn_nmax = 800 ! maximum of iterations for the SOR solver 614 nn_nmod = 10 ! frequency of test for the SOR solver 615 rn_resmax = 1.e-10 ! absolute precision for the SOR solver 616 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 617 / 618 !----------------------------------------------------------------------- 619 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) 620 !----------------------------------------------------------------------- 621 cn_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 607 622 ! buffer blocking send or immediate non-blocking sends, resp. 608 nn_buffer = 0! size in bytes of exported buffer ('B' case), 0 no exportation609 / 610 !----------------------------------------------------------------------- 611 &nam _mpp_dyndist ! Massively Parallel Distribution("key_agrif" && "key_mpp_dyndist")612 !----------------------------------------------------------------------- 613 jpni = 1! jpni number of processors following i614 jpnj = 1! jpnj number of processors following j615 jpnij = 1! jpnij number of local domains616 / 617 !----------------------------------------------------------------------- 618 &namctl ! Control prints & Benchmark619 !----------------------------------------------------------------------- 620 ln_ctl = .false.! trends control print (expensive!)621 n print = 0! level of print (0 no extra print)622 n ictls = 0! start i indice of control sum (use to compare mono versus623 n ictle = 0! end i indice of control sum multi processor runs624 n jctls = 0! start j indice of control over a subdomain)625 n jctle = 0 ! end j indice of control626 isplt = 1! number of processors in i-direction627 jsplt = 1! number of processors in j-direction628 n bench = 0! Bench mode (1/0): CAUTION use zero except for bench623 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 624 / 625 !----------------------------------------------------------------------- 626 &nammpp_dyndist ! Massively Parallel Distribution for AGRIF zoom ("key_agrif" && "key_mpp_dyndist") 627 !----------------------------------------------------------------------- 628 jpni = 1 ! jpni number of processors following i 629 jpnj = 1 ! jpnj number of processors following j 630 jpnij = 1 ! jpnij number of local domains 631 / 632 !----------------------------------------------------------------------- 633 &namctl ! Control prints & Benchmark 634 !----------------------------------------------------------------------- 635 ln_ctl = .false. ! trends control print (expensive!) 636 nn_print = 0 ! level of print (0 no extra print) 637 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 638 nn_ictle = 0 ! end i indice of control sum multi processor runs 639 nn_jctls = 0 ! start j indice of control over a subdomain) 640 nn_jctle = 0 ! end j indice of control 641 nn_isplt = 1 ! number of processors in i-direction 642 nn_jsplt = 1 ! number of processors in j-direction 643 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 629 644 ! (no physical validity of the results) 630 n bit_cmp = 0! bit comparison mode (1/0): CAUTION use zero except for test645 nn_bit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test 631 646 ! of comparison between single and multiple processor runs 632 647 / 648 633 649 !!====================================================================== 634 650 !! *** Diagnostics namelists *** … … 636 652 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 637 653 !! namgap level mean model-data gap ("key_diagap") 638 !! namspr surface pressure diagnosed in rigid-lid ("key_diaspr")639 654 !! namflo float parameters ("key_float") 640 655 !! namptr Poleward Transport Diagnostics … … 642 657 643 658 !----------------------------------------------------------------------- 644 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 645 ! or mixed-layer trends ('key_trdmld') 646 ! or barotropic vorticity ("key_trdvor") 647 !----------------------------------------------------------------------- 648 ntrd = 365 ! time step frequency dynamics and tracers trends 649 nctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 650 ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 651 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 652 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 653 ln_trdmld_restart = .false. ! restart for ML diagnostics 654 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 659 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 660 ! ! or mixed-layer trends or barotropic vorticity ('key_trdmld' or "key_trdvor") 661 !----------------------------------------------------------------------- 662 nn_trd = 365 ! time step frequency dynamics and tracers trends 663 nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 664 rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 665 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 666 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 667 ln_trdmld_restart = .false. ! restart for ML diagnostics 668 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 655 669 / 656 670 !----------------------------------------------------------------------- 657 671 &namgap ! level mean model-data gap ('key_diagap') 658 672 !----------------------------------------------------------------------- 659 ngap = 15 ! time-step frequency of model-data gap computation 660 nprg = 10 ! time-step frequency of gap print in model output 661 / 662 !----------------------------------------------------------------------- 663 &namspr ! surface pressure diagnostic 664 !----------------------------------------------------------------------- 665 nmaxp = 1000 ! maximum of iterations for the solver 666 epsp = 1.e-3 ! absolute precision of the solver 667 niterp = 400 ! number of iteration done by the solver 673 nn_gap = 15 ! time-step frequency of model-data gap computation 674 nn_prg = 10 ! time-step frequency of gap print in model output 668 675 / 669 676 !----------------------------------------------------------------------- … … 671 678 !----------------------------------------------------------------------- 672 679 ln_rstflo = .false. ! float restart (T) or not (F) 673 n writefl= 75 ! frequency of writing in float output file674 n stockfl= 5475 ! frequency of creation of the float restart file680 nn_writefl= 75 ! frequency of writing in float output file 681 nn_stockfl= 5475 ! frequency of creation of the float restart file 675 682 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 676 683 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) … … 680 687 &namptr ! Poleward Transport Diagnostic 681 688 !----------------------------------------------------------------------- 682 ln_diaptr = .false. 683 ln_diaznl = .false. 684 ln_subbas = .false. 689 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 690 ln_diaznl = .false. ! Add zonal means and meridional stream functions 691 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 685 692 ! (orca configuration only, need input basins mask file named "subbasins.nc" 686 693 nf_ptr = 1 ! Frequency of ptr computation [time step] -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/PARAM/namelist_pisces
r1379 r1631 114 114 &nampissed ! parameters for inputs deposition 115 115 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 116 bdustfer= .true. ! boolean for dust input from the atmosphere117 briver= .true. ! boolean for river input of nutrients118 bndepo= .true. ! boolean for atmospheric deposition of N119 bsedinput= .true. ! boolean for Fe input from sediments120 sedfeinput = 1E-9 ! Coastal release of Iron121 dustsolub = 0.014 ! Solubility of the dust116 ln_dustfer = .true. ! boolean for dust input from the atmosphere 117 ln_river = .true. ! boolean for river input of nutrients 118 ln_ndepo = .true. ! boolean for atmospheric deposition of N 119 ln_sedinput = .true. ! boolean for Fe input from sediments 120 sedfeinput = 1E-9 ! Coastal release of Iron 121 dustsolub = 0.014 ! Solubility of the dust 122 122 / 123 123 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' … … 165 165 pisdia3d(5) = 'PPPHY ' , 'Primary production of nanophyto ', 'molC/m3/s ' 166 166 pisdia3d(6) = 'PPPHY2 ' , 'Primary production of diatoms ', 'molC/m3/s ' 167 pisdia3d(7) = 'PP ZOO ' , 'Primary production of microzoo ', 'molC/m3/s '168 pisdia3d(8) = 'PP ZOO2 ' , 'Primary production of mesozoo', 'molC/m3/s '167 pisdia3d(7) = 'PPNEWN ' , 'New Primary production of nano ', 'molC/m3/s ' 168 pisdia3d(8) = 'PPNEWD ' , 'New Primary production of diat ', 'molC/m3/s ' 169 169 pisdia3d(9) = 'PBSi ' , 'Primary production of Si diatoms ', 'molSi/m3/s ' 170 170 pisdia3d(10) = 'PFeN ' , 'Primary production of nano iron ', 'molFe/m3/s ' -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/PARAM/namelist_top
r1379 r1631 13 13 ndttrc = 1 ! time step frequency for passive tracers 14 14 nwritetrc = 5475 ! time step frequency for tracer outputs 15 l rsttr= .false. ! start from a restart file (T) or not (F)15 ln_rsttr = .false. ! start from a restart file (T) or not (F) 16 16 nrsttr = 0 ! restart control = 0 initial time step is not compared to the restart file value 17 17 ! = 1 do not use the value in the restart file … … 23 23 ! ! ! ! ! from file ! or not ! 24 24 ! ! ! ! ! or not ! ! 25 tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', ' umol/L' , .false. , .true.26 tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', ' ueq/L ' , .false. , .true.27 tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', ' umol/L' , .false. , .true.28 tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', ' umol/L' , .false. , .true.29 tracer(5) = 'PO4 ' , 'Phosphate Concentration ', ' umol/L' , .false. , .true.30 tracer(6) = 'POC ' , 'Small organic carbon Concentration ', ' umol/L' , .false. , .true.31 tracer(7) = 'Si ' , 'Silicate Concentration ', ' umol/L' , .false. , .true.32 tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', ' umol/L' , .false. , .true.33 tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', ' umol/L' , .false. , .true.34 tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', ' umol/L' , .false. , .true.35 tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', ' umol/L' , .false. , .true.36 tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', ' umol/L' , .false. , .true.37 tracer(13) = 'BSi ' , 'Diatoms Silicate Concentration ', ' umol/L' , .false. , .true.38 tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', ' umol/L' , .false. , .true.39 tracer(15) = 'BFe ' , 'Big iron particles Concentration ', ' umol/L' , .false. , .true.40 tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', ' umol/L' , .false. , .true.41 tracer(17) = 'SFe ' , 'Small iron particles Concentration ', ' umol/L' , .false. , .true.42 tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', ' umol/L' , .false. , .true.43 tracer(19) = 'DSi ' , 'Sinking biogenic Silicate Concentration', ' umol/L' , .false. , .true.44 tracer(20) = 'NFe ' , 'Nano iron Concentration ', ' umol/L' , .false. , .true.45 tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', ' umol/L' , .false. , .true.46 tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', ' umol/L' , .false. , .true.47 tracer(23) = 'NO3 ' , 'Nitrates Concentration ', ' umol/L' , .false. , .true.48 tracer(24) = 'NH4 ' , 'Ammonium Concentration ', ' umol/L' , .false. , .true.25 tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'mol-C/L' , .true. , .true. 26 tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true. 27 tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'mol-C/L' , .true. , .true. 28 tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'mol-C/L' , .false. , .true. 29 tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'mol-C/L' , .true. , .true. 30 tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'mol-C/L' , .false. , .true. 31 tracer(7) = 'Si ' , 'Silicate Concentration ', 'mol-C/L' , .true. , .true. 32 tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'mol-C/L' , .false. , .true. 33 tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'mol-C/L' , .false. , .true. 34 tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'mol-C/L' , .true. , .true. 35 tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'mol-C/L' , .false. , .true. 36 tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'mol-C/L' , .false. , .true. 37 tracer(13) = 'BSi ' , 'Diatoms Silicate Concentration ', 'mol-C/L' , .false. , .true. 38 tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'mol-C/L' , .true. , .true. 39 tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'mol-C/L' , .false. , .true. 40 tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'mol-C/L' , .false. , .true. 41 tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'mol-C/L' , .false. , .true. 42 tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'mol-C/L' , .false. , .true. 43 tracer(19) = 'DSi ' , 'Sinking biogenic Silicate Concentration', 'mol-C/L' , .false. , .true. 44 tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'mol-C/L' , .false. , .true. 45 tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'mol-C/L' , .false. , .true. 46 tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'mol-C/L' , .false. , .true. 47 tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'mol-C/L' , .true. , .true. 48 tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'mol-C/L' , .false. , .true. 49 49 / 50 50 !----------------------------------------------------------------------- -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/README
r1379 r1631 4 4 trunk of 12 months 5 5 6 Wit hout monitoring andatlas6 Wit monitoring, without atlas 7 7 8 8 Tested on NEC SX-8R mercure at CCRT. -
branches/libIGCM/ORCA2_LIM_PISCES/IGCM00/config.card
r1379 r1631 16 16 #-- leap, noleap, 360d 17 17 CalendarType=360d 18 #-- Début et fin deJob18 #-- Begin and end of Job 19 19 #-- "YYYY-MM-DD" 20 20 DateBegin=2001-01-01 … … 39 39 # config.card configuration options : 40 40 #R_INIT=/home/mancip/PROG/IPSL/FORCAGE_OL/INIT 41 #R_BC=/home/mancip/PROG/IPSL/FORCAGE_OL/BC 42 R_BC=/dmnfs/cont003/p48ethe/IGCM/BC 41 R_BC=/dmnfs/cont003/p86ipsl/IGCM/BC 43 42 #======================================================================== 44 43 #D-- ListOfComponents - 45 44 [ListOfComponents] 46 45 #D- For each component, Name of component, Tag of component 47 OCE= (opa9, NEMO_v3 _1)46 OCE= (opa9, NEMO_v3) 48 47 ICE= (lim2, LIM_2) 49 TOP= (pisces, PISCES)48 MBG= (pisces, PISCES) 50 49 #======================================================================== 51 50 #D-- Executable - … … 55 54 OCE= (opa, nemo) 56 55 ICE= ("", "") 57 TOP= ("", "") 58 #D- TRENDS component required only for trends output 59 ###TRENDS= ("", "") 60 56 MBG= ("", "") 61 57 #======================================================================== 62 58 #D-- Restarts - … … 80 76 #D- If you want to produce time series, this flag determines 81 77 #D- frequency of post-processing submission 82 TimeSeriesFrequency=10Y 78 RebuildFromArchive=true 79 #D- If you want to produce time series, this flag determines 80 #D- frequency of post-processing submission 81 TimeSeriesFrequency=1M 83 82 #D- If you want to produce seasonal average, this flag determines 84 83 #D- the period of this average … … 114 113 115 114 #======================================================================== 116 #D-- TOP-117 [ TOP]115 #D-- MBG - 116 [MBG] 118 117 WriteFrequency="1M" 119 118 #WriteFrequency="5D"
Note: See TracChangeset
for help on using the changeset viewer.