Vous regardez une version antérieure (v. /display/EDDSDLTEL/GenPhysX) de cette page.

afficher les différences afficher l'historique de la page

« Afficher la version précédente Vous regardez la version actuelle de cette page. (v. 4) afficher la version suivante »

At UQAM GenPhysX is currently only working on alea.

To get access to GenPhysX (which uses SPI) you need to load:

      module load utils/SPI

As input GenPhysX needs to know the model grid in form of the tictacs of the full model grid. On our UQAM servers you can create them from your gem_settings.nml with the command:

     ~winger/ovbin/k.grille

The file you need as input for GenPhysX is called "tape1_model".

For the different available datasets check here.

When using USGS set:

    TOPO=USGS
    VEGE=USGS_R
    MASK=USGS_R
    SOIL=JPL

When using CCI-LC set:

    TOPO=GMTED75  # missing Antarctica and Greenland!!!
or
    TOPO=GMTED30  # lower resolution

    VEGE=CCI_LC
    MASK=CCI_LC
    SOIL=BNU

Example script to create geophysical fields:


#!/bin/bash

#  Set the following parameters
exp=domain_name    # Domain name (name of geophys files)
cd /Path/of/Geophys/directory/Geophys_${exp}   # Destination directory containing already the file 'tape1.model'!!!
                                           #   -> must NOT be in the home!!!
email="your.email@address"                # Email address to which a confirmation email will get send after completion

# ======================================================================================

module load utils/SPI
geophys=geophys_${exp}
GENPHYSX=GenPhysX
tictacs=tape1_model

# Set following parameters to dataset of your choice, following above examples:

TOPO=...
VEGE=...
MASK=...
SOIL=...

# Using Ayrton's new sharp low pass filter:
cat > parameter_file.txt << EOF
set GenX::Param(DBase)                      /BIG3/data/Geophys/Database/Dorval

set GenX::Settings(GRD_TYP_S)               LU
set GenX::Settings(TOPO_DGFMS_L)            False
set GenX::Settings(TOPO_DGFMX_L)            False
set GenX::Settings(TOPO_FILMX_L)            False
set GenX::Settings(TOPO_CLIP_ORO_L)         False
#                         
set Log::Param(Level)                       DEBUG
set GeoPhysX::Opt(LegacyMode)               False
set GenX::Param(Z0Topo)                     LEGACY
set GenX::Settings(TOPO_ZREF_ZV_RATIO_C)    True
set GenX::Settings(TOPO_RUGV_ZVG2)          True
set GenX::Settings(TOPO_RUGV_MG)            True
set GenX::Settings(TOPO_ZV_MIN_THRESHOLD)   0.0003
set GenX::Settings(TOPO_WATER_RUGV)         0.001

# ne pas filter le champs MENF a l entree du calcul de Z0, ZP
set GenX::Param(MEFilterForZ0)              None

# Ayrton's new sharp low pass filter
set GenX::Settings(LPASSFLT_RC_DELTAX)         r53.0
set GenX::Settings(LPASSFLT_P)                 20
set GenX::Settings(LPASSFLT_MASK_OPERATOR)     1
set GenX::Settings(LPASSFLT_MASK_THRESHOLD)    100.0
set GenX::Settings(LPASSFLT_APPLY_MINMAX)      True

# Use CCI LC's distribution
set GenX::Path(CCILC_LUT_CSV)              /home/winger/Scripts/GenPhysX/CCI_LC/CCI_LC_lut.csv
set GenX::Path(USGS_GLCC_LUT_CSV)          /home/winger/Scripts/GenPhysX/CCI_LC/USGS_GLCC_lut2.csv
set GenX::Path(NALCMS_LUT_CSV)             /home/winger/Scripts/GenPhysX/CCI_LC/NALCMS_lut2.csv
set GenX::Param(UseVegeLUT)                True
#set GenX::Param(AddHydroLakesToMask)       True

EOF

# Execute GenPhysX
$GENPHYSX -gridfile ${tictacs} -param parameter_file.txt -mask $MASK -vege $VEGE -topo $TOPO -soil $SOIL -check STD -cel
ldim 2 -z0topo LEGACY -mefilter LPASS -result ${geophys} -subgrid STD  -mail ${email}






  • Aucune étiquette