Ignore:
Timestamp:
02/18/09 12:30:42 (15 years ago)
Author:
aclsce
Message:

Modified IPSLCM5_WORK configuration to be coherent with IPSLCM5
configuration

Location:
CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00
Files:
3 added
10 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/COMP/lmdz.card

    r396 r561  
    1616        (${R_INIT}/ATM/${config_UserChoices_TagName}/Rugos.nc, .), \ 
    1717        (${R_INIT}/ATM/${config_UserChoices_TagName}/landiceref.nc, .), \ 
    18         (${R_INIT}/ATM/${config_UserChoices_TagName}/${RESOL}/o2a.nc, .)  
     18        (${R_INIT}/ATM/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/o2a.nc, .)  
    1919 
    2020[BoundaryFiles] 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/COMP/nemo.card

    r396 r561  
    44[UserChoices] 
    55OPA_NPDT_JOURS=15 
     6ORCA_version=ORCA2.1 
    67 
    78[InitialStateFiles] 
     
    1011[BoundaryFiles] 
    1112List=   () 
    12 ListNonDel= (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/ahmcoef, .), \ 
    13             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/bathy_level.nc, .), \ 
    14             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/bathy_meter.nc, .), \ 
    15             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/bathy_updated.nc, .), \ 
    16             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/coordinates.nc, .), \ 
    17             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/data_1m_potential_temperature_nomask.nc, .), \ 
    18             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/data_1m_salinity_nomask.nc, .), \ 
    19             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/geothermal_heating.nc, .),\ 
    20             (${R_BC}/OCE/${config_UserChoices_TagName}/${RESOL_OCE}/runoff_1m_nomask.nc, .),\ 
     13ListNonDel= (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/ahmcoef, .), \ 
     14            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/bathy_level.nc, .), \ 
     15            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/bathy_meter.nc, .), \ 
     16            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/coordinates.nc, .), \ 
     17            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/data_1m_potential_temperature_nomask.nc, .), \ 
     18            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/data_1m_salinity_nomask.nc, .), \ 
     19            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/geothermal_heating.nc, .),\ 
     20            (${R_BC}/OCE/${config_UserChoices_TagName}/\${${compname}_UserChoices_ORCA_version}/runoff_1m_nomask.nc, .),\ 
    2121 
    2222[ParametersFiles] 
     
    2727 
    2828[OutputText] 
    29 List=   (ocean.output, opa.xx.prt, solver.stat) 
     29List=   (ocean.output, oceanx.prt, solver.stat) 
    3030 
    3131[OutputFiles] 
     
    4545Patches = () 
    4646GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter) 
    47 TimeSeriesVars = (iowaflup, sohtc300, sohefldo, soicecov, somxl010, sorunoff, sosaline, sossheig, sosstsst, sowaflep, sowaflcd, sowaflup) 
     47TimeSeriesVars = (sohtc300, sohefldo, soicecov, somxl010, sosaline, sossheig, sosstsst, sowaflcd, sowaflup) 
    4848 
    4949[Post_1M_grid_U] 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/COMP/nemo.driver

    r396 r561  
    1414    fi 
    1515 
    16     ##--Transformation namelist fortran 77 --> fortran 90 
    17     #   ?       ? Is it really necessary ?       ? 
    18     #sed -e /:/d -e s/" \&END"/"\/"/ -e s/" \&"/"\&"/ ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} > ${SUBMIT_DIR}/PARAM/namelist 
    19  
     16# Local function to find namelists parameters 
     17        supergrep () { 
     18            grep "^ *$1 *=" $2 | sed -e "s% *\!.*%%" 
     19        } 
     20         
    2021    ##--Variables used by OPA -- 
    2122 
     
    3233    # nmsh  =1 create a mesh file (coordinates, scale factors, masks) 
    3334 
    34     PAT_NO=$(     grep 'no *='     ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v off | grep -v \!) 
    35     PAT_CEXPER=$( grep 'cexper *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    36     PAT_NIT000=$( grep 'nit000 *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    37     PAT_NITEND=$( grep 'nitend *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    38     PAT_NBISEX=$( grep 'nbisex *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE}| sed 's/ nbisex  =//' |sed '/^[ ]$/d'|sed 's/,//' | grep -v \!) 
    39     PAT_NWRITE=$( grep 'nwrite *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    40     PAT_RESTAR=$( grep 'ln_rstart *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    41     PAT_NSTOCK=$( grep 'nstock *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    42     PAT_NRSTAR=$( grep 'nrstdt *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    43     PAT_NPRG=$(   grep 'nprg *='   ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    44     PAT_NDATE0=$( grep 'ndate0 *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!) 
    45     PAT_NMSH=$(   grep 'nmsh *='   ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | grep -v \!)  
    46 #######    PAT_NWRIHF=$( grep 'nwrihf' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
    47  
    48     ORCA_NSTOCK=$( echo $PAT_NSTOCK | sed "s/[a-z,A-Z,=]//g" ) 
    49     ORCA_RDT=$( egrep 'rdt *=' ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | sed 's/ *rdt *=//' | sed 's/\.//' | grep -v \!) 
     35            PAT_CEXPER=$( supergrep cexper    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     36            PAT_NIT000=$( supergrep nit000    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     37            PAT_NITEND=$( supergrep nitend    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     38            PAT_NWRITE=$( supergrep nwrite    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     39            PAT_RESTAR=$( supergrep ln_rstart ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     40            PAT_NSTOCK=$( supergrep nstock    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     41            PAT_NRSTAR=$( supergrep nrstdt    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     42            PAT_NDATE0=$( supergrep ndate0    ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     43            PAT_NMSH=$(   supergrep nmsh      ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} ) 
     44         
     45            ORCA_NSTOCK=$( echo $PAT_NSTOCK | sed "s/[a-z,A-Z,=]//g" ) 
     46            ORCA_RDT=$( supergrep rdt  ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} | sed 's/ *rdt *=//' | sed 's/\. *//' ) 
     47 
     48    ##--Transformation namelist fortran 77 --> fortran 90 
     49    #   ?       ? Is it really necessary ?       ? 
     50    #sed -e /:/d -e s/" \&END"/"\/"/ -e s/" \&"/"\&"/ ${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE} > ${SUBMIT_DIR}/PARAM/namelist 
    5051 
    5152    # Period Length In Days between DateBegin and first day of calendar 0001 01 01 
     
    8889            (( ORCA_NWRITE = 0 ))                                                            ;; 
    8990    esac 
    90 ###            (( ORCA_NWRITE = 1 ))  ; 
    91 ###            PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}d ; 
    92 ###            R_OUT_OCE_NWRITE=${R_OUT_OCE_O_D} ; 
    93 ###            R_OUT_ICE_NWRITE=${R_OUT_ICE_O_D} ; 
    94 ###            WF1=${factor}D ; 
    95  
    9691 
    9792    ##--Write Frequency Purpose .... concerning the second opa file type. 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/COMP/oasis.card

    r396 r561  
    33          
    44[UserChoices] 
     5Output_mode=n 
    56 
    67[InitialStateFiles] 
    7 List=   (${R_INIT}/CPL/${config_UserChoices_TagName}/${RESOL}/flxat.nc, flxat.nc), \ 
    8         (${R_INIT}/CPL/${config_UserChoices_TagName}/${RESOL}/sstoc.nc, sstoc.nc) 
     8List=   (${R_INIT}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/flxat.nc, flxat.nc), \ 
     9        (${R_INIT}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/sstoc.nc, sstoc.nc) 
    910 
    1011[BoundaryFiles] 
    1112List=   () 
    12 ListNonDel= (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/grids.nc,                grids.nc),\ 
    13             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/masks.nc,                masks.nc),\ 
    14             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/areas.nc,                areas.nc),\ 
    15             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o,        wa2o.flx),\ 
    16             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wo2a,        wo2a.tsg),\ 
    17             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.runoff, wa2o.run),\ 
    18             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.calvin, wa2o.cal) 
     13ListNonDel= (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/grids.nc, grids.nc),\ 
     14            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/masks.nc, masks.nc),\ 
     15            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/areas.nc, areas.nc),\ 
     16            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/mozaic.wa2o, wa2o.flx),\ 
     17            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/mozaic.wo2a, wo2a.tsg),\ 
     18            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/mozaic.wa2o.runoff, wa2o.run),\ 
     19            (${R_BC}/CPL/${config_UserChoices_TagName}/\${nemo_UserChoices_ORCA_version}x${RESOL_ATM}/mozaic.wa2o.calvin, wa2o.cal) 
    1920 
    2021[ParametersFiles] 
     
    3031 
    3132[OutputFiles] 
    32 List=   (cpl_oce_tau.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_tau.nc,    NONE),\ 
    33         (cpl_oce_flx.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_flx.nc,    NONE),\ 
    34         (cpl_oce_sst.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_sst.nc,    NONE) 
     33List=   (SISUTESW_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_SISUTESW.nc,    NONE),\ 
     34        (SIICECOV_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_SIICECOV.nc,    NONE),\ 
     35        (SIICTEMW_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_SIICTEMW.nc,    NONE),\ 
     36        (SIICEALW_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_SIICEALW.nc,    NONE),\ 
     37        (COTAUXXU_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUXXU.nc,    NONE),\ 
     38        (COTAUXXV_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUXXV.nc,    NONE),\ 
     39        (COTAUYYU_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUYYU.nc,    NONE),\ 
     40        (COTAUYYV_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUYYV.nc,    NONE),\ 
     41        (COTAUZZU_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUZZU.nc,    NONE),\ 
     42        (COTAUZZV_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COTAUZZV.nc,    NONE),\ 
     43        (COWINDSP_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COWINDSP.nc,    NONE),\ 
     44        (COEMPSIC_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COEMPSIC.nc,    NONE),\ 
     45        (CONESOPR_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_CONESOPR.nc,    NONE),\ 
     46        (COEMPOCE_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COEMPOCE.nc,    NONE),\ 
     47        (COICEVAP_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COICEVAP.nc,    NONE),\ 
     48        (COSHFOCE_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COSHFOCE.nc,    NONE),\ 
     49        (CONSFOCE_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_CONSFOCE.nc,    NONE),\ 
     50        (COSHFICE_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COSHFICE.nc,    NONE),\ 
     51        (CONSFICE_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_CONSFICE.nc,    NONE),\ 
     52        (CODFLXDT_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_CODFLXDT.nc,    NONE),\ 
     53        (COCALVIN_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COCALVIN.nc,    NONE),\ 
     54        (COLIQRUN_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_COLIQRUN.nc,    NONE),\ 
     55        (O_SSTSST_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_SSTSST.nc,    NONE),\ 
     56        (OIceFrac_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OIceFrac.nc,    NONE),\ 
     57        (O_TepIce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_TepIce.nc,    NONE),\ 
     58        (O_AlbIce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_AlbIce.nc,    NONE),\ 
     59        (O_OTaux1_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTaux1.nc,    NONE),\ 
     60        (O_OTaux2_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTaux2.nc,    NONE),\ 
     61        (O_OTauy1_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTauy1.nc,    NONE),\ 
     62        (O_OTauy2_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTauy2.nc,    NONE),\ 
     63        (O_OTauz1_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTauz1.nc,    NONE),\ 
     64        (O_OTauz2_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_OTauz2.nc,    NONE),\ 
     65        (O_Wind10_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_Wind10.nc,    NONE),\ 
     66        (OSubMPre_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OSubMPre.nc,    NONE),\ 
     67        (OISubMSn_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OISubMSn.nc,    NONE),\ 
     68        (OOEvaMPr_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OOEvaMPr.nc,    NONE),\ 
     69        (OIceEvap_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OIceEvap.nc,    NONE),\ 
     70        (O_QsrOce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_QsrOce.nc,    NONE),\ 
     71        (O_QsrIce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_QsrIce.nc,    NONE),\ 
     72        (O_QnsOce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_QnsOce.nc,    NONE),\ 
     73        (O_QnsIce_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_QnsIce.nc,    NONE),\ 
     74        (O_dQnsdT_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_dQnsdT.nc,    NONE),\ 
     75        (O_Runoff_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_O_Runoff.nc,    NONE),\ 
     76        (OCalving_out.${DATE_OASIS}.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_OCalving.nc,    NONE) 
    3577 
    3678[Post_cpl_oce_tau] 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/COMP/oasis.driver

    r396 r561  
    2121CPL_Update () 
    2222{ 
    23  
    2423    IGCM_debug_PushStack "CPL_Update" 
    2524 
     
    2928 
    3029############ Update Other Parameters  ############## 
     30## Read Output_mode option in ${compname}.card --> Put "y" to have coupling variables in ouptut files 
     31IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices Output_mode  
     32eval Output_mode=\${${compname}_UserChoices_Output_mode} > /dev/null 2>&1 
     33echo ${Output_mode} 
     34Output_var='EXPORTED' 
     35if [ ${Output_mode} = y ] ; then 
     36    Output_var='EXPOUT' 
     37fi 
    3138 
    3239############ Update Parameter Files   ############## 
     
    3643    -e "s/<nproc_atm>/${NUM_PROC_ATM}/g" \ 
    3744    -e "s/<nproc_oce>/${NUM_PROC_OCE}/g" \ 
     45    -e "s/<output_mode>/${Output_var}/g" \ 
    3846        namcouple > namcouple.tmp 
    3947 
     
    5159CPL_Finalize () 
    5260{ 
     61    IGCM_debug_PushStack "CPL_Finalize" 
    5362 
    54     IGCM_debug_PushStack "CPL_Finalize" 
     63    DATE_OASIS=$(echo ${PeriodDateBegin} | sed -e "{ s/^\(.\{4\}\)/\1-/;s/\(-.\{2\}\)/\1-/; }"); 
     64    DATE_OASIS=${DATE_OASIS}"T00:00:00" 
    5565 
    5666    echo FINALIZE CPL !!! 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namcouple_ORCA2xLMD144142

    r409 r561  
    2929# 
    3030 $CHANNEL 
    31            <lib_mpi> 
     31<lib_mpi> 
    3232<nproc_atm>  <nproc_atm> 
    3333<nproc_oce>  <nproc_oce>         
     
    5959# 
    6060 $RUNTIME 
    61              <total_time> 
     61 <total_time> 
    6262 $END 
    6363############################################################################### 
     
    6767# 
    6868 $INIDATE 
    69              <initial_date> 
     69 <initial_date> 
    7070 $END 
    7171############################################################################### 
     
    108108# Field 1 : Weighted sea surface temperature (o->a 1) 
    109109# 
    110 O_SSTSST SISUTESW 1 86400  5  sstoc.nc  EXPORTED  
     110O_SSTSST SISUTESW 1 86400  5  sstoc.nc  <output_mode>  
    111111torc  tlmd  LAG=5760 
    112112P  2 P  0 
    113113LOCTRANS CHECKIN MOZAIC CHECKOUT REVERSE 
     114# LOCTRANS: AVERAGE to average value over coupling period 
    114115 AVERAGE 
    115116# CHECKIN: indicate computation of global, land and sea field integrals. 
     
    126127# Field 2 : Sea ice extent  (o->a 2) 
    127128# 
    128 OIceFrac SIICECOV 2 86400  5  sstoc.nc EXPORTED 
     129OIceFrac SIICECOV 2 86400  5  sstoc.nc <output_mode> 
    129130torc  tlmd   LAG=5760 
    130131P  2 P  0 
     
    144145# Field 3 : Weighted Sea Ice Temperature (o->a 3) 
    145146# 
    146 O_TepIce SIICTEMW 34 86400  5  sstoc.nc  EXPORTED 
     147O_TepIce SIICTEMW 34 86400  5  sstoc.nc  <output_mode> 
    147148torc  tlmd   LAG=5760 
    148149P  2 P  0 
     
    161162# Field 4 : Weighted Sea ice albedo (o->a 4) 
    162163# 
    163 O_AlbIce  SIICEALW 17 86400  5  sstoc.nc  EXPORTED 
     164O_AlbIce  SIICEALW 17 86400  5  sstoc.nc  <output_mode> 
    164165torc  tlmd  LAG=5760   
    165166P  2 P  0 
     
    182183# Field 5 : wind stress along X axis (a->o tau 1) 
    183184# 
    184 COTAUXXU O_OTaux1 23 86400  6  flxat.nc   EXPORTED 
     185COTAUXXU O_OTaux1 23 86400  6  flxat.nc   <output_mode> 
    185186tlmd    uorc  LAG=1800  
    186187P  0 P  2 
     
    202203# Field 6 : stress along Y axis (a->o tau 1) 
    203204# 
    204 COTAUYYU O_OTauy1 23 86400  6  flxat.nc  EXPORTED 
     205COTAUYYU O_OTauy1 23 86400  6  flxat.nc  <output_mode> 
    205206tlmd    uorc  LAG=1800  
    206207P  0 P  2 
     
    222223# Field 7 :wind stress along Z axis (a->o tau 1) 
    223224# 
    224 COTAUZZU O_OTauz1 23 86400  6  flxat.nc  EXPORTED 
     225COTAUZZU O_OTauz1 23 86400  6  flxat.nc  <output_mode> 
    225226tlmd    uorc  LAG=1800  
    226227P  0 P  2 
     
    242243# Field 8 : wind stress along X axis 2 (a->o tau 2) 
    243244# 
    244 COTAUXXV O_OTaux2 24 86400  6    flxat.nc  EXPORTED 
     245COTAUXXV O_OTaux2 24 86400  6    flxat.nc  <output_mode> 
    245246tlmd    vorc  LAG=1800  
    246247P  0 P  2 
     
    262263# Field 9 : wind stress along Y axis 2 (a->o tau 2)  
    263264# 
    264 COTAUYYV O_OTauy2 24 86400  6    flxat.nc  EXPORTED 
     265COTAUYYV O_OTauy2 24 86400  6    flxat.nc  <output_mode> 
    265266tlmd    vorc  LAG=1800 
    266267P  0 P  2 
     
    282283# Field 10 :  wind stress along Z axis 2 (a->o tau 2)  
    283284# 
    284 COTAUZZV O_OTauz2 24 86400  6    flxat.nc  EXPORTED 
     285COTAUZZV O_OTauz2 24 86400  6    flxat.nc  <output_mode> 
    285286tlmd    vorc  LAG=1800 
    286287P  0 P  2 
     
    299300 INT=1 
    300301## 
    301 ########################################################################## 
    302 # Field 11 : 10-m wind speed (a->o tau 2)  
    303 # 
    304 COWINDSP O_Wind10 36 86400  6    flxat.nc EXPORTED 
    305 tlmd  torc LAG=1800 
    306 P  0 P  2 
    307 INVERT CHECKIN MASK EXTRAP INTERP CHECKOUT 
    308 # Reverse 
    309 NORSUD WSTEST 
    310 # CHECKIN: indicate computation of global, land and sea field integrals. 
    311  INT=1 
    312 # Mask value 
    313 888888. 
    314 # Extrapolation method 
    315 NINENN 2 1 1 
    316 # Interpolation method 
    317 BILINEAR  Z   SCALAR 
    318 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    319  INT=1 
    320 # 
    321302######################################################################### 
    322 # Field 12 : tprecip = ( rain_ice + snow_ice ) - evap_ice (a->o flx 7) 
    323 # 
    324 COPEFWAT OIPr-Sub 25 86400  5   flxat.nc   EXPORTED 
    325 tlmd    torc  LAG=1800  
    326 P  0 P  2 
    327 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
     303# Field 11 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 
     304# 
     305COEMPSIC OSubMPre 25 86400  4   flxat.nc   <output_mode> 
     306tlmd    torc  LAG=1800  
     307P  0 P  2 
     308INVERT CHECKIN MOZAIC CHECKOUT 
    328309# Reverse 
    329310NORSUD WSTEST 
     
    332313# Interpolation method ou parametres mozaic 
    333314wa2o.flx 92  2   32 
    334 # change of sign 
    335 -1.0 0 
    336 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    337  INT=1 
    338 # 
    339 ########################################################################## 
    340 # Field 13 :  sprecip = snow_ice - evap_ice (o->a flx 6) 
    341 # 
    342 COPEFICE  OISn-Sub 25 86400  5   flxat.nc   EXPORTED 
    343 tlmd    torc  LAG=1800  
    344 P  0 P  2 
    345 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
     315# CHECKOUT: indicate computation of global, land and sea field integrals. 
     316 INT=1 
     317# 
     318########################################################################## 
     319# Field 12 :  sprecip = evap_ice - snow_ice (o->a flx 6) 
     320# 
     321CONESOPR  OISubMSn 25 86400  4   flxat.nc   <output_mode> 
     322tlmd    torc  LAG=1800  
     323P  0 P  2 
     324INVERT CHECKIN MOZAIC CHECKOUT 
    346325# Reverse 
    347326NORSUD WSTEST 
     
    350329# Interpolation method ou parametres mozaic 
    351330wa2o.flx 92  2    32 
    352 # change of sign 
    353 -1.0 0 
    354 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    355  INT=1 
    356 # 
    357 ########################################################################## 
    358 # Field 14: emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 
    359 # 
    360 COTOSPSU  OOEv-OPr 28 86400  4   flxat.nc   EXPORTED 
     331# CHECKOUT: indicate computation of global, land and sea field integrals. 
     332 INT=1 
     333# 
     334########################################################################## 
     335# Field 13 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 
     336# 
     337COEMPOCE  OOEvaMPr  28 86400  4   flxat.nc   <output_mode> 
    361338tlmd    torc  LAG=1800  
    362339P  0 P  2 
     
    368345# Interpolation method ou parametres mozaic 
    369346wa2o.flx 92  2    32 
    370 # change of sign 
    371 ###-1.0 0 
    372 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    373  INT=1 
     347# CHECKOUT: indicate computation of global, land and sea field integrals. 
     348 INT=1 
     349# 
     350########################################################################## 
     351# Field 14 : Evaporation Ice  (a->o flx 8) 
     352# 
     353COICEVAP OIceEvap 27 86400  4   flxat.nc   <output_mode> 
     354tlmd    torc  LAG=1800  
     355P  0 P  2 
     356INVERT CHECKIN MOZAIC CHECKOUT 
     357# Reverse 
     358NORSUD WSTEST 
     359# CHECKIN: indicate computation of global, land and sea field integrals. 
     360 INT=1 
     361# Interpolation method ou parametres mozaic 
     362wa2o.flx 92  2    32 
     363# CHECKOUT: indicate computation of global, land and sea field integrals. 
     364 INT=1 
     365# 
    374366########################################################################## 
    375367# Field 15 : solar heat flux on ocean (a->o flx 2) 
    376368# 
    377 COSWFLDO O_QsrOce  7 86400 4   flxat.nc  EXPORTED 
     369COSHFOCE O_QsrOce  7 86400 4   flxat.nc  <output_mode> 
    378370tlmd    torc  LAG=1800 
    379371P  0 P  2 
     
    387379# CHECKOUT: indicate computation of global, land and sea field integrals. 
    388380 INT=1 
     381# 
    389382######################################################################### 
    390383# Field 16 : Non solar heat flux on ocean (a->o flx 4) 
    391384# 
    392 CONSFLDO O_QnsOce 6 86400  4   flxat.nc  EXPORTED 
     385CONSFOCE O_QnsOce 6 86400  4   flxat.nc  <output_mode> 
    393386tlmd    torc   LAG=1800 
    394387P  0 P  2 
     
    402395# CHECKOUT: indicate computation of global, land and sea field integrals. 
    403396 INT=1 
     397# 
    404398########################################################################## 
    405399# Field 17 : solar heat flux on ice (a->o flx 1) 
    406400# 
    407 COSHFLIC O_QsrIce  7 86400 4   flxat.nc  EXPORTED 
     401COSHFICE O_QsrIce  7 86400 4   flxat.nc  <output_mode> 
    408402tlmd    torc  LAG=1800 
    409403P  0 P  2 
     
    421415# Field 18 : Non solar heat flux on ice (a->o flx 3) 
    422416# 
    423 CONSFLIC O_QnsIce 6 86400  4  flxat.nc  EXPORTED 
     417CONSFICE  O_QnsIce 6 86400  4  flxat.nc  <output_mode> 
    424418tlmd    torc  LAG=1800 
    425419P  0 P  2 
     
    437431# Field 19 : Non solar heat flux derivative (o->a flx 5) 
    438432# 
    439 CODFLXDT O_dQnsdT  35 86400  4   flxat.nc  EXPORTED 
     433CODFLXDT O_dQnsdT  35 86400  4   flxat.nc  <output_mode> 
    440434tlmd    torc  LAG=1800 
    441435P  0 P  2 
     
    449443# CHECKOUT: indicate computation of global, land and sea field integrals. 
    450444 INT=1 
     445# 
    451446########################################################################### 
    452447# Field 20 : iceberg calving (a->o flx 12) 
    453448# 
    454 CRWOCEIS OCalving  32 86400  5  flxat.nc  EXPORTED 
    455 tlmd    torc  LAG=1800  
    456 P  0 P  2 
    457 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
     449COCALVIN OCalving  32 86400  4  flxat.nc  <output_mode> 
     450tlmd    torc  LAG=1800  
     451P  0 P  2 
     452INVERT CHECKIN MOZAIC CHECKOUT 
    458453# Reverse 
    459454NORSUD WSTEST 
     
    462457# Interpolation method ou parametres mozaic 
    463458wa2o.cal  98   5   71 
    464 # change of sign 
    465 -1.0 0 
    466 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    467  INT=1 
     459# CHECKOUT: indicate computation of global, land and sea field integrals. 
     460 INT=1 
     461# 
    468462############################################################################ 
    469463# Field 21 : run-off (a->o flx 11) 
    470464# 
    471 CRWOCERD  O_Runoff 32 86400  5  flxat.nc   EXPORTED 
     465COLIQRUN  O_Runoff 32 86400  5  flxat.nc   <output_mode> 
    472466tlmd    torc  LAG=1800  
    473467P  0 P  2 
     
    481475# change units from m/s to kg/s, change of sign 
    482476# weights convert from kg/s to kg/m^2/s 
    483 -1000.0 0 
    484 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    485  INT=1 
    486 # 
    487 ############################################################################ 
    488 # 
    489  $END 
     4771000.0 0 
     478# CHECKOUT: indicate computation of global, land and sea field integrals. 
     479 INT=1 
     480# 
     481############################################################################ 
     482# 
     483 $END 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namcouple_ORCA2xLMD9671

    r396 r561  
    5555 $END 
    5656############################################################################### 
    57 #  
     57# 
    5858# RUNTIME: total simulated time for the actual run in seconds (<I8) 
    5959# 
     
    6464# 
    6565# INIDATE (I8) 
    66 #         initial date of the run (YYYYMMDD)          
     66#         initial date of the run (YYYYMMDD) 
    6767# 
    6868 $INIDATE 
     
    105105#                      OCEAN  --->>>  ATMOS 
    106106#                      -------------------- 
    107 # Field 1 : Sea surface temperature (o->a 1) 
    108 O_SSTSST SISUTESW 1 86400  5  sstoc.nc  EXPORTED  
     107# Field 1 : Weighted Sea surface temperature (o->a 1) 
     108# 
     109O_SSTSST SISUTESW 1 86400  5  sstoc.nc  <output_mode>  
    109110torc  tlmd  LAG=5760 
    110111P  2 P  0 
    111112LOCTRANS CHECKIN MOZAIC CHECKOUT REVERSE 
     113# LOCTRANS: AVERAGE to average value over coupling period 
    112114 AVERAGE 
    113115# CHECKIN: indicate computation of global, land and sea field integrals. 
     
    121123  NORSUD   WSTEST 
    122124 
    123 # 
    124 ############################################################################ 
    125 # 
    126125############################################################################ 
    127126# Field 2 : Sea ice extent  (o->a 2) 
    128127# 
    129 OIceFrac SIICECOV 2 86400  5  sstoc.nc EXPORTED 
     128OIceFrac SIICECOV 2 86400  5  sstoc.nc <output_mode> 
    130129torc  tlmd   LAG=5760 
    131130P  2 P  0 
     
    143142# 
    144143############################################################################ 
    145 # Field 3 : Sea Ice Temperature (o->a 4) 
    146 # 
    147 O_TepIce SIICTEMW 34 86400  5  sstoc.nc  EXPORTED 
     144# Field 3 : Weighted Sea Ice Temperature (o->a 3) 
     145# 
     146O_TepIce SIICTEMW 34 86400  5  sstoc.nc  <output_mode> 
    148147torc  tlmd   LAG=5760 
    149148P  2 P  0 
     
    160159# 
    161160############################################################################ 
    162 # Field 4 : Sea ice albedo (o->a 3) 
    163 # 
    164 O_AlbIce SIICEALW 17 86400  5  sstoc.nc  EXPORTED 
     161# Field 4 : Weighted Sea ice albedo (o->a 4) 
     162# 
     163O_AlbIce SIICEALW 17 86400  5  sstoc.nc  <output_mode> 
    165164torc  tlmd  LAG=5760   
    166165P  2 P  0 
     
    182181#                      ------------------------- 
    183182############################################################################    
    184 # Field 17 : wind stress along X axis (a->o tau 1) 
    185 # 
    186 COTAUXXU O_OTaux1 23 86400  6  flxat.nc   EXPORTED 
     183# 
     184# Field 5 : wind stress along X axis (a->o tau 1) 
     185# 
     186COTAUXXU O_OTaux1 23 86400  6  flxat.nc   <output_mode> 
    187187tlmd    uorc  LAG=1800  
    188188P  0 P  2 
     
    203203########################################################################### 
    204204# 
    205 # 
    206 # Field 18 : stress along Y axis (a->o tau 1) 
    207 # 
    208 COTAUYYU O_OTauy1 23 86400  6  flxat.nc  EXPORTED 
     205# Field 6 : stress along Y axis (a->o tau 1) 
     206# 
     207COTAUYYU O_OTauy1 23 86400  6  flxat.nc  <output_mode> 
    209208tlmd    uorc  LAG=1800  
    210209P  0 P  2 
     
    225224########################################################################### 
    226225# 
    227 # 
    228 # Field 19 :wind stress along Z axis (a->o tau 1) 
    229 # 
    230 COTAUZZU O_OTauz1 23 86400  6  flxat.nc  EXPORTED 
     226# Field 7 :wind stress along Z axis (a->o tau 1) 
     227# 
     228COTAUZZU O_OTauz1 23 86400  6  flxat.nc  <output_mode> 
    231229tlmd    uorc  LAG=1800  
    232230P  0 P  2 
     
    247245########################################################################## 
    248246# 
    249 # Field 20 : wind stress along X axis 2 (a->o tau 2) 
    250 # 
    251 COTAUXXV O_OTaux2 24 86400  6    flxat.nc  EXPORTED 
     247# Field 8 : wind stress along X axis 2 (a->o tau 2) 
     248# 
     249COTAUXXV O_OTaux2 24 86400  6    flxat.nc  <output_mode> 
    252250tlmd    vorc  LAG=1800  
    253251P  0 P  2 
     
    268266########################################################################## 
    269267# 
    270 # Field 21 : wind stress along Y axis 2 (a->o tau 2)  
    271 # 
    272 COTAUYYV O_OTauy2 24 86400  6    flxat.nc  EXPORTED 
     268# Field 9 : wind stress along Y axis 2 (a->o tau 2)  
     269# 
     270COTAUYYV O_OTauy2 24 86400  6    flxat.nc  <output_mode> 
    273271tlmd    vorc  LAG=1800 
    274272P  0 P  2 
     
    289287########################################################################## 
    290288# 
    291 # Field 22 :  wind stress along Z axis 2 (a->o tau 2)  
    292 # 
    293 COTAUZZV O_OTauz2 24 86400  6    flxat.nc  EXPORTED 
     289# Field 10 :  wind stress along Z axis 2 (a->o tau 2)  
     290# 
     291COTAUZZV O_OTauz2 24 86400  6    flxat.nc  <output_mode> 
    294292tlmd    vorc  LAG=1800 
    295293P  0 P  2 
     
    307305# CHECKOUT: indicate computation of global, land and sea field integrals. 
    308306 INT=1 
    309 ## 
     307# 
     308######################################################################### 
     309# 
     310# Field 11 : emp_ice = evap_ice - ( rain_ice + snow_ice ) (a->o flx 7) 
     311# 
     312COEMPSIC OSubMPre 25 86400  4   flxat.nc   <output_mode> 
     313tlmd    torc  LAG=1800  
     314P  0 P  2 
     315INVERT CHECKIN MOZAIC CHECKOUT 
     316# Reverse 
     317NORSUD WSTEST 
     318# CHECKIN: indicate computation of global, land and sea field integrals. 
     319 INT=1 
     320# Interpolation method ou parametres mozaic 
     321wa2o.flx 92  2    14 
     322# CHECKOUT: indicate computation of global, land and sea field integrals. 
     323 INT=1 
     324# 
    310325########################################################################## 
    311326# 
    312 # Field 11 : 10-m wind speed (a->o tau 2)  
    313 # 
    314 COWINDSP O_Wind10 36 86400  6    flxat.nc EXPORTED 
    315 tlmd  torc LAG=1800 
    316 P  0 P  2 
    317 INVERT CHECKIN MASK EXTRAP INTERP CHECKOUT 
    318 # Reverse 
    319 NORSUD WSTEST 
    320 # CHECKIN: indicate computation of global, land and sea field integrals. 
    321  INT=1 
    322 # Mask value 
    323 888888. 
    324 # Extrapolation method 
    325 NINENN 2 1 1 
    326 # Interpolation method 
    327 BILINEAR  Z   SCALAR 
     327# Field 12 :  sprecip = evap_ice - snow_ice (o->a flx 6) 
     328# 
     329CONESOPR  OISubMSn 25 86400  4   flxat.nc   <output_mode> 
     330tlmd    torc  LAG=1800  
     331P  0 P  2 
     332INVERT CHECKIN MOZAIC CHECKOUT 
     333# Reverse 
     334NORSUD WSTEST 
     335# CHECKIN: indicate computation of global, land and sea field integrals. 
     336 INT=1 
     337# Interpolation method ou parametres mozaic 
     338wa2o.flx 92  2    14 
    328339# CHECKOUT: indicate computation of global, land and sea field integrals. 
    329340 INT=1 
     
    331342######################################################################### 
    332343# 
    333 # Field 10 : tprecip = ( rain_ice + snow_ice ) - evap_ice  
    334 ###  (a->o flx 7) 
    335 # 
    336 COPEFWAT OIPr-Sub 25 86400  5   flxat.nc   EXPORTED 
     344# Field 13 : emp = emp_oce = evap_oce - ( rain_oce + snow_oce ) (a->o flx 9) 
     345# 
     346COEMPOCE  OOEvaMPr 29 86400  4   flxat.nc   <output_mode> 
    337347tlmd    torc  LAG=1800  
    338348P  0 P  2 
    339 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    340 # Reverse 
    341 NORSUD WSTEST 
    342 # CHECKIN: indicate computation of global, land and sea field integrals. 
    343  INT=1 
    344 # Interpolation method ou parametres mozaic 
    345 wa2o.flx 92  2    14 
    346 # change of sign 
    347 -1.0 0 
     349INVERT CHECKIN MOZAIC CHECKOUT 
     350# Reverse 
     351NORSUD WSTEST 
     352# CHECKIN: indicate computation of global, land and sea field integrals. 
     353 INT=1 
     354# Interpolation method ou parametres mozaic 
     355wa2o.flx 92  2    14 
    348356# CHECKOUT: indicate computation of global, land and sea field integrals. 
    349357 INT=1 
     
    351359########################################################################## 
    352360# 
    353 # Field 11 :  sprecip = snow_ice - evap_ice (o->a flx 6) 
    354 # 
    355 COPEFICE  OISn-Sub 25 86400  5   flxat.nc   EXPORTED 
     361# Field 14 : Evaporation Ice  (a->o flx 8) 
     362# 
     363COICEVAP OIceEvap 41 86400  4   flxat.nc   <output_mode> 
    356364tlmd    torc  LAG=1800  
    357365P  0 P  2 
    358 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    359 # Reverse 
    360 NORSUD WSTEST 
    361 # CHECKIN: indicate computation of global, land and sea field integrals. 
    362  INT=1 
    363 # Interpolation method ou parametres mozaic 
    364 wa2o.flx 92  2    14 
    365 # change of sign 
    366 -1.0 0 
    367 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    368  INT=1 
     366INVERT CHECKIN MOZAIC CHECKOUT 
     367# Reverse 
     368NORSUD WSTEST 
     369# CHECKIN: indicate computation of global, land and sea field integrals. 
     370 INT=1 
     371# Interpolation method ou parametres mozaic 
     372wa2o.flx 92  2    14 
     373# CHECKOUT: indicate computation of global, land and sea field integrals. 
     374 INT=1 
     375# 
     376########################################################################### 
     377# 
     378# Field 15 : solar heat flux on ocean (a->o flx 2) 
     379# 
     380COSHFOCE O_QsrOce  7 86400 4   flxat.nc  <output_mode> 
     381tlmd    torc  LAG=1800 
     382P  0 P  2 
     383INVERT CHECKIN MOZAIC CHECKOUT 
     384# Reverse 
     385NORSUD WSTEST 
     386# CHECKIN: indicate computation of global, land and sea field integrals. 
     387 INT=1 
     388# Interpolation method or mozaic parameters 
     389wa2o.flx 92  2    14 
     390# CHECKOUT: indicate computation of global, land and sea field integrals. 
     391 INT=1 
     392# 
    369393######################################################################### 
    370394# 
    371 # Field 12: emp = emp_oce = evap_oce -  
    372 ### ( rain_oce + snow_oce ) (a->o flx 9) 
    373 # 
    374 COTOSPSU  OOEv-OPr 28 86400  4   flxat.nc   EXPORTED 
    375 tlmd    torc  LAG=1800  
    376 P  0 P  2 
    377 INVERT CHECKIN MOZAIC CHECKOUT 
    378 # Reverse 
    379 NORSUD WSTEST 
    380 # CHECKIN: indicate computation of global, land and sea field integrals. 
    381  INT=1 
    382 # Interpolation method ou parametres mozaic 
    383 wa2o.flx 92  2    14 
    384 # change of sign 
    385 ###-1.0 0 
    386 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    387  INT=1 
    388 ########################################################################## 
    389 # 
    390 # Field 13 : Evaporation Ice  (a->o flx 8) 
    391 # 
    392 #####COICEVAP OIceEvap 27 86400  4   flxat.nc   EXPORTED 
    393 #####tlmd    torc  LAG=1800  
    394 #####P  0 P  2 
    395 ###INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    396 #####INVERT CHECKIN MOZAIC CHECKOUT 
    397 # Reverse 
    398 #####NORSUD WSTEST 
    399 # CHECKIN: indicate computation of global, land and sea field integrals. 
    400 ##### INT=1 
    401 # Interpolation method ou parametres mozaic 
    402 #####wa2o.flx 92  2    14 
    403 # change of sign 
    404 ###-1.0 0 
    405 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    406 ##### INT=1 
     395# Field 16 : Non solar heat flux on ocean (a->o flx 4) 
     396# 
     397CONSFOCE O_QnsOce 6 86400  4   flxat.nc  <output_mode> 
     398tlmd    torc   LAG=1800 
     399P  0 P  2 
     400INVERT CHECKIN MOZAIC CHECKOUT 
     401# Reverse 
     402NORSUD WSTEST 
     403# CHECKIN: indicate computation of global, land and sea field integrals. 
     404 INT=1 
     405# Interpolation method ou parametres mozaic 
     406wa2o.flx 92  2    14  
     407# CHECKOUT: indicate computation of global, land and sea field integrals. 
     408 INT=1 
    407409# 
    408410########################################################################### 
    409 # Field 6 : solar heat flux on ocean (a->o flx 2) 
    410 # 
    411 COSWFLDO O_QsrOce  7 86400 4   flxat.nc  EXPORTED 
    412 tlmd    torc  LAG=1800 
    413 P  0 P  2 
    414 INVERT CHECKIN MOZAIC CHECKOUT 
    415 # Reverse 
    416 NORSUD WSTEST 
    417 # CHECKIN: indicate computation of global, land and sea field integrals. 
    418  INT=1 
    419 # Interpolation method or mozaic parameters 
    420 wa2o.flx 92  2    14 
    421 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    422  INT=1 
    423 ######################################################################### 
    424 # 
    425 # Field 8 : Non solar heat flux on ocean (a->o flx 4) 
    426 # 
    427 CONSFLDO O_QnsOce 6 86400  4   flxat.nc  EXPORTED 
    428 tlmd    torc   LAG=1800 
    429 P  0 P  2 
    430 INVERT CHECKIN MOZAIC CHECKOUT 
    431 # Reverse 
    432 NORSUD WSTEST 
    433 # CHECKIN: indicate computation of global, land and sea field integrals. 
    434  INT=1 
    435 # Interpolation method ou parametres mozaic 
    436 wa2o.flx 92  2    14  
    437 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    438  INT=1 
    439 # 
    440 ########################################################################### 
    441 # Field 5 : solar heat flux on ice (a->o flx 1) 
    442 # 
    443 COSHFLIC O_QsrIce  7 86400 4   flxat.nc  EXPORTED 
     411# 
     412# Field 17 : solar heat flux on ice (a->o flx 1) 
     413# 
     414COSHFICE O_QsrIce  7 86400 4   flxat.nc  <output_mode> 
    444415tlmd    torc  LAG=1800 
    445416P  0 P  2 
     
    454425 INT=1 
    455426# 
    456 # 
    457427############################################################################ 
    458428# 
    459 # Field 7 : Non solar heat flux on ice (a->o flx 3) 
    460 # 
    461 CONSFLIC O_QnsIce 6 86400  4  flxat.nc  EXPORTED 
     429# Field 18 : Non solar heat flux on ice (a->o flx 3) 
     430# 
     431CONSFICE O_QnsIce 6 86400  4  flxat.nc  <output_mode> 
    462432tlmd    torc  LAG=1800 
    463433P  0 P  2 
     
    474444############################################################################ 
    475445# 
    476 # Field 9 : Non solar heat flux derivative (o->a flx 5) 
    477 # 
    478 CODFLXDT O_dQnsdT  35 86400  4   flxat.nc  EXPORTED 
     446# Field 19 : Non solar heat flux derivative (o->a flx 5) 
     447# 
     448CODFLXDT O_dQnsdT  35 86400  4   flxat.nc  <output_mode> 
    479449tlmd    torc  LAG=1800 
    480450P  0 P  2 
     
    491461######################################################################### 
    492462# 
    493 # Field 16 : iceberg calving (a->o flx 12) 
    494 # 
    495 CRWOCEIS OCalving  32 86400  5  flxat.nc  EXPORTED 
     463# Field 20 : iceberg calving (a->o flx 12) 
     464# 
     465COCALVIN OCalving  36 86400  4  flxat.nc  <output_mode> 
    496466tlmd    torc  LAG=1800  
    497467P  0 P  2 
     468INVERT CHECKIN MOZAIC CHECKOUT 
     469# Reverse 
     470NORSUD WSTEST 
     471# CHECKIN: indicate computation of global, land and sea field integrals. 
     472 INT=1 
     473# Interpolation method ou parametres mozaic 
     474wa2o.cal  98   5   36 
     475# CHECKOUT: indicate computation of global, land and sea field integrals. 
     476 INT=1 
     477# 
     478######################################################################### 
     479# 
     480# Field 21 : liquid run-off (river + direct)  (a->o flx 11) 
     481# 
     482COLIQRUN  O_Runoff 32 86400  5  flxat.nc   <output_mode> 
     483tlmd    torc  LAG=1800  
     484P  0 P  2 
    498485INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    499486# Reverse 
     
    502489 INT=1 
    503490# Interpolation method ou parametres mozaic 
    504 wa2o.cal  98   5   36 
    505 # change of sign 
    506 -1.0 0 
    507 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    508  INT=1 
    509 ######################################################################### 
    510 # 
    511 # Field 15 : river run-off (a->o flx 11) 
    512 # 
    513 CRWOCERD  O_Runoff 32 86400  5  flxat.nc   EXPORTED 
    514 tlmd    torc  LAG=1800  
    515 P  0 P  2 
    516 INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    517 # Reverse 
    518 NORSUD WSTEST 
    519 # CHECKIN: indicate computation of global, land and sea field integrals. 
    520  INT=1 
    521 # Interpolation method ou parametres mozaic 
    522491wa2o.run 94    4   158 
    523 # change units from m/s to kg/s, change of sign 
     492# change units from m/s to kg/s 
    524493# weights convert from kg/s to kg/m^2/s 
    525 -1000.0 0 
    526 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    527  INT=1 
    528 # 
    529 ######################################################################### 
    530 # 
    531 # Field 14 : direct run-off (a->o flx 10) 
    532 # 
    533 ##CRWOCECD  FRWOCECD 32 86400  5  flxat.nc   EXPORTED 
    534 ##tlmd    torc  LAG=1800  
    535 ##P  0 P  2 
    536 ##INVERT CHECKIN MOZAIC BLASNEW CHECKOUT 
    537 # Reverse 
    538 ##NORSUD WSTEST 
    539 # CHECKIN: indicate computation of global, land and sea field integrals. 
    540 ## INT=1 
    541 # Interpolation method ou parametres mozaic 
    542 ##wa2o.run 94    4   158 
    543 # change units from m/s to kg/s, change of sign 
    544 # weights convert from kg/s to kg/m^2/s 
    545 ##-1000.0 0 
    546 # CHECKOUT: indicate computation of global, land and sea field integrals. 
    547 ## INT=1 
    548 # 
    549 ########################################################################## 
     4941000.0 0 
     495# CHECKOUT: indicate computation of global, land and sea field integrals. 
     496 INT=1 
     497# 
    550498########################################################################### 
    551499# 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namelist_ORCA2

    r396 r561  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
    33!! namelists    2 - Domain           (nam_zgr, nam_zgr_sco, namdom) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
    55!!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 
    6 !!              4 - lateral boundary (namlbc, namcla, namobc) 
     6!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, namtide) 
    77!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    88!!              6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp) 
     
    1010!!              8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm) 
    1111!!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
    12 !!              9 - miscellaneous    (namsol) 
     12!!              9 - miscellaneous    (namsol, nam_mpp, nam_mpp_dyndist, namctl) 
    1313!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    14 !  CAUTION: some scripts does not support CAPITALs for logical definition  
    15 !  *******       use .true. or .false. and NOT .TRUE. or .FALSE. 
     14!  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. 
    1615 
    1716!!====================================================================== 
     
    1918!!====================================================================== 
    2019!!   namrun            parameters of the run 
    21 !!   nam_mpp           Massively Parallel Processing 
    22 !!   nam_mpp_dyndist    
    23 !!   nam_ctl           Control prints & Benchmark 
    24 !!====================================================================== 
    25  
    26 !----------------------------------------------------------------------- 
    27 !       namrun   parameters of the run 
    28 !----------------------------------------------------------------------- 
    29 !  no         job number 
    30 !  cexper     experience name for vairmer format 
    31 !  ln_rstart  boolean term for restart (true or false) 
    32 !  nrstdt     restart control = 0 restart, do not control nit000 in the restart file. 
    33 !             !               = 1 restart, control nit000 in the restart file. Do not 
    34 !             !                   use the date in the restart file (use ndate0 in namelist) 
    35 !             !               = 2 restart, control nit000 in the restart file, use the date 
    36 !             !                   in the restart file. ndate0 in the namelist is ignored. 
    37 !  nit000     number of the first time step 
    38 !  nitend     number of the last time step 
    39 !  ndate0     initial calendar date aammjj 
    40 !  nleapy     Leap year calendar (0/1) 
    41 !  ninist     initial state output flag (0/1) 
    42 !  nstock     frequency of restart file 
    43 !  nwrite     frequency of OUTPUT file 
    44 !  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 
    45 &namrun 
    46    no         =       0 
    47    cexper     =  "ORCA2" 
    48    ln_rstart  =  .false. 
    49    nrstdt     =       1 
    50    nit000     =      76 
    51    nitend     =     150 
    52    ndate0     =18600106 
    53    nleapy     =      30 
    54    ninist     =       0 
    55    nstock     =    5400  
    56    nwrite     =       1 
    57    ln_dimgnnn = .false. 
    58 / 
    59 !----------------------------------------------------------------------- 
    60 !       nam_mpp      Massively Parallel Processing 
    61 !----------------------------------------------------------------------- 
    62 !  c_mpi_send         mpi send/recieve type 
    63 !                      = 'S'  : standard blocking send 
    64 !                      = 'B'  : buffer blocking send 
    65 !                      = 'I'  : immediate non-blocking send 
    66 &nam_mpp 
    67    c_mpi_send =  'S' 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 !       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist") 
    71 !----------------------------------------------------------------------- 
    72 !  jpni   number of processors following i 
    73 !  jpnj   number of processors following j 
    74 !  jpnij  number of local domains 
    75 &nam_mpp_dyndist 
    76    jpni  = 1 
    77    jpnj  = 1 
    78    jpnij = 1 
    79 / 
    80 !----------------------------------------------------------------------- 
    81 !       nam_ctl      Control prints & Benchmark 
    82 !----------------------------------------------------------------------- 
    83 !  ln_ctl     trends control print (expensive!) 
    84 !  nprint     level of print (0 no print) 
    85 !  nictls     start i indice to make the control SUM (very usefull to compare mono- 
    86 !  nictle     end   i indice to make the control SUM (-versus multi processor runs) 
    87 !  njctls     start j indice to make the control SUM (very usefull to compare mono- 
    88 !  njctle     end   j indice to make the control SUM (-versus multi processor runs) 
    89 !  nisplt     number of processors following i 
    90 !  njsplt     number of processors following j 
    91 !  nbench     Bench parameter (0/1): CAUTION it must be zero except for bench 
    92 !             for which we don't care about physical meaning of the results 
    93 !  nbit_cmp   bit comparison mode parameter (0/1): enables bit comparison between  
    94 !             single and multiple processor runs. 
    95 &namctl 
    96    ln_ctl = .false. 
    97    nprint =       0 
    98    nictls =       0 
    99    nictle =       0 
    100    njctls =       0 
    101    njctle =       0 
    102    isplt  =       1 
    103    jsplt  =       1 
    104    nbench =       0 
    105    nbit_cmp =     0 
    106 / 
    107  
     20!!====================================================================== 
     21 
     22!----------------------------------------------------------------------- 
     23&namrun        !   parameters of the run 
     24!----------------------------------------------------------------------- 
     25   no          =       0       !  job number 
     26   cexper      =  "ORCA2"      !  experience name  
     27   cn_ocerst_in  = "restartopa"!  suffix of ocean restart name (input) 
     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) 
     41/ 
    10842!!====================================================================== 
    10943!!                      ***  Domain namelists  *** 
     
    11549 
    11650!----------------------------------------------------------------------- 
    117 !       nam_zgr       vertical coordinate 
    118 !----------------------------------------------------------------------- 
    119 !  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined) 
    120 !  ln_zps     z-coordinate - partial steps   (T/F) 
    121 !  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
    122 &nam_zgr 
    123    ln_zco   =  .false. 
    124    ln_zps   =  .true. 
    125    ln_sco   =  .false. 
    126 / 
    127 !----------------------------------------------------------------------- 
    128 !       nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
    129 !----------------------------------------------------------------------- 
    130 !  sbot_min   minimum depth of s-bottom surface (>0) (m) 
    131 !  sbot_max   maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    132 !  theta      surface control parameter (0<=theta<=20) 
    133 !  thetb      bottom control parameter  (0<=thetb<= 1) 
    134 !  r_max      maximum cut-off r-value allowed (0<r_max<1) 
    135 &nam_zgr_sco 
    136    sbot_min =  300. 
    137    sbot_max = 5250. 
    138    theta    =    6.0 
    139    thetb    =    0.75 
    140    r_max    =    0.15 
    141 / 
    142 !----------------------------------------------------------------------- 
    143 !       namdom   space and time domain (bathymetry, mesh, timestep) 
    144 !----------------------------------------------------------------------- 
    145 !  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 
    146 !  e3zps_min  the thickness of the partial step is set larger than the 
    147 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 
    148 !  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    149 !  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra 
    150 !                                                  = 1,     used, rdt < rdttra(k) 
    151 !  atfp       asselin time filter parameter 
    152 !  rdt        time step for the dynamics (and tracer if nacc=0) 
    153 !  rdtmin     minimum time step on tracers 
    154 !  rdtmax     maximum time step on tracers 
    155 !  rdth       depth variation of tracer time step 
    156 !  rdtbt      barotropic time step (for the time splitting algorithm) 
    157 !  nclosea    =0 remove the closed sea from the global domain (orca configuration) 
    158 !             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
    159 &namdom 
    160    ntopo     =     1 
    161    e3zps_min =     5. 
    162    e3zps_rat =     0.1 
    163    nmsh      =     0 
    164    nacc      =     0 
    165    atfp      =     0.1 
    166    rdt       =  5760. 
    167    rdtmin    =  5760. 
    168    rdtmax    =  5760. 
    169    rdth      =   800. 
    170    rdtbt     =    90. 
    171    nclosea   =     0 
    172 / 
    173  
     51&nam_zgr       !   vertical coordinate 
     52!----------------------------------------------------------------------- 
     53   ln_zco      = .true.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
     54   ln_zps      = .false.    !  z-coordinate - partial steps   (T/F) 
     55   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     56/ 
     57!----------------------------------------------------------------------- 
     58&nam_zgr_sco   !   s-coordinate or hybrid z-s-coordinate 
     59!----------------------------------------------------------------------- 
     60   sbot_min    =  300.     !  minimum depth of s-bottom surface (>0) (m) 
     61   sbot_max    = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     62   theta       =    6.0    !  surface control parameter (0<=theta<=20) 
     63   thetb       =    0.75   !  bottom control parameter  (0<=thetb<= 1) 
     64   r_max       =    0.15   !  maximum cut-off r-value allowed (0<r_max<1) 
     65/ 
     66!----------------------------------------------------------------------- 
     67&namdom        !   space and time domain (bathymetry, mesh, timestep) 
     68!----------------------------------------------------------------------- 
     69   ntopo       =    1      !  compute (=0) or read(=1) the bathymetry file 
     70   e3zps_min   =    5.     !  the thickness of the partial step is set larger than the minimum 
     71   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   nacc        =    0      !  =1 acceleration of convergence method used, rdt < rdttra(k) 
     74                           !  =0, no acceleration, rdt = rdttra 
     75   atfp        =    0.1    !  asselin time filter parameter 
     76   rdt         = 5760.     !  time step for the dynamics (and tracer if nacc=0) 
     77   rdtmin      = 5760.     !  minimum time step on tracers (used if nacc=1) 
     78   rdtmax      = 5760.     !  maximum time step on tracers (used if nacc=1) 
     79   rdth        =  800.     !  depth variation of tracer time step  (used if nacc=1) 
     80   nn_baro     =   64      !  number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 
     81   nclosea     =    1      !  = 0 no closed sea in the model domain 
     82                           !  = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
     83/ 
    17484!!====================================================================== 
    17585!!            ***  Surface Boundary Condition namelists  *** 
     
    18898 
    18999!----------------------------------------------------------------------- 
    190 !       namsbc       Surface Boundary Condition (surface module) 
    191 !----------------------------------------------------------------------- 
    192 !  nn_fsbc       frequency of surface boundary condition computation  
    193 !                (= the frequency of sea-ice model call) 
    194 !  ln_ana        analytical formulation (fill namsbc_ana )  
    195 !  ln_flx        flux formulation       (fill namsbc_flx )  
    196 !  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)  
    197 !  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)  
    198 !  ln_cpl        Coupled formulation    (fill namsbc_cpl )  
    199 !  nn_ice        =0 no ice boundary condition   , 
    200 !                =1 observed ice-cover          , 
    201 !                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
    202 !  nn_ico_cpl    ice-ocean coupling : =0 LIM-3 old case  
    203 !                                     =1 stresses computed using now ocean velocity 
    204 !                                     =2 combination of 0 and 1 cases 
    205 !  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    206 !  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
    207 !  nn_fwb       FreshWater Budget: =0 no control                              ,  
    208 !                                  =1 annual global mean of e-p-r set to zero   , 
    209 !                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step 
    210 &namsbc 
    211    nn_fsbc     = 5 
    212    ln_ana      = .false. 
    213    ln_flx      = .false. 
    214    ln_blk_clio = .false. 
    215    ln_blk_core = .false. 
    216    ln_cpl      = .true. 
    217    nn_ice      = 2 
    218    nn_ico_cpl  = 0 
    219    ln_dm2dc    = .false. 
    220    ln_rnf      = .false. 
    221    ln_ssr      = .false. 
    222    nn_fwb      = 0 
    223 / 
    224 !----------------------------------------------------------------------- 
    225 !       namsbc_ana   analytical surface boundary condition 
    226 !----------------------------------------------------------------------- 
    227 !  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
    228 !  rn_utau0    uniform value for the i-stress 
    229 !  rn_vtau0    uniform value for the j-stress 
    230 !  rn_q0       uniform value for the total heat flux 
    231 !  rn_qsr0     uniform value for the solar radiation 
    232 !  rn_emp0     uniform value for the freswater budget (E-P) 
    233 &namsbc_ana 
    234    nn_tau000  =   0 
    235    rn_utau0   =   0.5 
    236    rn_vtau0   =   0.e0 
    237    rn_q0      =   0.e0 
    238    rn_qsr0    =   0.e0 
    239    rn_emp0    =   0.e0 
    240 / 
    241 !----------------------------------------------------------------------- 
    242 !       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
    243 !----------------------------------------------------------------------- 
    244 !  cn_dir   directory in which the model is executed 
     100&namsbc        !   Surface Boundary Condition (surface module) 
     101!----------------------------------------------------------------------- 
     102   nn_fsbc     = 5         !  frequency of surface boundary condition computation  
     103                           !               (= the frequency of sea-ice model call) 
     104   ln_ana      = .false.   !  analytical formulation (T => fill namsbc_ana )  
     105   ln_flx      = .false.   !  flux formulation       (T => fill namsbc_flx ) 
     106   ln_blk_clio = .false.    !  CLIO bulk formulation  (T => fill namsbc_clio)  
     107   ln_blk_core = .false.   !  CORE bulk formulation  (T => fill namsbc_core)  
     108   ln_cpl      = .true.   !  Coupled formulation    (T => fill namsbc_cpl ) 
     109   nn_ice      = 2         !  =0 no ice boundary condition   , 
     110                           !  =1 use observed ice-cover      , 
     111                           !  =2 ice-model used                             ("key_lim3" or "key_lim2) 
     112   nn_ico_cpl  = 0         !  ice-ocean coupling : =0 each nn_fsbc  
     113                           !                       =1 stresses recomputed each ocean time step ("key_lim3" only) 
     114                           !                       =2 combination of 0 and 1 cases             ("key_lim3" only) 
     115   ln_dm2dc    = .false.   !  daily mean to diurnal cycle short wave (qsr) 
     116   ln_rnf      = .true.    !  runoffs (T => fill namsbc_rnf) 
     117   ln_ssr      = .false.    !  Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 
     118   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked                              ,  
     119                           !                     =1 global mean of e-p-r set to zero at each nn_fsbc time step   , 
     120                           !                     =2 annual global mean of e-p-r set to zero 
     121/ 
     122!----------------------------------------------------------------------- 
     123&namsbc_ana    !   analytical surface boundary condition 
     124!----------------------------------------------------------------------- 
     125   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps 
     126   rn_utau0    =   0.5     !  uniform value for the i-stress 
     127   rn_vtau0    =   0.e0    !  uniform value for the j-stress 
     128   rn_q0       =   0.e0    !  uniform value for the total heat flux 
     129   rn_qsr0     =   0.e0    !  uniform value for the solar radiation 
     130   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P) 
     131/ 
     132!----------------------------------------------------------------------- 
     133&namsbc_flx    !   surface boundary condition : flux formulation 
     134!----------------------------------------------------------------------- 
     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' 
    245142! 
    246 !  THE ORDER OF THE FILES MATTER: 
    247 !               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    248 !               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    249 &namsbc_flx 
    250    cn_dir       = './' 
    251    sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0 
    252    sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0 
    253    sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0 
    254    sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0 
    255    sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0 
     143   cn_dir      = './'      !  root directory for the location of the flux files 
    256144/       
    257145!----------------------------------------------------------------------- 
    258 !       namsbc_clio  CLIO bulk formulea 
    259 !----------------------------------------------------------------------- 
    260 !  cn_dir   directory in which the model is executed 
     146&namsbc_clio   !   namsbc_clio  CLIO bulk formulea 
     147!----------------------------------------------------------------------- 
     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'        
    261157! 
    262 !  THE ORDER OF THE FILES MATTER: 
    263 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    264 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    265 &namsbc_clio 
    266    cn_dir       = './' 
    267    sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0     
    268    sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0     
    269    sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0       
    270    sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0      
    271    sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0     
    272    sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0       
    273    sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0        
    274 / 
    275 !----------------------------------------------------------------------- 
    276 !       namsbc_core  CORE bulk formulea 
    277 !----------------------------------------------------------------------- 
    278 !  cn_dir         directory in which the model is executed 
    279 !  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
    280 !  alpha_precip   multiplication factor for precipitation (total & snow) 
     158   cn_dir      = './'      !  root directory for the location of the bulk files are 
     159/ 
     160!----------------------------------------------------------------------- 
     161&namsbc_core   !   namsbc_core  CORE bulk formulea 
     162!----------------------------------------------------------------------- 
     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     =    't10'           ,       24.         , 'T_10_MOD' ,    .false.     , .true.  ,   'yearly'      
     170   sn_humi     =    'q10'           ,       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'        
    281173! 
    282 !  THE ORDER OF THE FILES MATTER: 
    283 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    284 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    285 &namsbc_core 
    286    cn_dir       = './' 
    287    ln_2m        = .false. 
    288    alpha_precip = 1. 
    289    sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0     
    290    sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0     
    291    sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0       
    292    sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0       
    293    sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0      
    294    sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0     
    295    sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0       
    296    sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0        
    297 / 
    298 !----------------------------------------------------------------------- 
    299 !       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled") 
    300 !----------------------------------------------------------------------- 
    301 &namsbc_cpl 
     174   cn_dir      = './'      !  root directory for the location of the bulk files 
     175   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     176   alpha_precip= 1.        !  multiplicative factor for precipitation (total & snow) 
     177/ 
     178!----------------------------------------------------------------------- 
     179&namsbc_cpl    !   coupled ocean/atmosphere model                        ("key_coupled") 
     180!----------------------------------------------------------------------- 
    302181! SEND 
    303182cn_snd_temperature= 'weighted oce and ice'  ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 
    304183cn_snd_albedo     = 'weighted ice'          ! 'none' 'weighted ice' 'mixed oce-ice' 
    305184cn_snd_thickness  = 'none'                  ! 'none' 'weighted ice and snow' 
    306 cn_snd_current_1 = 'none'                  ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 
    307 cn_snd_current_2 = 'spherical'             ! 'spherical' 'cartesian' 
    308 cn_snd_current_3 = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
    309 cn_snd_current_4  = 'T'                     ! 'T''U,V' 
     185cn_snd_crt_nature = 'none'                  ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 
     186cn_snd_crt_refere = 'spherical'             ! 'spherical' 'cartesian' 
     187cn_snd_crt_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
     188cn_snd_crt_grid   = 'T'                     ! 'T' 
    310189! RECEIVE 
    311 cn_rcv_w10m       = 'coupled'               ! 'none' 'coupled' 
    312 cn_rcv_stress_1   = 'oce only'              ! 'oce only' 'oce and ice' 'mixed oce-ice' 
    313 cn_rcv_stress_2  = 'cartesian'             ! 'spherical' 'cartesian' 
    314 cn_rcv_stress_3  = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
    315 cn_rcv_stress_4   = 'U,V'                   ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 
     190cn_rcv_w10m       = 'none'                  ! 'none' 'coupled' 
     191cn_rcv_tau_nature = 'mixed oce-ice'         ! 'oce only' 'oce and ice' 'mixed oce-ice' 
     192cn_rcv_tau_refere = 'cartesian'             ! 'spherical' 'cartesian' 
     193cn_rcv_tau_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
     194cn_rcv_tau_grid   = 'U,V'                   ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 
    316195cn_rcv_dqnsdt     = 'coupled'               ! 'none' 'coupled' 
    317196cn_rcv_qsr        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    318197cn_rcv_qns        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    319 cn_rcv_emp        = 'oce and ice'          ! 'conservative' 'oce and ice' 'mixed oce-ice' 
     198cn_rcv_emp        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    320199cn_rcv_rnf        = 'coupled'               ! 'coupled' 'climato' 'mixed' 
    321200cn_rcv_cal        = 'coupled'               ! 'none' 'coupled' 
    322201/ 
    323202!----------------------------------------------------------------------- 
    324 !       namqsr   penetrative solar radiation 
    325 !----------------------------------------------------------------------- 
    326 !  ln_traqsr : penetrative solar radiation (T) or not (F)     (Default=T) 
    327 !  rabs       fraction of qsr associated with xsi1 
    328 !  xsi1       first depth of extinction 
    329 !  xsi2       second depth of extinction 
    330 &namqsr 
    331    ln_traqsr = .true. 
    332    rabs     =   0.58 
    333    xsi1     =   0.35 
    334    xsi2     =   23.0 
    335 / 
    336 !----------------------------------------------------------------------- 
    337 !       namsbc_rnf   runoffs namelist surface boundary condition 
    338 !----------------------------------------------------------------------- 
    339 !cn_dir      Root directory for location of ssr files 
    340 !nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
    341 !rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
    342 !rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
    343 !sn_rnf      informations about the runoff file to be read 
    344 !sn_cnf      informations about the runoff mouth file to be read 
    345 !  THE ORDER OF THE FILES MATTER: 
    346 !               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
    347 !               !                    !   (if <0  months) !   name     !     (logical)      ! 
    348 &namsbc_rnf 
    349    cn_dir     = './' 
    350    nn_runoff  = 2 
    351    rn_hrnf    = 0.e0 
    352    rn_avt_rnf = 1.e-3 
    353    sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0      
    354    sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0               
    355 / 
    356 !----------------------------------------------------------------------- 
    357 !       namsbc_ssr   surface boundary condition : sea surface restoring 
    358 !----------------------------------------------------------------------- 
    359 !  cn_dir    directory in which the model is executed 
    360 !  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
    361 !  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
    362 !             (=1, salt flux, concentration/dillution emps only) 
    363 !             (=2, volume flux, both emp and emps are updated)  
    364 !  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
    365 !  deds       magnitude of the damping on salinity [mm/day] 
    366 !  THE ORDER OF THE FILES MATTER: 
    367 !               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    368 !               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    369 &namsbc_ssr 
    370    cn_dir       = './' 
    371    nn_sstr      = 0 
    372    nn_sssr      = 0 
    373    dqdt         = -40. 
    374    deds         =  27.7 
    375    sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0 
    376    sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0 
     203&namqsr        !   penetrative solar radiation 
     204!----------------------------------------------------------------------- 
     205   ln_traqsr   = .true.    !  penetrative solar radiation (T) or not (F) 
     206   rabs        =   0.58    !  fraction of qsr associated with xsi1 
     207   xsi1        =   0.35    !  first depth of extinction 
     208   xsi2        =   23.0    !  second depth of extinction 
     209/ 
     210!----------------------------------------------------------------------- 
     211&namsbc_rnf    !   runoffs namelist surface boundary condition 
     212!----------------------------------------------------------------------- 
     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'               
     217! 
     218   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 mouths 
     221   rn_hrnf      =   0.e0    !  depth over which enhanced vertical mixing is used 
     222   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] 
     223/ 
     224!----------------------------------------------------------------------- 
     225&namsbc_ssr    !   surface boundary condition : sea surface restoring 
     226!----------------------------------------------------------------------- 
     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' 
     231!    
     232   cn_dir      = './'      !  root directory for the location of the runoff files 
     233   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0) 
     234   nn_sssr     =     1     !  add a damping     term in the surface freshwater flux (=1) or not (=0) 
     235   dqdt        =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K] 
     236   deds        =   -27.7   !  magnitude of the damping on salinity   [mm/day/psu] 
    377237/       
    378238!----------------------------------------------------------------------- 
    379 !       namalb   albedo parameters 
    380 !----------------------------------------------------------------------- 
    381 !  cgren    correction of the snow or ice albedo to take into account 
    382 !  albice   albedo of melting ice in the arctic and antarctic 
    383 !  alphd    coefficients for linear interpolation used to compute albedo 
    384 !           between two extremes values (Pyane, 1972) 
    385 !  alphc     "                                         " 
    386 !  alphdi    "                                         " 
    387 &namalb 
    388    cgren    =      0.06 
    389    albice   =      0.5 
    390    alphd    =      0.80 
    391    alphc    =      0.65 
    392    alphdi   =      0.72 
    393 / 
    394  
     239&namalb        !   albedo parameters 
     240!----------------------------------------------------------------------- 
     241   cgren       =    0.0    !  correction of the snow or ice albedo to take into account the 
     242   albice      =    0.5    !  albedo of melting ice in the arctic and antarctic 
     243   alphd       =    0.80   !  coefficients for linear interpolation used to 
     244   alphc       =    0.65   !  compute albedo between two extremes values  
     245   alphdi      =    0.72   !  (Pyane, 1972) 
     246/ 
    395247!!====================================================================== 
    396248!!               ***  Lateral boundary condition  *** 
     
    399251!!   namcla        cross land advection 
    400252!!   namobc        open boundaries parameters                           ("key_obc") 
    401 !!====================================================================== 
    402  
    403 !----------------------------------------------------------------------- 
    404 !       namlbc   lateral momentum boundary condition 
    405 !----------------------------------------------------------------------- 
    406 !  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip 
    407 !                   shlat = 2 , no slip     ;   2 < shlat     , strong slip 
    408 &namlbc 
    409    shlat  =      2. 
    410 / 
    411 !----------------------------------------------------------------------- 
    412 !       namcla   cross land advection 
    413 !----------------------------------------------------------------------- 
    414 !  n_cla   advection between 2 ocean pts separates by land  
    415 &namcla 
    416    n_cla   = 1 
    417 / 
    418 !----------------------------------------------------------------------- 
    419 !       namobc    open boundaries parameters (#ifdef key_obc) 
    420 !----------------------------------------------------------------------- 
    421 !  nobc_dta   = 0 the obc data are equal to the initial state 
    422 !             = 1 the obc data are read in 'obc   .dta' files 
    423 !  rdpeob  time relaxation (days) for the east open boundary 
    424 !  rdpwob  time relaxation (days) for the west open boundary 
    425 !  rdpnob  time relaxation (days) for the north open boundary 
    426 !  rdpsob  time relaxation (days) for the south open boundary 
    427 !  zbsic1  barotropic stream function on isolated coastline 1 
    428 !  zbsic2  barotropic stream function on isolated coastline 2 
    429 !  zbsic3  barotropic stream function on isolated coastline 3 
    430 !  ln_obc_clim  climatological obc data files (default T) 
    431 !  ln_vol_cst   total volume conserved 
    432 &namobc 
    433     nobc_dta =    0 
    434     rdpein   =    1. 
    435     rdpwin   =    1. 
    436     rdpnin   =   30. 
    437     rdpsin   =    1. 
    438     rdpeob   = 1500. 
    439     rdpwob   =   15. 
    440     rdpnob   =  150. 
    441     rdpsob   =   15. 
    442     zbsic1   =  140.e+6 
    443     zbsic2   =    1.e+6 
    444     zbsic3   =    0. 
    445     ln_obc_clim = .true. 
    446     ln_vol_cst  = .false. 
    447 / 
    448  
     253!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
     254!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
     255!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides") 
     256!!====================================================================== 
     257 
     258!----------------------------------------------------------------------- 
     259&namlbc        !   lateral momentum boundary condition 
     260!----------------------------------------------------------------------- 
     261   shlat       =    2.     !      shlat = 0 : free slip 
     262                           !  0 < shlat < 2 : partial slip 
     263                           !      shlat = 2 : no slip 
     264                           !  2 < shlat     : strong slip 
     265/ 
     266!----------------------------------------------------------------------- 
     267&namcla        !   cross land advection 
     268!----------------------------------------------------------------------- 
     269   n_cla       =    0      !  advection between 2 ocean pts separates by land 
     270/ 
     271!----------------------------------------------------------------------- 
     272&namobc        !   open boundaries parameters                           ("key_obc") 
     273!----------------------------------------------------------------------- 
     274    nobc_dta   =    1      !  = 0 the obc data are equal to the initial state 
     275                           !  = 1 the obc data are read in 'obc.dta' files 
     276    cffile     = 'annual'  !  set to annual if obc datafile hold 1 year of data 
     277                           !  set to monthly if obc datafile hold 1 month of data 
     278    rdpein     =    1.     !  ??? 
     279    rdpwin     =    1.     !  ??? 
     280    rdpnin     =    1.     !  ??? 
     281    rdpsin     =    1.     !  ??? 
     282    rdpeob     = 3000.     !  time relaxation (days) for the east  open boundary 
     283    rdpwob     =   15.     !    "        "             "     west         " 
     284    rdpnob     = 3000.     !    "        "             "     north        " 
     285    rdpsob     =   15.     !    "        "             "     south        " 
     286    zbsic1     =  140.e+6  !   barotropic stream function on first  isolated coastline 
     287    zbsic2     =    1.e+6  !    "                   "        second       " 
     288    zbsic3     =    0.     !    "                   "        thrid        " 
     289    ln_obc_clim= .false.   !  climatological obc data files (T) or not (F) 
     290    ln_vol_cst = .true.    !  impose the total volume conservation (T) or not (F) 
     291/ 
     292!----------------------------------------------------------------------- 
     293&namagrif      !                                                        ("key_agrif") 
     294!----------------------------------------------------------------------- 
     295    nbclineupdate = 3      !  baroclinic update frequency 
     296    ln_spc_dyn    = .true. !  use 0 as special value for dynamics 
     297    visc_tra      = 2880.  !  viscosity coeeficient for tracers sponge layer 
     298    visc_dyn      = 2880.  !  viscosity coeeficient for dynamics sponge layer 
     299/ 
     300!----------------------------------------------------------------------- 
     301&nambdy        !  unstructured open boundaries parameters               ("key_bdy") 
     302!----------------------------------------------------------------------- 
     303    filbdy_mask    =  ''                  !  name of mask file (if ln_bdy_mask=.TRUE.) 
     304    filbdy_data_T  = 'bdydata_grid_T.nc'  !  name of data file (T-points) 
     305    filbdy_data_U  = 'bdydata_grid_U.nc'  !  name of data file (U-points) 
     306    filbdy_data_V  = 'bdydata_grid_V.nc'  !  name of data file (V-points) 
     307    ln_bdy_clim    = .false.              !  contain 1 (T) or 12 (F) time dumps and be cyclic 
     308    ln_bdy_vol     = .true.               !  total volume correction (see volbdy parameter) 
     309    ln_bdy_mask    = .false.              !  boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) 
     310    ln_bdy_tides   = .true.               !  Apply tidal harmonic forcing with Flather condition 
     311    ln_bdy_dyn_fla = .true.               !  Apply Flather condition to velocities 
     312    ln_bdy_tra_frs = .false.              !  Apply FRS condition to temperature and salinity  
     313    ln_bdy_dyn_frs = .false.              !  Apply FRS condition to velocities 
     314    nbdy_dta       =  1                   !  = 0, bdy data are equal to the initial state 
     315                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     316    nb_rimwidth    = 9                    !  width of the relaxation zone 
     317    volbdy         = 0                    !  = 0, the total water flux across open boundaries is zero 
     318                                          !  = 1, the total volume of the system is conserved 
     319/ 
     320!----------------------------------------------------------------------- 
     321&namtide        ! tidal forcing at unstructured boundaries               
     322!----------------------------------------------------------------------- 
     323    filtide      = 'bdytide_'           !  file name root of tidal forcing files 
     324    tide_cpt     = 'M2','S1'            !  names of tidal components used 
     325    tide_speed   = 28.984106, 15.000001 !  phase speeds of tidal components (deg/hour) 
     326    ln_tide_date = .false.              !  adjust tidal harmonics for start date of run 
     327/ 
    449328!!====================================================================== 
    450329!!                 ***  Bottom boundary condition  *** 
     
    454333!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
    455334!!====================================================================== 
    456 !----------------------------------------------------------------------- 
    457 !       nambfr   bottom friction 
    458 !----------------------------------------------------------------------- 
    459 !  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction 
    460 !                                     = 3 : free slip  ;  = 2 : non linear friction 
    461 !  bfri1    bottom drag coefficient (linear case) 
    462 !  bfri2    bottom drag coefficient (non linear case) 
    463 !  bfeb2    bottom turbulent kinetic energy  (m2/s2) 
    464 &nambfr 
    465    nbotfr =       1 
    466    bfri1  =   4.e-4 
    467    bfri2  =   1.e-3 
    468    bfeb2  =  2.5e-3 
    469 / 
    470 !----------------------------------------------------------------------- 
    471 !       nambbc   bottom temperature boundary condition 
    472 !----------------------------------------------------------------------- 
    473 !  ngeo_flux  = 0 no geothermal heat flux 
    474 !             = 1 constant geothermal heat flux 
    475 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
    476 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    477 &nambbc 
    478    ngeo_flux =  2 
    479    ngeo_flux_const = 86.4e-3 
    480 / 
    481 !----------------------------------------------------------------------- 
    482 !       nambbl   bottom boundary layer scheme 
    483 !----------------------------------------------------------------------- 
    484 !                                !  diffusive bbl                       ("key_trabbl") 
    485 !                                !  advective bbl                       ("key_trabbl_adv") 
    486 !  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s) 
    487 &nambbl 
    488    atrbbl = 10000. 
    489 / 
    490  
     335 
     336!----------------------------------------------------------------------- 
     337&nambfr        !   bottom friction 
     338!----------------------------------------------------------------------- 
     339   nbotfr      =    1      !  type of bottom friction :   = 0 : no   slip,  = 2 : nonlinear friction 
     340                           !                              = 3 : free slip,  = 1 :    linear friction 
     341   bfri1       =    4.e-4  !  bottom drag coefficient (linear case) 
     342   bfri2       =    1.e-3  !  bottom drag coefficient (non linear case) 
     343   bfeb2       =    2.5e-3 !  bottom turbulent kinetic energy background  (m^2/s^2) 
     344/ 
     345!----------------------------------------------------------------------- 
     346&nambbc        !   bottom temperature boundary condition 
     347!----------------------------------------------------------------------- 
     348   ngeo_flux   =    2      !  geothermal heat flux = 0 no flux considered  
     349                           !                       = 1 constant flux 
     350                           !                       = 2 variable flux (read in geothermal_heating.nc in mW/m2)  
     351   ngeo_flux_const = 86.4e-3  !  Constant value of geothermal heat flux [W/m2] 
     352/ 
     353!----------------------------------------------------------------------- 
     354&nambbl        !   bottom boundary layer scheme 
     355!----------------------------------------------------------------------- 
     356!                          !  diffusive bbl                             ("key_trabbl") 
     357!                          !  advective bbl                             ("key_trabbl_adv") 
     358   atrbbl      =  10000.   !  lateral mixing coefficient in the bbl  [m2/s] 
     359/ 
    491360!!====================================================================== 
    492361!!                        Tracer (T & S ) namelists 
     
    499368 
    500369!----------------------------------------------------------------------- 
    501 !       nameos   ocean physical parameters 
    502 !----------------------------------------------------------------------- 
    503 !  neos    type of equation of state and Brunt-Vaisala frequency 
    504 !          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
    505 !          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    506 !          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    507 !  ralpha  thermal expension coefficient (neos= 1 or 2) 
    508 !  rbeta   saline  expension coefficient (neos= 2) 
    509 &nameos 
    510    neos   =      0 
    511    ralpha =  2.e-4 
    512    rbeta  =  0.001 
    513 / 
    514 !----------------------------------------------------------------------- 
    515 !       nam_traadv   advection scheme for tracer 
    516 !----------------------------------------------------------------------- 
    517 !  ln_traadv_cen2     2nd order centered scheme          (default T) 
    518 !  ln_traadv_tvd      TVD scheme                         (default F) 
    519 !  ln_traadv_muscl    MUSCL scheme                       (default F) 
    520 !  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F) 
    521 !  ln_traadv_ubs      UBS scheme                         (default F) 
    522 !  ln_traadv_qck      QUICKEST scheme                    (default F) 
    523 &nam_traadv 
    524    ln_traadv_cen2   =  .true. 
    525    ln_traadv_tvd    =  .false. 
    526    ln_traadv_muscl  =  .false. 
    527    ln_traadv_muscl2 =  .false. 
    528    ln_traadv_ubs    =  .false. 
    529    ln_traadv_qck    =  .false. 
    530 / 
    531 !----------------------------------------------------------------------- 
    532 !       nam_traldf   lateral diffusion scheme for tracer 
    533 !----------------------------------------------------------------------- 
    534 !  Type of the operator :  
    535 !     ln_traldf_lap    laplacian operator          (default T) 
    536 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    537 !  Direction of action  :  
    538 !     ln_traldf_level  iso-level                   (default F) 
    539 !     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T) 
    540 !     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
    541 !  Coefficient 
    542 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    543 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    544 !     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv") 
    545 &nam_traldf 
    546    ln_traldf_lap    =  .true. 
    547    ln_traldf_bilap  =  .false. 
    548    ln_traldf_level  =  .false. 
    549    ln_traldf_hor    =  .false. 
    550    ln_traldf_iso    =  .true. 
    551    aht0    =  2000. 
    552    ahtb0   =     0. 
    553    aeiv0   =  2000. 
    554 / 
    555 !----------------------------------------------------------------------- 
    556 !       namtdp   tracer newtonian damping ('key_tradmp') 
    557 !----------------------------------------------------------------------- 
    558 !  ndmp    type of damping in temperature and salinity  
    559 !          !   ='latitude', damping poleward of 'ndmp' degrees and function  
    560 !          !                of the distance-to-coast. Red and Med Seas as ndmp=-1 
    561 !          !   =-1          damping only in Med and Red Seas 
    562 !  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    563 !  nmldmp  type of damping in the mixed layer 
    564 !          !   =0 damping throughout the water column 
    565 !               !   =1 no damping in the mixed layer defined by avt >5cm2/s  
    566 !               !   =2 no damping in the mixed layer defined rho<rho(surf)+.01  
    567 !  sdmp    surface time scale for internal damping (days) 
    568 !  bdmp    bottom time scale for internal damping (days) 
    569 !  hdmp    depth of transition between sdmp and bdmp (meters) 
    570 &namtdp 
    571    ndmp   =   -1 
    572    ndmpf  =    1 
    573    nmldmp =    1 
    574    sdmp   =  50. 
    575    bdmp   = 360. 
    576    hdmp   = 800. 
    577 / 
    578  
     370&nameos        !   ocean physical parameters 
     371!----------------------------------------------------------------------- 
     372   neos        =    0      !  type of equation of state and Brunt-Vaisala frequency 
     373                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     374                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     375                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     376   ralpha      =    2.e-4  !  thermal expension coefficient (neos= 1 or 2) 
     377   rbeta       =    0.001  !  saline  expension coefficient (neos= 2) 
     378/ 
     379!----------------------------------------------------------------------- 
     380&nam_traadv    !   advection scheme for tracer  
     381!----------------------------------------------------------------------- 
     382   ln_traadv_cen2   =  .true.   !  2nd order centered scheme    
     383   ln_traadv_tvd    =  .false.  !  TVD scheme                 
     384   ln_traadv_muscl  =  .false.  !  MUSCL scheme              
     385   ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
     386   ln_traadv_ubs    =  .false.  !  UBS scheme                  
     387/ 
     388!----------------------------------------------------------------------- 
     389&nam_traldf    !   lateral diffusion scheme for tracer  
     390!----------------------------------------------------------------------- 
     391!                               !  Type of the operator :  
     392   ln_traldf_lap    =  .true.   !     laplacian operator        
     393   ln_traldf_bilap  =  .false.  !     bilaplacian operator      
     394                                !  Direction of action  : 
     395   ln_traldf_level  =  .false.  !     iso-level                 
     396   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
     397   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     398!                               !  Coefficient 
     399   aht0        =  2000.         !     horizontal eddy diffusivity for tracers [m2/s] 
     400   ahtb0       =     0.         !     background eddy diffusivity for ldf_iso [m2/s] 
     401   aeiv0       =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     402/ 
     403!----------------------------------------------------------------------- 
     404&namtdp        !   tracer newtonian damping                             ('key_tradmp') 
     405!----------------------------------------------------------------------- 
     406   ndmp        =   -1      !  type of damping in temperature and salinity  
     407                           !     ='latitude', damping poleward of 'ndmp' degrees and function  
     408                           !                  of the distance-to-coast. Red and Med Seas as ndmp=-1 
     409                           !     =-1 damping only in Med and Red Seas 
     410   ndmpf       =    1      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     411   nmldmp      =    1      !  type of damping: =0 damping throughout the water column 
     412                           !                   =1 no damping in the mixed layer defined by avt >5cm2/s ) 
     413                           !                   =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
     414   sdmp        =   50.     !  surface time scale for internal damping (days) 
     415   bdmp        =  360.     !  bottom  time scale for internal damping (days) 
     416   hdmp        =  800.     !  depth of transition between sdmp and bdmp (meters) 
     417/ 
    579418!!====================================================================== 
    580419!!                      ***  Dynamics namelists  *** 
     
    589428 
    590429!----------------------------------------------------------------------- 
    591 !       nam_dynadv   formulation of the momentum advection 
    592 !----------------------------------------------------------------------- 
    593 !  ln_dynadv_vec      vector form (T) or flux form (F)         (default T) 
    594 !  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
    595 !  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
    596 &nam_dynadv    
    597    ln_dynadv_vec   = .true. 
    598    ln_dynadv_cen2  = .false. 
    599    ln_dynadv_ubs   = .false. 
     430&nam_dynadv    !   formulation of the momentum advection 
     431!----------------------------------------------------------------------- 
     432   ln_dynadv_vec = .false.  !  vector form (T) or flux form (F)   
     433   ln_dynadv_cen2= .true. !  flux form - 2nd order centered scheme 
     434   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme  
    600435 
    601436!----------------------------------------------------------------------- 
    602 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
    603 !----------------------------------------------------------------------- 
    604 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    605 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    606 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    607 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    608 &nam_dynvor 
    609    ln_dynvor_ene = .false. 
    610    ln_dynvor_ens = .true. 
    611    ln_dynvor_mix = .false. 
    612    ln_dynvor_een = .false. 
    613 / 
    614 !----------------------------------------------------------------------- 
    615 !       nam_dynhpg   Hydrostatic pressure gradient option 
    616 !----------------------------------------------------------------------- 
    617 !  type of pressure gradient scheme (choose one only!) 
    618 !     ln_hpg_zco    z-coordinate - full steps                   (default T) 
    619 !     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
    620 !     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
    621 !     ln_hpg_hel    s-coordinate (helsinki modification) 
    622 !     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
    623 !     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
    624 !     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
    625 !  parameters 
    626 !     gamm          weighting coefficient (wdj scheme) 
    627 &nam_dynhpg 
    628    ln_hpg_zco = .false. 
    629    ln_hpg_zps = .true. 
    630    ln_hpg_sco = .false. 
    631    ln_hpg_hel = .false. 
    632    ln_hpg_wdj = .false. 
    633    ln_hpg_djc = .false. 
    634    ln_hpg_rot = .false. 
    635    gamm       = 0.e0 
    636 / 
    637 !----------------------------------------------------------------------- 
    638 !       namflg   algorithm flags (algorithm not control by CPP keys) 
    639 !----------------------------------------------------------------------- 
    640 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    641 !                                                  centered     time scheme  (F) 
    642 !   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
    643 &namflg 
    644    ln_dynhpg_imp   =  .false. 
    645    nn_dynhpg_rst   =  0 
    646 / 
    647 !----------------------------------------------------------------------- 
    648 !       nam_dynspg   surface pressure gradient   (CPP key only) 
    649 !----------------------------------------------------------------------- 
    650 !                                !  explicit free surface               ("key_dynspg_exp") 
    651 !                                !  filtered free surface               ("key_dynspg_flt") 
    652 !                                !  split-explicit free surface         ("key_dynspg_ts") 
    653 !                                !  rigid-lid                           ("key_dynspg_rl") 
    654  
    655 !----------------------------------------------------------------------- 
    656 !       nam_dynldf   lateral diffusion on momentum 
    657 !----------------------------------------------------------------------- 
    658 !  Type of the operator :  
    659 !     ln_dynldf_lap    laplacian operator          (default T) 
    660 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    661 !  Direction of action  :  
    662 !     ln_dynldf_level  iso-level                   (default F) 
    663 !     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T) 
    664 !     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
    665 !  Coefficient 
    666 !     ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    667 !     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    668 &nam_dynldf 
    669    ln_dynldf_lap    =  .true. 
    670    ln_dynldf_bilap  =  .false. 
    671    ln_dynldf_level  =  .false. 
    672    ln_dynldf_hor    =  .true. 
    673    ln_dynldf_iso    =  .false. 
    674    ahm0    = 40000. 
    675    ahmb0   =     0. 
    676 / 
    677  
     437&nam_dynvor    !   option of physics/algorithm (not control by CPP keys) 
     438!----------------------------------------------------------------------- 
     439   ln_dynvor_ene = .false. !  enstrophy conserving scheme   
     440   ln_dynvor_ens = .true.  !  energy conserving scheme     
     441   ln_dynvor_mix = .false. !  mixed scheme                
     442   ln_dynvor_een = .false. !  energy & enstrophy scheme   
     443/ 
     444!----------------------------------------------------------------------- 
     445&nam_dynhpg    !   Hydrostatic pressure gradient option 
     446!----------------------------------------------------------------------- 
     447   ln_hpg_zco  = .true.   !  z-coordinate - full steps                    
     448   ln_hpg_zps  = .false.    !  z-coordinate - partial steps (interpolation) 
     449   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
     450   ln_hpg_hel  = .false.   !  s-coordinate (helsinki modification) 
     451   ln_hpg_wdj  = .false.   !  s-coordinate (weighted density jacobian) 
     452   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
     453   ln_hpg_rot  = .false.   !  s-coordinate (ROTated axes scheme) 
     454   gamm        = 0.e0      !  weighting coefficient (wdj scheme) 
     455/ 
     456!----------------------------------------------------------------------- 
     457&namflg        !   algorithm flags (algorithm not control by CPP keys) 
     458!----------------------------------------------------------------------- 
     459   ln_dynhpg_imp = .false. !  hydrostatic pressure gradient: semi-implicit time scheme  (T) 
     460                           !                                 centered      time scheme  (F) 
     461   nn_dynhpg_rst =  0      !  add dynhpg implicit variables in restart ot not (1/0) 
     462/ 
     463!----------------------------------------------------------------------- 
     464!nam_dynspg    !   surface pressure gradient   (CPP key only) 
     465!----------------------------------------------------------------------- 
     466!                          !  explicit free surface                     ("key_dynspg_exp") 
     467!                          !  filtered free surface                     ("key_dynspg_flt") 
     468!                          !  split-explicit free surface               ("key_dynspg_ts") 
     469!                          !  rigid-lid                                 ("key_dynspg_rl") 
     470 
     471!----------------------------------------------------------------------- 
     472&nam_dynldf    !   lateral diffusion on momentum 
     473!----------------------------------------------------------------------- 
     474!                               !  Type of the operator :  
     475   ln_dynldf_lap    =  .true.   !     laplacian operator          
     476   ln_dynldf_bilap  =  .false.  !     bilaplacian operator     
     477!                               !  Direction of action  :  
     478   ln_dynldf_level  =  .false.  !     iso-level                
     479   ln_dynldf_hor    =  .true.   !     horizontal (geopotential)        (require "key_ldfslp" in s-coord.) 
     480   ln_dynldf_iso    =  .false.  !     iso-neutral                      (require "key_ldfslp") 
     481                                !  Coefficient 
     482   ahm0    = 40000.             !     horizontal eddy viscosity   [m2/s] 
     483   ahmb0   =     0.             !     background eddy viscosity for ldf_iso [m2/s] 
     484/ 
    678485!!====================================================================== 
    679486!!             Tracers & Dynamics vertical physics namelists 
     
    688495 
    689496!----------------------------------------------------------------------- 
    690 !       namzdf   vertical physics 
    691 !----------------------------------------------------------------------- 
    692 !  ln_zdfevd  enhanced vertical diffusion         (default T) 
    693 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    694 !  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst") 
    695 !  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst") 
    696 !  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
    697 !  n_evdm     = 0  apply enhanced mixing on tracer only 
    698 !             = 1  apply enhanced mixing on both tracer and momentum 
    699 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    700 !                               (=F)  euler backward (F) 
    701 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    702 &namzdf 
    703    ln_zdfevd = .true. 
    704    ln_zdfnpc = .false. 
    705    avm0      = 1.2e-4 
    706    avt0      = 1.2e-5 
    707    avevd     =   100. 
    708    n_evdm    =     1 
    709    ln_zdfexp =  .false. 
    710    n_zdfexp  =     3 
    711 / 
    712 !----------------------------------------------------------------------- 
    713 !       namnpc   non penetrative convection 
    714 !----------------------------------------------------------------------- 
    715 !  nnpc1   non penetrative convective scheme frequency    
    716 !  nnpc2   non penetrative convective scheme print frequency 
    717 &namnpc 
    718    nnpc1  =      1 
    719    nnpc2  =    365 
    720 / 
    721 !----------------------------------------------------------------------- 
    722 !       namric   richardson number dependent vertical diffusion         ("key_zdfric" ) 
    723 !----------------------------------------------------------------------- 
    724 !  avmri   maximum value of the vertical viscosity 
    725 !  alp     coefficient of the parameterization 
    726 !  nric    coefficient of the parameterization 
    727 &namric 
    728    avmri = 100.e-4 
    729    alp   =      5. 
    730    nric  =       2 
    731 / 
    732 !----------------------------------------------------------------------- 
    733 !       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
    734 !----------------------------------------------------------------------- 
    735 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    736 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    737 !  ediss    coef. of the Kolmogoroff dissipation   
    738 !  ebb      coef. of the surface input of tke 
    739 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    740 !  emin     minimum value of tke (m^2/s^2) 
    741 !  emin0    surface minimum value of tke (m^2/s^2) 
    742 !  nitke    number of restart iterative loops 
    743 !  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
    744 !  nmxl     length used   = 0 bounded by the distance to surface and bottom 
    745 !           !             = 1 bounded by the local vertical scale factor 
    746 !           !             = 2 first vertical derivative of mixing length bounded by 1 
    747 !  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
    748 !           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
    749 !           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
    750 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    751 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    752 &namtke 
    753    ln_rstke = .false. 
    754    ediff =       0.1 
    755    ediss =       0.7 
    756    ebb   =      3.75 
    757    efave =        1. 
    758    emin  =     1.e-6 
    759    emin0 =     1.e-4 
    760    nitke =        50 
    761    nmxl  =         2 
    762    npdl  =         1 
    763    navb  =         0 
    764 / 
    765 !----------------------------------------------------------------------- 
    766 !       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" ) 
    767 !----------------------------------------------------------------------- 
    768 !  ln_kpprimix   shear instability mixing  (default T) 
    769 !  difmiw        constant internal wave viscosity (m2/s) 
    770 !  difsiw        constant internal wave diffusivity (m2/s) 
    771 !  Riinfty       local Richardson Number limit for shear instability 
    772 !  difri         maximum shear mixing at Rig = 0    (m2/s) 
    773 !  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
    774 !  difcon        maximum mixing in interior convection (m2/s) 
    775 !  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
    776 !  navb          = 0/1 flag for constant or profile background avt 
    777 &namkpp 
    778    ln_kpprimix  = .true. 
    779    difmiw       =  1.e-04 
    780    difsiw       =  0.1e-04 
    781    Riinfty      =  0.8 
    782    difri        =  0.0050 
    783    bvsqcon      = -0.01e-07 
    784    difcon       =  1. 
    785    navb         =  0 
    786    nave         =  1 
    787 / 
    788 !----------------------------------------------------------------------- 
    789 !       namddm   double diffusive mixing parameterization               ("key_zdfddm") 
    790 !----------------------------------------------------------------------- 
    791 !   avts    maximum avs for dd mixing  
    792 !   hsbfr   heat/salt buoyancy flux ratio 
    793 &namddm 
    794       avts  = 1.e-4 
    795       hsbfr = 1.6 
    796 / 
    797  
     497&namzdf        !   vertical physics 
     498!----------------------------------------------------------------------- 
     499   avm0        =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
     500   avt0        =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
     501   ln_zdfnpc   = .false.   !  convection: Non-Penetrative algorithm (T) or not (F) 
     502   ln_zdfevd   = .true.    !  convection: enhanced vertical diffusion (T) or not (F)     
     503   avevd       = 100.      !  vertical coefficient for enhanced diffusion scheme [m2/s] 
     504   n_evdm      =   1       !  enhanced mixing apply on tracer (=0) or on tracer and momentum (=1) 
     505   ln_zdfexp   =  .false.  !  split explicit (T) or implicit (F) time stepping 
     506   n_zdfexp    =   3       !  number of sub-timestep for ln_zdfexp=T 
     507/ 
     508!----------------------------------------------------------------------- 
     509&namnpc        !   non penetrative convection 
     510!----------------------------------------------------------------------- 
     511   nnpc1       =    1      !  non penetrative convective scheme computation frequency 
     512   nnpc2       =  365      !  non penetrative convective scheme print frequency 
     513/ 
     514!----------------------------------------------------------------------- 
     515&namric        !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
     516!----------------------------------------------------------------------- 
     517   avmri       = 100.e-4   !  maximum value of the vertical viscosity 
     518   alp         =   5.      !  coefficient of the parameterization 
     519   nric        =   2       !  coefficient of the parameterization 
     520/ 
     521!----------------------------------------------------------------------- 
     522&namtke        !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
     523!----------------------------------------------------------------------- 
     524   ln_rstke    = .false.   !  restart with tke from a run without tke (T) or not (F) 
     525   nn_itke     =  50       !  number of iterative loops if ln_rstke=T  
     526   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     527   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation 
     528   rn_ebb      =   3.75    !  coef. of the surface input of tke 
     529   rn_efave    =   1.      !  boost of the tke diffusion ( avtke=rn_efave*avm ) 
     530   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2] 
     531   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2] 
     532   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
     533                           !                 = 1 bounded by the local vertical scale factor 
     534                           !                 = 2 first vertical derivative of mixing length bounded by 1 
     535                           !                 = 3 same criteria as case 2 but applied in a different way 
     536   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
     537   nn_avb      =   0       !  profile for constant background used on avt & avm (=1) or not (=0) 
     538   nn_ave      =   1       !  horizontal averaged on avt (=1) or not (=0)  
     539   ln_mxl0     = .false.   !  mixing length scale surface value as function of wind stress (T) or not (F) 
     540   rn_lmin     =   0.4     !  interior buoyancy lenght scale minimum value 
     541   rn_lmin0    =   0.4     !  surface  buoyancy lenght scale minimum value 
     542   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
     543                           !        = 0 no penetration ( O(2 km) resolution) 
     544                           !        = 1 additional tke source  
     545                           !        = 2 additional tke source applied only at the base of the mixed layer 
     546   nn_htau     =   2       !  type of exponential decrease of tke penetration 
     547                           !        = 0  constant 10 m length scale 
     548                           !        = 1  ??? 
     549                           !        = 2  ??? 
     550   rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     551   ln_lc       = .false.   !  Langmuir cell effect 
     552   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
     553   nn_havtb    =   0       !  horizontal shape for avtb (=1) or not (=0) 
     554/ 
     555!------------------------------------------------------------------------ 
     556&namkpp        !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionnally: 
     557!                                                                         "key_kppcustom" or "key_kpplktb") 
     558!------------------------------------------------------------------------ 
     559   ln_kpprimix = .true.    !  shear instability mixing  
     560   difmiw      =  1.0e-04  !  constant internal wave viscosity [m2/s] 
     561   difsiw      =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
     562   Riinfty     =  0.8      !  local Richardson Number limit for shear instability 
     563   difri       =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
     564   bvsqcon     = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]  
     565   difcon      =  1.       !  maximum mixing in interior convection [m2/s] 
     566   navb        =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
     567   nave        =  1        !  constant (=0) or profile (=1) background on avt 
     568/ 
     569!----------------------------------------------------------------------- 
     570&namddm        !   double diffusive mixing parameterization                  ("key_zdfddm") 
     571!----------------------------------------------------------------------- 
     572      avts     = 1.e-4     !  maximum avs (vertical mixing on salinity) 
     573      hsbfr    = 1.6       !  heat/salt buoyancy flux ratio 
     574/ 
    798575!!====================================================================== 
    799576!!                  ***  Miscelaneous namelists  *** 
    800577!!====================================================================== 
    801 !!   namsol      elliptic solver / island / free surface  
    802 !!====================================================================== 
    803  
    804 !----------------------------------------------------------------------- 
    805 !       namsol   elliptic solver / island / free surface  
    806 !----------------------------------------------------------------------- 
    807 !  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
    808 !                            (=2 successive-over-relaxation: sor) 
    809 !                            (=3 FETI currently it does not work!       ("key_feti") 
    810 !  nsol_arp  absolute/relative (0/1) precision convergence test 
    811 !  nmin      minimum of iterations for the SOR solver 
    812 !  nmax      maximum of iterations for the SOR solver 
    813 !  nmod      frequency of test for the SOR solver 
    814 !  eps       absolute precision of the solver 
    815 !  resmax    absolute precision for the SOR solver 
    816 !  sor       optimal coefficient for SOR solver 
    817 !  epsisl    absolute precision on stream function solver 
    818 !  nmisl     maximum pcg iterations for island 
    819 !  rnu       strength of the additional force used in free surface b.c. 
    820 &namsol 
    821    nsolv     =      1 
    822    nsol_arp  =      0 
    823    nmin      =    300 
    824    nmax      =    800 
    825    nmod      =     10 
    826    eps       =  1.E-6 
    827    resmax    = 1.E-10 
    828    sor       =   1.92 
    829    epsisl    = 1.e-10 
    830    nmisl     =   4000 
    831    rnu       =     1. 
     578!!   nam_mpp           Massively Parallel Processing                    ("key_mpp_mpi) 
     579!!   nam_mpp_dyndist   Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist") 
     580!!   namctl            Control prints & Benchmark 
     581!!   namsol            elliptic solver / island / free surface  
     582!!====================================================================== 
     583 
     584!----------------------------------------------------------------------- 
     585&namsol        !   elliptic solver / island / free surface  
     586!----------------------------------------------------------------------- 
     587   nsolv       =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
     588                           !                   =2 successive-over-relaxation (sor) 
     589                           !                   =3 FETI (fet)                               ("key_feti") 
     590                           !                   =4 sor with extra outer halo 
     591   nsol_arp    =      0    !  absolute/relative (0/1) precision convergence test 
     592   nmin        =    300    !  minimum of iterations for the SOR solver 
     593   nmax        =    800    !  maximum of iterations for the SOR solver 
     594   nmod        =     10    !  frequency of test for the SOR solver 
     595   eps         =  1.e-6    !  absolute precision of the solver 
     596   resmax      =  1.e-10   !  absolute precision for the SOR solver 
     597   sor         =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
     598   epsisl      =  1.e-10   !  absolute precision on stream function solver 
     599   nmisl       =   4000    !  maximum pcg iterations for island                            ("key_islands") 
     600   rnu         =      1.   !  strength of the additional force used in filtered free surface 
     601/ 
     602!----------------------------------------------------------------------- 
     603&nam_mpp      !   Massively Parallel Processing                         ("key_mpp_mpi) 
     604!----------------------------------------------------------------------- 
     605   c_mpi_send =  'S'       !  mpi send/recieve type   ='S', 'B', or 'I' for standard send, 
     606                           !  buffer blocking send or immediate non-blocking sends, resp. 
     607   nn_buffer  =   0        !  size in bytes of exported buffer ('B' case), 0 no exportation 
     608/ 
     609!----------------------------------------------------------------------- 
     610&nam_mpp_dyndist !   Massively Parallel Distribution                    ("key_agrif" && "key_mpp_dyndist") 
     611!----------------------------------------------------------------------- 
     612   jpni       =    1       !  jpni   number of processors following i 
     613   jpnj       =    1       !  jpnj   number of processors following j 
     614   jpnij      =    1       !  jpnij  number of local domains 
     615/ 
     616!----------------------------------------------------------------------- 
     617&namctl       !   Control prints & Benchmark 
     618!----------------------------------------------------------------------- 
     619   ln_ctl     = .false.    !  trends control print (expensive!) 
     620   nprint     =    0       !  level of print (0 no extra print) 
     621   nictls     =    0       !  start i indice of control sum (use to compare mono versus 
     622   nictle     =    0       !  end   i indice of control sum        multi processor runs 
     623   njctls     =    0       !  start j indice of control               over a subdomain) 
     624   njctle     =    0       !  end   j indice of control  
     625   isplt      =    1       !  number of processors in i-direction 
     626   jsplt      =    1       !  number of processors in j-direction 
     627   nbench     =    0       !  Bench mode (1/0): CAUTION use zero except for bench 
     628                           !     (no physical validity of the results) 
     629   nbit_cmp   =    0       !  bit comparison mode (1/0): CAUTION use zero except for test 
     630                           !     of comparison between single and multiple processor runs 
    832631/ 
    833632!!====================================================================== 
     
    842641 
    843642!----------------------------------------------------------------------- 
    844 !       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra") 
     643&namtrd       !   diagnostics on dynamics and/or tracer trends             ("key_trddyn" and/or "key_trdtra") 
    845644!                          or mixed-layer trends                        ('key_trdmld') 
    846645!                          or barotropic vorticity                      ("key_trdvor") 
    847646!----------------------------------------------------------------------- 
    848 !  ntrd                time step frequency dynamics and tracers trends 
    849 !  nctls               control surface type in mixed-layer trends (0,1 or n<jpk) 
    850 !  ln_trdmld_restart   restart for ML diagnostics 
    851 !  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day) 
    852 !  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S 
    853 &namtrd 
    854    ntrd              = 365 
    855    nctls             =   0 
    856    ln_trdmld_restart = .false. 
    857    ucf               =  1. 
    858    ln_trdmld_instant = .false. 
    859 / 
    860 !----------------------------------------------------------------------- 
    861 !       namgap    level mean model-data gap ('key_diagap') 
    862 !----------------------------------------------------------------------- 
    863 !  ngap    time-step frequency of model-data gap computation 
    864 !  nprg    time-step frequency of gap print in model output 
    865 &namgap 
    866    ngap =  15 
    867    nprg =  10 
    868 / 
    869 !----------------------------------------------------------------------- 
    870 !       namspr  surface pressure diagnostic 
    871 !----------------------------------------------------------------------- 
    872 !  nmaxp   maximum of iterations for the solver 
    873 !  epsp    absolute precision of the solver 
    874 !  niterp  number of iteration done by the solver 
    875 &namspr 
    876    nmaxp   =   1000 
    877    epsp    =  1.e-3 
    878    niterp  =    400 
    879 / 
    880 !----------------------------------------------------------------------- 
    881 !       namflo    float parameters                                      ("key_float") 
    882 !----------------------------------------------------------------------- 
    883 !  ln_rstflo   boolean term for float restart (true or false) 
    884 !  nwritefl   frequency of float output file  
    885 !  nstockfl   frequency of float restart file  
    886 !  ln_argo    Argo type floats (stay at the surface each 10 days) 
    887 !  ln_flork4  = T trajectories computed with a 4th order Runge-Kutta 
    888 !             = F  (default)   computed with Blanke' scheme 
    889 &namflo 
    890     ln_rstflo = .false. 
    891     nwritefl  =      75 
    892     nstockfl  =    5475 
    893     ln_argo   = .false. 
    894     ln_flork4 = .false. 
    895 / 
    896 !----------------------------------------------------------------------- 
    897 !       namptr   Poleward Transport Diagnostic 
    898 !----------------------------------------------------------------------- 
    899 !  ln_diaptr  logical flag for Poleward transport computation 
    900 !  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
    901 !             need input basins mask file named "subbasins.nc" 
    902 !  nf_ptr     Frequency of computation 
    903 &namptr 
    904    ln_diaptr = .false. 
    905    ln_subbas = .false. 
    906    nf_ptr    =  15 
    907 / 
     647   ntrd       = 365        !  time step frequency dynamics and tracers trends 
     648   nctls      =   0        !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     649   ucf        =   1.       !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     650   cn_trdrst_in  = "restart_mld"   ! suffix of ocean restart name (input) 
     651   cn_trdrst_out = "restart_mld"   ! suffix of ocean restart name (output) 
     652   ln_trdmld_restart = .false.     !  restart for ML diagnostics 
     653   ln_trdmld_instant = .false.     !  flag to diagnose trends of instantantaneous or mean ML T/S 
     654/ 
     655!----------------------------------------------------------------------- 
     656&namgap       !   level mean model-data gap                             ('key_diagap') 
     657!----------------------------------------------------------------------- 
     658   ngap       =  15        !  time-step frequency of model-data gap computation 
     659   nprg       =  10        !  time-step frequency of gap print in model output 
     660/ 
     661!----------------------------------------------------------------------- 
     662&namspr       !   surface pressure diagnostic 
     663!----------------------------------------------------------------------- 
     664   nmaxp      =   1000     !  maximum of iterations for the solver 
     665   epsp       =  1.e-3     !  absolute precision of the solver 
     666   niterp     =    400     !  number of iteration done by the solver 
     667/ 
     668!----------------------------------------------------------------------- 
     669&namflo       !   float parameters                                      ("key_float") 
     670!----------------------------------------------------------------------- 
     671    ln_rstflo = .false.    !  float restart (T) or not (F) 
     672    nwritefl  =      75    !  frequency of writing in float output file  
     673    nstockfl  =    5475    !  frequency of creation of the float restart file  
     674    ln_argo   = .false.    !  Argo type floats (stay at the surface each 10 days) 
     675    ln_flork4 = .false.    !  trajectories computed with a 4th order Runge-Kutta (T) 
     676                           !  or computed with Blanke' scheme (F) 
     677/ 
     678!----------------------------------------------------------------------- 
     679&namptr       !   Poleward Transport Diagnostic 
     680!----------------------------------------------------------------------- 
     681   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
     682   ln_subbas  = .false.    !  Atlantic/Pacific/Indian basins computation (T) or not  
     683                           !  (orca configuration only, need input basins mask file named "subbasins.nc" 
     684   nf_ptr     =  15        !  Frequency of ptr computation [time step] 
     685/ 
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/PARAM/namelist_ice_ORCA2

    r396 r561  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/LIM2 :  1 - dynamics/advection/thermo          (namicerun) 
     3!! namelists    2 - ice intialisation                  (namiceini) 
     4!!              3 - ice dynamic                        (namicedyn) 
     5!!              5 - ice advection                      (namicetrp) 
     6!!              6 - thermodynamic                      (namicethd) 
     7!!              7 - ice salinity                       (namicesal) 
     8!!              8 - mechanical redistribution of ice   (namiceitdme) 
     9!!              3 - ice diagnostics                    (namicedia) 
     10!!              9 - ice outputs                        (namiceout) 
    111!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! 
    3 ! ICE namelist :  ice model option and parameter input 
    4 ! ------------- 
    5 ! 
    6 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    7 ! 
    8 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    9 !       namicerun   Share parameters for dynamics/advection/thermo 
    10 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    11 ! 
    12 !  ln_limdyn   : switch for ice dynamics (true) or not (false) 
    13 !  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere 
    14 !  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0) 
    15 !  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0) 
    16 ! 
    17 &namicerun 
    18    ln_limdyn = .true.  
    19    acrit     = 1.0e-06 , 1.0e-06 
    20    hsndif    = 0.0  
    21    hicdif    = 0.0 
     12 
     13!----------------------------------------------------------------------- 
     14&namicerun     !   Share parameters for dynamics/advection/thermo 
     15!----------------------------------------------------------------------- 
     16   cn_icerst_in  = "restart_ice_in"   !  suffix of ice restart name (input) 
     17   cn_icerst_out = "restart_ice"      !  suffix of ice restart name (output) 
     18   ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F) 
     19   ln_limdmp   = .false.   !  restoring ice thickness and fraction leads (T) or not (F) 
     20   acrit       = 1.0e-06   , 1.0e-06   ! minimum fraction for leads in the Northern (Southern) Hemisphere 
     21   hsndif      =  0.0      !  computation of temperature in snow (=0.0) or not 
     22   hicdif      =  0.0      !  computation of temperature in ice  (=0.0) or not (=9999.0) 
    2223/ 
    23 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    24 !       namiceini   parameters for ice initialisation 
    25 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    26 ! 
    27 !  ttest : threshold water temperature for initial sea ice 
    28 !  hninn : initial snow thickness in the north 
    29 !  hginn : initial ice thickness in the north 
    30 !  alinn : initial leads area in the north 
    31 !  hnins : initial snow thickness in the south 
    32 !  hgins : initial ice thickness in the south 
    33 !  alins : initial leads area in the south 
    34 ! 
    35 &namiceini 
    36    ttest = 2.0 
    37    hninn = 0.5 
    38    hginn = 3.0 
    39    alinn = 0.05 
    40    hnins = 0.1 
    41    hgins = 1.0 
    42    alins = 0.1 
     24!----------------------------------------------------------------------- 
     25&namiceini     !   ice initialisation 
     26!----------------------------------------------------------------------- 
     27   ln_limini   = .false.   !  read the ice initial state in the file 'Ice_initialization.nc' (T) or not (F) 
     28   ttest       =  2.0      !  threshold water temperature for initial sea ice 
     29   hninn       =  0.5      !  initial snow thickness in the north 
     30   hginn       =  3.0      !  initial ice  thickness in the north 
     31   alinn       =  0.05     !  initial leads area     in the north 
     32   hnins       =  0.1      !  same  three parameter  in the south 
     33   hgins       =  1.0      !        "                 "     south    
     34   alins       =  0.1      !        "                 "     south 
    4335/ 
    44 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    45 !       namicedia   parameters for ice diagnostics 
    46 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    47 ! 
    48 !  fmtinf : format of the output values  
    49 !  nfrinf : number of variables written in one line 
    50 !  ntmoy  : instantaneous values of ice evolution or averaging 
    51 !  ninfo  : frequency of ouputs on file ice_evolu in case of averaging 
    52 ! 
    53 &namicedia 
    54    fmtinf = '1PE13.5 ' 
    55    nfrinf = 4 
    56    ntmoy  = 1 
    57    ninfo  = 1 
     36!----------------------------------------------------------------------- 
     37&namicedyn     !   ice dynamic 
     38!----------------------------------------------------------------------- 
     39   epsd        =   1.0e-20 !  tolerance parameter 
     40   alpha       =   0.5     !  coefficient for semi-implicit coriolis 
     41   dm          =   0.6e+03 !  diffusion constant for dynamics 
     42   nbiter      =   1       !  number of sub-time steps for relaxation 
     43   nbitdr      = 100       !  maximum number of iterations for relaxation 
     44   om          =   0.5     !  relaxation constant  
     45   resl        =   5.0e-05 !  maximum value for the residual of relaxation 
     46   cw          =   5.0e-03 !  drag coefficient for oceanic stress 
     47   angvg       =   0.0     !  turning angle for oceanic stress 
     48   pstar       =   1.0e+04 !  1st bulk-rheology parameter 
     49   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
     50   etamn       =   0.0e+07 !  minimun value for viscosity 
     51   creepl      =   2.0e-08 !  creep limit 
     52   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
     53   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    5854/ 
    59 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    60 !       namicedyn   parameters for ice dynamic 
    61 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    62 ! 
    63 !  epsd   : tolerance parameter 
    64 !  alpha  : coefficient for semi-implicit coriolis 
    65 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    66 !  dm     : diffusion constant for dynamics. 
    67 !  nbiter : number of sub-time steps for relaxation 
    68 !  nbitdr : maximum number of iterations for relaxation 
    69 !  om     : relaxation constant  
    70 !  resl   : maximum value for the residual of relaxation  
    71 !  cw     : drag coefficient for oceanic stress 
    72 !  angvg  : turning angle for oceanic stress 
    73 !  pstar  : first bulk-rheology parameter 
    74 !  c_rhg  : second bulk-rhelogy parameter 
    75 !  etamn  : minimun value for viscosity 
    76 !  creepl : creep limit 
    77 !  ecc    : eccentricity of the elliptical yield curve 
    78 !  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s) 
    79 ! 
    80 &namicedyn  
    81    epsd   = 1.0e-20 
    82    alpha  = 0.5 
    83    dm     = 0.6e+03 
    84    nbiter = 1 
    85    nbitdr = 100 
    86    om     = 0.5 
    87    resl   = 5.0e-05 
    88    cw     = 5.0e-03 
    89    angvg  = 0.0 
    90    pstar  = 1.0e+04  
    91    c_rhg  = 20.0 
    92    etamn  = 0.0e+07 
    93    creepl = 2.0e-08 
    94    ecc    = 2.0 
    95    ahi0   = 350.e0 
     55!----------------------------------------------------------------------- 
     56&namicetrp     !   ice transport 
     57!----------------------------------------------------------------------- 
     58   bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    9659/ 
    97 ! 
    98 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    99 !       namicetrp   parameters for ice advection 
    100 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    101 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    102 &namicetrp  
    103    bound  = 0. 
    104 / 
    105  
    106 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    107 !       namicethd   parameters for thermodynamic computation 
    108 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    109 ! 
     60!----------------------------------------------------------------------- 
     61&namicethd     !   ice thermodynamic 
     62!----------------------------------------------------------------------- 
    11063!  hmelt       : maximum melting at the bottom 
    11164!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     
    12780!  alphs       : coefficient for snow density when snow ice formation 
    12881!  
    129 &namicethd 
    130    hmelt   = -0.15 
    131    hiccrit = 0.3 , 0.3 
    132    hicmin  = 0.2 
    133    hiclim  = 0.05 
    134    amax    = 0.999  
    135    swiqst  = 1. 
    136    sbeta   = 1. 
    137    parlat  = 0.0 
    138    hakspl  = 0.5 
    139    hibspl  = 0.5 
    140    exld    = 2.0 
    141    hakdif  = 1.0 
    142    thth    = 0.2 
    143    hnzst   = 0.1 
    144    parsub  = 0.0 
    145    alphs   = 1.0 
     82   hmelt       = -0.15     !  maximum melting at the bottom 
     83   hiccrit     = 0.3 , 0.3 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     84   hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket 
     85   hiclim      = 0.05      !  minimum ice thickness 
     86   amax        = 0.999     !  maximum lead fraction 
     87   swiqst      = 1.        !  energy stored in brine pocket (=1) or not (=0)     
     88   sbeta       = 1.        !  numerical caracteritic of the scheme for diffusion in ice 
     89   parlat      = 0.0       !  percentage of energy used for lateral ablation 
     90   hakspl      = 0.5       !  slope of distr. for Hakkinen-Mellor's lateral melting 
     91   hibspl      = 0.5       !  slope of distribution for Hibler's lateral melting 
     92   exld        = 2.0       !  exponent for leads-closure rate 
     93   hakdif      = 1.0       !  coefficient for diffusions of ice and snow 
     94   thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity 
     95   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
     96   parsub      = 0.0       !  switch for snow sublimation or not 
     97   alphs       = 1.0       !  coefficient for snow density when snow ice formation 
    14698/ 
     99!----------------------------------------------------------------------- 
     100&namicedia     !   ice diagnostics 
     101!----------------------------------------------------------------------- 
     102   fmtinf      ='1PE13.5 ' !  format of the output values 
     103   nfrinf      = 4         !  number of variables written in one line 
     104   ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
     105   ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     106/ 
     107!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     108&namiceout     !   parameters for outputs 
     109!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     110   noumef      =   19      !  number of fields 
    147111! 
    148 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    149 !       namiceout   parameters for outputs 
    150 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    151 ! 
    152 !   noumef      : number of fields 
    153 !   field 
    154 !        title  : title of the field 
    155 !        name   : name of the field 
    156 !        unit   : unit of the field 
    157 !        nc     : switch for saving field ( = 1 ) or not ( = 0 ) 
    158 !        cmulti : multiplicative constant 
    159 !        cadd   : additive constant 
    160 ! 
    161 &namiceout 
    162    noumef   = 19 
    163    field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0 
    164    field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0 
    165    field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ', 1, 1.0, 0.0 
    166    field_4  = 'Ice concentration                  ', 'ileadfra', '%       ', 1, -1.0, 1.0 
    167    field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ', 1, 1.0, -273.15 
    168    field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0 
    169    field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0 
    170    field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0 
    171    field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15 
    172    field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0 
    173    field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0 
    174    field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0 
    175    field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0 
    176    field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt', 1, 1.0, 0.0 
    177    field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0 
    178    field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0  
    179    field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0 
    180    field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0 
    181    field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0 
     112!           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
     113!           !                                      !           !           ! or not !  factor  !  factor  ! 
     114   field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
     115   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
     116   field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ',    1   ,  1.0     ,    0.0 
     117   field_4  = 'Ice concentration                  ', 'ileadfra', '%       ',    1   , -1.0     ,    1.0 
     118   field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
     119   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
     120   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
     121   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
     122   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
     123   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
     124   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     125   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     126   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
     127   field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt',    1   ,  1.0     ,    0.0 
     128   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
     129   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
     130   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     131   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     132   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ',    1   ,  1.0     ,    0.0 
    182133/       
  • CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_WORK/EXP00/config.card

    r396 r561  
    8080#D- frequency of rebuild submission 
    8181RebuildFrequency=NONE 
     82#D- Do we rebuild parallel output from archive 
     83RebuildFromArchive=true 
    8284#D- If you want to produce time series, this flag determines 
    8385#D- frequency of post-processing submission 
Note: See TracChangeset for help on using the changeset viewer.