Changeset 5586 for CONFIG/UNIFORM


Ignore:
Timestamp:
02/11/21 16:46:40 (3 years ago)
Author:
acosce
Message:

new configuration for LMDZORINCA coupled with Dynamico

Location:
CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7
Files:
70 added
58 deleted
10 edited

Legend:

Unmodified
Added
Removed

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/DRIVER/inca.driver

    r5543 r5586  
    419419    #  <context id="inca" src="./context_inca.xml"/> 
    420420    echo '<context id="inca" src="./context_inca.xml"/>' > add.tmp 
    421     echo '<context id="inca" src="./context_input_inca.xml"/>' >> add.tmp 
    422     echo '<context id="inca" src="./context_start_inca.xml"/>' >> add.tmp 
    423     echo '<context id="inca" src="./context_sflx_inca.xml"/>' >> add.tmp 
    424     echo '<context id="inca" src="./context_winds_inca.xml"/>' >> add.tmp 
    425     if [ X${CHEM} =  XGES ] ; then  
    426         echo '<context id="inca" src="./context_co2_bis_inca.xml"/>' >> add.tmp 
    427     fi 
    428     echo '<context id="inca" src="./context_oxydants_inca.xml"/>' >> add.tmp 
    429  
     421    if [ -f context_input_inca.xml ] ; then 
     422        echo '<context id="inca" src="./context_input_inca.xml"/>' >> add.tmp 
     423  #      echo '<context id="inca" src="./context_start_inca.xml"/>' >> add.tmp 
     424        echo '<context id="inca" src="./context_sflx_inca.xml"/>' >> add.tmp 
     425        echo '<context id="inca" src="./context_co2_bis_inca.xml"/>' >> add.tmp 
     426        echo '<context id="inca" src="./context_winds_inca.xml"/>' >> add.tmp 
     427        echo '<context id="inca" src="./context_oxydants_inca.xml"/>' >> add.tmp 
     428    fi 
    430429    cp iodef.xml iodef.xml.tmp 
    431430    sed -e "/COMPONENT CONTEXT/r add.tmp" \ 
    … …  
    449448 
    450449 
    451     echo "retire UNLIMITED sur axe temps " 
    452     nccopy -u u10xios.nc tmp.nc 
    453     mv tmp.nc u10xios.nc 
    454     nccopy -u v10xios.nc tmp.nc 
    455     mv tmp.nc v10xios.nc 
     450 
     451#    if [ X${inca_UserChoices_LMDZ_10m_winds} = Xn ] ; then  
     452        echo "retire UNLIMITED sur axe temps " 
     453        nccopy -u u10xios.nc tmp.nc 
     454        mv tmp.nc u10xios.nc 
     455        nccopy -u v10xios.nc tmp.nc 
     456        mv tmp.nc v10xios.nc 
     457#    fi 
     458 
    456459    nccopy -u oxydants.nc tmp.nc 
    457460    mv tmp.nc oxydants.nc 
     461 
     462 
    458463 
    459464 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/DRIVER/lmdz.driver

    r5542 r5586  
    55    IGCM_debug_PushStack "ATM_Initialize" 
    66 
    7     # Read the resolution from the variables ResolAtm set in config.card UserChoices section or from .resol file 
    8     if [ X$ResolAtm != X ] ; then 
    9         # ResolAtm is set in config.card 
    10         RESOL_ATM_3D=${ResolAtm} 
    11     elif [ -f ${SUBMIT_DIR}/../.resol ] ; then 
    12         # ResolAtm was not set in config.card, try to read .resol file 
    13         eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x39 
    14     else 
    15         IGCM_debug_Exit "ResolAtm is not set in config.card and the .resol file does not exist." 
    16         IGCM_debug_Verif_Exit 
    17     fi 
    18      
    19     IGCM_debug_Print 1 "The atomsopheric resolution is ${RESOL_ATM_3D} "  
    20  
    21     RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' ) 
    22     RESOL_ATM_X=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $1}' ) 
    23     RESOL_ATM_Y=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $2}' ) 
    24     RESOL_ATM_XY="${RESOL_ATM_X}x${RESOL_ATM_Y}" 
    25     RESOL_ATM=LMD${RESOL_ATM_X}${RESOL_ATM_Y} 
     7 
     8    # Define the resolution only if the composant ICO for DYNAMICO is not set in config.card 
     9    if [ X${config_ListOfComponents_ICO} = X ] ; then 
     10 
     11        # Read the resolution from the variables ResolAtm set in config.card UserChoices section or from .resol file 
     12        if [ X$ResolAtm != X ] ; then 
     13            # ResolAtm is set in config.card 
     14            RESOL_ATM_3D=${ResolAtm} 
     15        elif [ -f ${SUBMIT_DIR}/../.resol ] ; then 
     16            # ResolAtm was not set in config.card, try to read .resol file 
     17            eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x39 
     18        else 
     19            IGCM_debug_Exit "ResolAtm is not set in config.card and the .resol file does not exist." 
     20            IGCM_debug_Verif_Exit 
     21        fi 
     22         
     23        RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' ) 
     24        RESOL_ATM_X=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $1}' ) 
     25        RESOL_ATM_Y=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $2}' ) 
     26        RESOL_ATM_XY="${RESOL_ATM_X}x${RESOL_ATM_Y}" 
     27        RESOL_ATM=LMD${RESOL_ATM_X}${RESOL_ATM_Y} 
     28 
     29        IGCM_debug_Print 1 "The atomsopheric resolution is ${RESOL_ATM_3D} "  
     30    else 
     31        RESOL_ATM_Z=${dynamico_UserChoices_RESOL_ATM_Z} 
     32        RESOL_NBP=${dynamico_UserChoices_RESOL_NBP} 
     33    fi 
    2634 
    2735    ##- LMDZ physics version 
    … …  
    4856    ConfType=${lmdz_UserChoices_ConfType} 
    4957 
     58 
    5059    IGCM_debug_PopStack "ATM_Initialize" 
    5160} 
    … …  
    193202 
    194203 
    195  
    196204    ## Activate diagnostic output files and set output_level for each file 
    197205 
    … …  
    261269    # histdaystrataer 
    262270    if [ X${lmdz_UserChoices_output_level_histdaystrataer} = X ] || [ X${lmdz_UserChoices_output_level_histdaystrataer} = XNONE ] ; then 
    263         IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer enabled FALSE 
    264         IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer output_level 0 
    265     else 
    266         IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer enabled TRUE 
    267         IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer output_level ${lmdz_UserChoices_output_level_histstrataer} 
     271    IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer enabled FALSE 
     272    IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer output_level 0 
     273    else 
     274    IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer enabled TRUE 
     275    IGCM_comp_modifyXmlFile nonblocker file_def_histdaystrataer_lmdz.xml histdaystrataer output_level ${lmdz_UserChoices_output_level_histstrataer} 
    268276    fi 
    269277 
    … …  
    312320 
    313321    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def 
     322 
    314323    if [ ${CumulPeriod} -eq 1 ] ; then 
    315324        IGCM_comp_modifyDefFile blocker run.def raz_date  1 
    … …  
    319328 
    320329    ##-- Calendar type for LMDZ and create_etat0_limit 
     330    #    calend is read from both LMDZ dynamics and physics parts of the model 
    321331    case ${config_UserChoices_CalendarType} in 
    322332        leap|gregorian) 
    … …  
    330340    esac 
    331341 
     342 
     343    ## -- Modify more calendar related variables. 
    332344    IGCM_comp_modifyDefFile blocker run.def dayref    ${InitDay} 
     345    # Following variables are read by LMDZ dynamics  
    333346    IGCM_comp_modifyDefFile blocker run.def nday      ${PeriodLengthInDays} 
    334347 
    335     # Set anneeref different for gcm and ce0l 
    336     # Test if executable create_etat0_limit is present and lmdz.x is not present 
    337     if [ -f create_etat0_limit* ] || [ -f ce0l* ] && [ ! -f lmdz.x ] ; then 
    338         # for case ce0l : always take current year 
    339         IGCM_comp_modifyDefFile blocker run.def anneeref ${year} 
    340  
    341         # Temporary, set use_filtre_fft=n because we now use dyn3d for ce0l version where fft is not implemented.  
    342         # In more recent versions of LMDZ, this will not be needed.  
    343         IGCM_comp_modifyDefFile force   run.def use_filtre_fft n 
     348    # Set anneeref different for gcm and create_etat0_limit mode 
     349    # Test if keyword create_etat0_limit=y is set in dynamico.card or in lmdz.card 
     350    if [ X${lmdz_UserChoices_create_etat0_limit} = Xy ] || [ X${dynamico_UserChoices_create_etat0_limit} = Xy ] ; then 
     351       # for case create_etat0_limit : always take current year 
     352       IGCM_comp_modifyDefFile blocker run.def anneeref ${year} 
     353         
     354       # Temporary, set use_filtre_fft=n because we now use dyn3d for ce0l version where fft is not implemented.  
     355       # In more recent versions of LMDZ, this will not be needed.  
     356           IGCM_comp_modifyDefFile force   run.def use_filtre_fft n 
    344357    else 
    345358        # for case gcm : take first year of simulation 
    346359        IGCM_comp_modifyDefFile blocker run.def anneeref ${InitYear} 
    347360    fi 
     361     
    348362 
    349363    ## Determine from the variable ListOfComponents in config.card coupling to external models  
    … …  
    370384        echo "Activate coupling to INCA, set type_trac=inca in run.def" 
    371385        IGCM_comp_modifyDefFile blocker run.def type_trac inca 
    372         IGCM_comp_modifyDefFile blocker run.def iflag_phytrac 1   
     386        IGCM_comp_modifyDefFile blocker run.def iflag_phytrac 1 
    373387    elif [ X${config_ListOfComponents_CHM} = Xreprobus ] ; then 
    374388        echo "Activate coupling to REPROBUS, set type_trac=repr in run.def" 
    375389        IGCM_comp_modifyDefFile blocker run.def type_trac repr 
    376390        IGCM_comp_modifyDefFile blocker run.def config_inca none 
    377         IGCM_comp_modifyDefFile blocker run.def iflag_phytrac 1   
     391        IGCM_comp_modifyDefFile blocker run.def iflag_phytrac 1 
    378392    elif [ X${config_UserChoices_ConfigType} = XESMCO2 ] ; then  
    379393        echo "ESM config, set type_trac=co2i in run.def" 
    … …  
    503517    fi 
    504518 
     519 
    505520    # physiq.def : Modify parameter tau_gl if it is set in lmdz.card 
    506521    if [ ! X${lmdz_UserChoices_tau_gl} = X ]; then  
    … …  
    520535 
    521536 
    522     ## Read ByPass_hgardfou_teta option in lmdz.card --> divide teta* by 2 if [ $ByPass_hgardfou_teta = 1 ] 
    523     ByPass_hgardfou_teta=${lmdz_UserChoices_ByPass_hgardfou_teta} 
    524     if [ X"${ByPass_hgardfou_teta}" = X"y" ] ; then 
    525         awk '{ if ($0 ~ /^teta.*=/) {split($0,a,"=") ; print a[1]"="a[2]/2"."} else print $0}' gcm.def > gcm.def.tmp 
    526         IGCM_sys_Mv gcm.def.tmp gcm.def 
    527         echo 
    528         IGCM_debug_Print 1 "ByPass_hgardfou_teta : ^teta*/2 in gcm.def" 
    529         echo 
    530         cat gcm.def 
    531         ByPass_hgardfou_teta=n 
    532         IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_teta "${ByPass_hgardfou_teta}" 
    533     fi 
    534  
    535     ## Read ByPass_hgardfou_mats option in lmdz.card --> purmats=y 2 if [ $ByPass_hgardfou_mats = 1 ] 
    536     ByPass_hgardfou_mats=${lmdz_UserChoices_ByPass_hgardfou_mats} 
    537     if [ X"${ByPass_hgardfou_mats}" = X"y" ] ; then 
    538         sed -e "s/^purmats=.*/purmats=y/" gcm.def > gcm.def.tmp 
    539         IGCM_sys_Mv gcm.def.tmp gcm.def 
    540         echo 
    541         IGCM_debug_Print 1 "ByPass_hgardfou_mats : purmats=y in gcm.def" 
    542         echo 
    543         cat gcm.def 
    544         ByPass_hgardfou_mats=n 
    545         IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_mats "${ByPass_hgardfou_mats}" 
    546     fi 
     537    # Specific ByPass option only used with LMDZ dynamcis (without DYNMAICO) 
     538    if [ X${config_ListOfComponents_ICO} = X ] ; then 
     539 
     540        ## Read ByPass_hgardfou_teta option in lmdz.card --> divide teta* by 2 if [ $ByPass_hgardfou_teta = 1 ] 
     541        ByPass_hgardfou_teta=${lmdz_UserChoices_ByPass_hgardfou_teta} 
     542        if [ X"${ByPass_hgardfou_teta}" = X"y" ] ; then 
     543            awk '{ if ($0 ~ /^teta.*=/) {split($0,a,"=") ; print a[1]"="a[2]/2"."} else print $0}' gcm.def > gcm.def.tmp 
     544            IGCM_sys_Mv gcm.def.tmp gcm.def 
     545            echo 
     546            IGCM_debug_Print 1 "ByPass_hgardfou_teta : ^teta*/2 in gcm.def" 
     547            echo 
     548            cat gcm.def 
     549            ByPass_hgardfou_teta=n 
     550            IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_teta "${ByPass_hgardfou_teta}" 
     551        fi 
     552         
     553        ## Read ByPass_hgardfou_mats option in lmdz.card --> purmats=y 2 if [ $ByPass_hgardfou_mats = 1 ] 
     554        ByPass_hgardfou_mats=${lmdz_UserChoices_ByPass_hgardfou_mats} 
     555        if [ X"${ByPass_hgardfou_mats}" = X"y" ] ; then 
     556            sed -e "s/^purmats=.*/purmats=y/" gcm.def > gcm.def.tmp 
     557            IGCM_sys_Mv gcm.def.tmp gcm.def 
     558            echo 
     559            IGCM_debug_Print 1 "ByPass_hgardfou_mats : purmats=y in gcm.def" 
     560            echo 
     561            cat gcm.def 
     562            ByPass_hgardfou_mats=n 
     563            IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_mats "${ByPass_hgardfou_mats}" 
     564        fi 
     565    fi 
     566 
    547567 
    548568    # Add include of LMDZ context in iodef.xml 
    … …  
    575595        rm add.tmp 
    576596    fi 
     597             
    577598 
    578599    # Compression level (if CompressionLevel is not empty) 
    579600    if [ ! "X${config_UserChoices_CompressionLevel}" = "X" ] ; then 
    580     echo "NetCDF output files compression level is " ${config_UserChoices_CompressionLevel} 
    581     listfile=$(ls file_def*lmdz.xml) 
    582     for file in ${listfile} ; do 
    583         sed -i -e "s/\(compression_level=\"\)[^\"]*\(\"\)/\1${config_UserChoices_CompressionLevel}\2/" ${file} 
    584     done 
    585     fi 
    586  
    587              
     601        echo "NetCDF output files compression level is " ${config_UserChoices_CompressionLevel} 
     602        listfile=$(ls file_def*lmdz.xml) 
     603        for file in ${listfile} ; do 
     604            sed -i -e "s/\(compression_level=\"\)[^\"]*\(\"\)/\1${config_UserChoices_CompressionLevel}\2/" ${file} 
     605        done 
     606    fi 
     607 
     608 
     609    # Set interpolated output grid depending on resoultion, only if ICO component is set in config.card 
     610    if [ ! "X${config_ListOfComponents_ICO}" = "X" ] ; then    
     611        if [ ${RESOL_NBP} == 40  ] ; then 
     612            valuei=144 
     613            valuej=143 
     614        elif [ ${RESOL_NBP} == 80  ] ; then 
     615            valuei=360 
     616            valuej=180 
     617        elif [ ${RESOL_NBP} == 160  ] ; then 
     618            valuei=720 
     619            valuej=360 
     620    elif [ ${RESOL_NBP} == 320  ] ; then 
     621        valuei=1440 
     622        valuej=720 
     623        else 
     624            IGCM_debug_Exit "RESOL_NBP=${RESOL_NBP} is not yet implemented in drivers" 
     625            IGCM_debug_Verif_Exit 
     626        fi 
     627 
     628        IGCM_debug_Print 1 "Modify output grid parameters for current resolution in context_lmdz.xml" 
     629        cp context_lmdz.xml context_lmdz.xml.init 
     630        sed -e "s/ni_glo.*=.*\"144\"/ni_glo=\"${valuei}\"/" \ 
     631            -e "s/nj_glo.*=.*\"143\"/nj_glo=\"${valuej}\"/" \ 
     632            context_lmdz.xml.init > context_lmdz.xml 
     633    fi 
     634 
     635 
    588636    IGCM_debug_PopStack "ATM_Update" 
    589637} 
    … …  
    599647} 
    600648 
     649 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/DRIVER/orchidee.driver

    r5542 r5586  
    3636        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate1 output_level 0 
    3737        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate1 output_freq 1mo 
     38 
    3839        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate2 enabled .FALSE. 
    3940        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate2 output_level 0 
    4041        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate2 output_freq 1mo 
    41         IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 enabled .FALSE. 
    42         IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 output_level 0 
    43         IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 output_freq 1mo 
     42 
     43        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 enabled .FALSE. 
     44        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 output_level 0 
     45        IGCM_comp_modifyXmlFile nonblocker file_def_orchidee.xml stomate3 output_freq 1mo 
    4446    fi 
    4547 
    … …  
    138140    fi 
    139141         
     142    # Desactivate orchidee variables which are not correct in ICOLMDZOR configuration 
     143    IGCM_comp_modifyXmlFile force field_def_orchidee.xml RESOLUTION_X enabled .FALSE. 
     144    IGCM_comp_modifyXmlFile force field_def_orchidee.xml RESOLUTION_Y enabled .FALSE. 
     145    IGCM_comp_modifyXmlFile force field_def_orchidee.xml CONTFRAC_STOMATE enabled .FALSE. 
     146 
    140147 
    141148    # Add include of orchidee context in iodef.xml 
    … …  
    169176    fi 
    170177 
     178 
    171179    # Compression level (if CompressionLevel is not empty) 
    172180    if [ ! "X${config_UserChoices_CompressionLevel}" = "X" ] ; then 
    173181        echo "NetCDF output files compression level is " ${config_UserChoices_CompressionLevel} 
    174         listfile=$(ls file_def*orchidee.xml) 
    175         for file in ${listfile} ; do 
    176             sed -i -e "s/\(compression_level=\"\)[^\"]*\(\"\)/\1${config_UserChoices_CompressionLevel}\2/" ${file} 
    177         done 
    178     fi 
    179  
    180  
    181  
     182        listfile=$(ls file_def*orchidee.xml) 
     183        for file in ${listfile} ; do 
     184            sed -i -e "s/\(compression_level=\"\)[^\"]*\(\"\)/\1${config_UserChoices_CompressionLevel}\2/" ${file} 
     185        done 
     186    fi 
     187 
     188 
     189    # Set interpolated output grid depending on resoultion, only if ICO component is set in config.card 
     190    if [ ! "X${config_ListOfComponents_ICO}" = "X" ] ; then    
     191        if [ ${RESOL_NBP} == 40  ] ; then 
     192            valuei=144 
     193            valuej=143 
     194        elif [ ${RESOL_NBP} == 80  ] ; then 
     195            valuei=360 
     196            valuej=180 
     197        elif [ ${RESOL_NBP} == 160  ] ; then 
     198            valuei=720 
     199            valuej=360 
     200    elif [ ${RESOL_NBP} == 320  ] ; then 
     201        valuei=1440 
     202        valuej=720 
     203        else 
     204            IGCM_debug_Exit "RESOL_NBP=${RESOL_NBP} is not yet implemented in drivers" 
     205            IGCM_debug_Verif_Exit 
     206        fi 
     207         
     208        IGCM_debug_Print 1 "Modify output grid parameters for current resolution in context_orchidee.xml" 
     209        cp context_orchidee.xml context_orchidee.xml.init 
     210        sed -e "s/ni_glo.*=.*\"144\"/ni_glo=\"${valuei}\"/" \ 
     211            -e "s/nj_glo.*=.*\"143\"/nj_glo=\"${valuej}\"/" \ 
     212            context_orchidee.xml.init > context_orchidee.xml 
     213         
     214    fi 
     215     
    182216    IGCM_debug_PopStack "SRF_Update" 
    183217} 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/DRIVER/xios.driver

    r5542 r5586  
    2525 
    2626    if [ -f namcouple ] ; then 
    27         sed -e "s/2  LMDZ oceanx/3 LMDZ oceanx ${config_Executable_IOS[1]}/" namcouple > namcouple.tmp  
     27        sed -e "s/2  icosa oceanx/3 icosa oceanx ${config_Executable_IOS[1]}/" namcouple > namcouple.tmp  
    2828        IGCM_sys_Mv namcouple.tmp namcouple 
    2929    fi 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/PARAM/iodef.xml

    r5542 r5586  
    3030          <variable id="print_file" type="bool">false</variable> 
    3131          <variable id="using_oasis" type="bool">false</variable> 
    32           <variable id="oasis_codes_id" type="string" >LMDZ,oceanx</variable> 
     32          <variable id="oasis_codes_id" type="string" >icosa,oceanx</variable> 
    3333        </variable_group> 
    3434      <variable id="activate_non_distributed_transformation" type="bool">true</variable> 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/PARAM/ping_orchidee.xml

    r5542 r5586  
    184184   <field id="CMIP6_irrLut"        field_ref="dummy_XY_landUse_not_provided" /> <!-- P1 (kg s-1) surface_downward_water_flux_due_to_irrigation : Irrigation flux including any irrigation for crops, trees, pasture, or urban lawns --> 
    185185   <field id="CMIP6_ksat"          field_ref="ksat"       /> <!-- P1 (1e-6 m s-1)  : Saturated Hydraulic Conductivity --> 
    186    <field id="CMIP6_lai"           field_ref="lai_ipcc"         /> <!-- P1 (1) leaf_area_index : a ratio obtained by dividing the total upper leaf surface area of vegetation by the (horizontal) surface area of the land on which it grows. Warning: here we take the average over the entire land fraction; not sure it's ok with the DataRequest comment "a ratio obtained by dividing the total upper leaf surface area of vegetation by the (horizontal) surface area of the land on which it grows." PC+JG : lai_ipcc and laimean are related to the same land fraction(note veget_max in stomate_lpj is not the same veget_max in sechiba!). We see other differencies between these output variables coming mainly from pft3 and some from 6,8,9,12. --> 
     186   <field id="CMIP6_lai"           field_ref="LAImean"         /> <!-- P1 (1) leaf_area_index : a ratio obtained by dividing the total upper leaf surface area of vegetation by the (horizontal) surface area of the land on which it grows. Warning: here we take the average over the entire land fraction; not sure it's ok with the DataRequest comment "a ratio obtained by dividing the total upper leaf surface area of vegetation by the (horizontal) surface area of the land on which it grows." PC+JG : lai_ipcc and laimean are related to the same land fraction(note veget_max in stomate_lpj is not the same veget_max in sechiba!). We see other differencies between these output variables coming mainly from pft3 and some from 6,8,9,12. --> 
    187187   <field id="CMIP6_laiLut"        field_ref="lailut" /> <!-- P1 (1)  : Leaf Area Index on Land Use Tile --> 
    188    <field id="CMIP6_landCoverFrac" field_ref="landCoverFrac" /> <!-- P1 (%) area_fraction : The categories may differ from model to model, depending on their  PFT definitions.   This may include natural PFTs, anthropogenic PFTs, bare soil, lakes, urban areas, etc.   Sum of all should equal the fraction of the grid-cell that is land. Warning: Note: here not multiplying by contfrac according to DR --> 
     188   <field id="CMIP6_landCoverFrac" field_ref="maxvegetfrac" > maxvegetfrac*100. </field> <!-- P1 (%) area_fraction : The categories may differ from model to model, depending on their  PFT definitions.   This may include natural PFTs, anthropogenic PFTs, bare soil, lakes, urban areas, etc.   Sum of all should equal the fraction of the grid-cell that is land. Warning: Note: here not multiplying by contfrac according to DR --> 
    189189   <field id="CMIP6_lwsnl"         field_ref="snowliqtot"         /> <!-- P1 (kg m-2) liquid_water_content_of_snow_layer : where land over land: this is computed as the total mass of liquid water contained interstitially within the snow layer of the land portion of a grid cell divided by the area of the land portion of the cell. Warning: Computed only for the layered scheme on vegetated and bare soil area. --> 
    190190   <field id="CMIP6_mrfso"         field_ref="humtot_frozen"         /> <!-- P1 (kg m-2) soil_frozen_water_content : the mass (summed over all all layers) of frozen water. Warning: To be calculated. --> 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/PARAM/run.def

    r5542 r5586  
    11# $Id$ 
    22# 
     3INCLUDEDEF=run_dynamico.def 
     4INCLUDEDEF=earth_const.def 
    35INCLUDEDEF=physiq.def 
    46INCLUDEDEF=gcm.def 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/PARAM/vert.def_L79

    r5542 r5586  
    66## The new version can be controled by following parameters 
    77## Al the parameters are given in km assuming a given scalehigh 
     8##  
     9## This file is read by both LMDZ and DYNAMICO.  
     10## Parameters in section 1 are only read by LMDZ when running without DYNAMICO 
     11## Parameters in section 2 are only read by DYNAMCIO 
     12## Parameters in section 3 are shared between LMDZ/dynamics and DYNAMCIO 
     13 
     14## Section 1: Following paramereters are read only when running LMDZ without DYNAMCIO 
     15##################################################################################### 
    816vert_sampling=strato_custom 
     17#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche 
     18iflag_top_bound=2 
     19#  Coefficient pour la couche eponge (valeur derniere couche) 
     20tau_top_bound=5.e-4 
     21#  mode_top_bound (defaults: 3)  
     22mode_top_bound=3 
     23 
     24## Section 2: Following paramereters are read only when running with DYNAMCIO 
     25############################################################################# 
     26# Number of vertical layer (llm) : integer (default=19) 
     27llm=79 
     28# disvert : vertical discretisation : string (default='std') : std, ncar, ncar30l 
     29disvert=strato_custom 
     30# sponge parmeters 
     31iflag_sponge=2 
     32tau_sponge=5.e-4 
     33mode_sponge=1 
     34 
     35## Section 3: Following paramereters are shared between LMDZ/dynamics and DYNAMICO 
     36##            They will be read either by LMDZ or DYNAMCIO 
     37################################################################################## 
    938## scale hight 
    1039vert_scale=7. 
    … …  
    3261dissip_zref = 50 
    3362##  Avec ou sans strato 
     63#   ok_strato est lu dans LMDZ/dynamique et phylmd (necessaire avec ou sans DYNAMICO) 
    3464ok_strato=y 
    3565##  Activation ou non de la parametrisation de Hines pour la strato 
     66# ok_hines est lu par LMDZ/phylmd (avec ou sans DYNAMICO) 
    3667ok_hines=n 
    37 #  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche 
    38 iflag_top_bound=2 
    39 #  Coefficient pour la couche eponge (valeur derniere couche) 
    40 tau_top_bound=5.e-4 
    41 # 
     68 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/POST/monitoring01_lmdz.cfg

    r5542 r5586  
    3939#  field | files patterns | files additionnal | operations | title | units | calcul of area 
    4040#----------------------------------------------------------------------------------------------------------------- 
    41  bils_global    | "bils"                       | "" | "bils[d=1]"               | "Surf. total heat flux (GLOBAL)"         | "W/m^2"     | "aire[d=1]"    
    42  nettop_global  | "tops topl"                  | "" | "(tops[d=1]-topl[d=2])"   | "TOA. total heat flux (GLOBAL)"          | "W/m^2"     | "aire[d=1]" 
    43  nettop0_global | "tops0 topl0"                | "" | "(tops0[d=1]-topl0[d=2])" | "TOA. total heat flux ClearSky (GLOBAL)" | "W/m^2"     | "aire[d=1]" 
    44  t2m_global_prio     | "t2m"                        | "" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (GLOBAL)"       | "degrees C" | "aire[d=1]" 
    45  t2m_land       | "t2m pourc_ter pourc_lic"    | "" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (LAND)"         | "degrees C" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
    46  tsol_global    | "tsol"                       | "" | "(tsol[d=1]-273.15)"      | "Surface Temperature (GLOBAL)"           | "degrees C" | "aire[d=1]" 
    47  tsol_land      | "tsol pourc_ter pourc_lic"   | "" | "(tsol[d=1]-273.15)"      | "Surface Temperature (LAND)"             | "degrees C" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
    48  precip_global  | "precip"                     | "" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (GLOBAL)"  | "mm/d"      | "aire[d=1]" 
    49  precip_land    | "precip pourc_ter pourc_lic" | "" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (LAND)"    | "mm/d"      | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
    50  evap_global    | "evap"                       | "" | "evap[d=1]"               | "Evaporation (GLOBAL)"                   | "kg/(s*m2)" | "aire[d=1]" 
    51  evap_land      | "evap pourc_ter pourc_lic"   | "" | "evap[d=1]"               | "Evaporation (LAND)"                     | "kg/(s*m2)" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
    52  flat_global    | "flat"                       | "" | "flat[d=1]"               | "Surf. latent heat flux (GLOBAL)"         | "W/m^2"     | "aire[d=1]" 
    53  prw_global     | "prw"                        | "" | "prw[d=1]"                | "Precipitable water (GLOBAL)"            | "kg/m^2"    | "aire[d=1]" 
    54  ener_global    | "nettop precip snow sols soll sens" | "" | "nettop[d=1,x=@ave,y=@ave]+(precip[d=2,x=@ave,y=@ave]-snow[d=3,x=@ave,y=@ave])*2.5008E+6+snow[d=3,x=@ave,y=@ave]*2.8345E+6-sols[d=4,x=@ave,y=@ave]-soll[d=5,x=@ave,y=@ave]-sens[d=6,x=@ave,y=@ave]" | "Destruction d energie par l atmosphere" | "W/m^2" | "1" | 
    55  bils_corrected_global  | "bils precip snow flat" | "" | "bils[d=1,x=@ave,y=@ave]-flat[d=4,x=@ave,y=@ave]-(precip[d=2,x=@ave,y=@ave]-snow[d=3,x=@ave,y=@ave])*2.5008E+6-snow[d=3,x=@ave,y=@ave]*2.8345E+6" | "Corrected energy budget at surface" | "W/m^2" | "1" | 
    56  vitu_qbo               | "vitu"                  | "" | "vitu[y=-5:5@ave,x=-180:180@ave,z=5000,d=1]"                | "vitu at 5000 Pa for QBO +5deg/-5deg"            | "m/s"    | "1" 
    57  od550aer_tropo_forcing | "od550aer"              | "" | "od550aer[d=1]"           | "Total aerosol optical depth at 550nm"   | "-"         | "aire[d=1]"    
    58  od550aer_strat_forcing | "od550_STRAT"           | "" | "od550_STRAT[d=1]"        | "Aerosol optical depth at 550nm"         | "-"         | "aire[d=1]"    
    59  od550lt1aer_forcing    | "od550lt1aer"           | "" | "od550lt1aer[d=1]"        | "Fine mode optical depth"                | "-"         | "aire[d=1]"    
     41 bils_global    | "bils"                       | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "bils[d=1]"               | "Surf. total heat flux (GLOBAL)"         | "W/m^2"     | "aire[d=2]"    
     42 nettop_global  | "tops topl"                  | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(tops[d=1]-topl[d=2])"   | "TOA. total heat flux (GLOBAL)"          | "W/m^2"     | "aire[d=3]" 
     43 nettop0_global | "tops0 topl0"                | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(tops0[d=1]-topl0[d=2])" | "TOA. total heat flux ClearSky (GLOBAL)" | "W/m^2"     | "aire[d=3]" 
     44 t2m_global_prio     | "t2m"                        | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (GLOBAL)"       | "degrees C" | "aire[d=2]" 
     45 t2m_land       | "t2m pourc_ter pourc_lic"    | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (LAND)"         | "degrees C" | "aire[d=4]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
     46 tsol_global    | "tsol"                       | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(tsol[d=1]-273.15)"      | "Surface Temperature (GLOBAL)"           | "degrees C" | "aire[d=2]" 
     47 tsol_land      | "tsol pourc_ter pourc_lic"   | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "(tsol[d=1]-273.15)"      | "Surface Temperature (LAND)"             | "degrees C" | "aire[d=4]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
     48 precip_global  | "precip"                     | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (GLOBAL)"  | "mm/d"      | "aire[d=2]" 
     49 precip_land    | "precip pourc_ter pourc_lic" | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (LAND)"    | "mm/d"      | "aire[d=4]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
     50 evap_global    | "evap"                       | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "evap[d=1]"               | "Evaporation (GLOBAL)"                   | "kg/(s*m2)" | "aire[d=2]" 
     51 evap_land      | "evap pourc_ter pourc_lic"   | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "evap[d=1]"               | "Evaporation (LAND)"                     | "kg/(s*m2)" | "aire[d=4]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)" 
     52 flat_global    | "flat"                       | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "flat[d=1]"               | "Surf. latent heat flux (GLOBAL)"         | "W/m^2"     | "aire[d=2]" 
     53 prw_global     | "prw"                        | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "prw[d=1]"                | "Precipitable water (GLOBAL)"            | "kg/m^2"    | "aire[d=2]" 
     54 ener_global    | "nettop precip snow sols soll sens" | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "nettop[d=1,x=@ave,y=@ave]+(precip[d=2,x=@ave,y=@ave]-snow[d=3,x=@ave,y=@ave])*2.5008E+6+snow[d=3,x=@ave,y=@ave]*2.8345E+6-sols[d=4,x=@ave,y=@ave]-soll[d=5,x=@ave,y=@ave]-sens[d=6,x=@ave,y=@ave]" | "Destruction d energie par l atmosphere" | "W/m^2" | "1" | 
     55 bils_corrected_global  | "bils precip snow flat" | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "bils[d=1,x=@ave,y=@ave]-flat[d=4,x=@ave,y=@ave]-(precip[d=2,x=@ave,y=@ave]-snow[d=3,x=@ave,y=@ave])*2.5008E+6-snow[d=3,x=@ave,y=@ave]*2.8345E+6" | "Corrected energy budget at surface" | "W/m^2" | "1" | 
     56 vitu_qbo               | "vitu"                  | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "vitu[y=-5:5@ave,x=-180:180@ave,z=5000,d=1]"                | "vitu at 5000 Pa for QBO +5deg/-5deg"            | "m/s"    | "1" 
     57 od550aer_tropo_forcing | "od550aer"              | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "od550aer[d=1]"           | "Total aerosol optical depth at 550nm"   | "-"         | "aire[d=2]"    
     58 od550aer_strat_forcing | "od550_STRAT"           | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "od550_STRAT[d=1]"        | "Aerosol optical depth at 550nm"         | "-"         | "aire[d=2]"    
     59 od550lt1aer_forcing    | "od550lt1aer"           | "/ccc/work/cont003/igcmg/igcmg/IGCM/ATM/GRID/aire_ICO40_to_144x143.nc" | "od550lt1aer[d=1]"        | "Fine mode optical depth"                | "-"         | "aire[d=2]"    
    6060 tro3_strato_forcing_DU | "colO3_strat"           | "" | "colO3_strat[x=@ave,y=@ave,d=1]"     | "colonne ozone strato" | "DU" | "" 
    6161 tro3_tropo_forcing_DU  | "colO3_trop"            | "" | "colO3_trop[x=@ave,y=@ave,d=1]"      | "colonne ozone trop"   | "DU" | "" 

  • CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/POST/monitoring01_orchidee.cfg

    r5542 r5586  
    3939#  field | files patterns | files additionnal | operations | title | units | calcul of area 
    4040#------------------------------------------------------------------------------------------------------------------------------------------------------ 
    41 delta_water_stock | "nobiofrac DelSoilMoist_daily DelIntercept_daily DelSWE_daily delstock_routing_daily" |  "" | "((1-nobiofrac[d=1,l=2:100000])*DelSoilMoist_daily[d=2,l=2:100000] + DelIntercept_daily[d=3,l=2:100000] + DelSWE_daily[d=4,l=2:100000] + delstock_routing_daily[d=5,l=2:100000])" | "Change in total water stock (LANDS) " | "kg/day/m2 of land" | "Areas[d=2]*Contfrac[d=2]" 
    42 water_budget_closure | "nobiofrac DelSoilMoist_daily DelIntercept_daily DelSWE_daily delstock_routing_daily rain snowf evap riverflow coastalflow" | "" | "(((1-nobiofrac[d=1,l=2:100000])*DelSoilMoist_daily[d=2,l=2:100000] + DelIntercept_daily[d=3,l=2:100000] + DelSWE_daily[d=4,l=2:100000] +  delstock_routing_daily[d=5,l=2:100000])*Areas[d=1]*Contfrac[d=1] -  ((rain[d=6,l=2:100000] + snowf[d=7,l=2:100000] - evap[d=8,l=2:100000])*Areas[d=1]*Contfrac[d=1]  - (riverflow[d=9,l=2:100000] + coastalflow[d=10,l=2:100000])*1000*86400 ))/(4*3.14*6370000*6370000*0.26)" | "Water budget closure (LANDS) " | "kg/day/m2 of land" | "2" 
    43 twbr_lands       | "TWBR"             |  "" | "twbr[d=1,l=2:100000]*86400"       | "Hydrol water buget residu (LANDS)" | "kg/m^2/d" | "Areas[d=1]*Contfrac[d=1]" 
    44 temp_sol_lands   | "temp_sol"         |  "" | "temp_sol[d=1]"                    | "TSOL (LANDS)"                      | "C"        | "Areas[d=1]*Contfrac[d=1]" 
    45 fluxsens_lands   | "fluxsens"         |  "" | "fluxsens[d=1]"                    | "Sensible Flux (LANDS)"             | "W/m^2"    | "Areas[d=1]*Contfrac[d=1]" 
    46 evap_lands       | "evap"             |  "" | "evap[d=1]"                        | "Evaporation (LANDS)"               | "mm/d"     | "Areas[d=1]*Contfrac[d=1]" 
    47 evapnu_lands     | "evapnu"           |  "" | "evapnu[d=1]"                      | "Soil evaporation (LANDS)"          | "mm/d"     | "Areas[d=1]*Contfrac[d=1]" 
    48 transpir_lands   | "transpir"         |  "" | "transpir[d=1,k=@sum]"             | "Transpiration (sum over all PFT, LANDS)" | "mm/d" | "Areas[d=1]*Contfrac[d=1]" 
    49 subli_lands      | "subli"            |  "" | "subli[d=1]"                       | "Sublimation (LANDS)"               | "mm/d"     | "Areas[d=1]*Contfrac[d=1]" 
    50 snow_swe_lands   | "frac_snow snow"   |  "" | "snow[d=2]"                        | "Snow Mass (LANDS)"                 | "kg/m^2"   | "Areas[d=2]*Contfrac[d=2]" 
    51 snownobio_lands  | "snownobio"        |  "" | "snownobio[d=1]"                   | "Snow Other Surfaces (LANDS)"       | "kg/m^2"   | "Areas[d=1]*Contfrac[d=1]" 
    52 snowf_lands      | "snowf"            |  "" | "snowf[d=1]"                       | "Snowfall (LANDS)"                  | "mm/d"     | "Areas[d=1]*Contfrac[d=1]" 
    53 precip_lands     | "rain snowf"       |  "" | "rain[d=1] + snowf[d=2]"           | "Total precipitaion (LANDS)"        | "mm/d"     | "Areas[d=1]*Contfrac[d=1]" 
    54 runoff_lands     | "runoff"           |  "" | "runoff[d=1]"                      | "Surface runoff  (LANDS)"           | "kg/m^2"   | "Areas[d=1]*Contfrac[d=1]" 
    55 drainage_lands   | "drainage"         |  "" | "drainage[d=1]"                    | "Drainage  (LANDS)"                 | "kg/m^2"   | "Areas[d=1]*Contfrac[d=1]" 
    56 humidity_lands   | "mrso"             |  "" | "mrso[d=1]"                        | "Total Soil Moisture  (LANDS)"      | "kg/m^2"   | "Areas[d=1]*Contfrac[d=1]" 
    57 surface_humidity_lands | "mrsos"      |  "" | "mrsos[d=1]"                       | "Top 10cm soil moisture  (LANDS)"   | "kg/m^2"   | "Areas[d=1]*Contfrac[d=1]" 
    58 albedo_lands     | "alb_vis alb_nir"  |  "" | "((alb_vis[d=1]+alb_nir[d=2])*0.5)"| "Mean Albedo (LANDS)"               | "1"        | "Areas[d=1]*Contfrac[d=1]" 
    59 lai_lands        | "LAImean"          |  "" | "LAImean[d=1]"                     | "LAI (mean excluding PFT=1) (LANDS)"| "1"        | "Areas[d=1]*Contfrac[d=1]" 
    60 nee_lands        | "nee maxvegetfrac" |  "" | "(nee[d=1]*365*86400*maxvegetfrac[d=2]*Areas[d=1]*Contfrac[d=1]/1e15)"   | "Net Ecosystem Exchange" | "PgC/yr" | "2" 
    61 maxvegetagri_forcing | "maxvegetfrac" |  "" | "maxvegetfrac[d=1,K=12:13@sum]"    | "veget_max fraction of agriculture [sum of pft=12,13] "| "1"  | "Areas[d=1]*Contfrac[d=1]" 
    62 maxvegetbare_forcing | "maxvegetfrac" |  "" | "maxvegetfrac[d=1,K=1]"            | "veget_max fraction on bare soil [pft=1] "             | "1"  | "Areas[d=1]*Contfrac[d=1]" 
    63 vegetn_lands     | "vegetfrac"        |  "" | "vegetfrac[d=1,K=2:9@sum]"         | "veget trees[sum of pft=2:9] (LANDS)"                  | "1"  | "Areas[d=1]*Contfrac[d=1]" 
    64 vegetg_lands     | "vegetfrac"        |  "" | "vegetfrac[d=1,K=10:11@sum]+vegetfrac[d=1,K=14:15@sum]"       | "veget natural grasses[sum of pft=10,11,14,15] (LANDS)"      | "1"  | "Areas[d=1]*Contfrac[d=1]" 
    65 vegeta_lands     | "vegetfrac"        |  "" | "vegetfrac[d=1,K=12:13@sum]"       | "veget agriculture [sum of pft=12,13] (LANDS)"         | "1"  | "Areas[d=1]*Contfrac[d=1]" 
     41delta_water_stock | "nobiofrac DelSoilMoist_daily DelIntercept_daily DelSWE_daily delstock_routing_daily" |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "((1-nobiofrac[d=1,l=2:100000])*DelSoilMoist_daily[d=2,l=2:100000] + DelIntercept_daily[d=3,l=2:100000] + DelSWE_daily[d=4,l=2:100000] + delstock_routing_daily[d=5,l=2:100000])" | "Change in total water stock (LANDS) " | "kg/day/m2 of land" | "Areas[d=6]*Contfrac[d=2]" 
     42water_budget_closure | "nobiofrac DelSoilMoist_daily DelIntercept_daily DelSWE_daily delstock_routing_daily rain snowf evap riverflow coastalflow" | "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "(((1-nobiofrac[d=1,l=2:100000])*DelSoilMoist_daily[d=2,l=2:100000] + DelIntercept_daily[d=3,l=2:100000] + DelSWE_daily[d=4,l=2:100000] +  delstock_routing_daily[d=5,l=2:100000])*Areas[d=11]*Contfrac[d=1] -  ((rain[d=6,l=2:100000] + snowf[d=7,l=2:100000] - evap[d=8,l=2:100000])*Areas[d=11]*Contfrac[d=1]  - (riverflow[d=9,l=2:100000] + coastalflow[d=10,l=2:100000])*1000*86400 ))/(4*3.14*6370000*6370000*0.26)" | "Water budget closure (LANDS) " | "kg/day/m2 of land" | "2" 
     43twbr_lands       | "TWBR"             |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "twbr[d=1,l=2:100000]*86400"       | "Hydrol water buget residu (LANDS)" | "kg/m^2/d" | "Areas[d=2]*Contfrac[d=1]" 
     44temp_sol_lands   | "temp_sol"         |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "temp_sol[d=1]"                    | "TSOL (LANDS)"                      | "C"        | "Areas[d=2]*Contfrac[d=1]" 
     45fluxsens_lands   | "fluxsens"         |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "fluxsens[d=1]"                    | "Sensible Flux (LANDS)"             | "W/m^2"    | "Areas[d=2]*Contfrac[d=1]" 
     46evap_lands       | "evap"             |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "evap[d=1]"                        | "Evaporation (LANDS)"               | "mm/d"     | "Areas[d=2]*Contfrac[d=1]" 
     47evapnu_lands     | "evapnu"           |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "evapnu[d=1]"                      | "Soil evaporation (LANDS)"          | "mm/d"     | "Areas[d=2]*Contfrac[d=1]" 
     48transpir_lands   | "transpir"         |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "transpir[d=1,k=@sum]"             | "Transpiration (sum over all PFT, LANDS)" | "mm/d" | "Areas[d=2]*Contfrac[d=1]" 
     49subli_lands      | "subli"            |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "subli[d=1]"                       | "Sublimation (LANDS)"               | "mm/d"     | "Areas[d=2]*Contfrac[d=1]" 
     50snow_swe_lands   | "frac_snow snow"   |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "snow[d=2]"                        | "Snow Mass (LANDS)"                 | "kg/m^2"   | "Areas[d=3]*Contfrac[d=2]" 
     51snownobio_lands  | "snownobio"        |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "snownobio[d=1]"                   | "Snow Other Surfaces (LANDS)"       | "kg/m^2"   | "Areas[d=2]*Contfrac[d=1]" 
     52snowf_lands      | "snowf"            |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "snowf[d=1]"                       | "Snowfall (LANDS)"                  | "mm/d"     | "Areas[d=2]*Contfrac[d=1]" 
     53precip_lands     | "rain snowf"       |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "rain[d=1] + snowf[d=2]"           | "Total precipitaion (LANDS)"        | "mm/d"     | "Areas[d=3]*Contfrac[d=1]" 
     54runoff_lands     | "runoff"           |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "runoff[d=1]"                      | "Surface runoff  (LANDS)"           | "kg/m^2"   | "Areas[d=2]*Contfrac[d=1]" 
     55drainage_lands   | "drainage"         |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "drainage[d=1]"                    | "Drainage  (LANDS)"                 | "kg/m^2"   | "Areas[d=2]*Contfrac[d=1]" 
     56humidity_lands   | "mrso"             |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "mrso[d=1]"                        | "Total Soil Moisture  (LANDS)"      | "kg/m^2"   | "Areas[d=2]*Contfrac[d=1]" 
     57surface_humidity_lands | "mrsos"      |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "mrsos[d=1]"                       | "Top 10cm soil moisture  (LANDS)"   | "kg/m^2"   | "Areas[d=2]*Contfrac[d=1]" 
     58albedo_lands     | "alb_vis alb_nir"  |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "((alb_vis[d=1]+alb_nir[d=2])*0.5)"| "Mean Albedo (LANDS)"               | "1"        | "Areas[d=3]*Contfrac[d=1]" 
     59lai_lands        | "LAImean"          |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "LAImean[d=1]"                     | "LAI (mean excluding PFT=1) (LANDS)"| "1"        | "Areas[d=2]*Contfrac[d=1]" 
     60nee_lands        | "nee maxvegetfrac" |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "(nee[d=1]*365*86400*maxvegetfrac[d=2]*Areas[d=3]*Contfrac[d=1]/1e15)"   | "Net Ecosystem Exchange" | "PgC/yr" | "2" 
     61maxvegetagri_forcing | "maxvegetfrac" |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "maxvegetfrac[d=1,K=12:13@sum]"    | "veget_max fraction of agriculture [sum of pft=12,13] "| "1"  | "Areas[d=2]*Contfrac[d=1]" 
     62maxvegetbare_forcing | "maxvegetfrac" |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "maxvegetfrac[d=1,K=1]"            | "veget_max fraction on bare soil [pft=1] "             | "1"  | "Areas[d=2]*Contfrac[d=1]" 
     63vegetn_lands     | "vegetfrac"        |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "vegetfrac[d=1,K=2:9@sum]"         | "veget trees[sum of pft=2:9] (LANDS)"                  | "1"  | "Areas[d=2]*Contfrac[d=1]" 
     64vegetg_lands     | "vegetfrac"        |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "vegetfrac[d=1,K=10:11@sum]+vegetfrac[d=1,K=14:15@sum]"       | "veget natural grasses[sum of pft=10,11,14,15] (LANDS)"      | "1"  | "Areas[d=2]*Contfrac[d=1]" 
     65vegeta_lands     | "vegetfrac"        |  "/ccc/work/cont003/igcmg/igcmg/IGCM/SRF/GRID/Areas_ICO40_to_144x143.nc" | "vegetfrac[d=1,K=12:13@sum]"       | "veget agriculture [sum of pft=12,13] (LANDS)"         | "1"  | "Areas[d=2]*Contfrac[d=1]" 
    6666 
Note: See TracChangeset for help on using the changeset viewer.