Changeset 568 for CONFIG/IPSLCM/IPSLCM5/trunk
- Timestamp:
- 02/27/09 11:19:52 (16 years ago)
- Location:
- CONFIG/IPSLCM/IPSLCM5/trunk/EXP00
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/IPSLCM/IPSLCM5/trunk/EXP00/COMP/nemo.driver
r533 r568 8 8 RESOL_OCE=$( echo $RESOL | awk "-Fx" '{print $1}' ) 9 9 10 if [ -z "${opa_UserChoices_OPA_NPDT_JOURS}" ] ; then 11 ORCA_NPDT_JOUR=15 12 else 13 ORCA_NPDT_JOUR=${opa_UserChoices_OPA_NPDT_JOURS} 14 fi 10 ORCA_NPDT_JOUR=${nemo_UserChoices_OPA_NPDT_JOURS:=15} 15 11 16 12 # Local function to find namelists parameters … … 22 18 23 19 # cexper experience name for vairmer format 20 # ln_rstart boolean term for restart (true or false) 21 # nrstdt control of the time step (0, 1 or 2) 24 22 # nit000 number of the first time step 25 23 # nitend number of the last time step 26 # nbisex Leap year calendar (0/1) (30 for 360d) 24 # ndate0 initial calendar date aammjj 25 # nleapy Leap year calendar (1) or not (0), or 360 days calendar (30) 26 # nstock frequency of restart file 27 27 # nwrite frequency of OUTPUT file 28 # nmsh =1 create a mesh file (coordinates, scale factors, masks) 28 29 # nwrihf frequency of HF OUTPUT file 29 # lrstar boolean term for restart (true or false) 30 # nstock frequency of restart file 31 # nrstdt control of the time step (0, 1 or 2) 32 # ndate0 initial calendar date aammjj 33 # nmsh =1 create a mesh file (coordinates, scale factors, masks) 30 # rdt time step in seconds (coming from namelist) 34 31 35 32 PAT_CEXPER=$( supergrep cexper ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) … … 41 38 PAT_NRSTAR=$( supergrep nrstdt ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 42 39 PAT_NDATE0=$( supergrep ndate0 ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 40 PAT_NLEAPY=$( supergrep nleapy ${SUBMIT_DIR}/PARAM/namelist ) 43 41 PAT_NMSH=$( supergrep nmsh ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 44 42 45 ORCA_NSTOCK=$( echo $PAT_NSTOCK | sed "s/[a-z,A-Z,=]//g" )46 43 ORCA_RDT=$( supergrep rdt ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | sed 's/ *rdt *=//' | sed 's/\. *//' ) 47 48 ##--Transformation namelist fortran 77 --> fortran 9049 # ? ? Is it really necessary ? ?50 #sed -e /:/d -e s/" \&END"/"\/"/ -e s/" \&"/"\&"/ ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} > ${SUBMIT_DIR}/PARAM/namelist51 44 52 45 # Period Length In Days between DateBegin and first day of calendar 0001 01 01 … … 126 119 (( ORCA_NITEND = ORCA_NIT000 + ORCA_NPDT - 1)) 127 120 128 if [ "${ORCA_NSTOCK}" -gt "${ORCA_NPDT}" ] 129 then 130 ORCA_NSTOCK="${ORCA_NPDT}" 131 fi 121 ORCA_NSTOCK="${ORCA_NPDT}" 132 122 133 123 ##-- verification du nb de pas de temps par jour … … 137 127 if [ ${NB_SEC_JOUR_MODEL} -ne ${NB_SEC_JOUR} ] 138 128 then 139 echo " VERIFIER ORCA_NPDT_JOUR dans le job ${JOB}"129 echo " VERIFIER ORCA_NPDT_JOUR dans nemo.card " 140 130 exit 141 131 fi … … 147 137 ORCA_NRSTDT=0 148 138 ORCA_NMSH=1 149 # ??? a garder ???150 # cat <<EOF > EMPave_old.dat151 # Old global EMP coeff = 1.000000000152 # EOF153 139 154 140 #echo "PAS DE RESTART OPA" … … 163 149 164 150 ORCA_LRSTAR=.TRUE. 165 ORCA_NRSTDT= 1151 ORCA_NRSTDT=2 166 152 ORCA_NMSH=0 167 153 168 154 #echo "RESTART OPA" 169 155 fi 170 typeset -r PRECIS=8 171 NEMO_END=$( echo $( awk "BEGIN { printf \"%0${PRECIS}d\",${ORCA_NITEND} }" ) ) 156 157 158 # nleapy configuration 159 case ${config_UserChoices_CalendarType} in 160 leap) 161 ORCA_NLEAPY=1;; 162 noleap) 163 ORCA_NLEAPY=0;; 164 360d) 165 ORCA_NLEAPY=30;; 166 *) 167 ORCA_NLEAPY=30 168 esac 169 170 171 typeset -r PRECIS=8 172 NEMO_END=$( echo $( awk "BEGIN { printf \"%0${PRECIS}d\",${ORCA_NITEND} }" ) ) 173 172 174 173 175 sed -e "s/${PAT_CEXPER}/ cexper=\"${config_UserChoices_JobName}\"/" \ 174 -e "s/${PAT_NIT000}/ nit000=${ORCA_NIT000}/" \ 175 -e "s/${PAT_NITEND}/ nitend=${ORCA_NITEND}/" \ 176 -e "s/${PAT_NWRITE}/ nwrite=${ORCA_NWRITE}/" \ 177 -e "s/${PAT_RESTAR}/ ln_rstart=${ORCA_LRSTAR}/" \ 178 -e "s/${PAT_NPRG}/ nprg=${ORCA_NWRITE}/" \ 179 -e "s/${PAT_NSTOCK}/ nstock=${ORCA_NSTOCK}/" \ 180 -e "s/${PAT_NRSTAR}/ nrstdt=${ORCA_NRSTDT}/" \ 181 -e "s/${PAT_NDATE0}/ ndate0=${PeriodDateBegin}/" \ 182 -e "s/${PAT_NMSH}/ nmsh=${ORCA_NMSH}/" \ 176 -e "s/${PAT_NIT000}/ nit000=${ORCA_NIT000}/" \ 177 -e "s/${PAT_NITEND}/ nitend=${ORCA_NITEND}/" \ 178 -e "s/${PAT_NWRITE}/ nwrite=${ORCA_NWRITE}/" \ 179 -e "s/${PAT_RESTAR}/ ln_rstart=${ORCA_LRSTAR}/" \ 180 -e "s/${PAT_NPRG}/ nprg=${ORCA_NWRITE}/" \ 181 -e "s/${PAT_NSTOCK}/ nstock=${ORCA_NSTOCK}/" \ 182 -e "s/${PAT_NRSTAR}/ nrstdt=${ORCA_NRSTDT}/" \ 183 -e "s/${PAT_NDATE0}/ ndate0=${PeriodDateBegin}/" \ 184 -e "s%${PAT_NLEAPY}% nleapy=${ORCA_NLEAPY}%" \ 185 -e "s/${PAT_NMSH}/ nmsh=${ORCA_NMSH}/" \ 183 186 namelist > namelist.tmp 184 187 185 188 IGCM_sys_Mv namelist.tmp namelist 189 190 echo 'Variables automatically updated in namelist' 191 grep AUTO namelist 186 192 187 193 IGCM_debug_PopStack "OCE_Update" … … 197 203 DATE_OPA=$( cat date.file | \ 198 204 sed "s/\ ${config_UserChoices_JobName}_[0-9]*[a-z]_\([0-9]*_[0-9]*\)_\ */\1/g" ) 199 ###DATE_OPA=$( $DATE_OPABF | cut -c1-30 )200 205 MainPrefix=${config_UserChoices_JobName}_1d_${DATE_OPA} 201 206 SecondPrefix=${config_UserChoices_JobName}_5d_${DATE_OPA} -
CONFIG/IPSLCM/IPSLCM5/trunk/EXP00/PARAM/namelist_ORCA2
r555 r568 24 24 !----------------------------------------------------------------------- 25 25 no = 0 ! job number 26 cexper = "ORCA2" ! experience name26 cexper = "ORCA2" ! AUTO - experience name 27 27 cn_ocerst_in = "restartopa"! suffix of ocean restart name (input) 28 28 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 = 1 ! restart control = 0 nit000 is not compared to the restart file value 31 ! = 1 use ndate0 in namelist (not the value in the restart file) 32 ! = 2 calendar parameters read in the restart file 33 nit000 = 76 ! first time step 34 nitend = 150 ! last time step 35 ndate0 =18600106 ! initial calendar date yymmdd (used if nrstdt=1) 36 nleapy = 30 ! Leap year calendar (1) or not (0) 37 ninist = 0 ! output the initial state (1) or not (0) 38 nstock = 5400 ! frequency of creation of a restart file 39 nwrite = 1 ! frequency of write in the output file 40 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 29 ln_rstart = .false. ! AUTO - start from rest (F) or from a restart file (T) 30 nrstdt = 0 ! AUTO - restart control = 0 nit000 is not compared to the restart file value 31 ! = 1 use ndate0 in namelist (not the value in the restart file) 32 ! = 2 calendar parameters read in the restart file 33 nit000 = 1 ! AUTO - first time step 34 nitend = 5475 ! AUTO - last time step 35 ndate0 = 010101 ! AUTO - initial calendar date yymmdd (used if nrstdt=1) 36 nleapy = 0 ! AUTO - Leap year calendar (1) or not (0) 37 ninist = 0 ! output the initial state (1) or not (0) 38 nstock = 5475 ! AUTO - frequency of creation of a restart file (modulo referenced to 1) 39 nwrite = 5475 ! AUTO - 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 = .true. ! mask land points in NetCDF outputs (costly: + ~15%) 41 42 / 42 43 !!====================================================================== … … 51 52 &nam_zgr ! vertical coordinate 52 53 !----------------------------------------------------------------------- 53 ln_zco = .true. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined)54 ln_zps = .false. 54 ln_zco = .true. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 55 ln_zps = .false. ! z-coordinate - partial steps (T/F) 55 56 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 56 57 / … … 70 71 e3zps_min = 5. ! the thickness of the partial step is set larger than the minimum 71 72 e3zps_rat = 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 72 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0)73 nmsh = 0 ! AUTO - create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 73 74 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 74 75 ! =0, no acceleration, rdt = rdttra … … 104 105 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 105 106 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 106 ln_blk_clio = .false. 107 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 107 108 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 108 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl )109 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 109 110 nn_ice = 2 ! =0 no ice boundary condition , 110 111 ! =1 use observed ice-cover , … … 115 116 ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) 116 117 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 117 ln_ssr = .false. 118 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 118 119 nn_fwb = 0 ! FreshWater Budget: =0 unchecked , 119 120 ! =1 global mean of e-p-r set to zero at each nn_fsbc time step , … … 133 134 &namsbc_flx ! surface boundary condition : flux formulation 134 135 !----------------------------------------------------------------------- 135 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 136 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 137 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' 138 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' 139 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' 140 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' 141 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' 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 143 ! 143 144 cn_dir = './' ! root directory for the location of the flux files … … 146 147 &namsbc_clio ! namsbc_clio CLIO bulk formulea 147 148 !----------------------------------------------------------------------- 148 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 149 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 150 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. , 'yearly' 151 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. , 'yearly' 152 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. , 'yearly' 153 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. , 'yearly' 154 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. , 'yearly' 155 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' 156 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' 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' , '' , '' 157 158 ! 158 159 cn_dir = './' ! root directory for the location of the bulk files are … … 161 162 &namsbc_core ! namsbc_core CORE bulk formulea 162 163 !----------------------------------------------------------------------- 163 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 164 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 165 sn_wndi = 'u10 ' , 24. , 'U_10_MOD' , .false. , .true. , 'yearly'166 sn_wndj = 'v10 ' , 24. , 'V_10_MOD' , .false. , .true. , 'yearly'167 sn_qsr = ' rad' , 24. , 'SWDN_MOD' , .false. , .true. , 'yearly'168 sn_qlw = ' rad' , 24. , 'LWDN_MOD' , .false. , .true. , 'yearly'169 sn_tair = 't 10' , 24. , 'T_10_MOD' , .false. , .true. , 'yearly'170 sn_humi = 'q 10' , 24. , 'Q_10_MOD' , .false. , .true. , 'yearly'171 sn_prec = 'precip ' , -1. , 'RAIN' , .false. , .true. , 'yearly'172 sn_snow = ' precip' , -1. , 'SNOW' , .false. , .true. , 'yearly'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', '' 173 174 ! 174 175 cn_dir = './' ! root directory for the location of the bulk files … … 211 212 &namsbc_rnf ! runoffs namelist surface boundary condition 212 213 !----------------------------------------------------------------------- 213 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 214 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 215 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' 216 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' 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' , '' , '' 217 218 ! 218 219 cn_dir = './' ! root directory for the location of the runoff files 219 ln_rnf_emp = .true. ! runoffs included into precipitation field (T) or into a file (F)220 ln_rnf_mouth = .true. ! specific treatment at rivers mouths220 ln_rnf_emp = .true. ! runoffs included into precipitation field (T) or into a file (F) 221 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 221 222 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used 222 223 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] … … 225 226 &namsbc_ssr ! surface boundary condition : sea surface restoring 226 227 !----------------------------------------------------------------------- 227 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! 228 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! 229 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' 230 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' 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' , '' , '' 231 232 ! 232 233 cn_dir = './' ! root directory for the location of the runoff files … … 430 431 &nam_dynadv ! formulation of the momentum advection 431 432 !----------------------------------------------------------------------- 432 ln_dynadv_vec = .false. 433 ln_dynadv_cen2= .true. ! flux form - 2nd order centered scheme433 ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 434 ln_dynadv_cen2= .true. ! flux form - 2nd order centered scheme 434 435 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 435 436 / … … 445 446 &nam_dynhpg ! Hydrostatic pressure gradient option 446 447 !----------------------------------------------------------------------- 447 ln_hpg_zco = .true. ! z-coordinate - full steps448 ln_hpg_zps = .false. 448 ln_hpg_zco = .true. ! z-coordinate - full steps 449 ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) 449 450 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 450 451 ln_hpg_hel = .false. ! s-coordinate (helsinki modification) … … 657 658 !----------------------------------------------------------------------- 658 659 ngap = 15 ! time-step frequency of model-data gap computation 659 nprg = 10 ! time-step frequency of gap print in model output660 nprg = 10 ! AUTO - time-step frequency of gap print in model output 660 661 / 661 662 !-----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.