CROCO-WW3

19.4.5.3. CROCO-WW3#

In your ${CHOME} repository you should have already filled in myenv_mypath.sh.

To prepare the run you need to modify the files myjob.sh and mynamelist.sh.

  • In myjob.sh , you will have to fill in information about dates, job sequence:

    # Real job duration in sec (converted to MACHINE format in submit job)
    export TIMEJOB=1800
    
    #-------------------------------------------------------------------------------
    # Run date settings
    #-------------------------------------------------------------------------------
    # Your run can be divided into several jobs (e.g.: 1 year run into 12 jobs of 1 month)
    
    # Start date of the first Job
    export YEAR_BEGIN_JOB=2005
    export MONTH_BEGIN_JOB=1
    export DAY_BEGIN_JOB=1
    
    # Duration of each Job
    export JOB_DUR_MTH=1
    export JOB_DUR_DAY=0
    
    # How many jobs do you want to launch?
    export NBJOB=1
    
    # Do we start from a restart?
    export RESTART_FLAG="FALSE"
    

    Along with the number of cpu you will use for each model

    # nb of CPUs for each model
    export NP_OCEX=2
    export NP_OCEY=2
    export NP_WAV=14
    
  • In mynamelist.sh, first specify the run type, and the name of the experiment:

    # Run type (o/a/w, w.Afrc, oa, 2o1a, owa, owa.full...)
    #  - Will select the models to use reading letters o/w/a/toywav/toyoce/toyatm
    #  - Will select the executables, and some options (see in the following sections)
    #  - In coupled mode corresponds to the suffix of the OASIS_IN/namcouple.base.$RUNtype to use
    export RUNtype=ow
    export MOD=`echo $RUNtype | cut -d . -f 1`
    
    #  Name of the experiment you are about to launch (max 30. CHAR)
    export CEXPER=BENGUELA_example_${RUNtype}
    

    Then, there is a section indicating where the run will be executed and where the outputs and restarts will be stored:

    #-------------------------------------------------------------------------------
    # RUN_DIR
    #-------------------------------------------------------------------------------
    
    export EXEDIR_ROOT="$CWORK/rundir/${CEXPER}_execute"
    export OUTPUTDIR_ROOT="$CWORK/rundir/${CEXPER}_outputs"
    export RESTDIR_ROOT="$CWORK/rundir/${CEXPER}_restarts"
    
    export  JOBDIR_ROOT=${CHOME}/jobs_${CEXPER}
    

    Then, there are sections for the different components.

    For the coupler settings:

    #-------------------------------------------------------------------------------
    # CPL
    #-------------------------------------------------------------------------------
    
    # Namelist
    #---------
    # Note: namelist example files are provided in OASIS_IN/
    # if you want to use a pre-built weight file for grid interpolations, point to
    # e.g. namcouple.base.oa.smtho2a
    export namcouplename=namcouple.base.${RUNtype}
    
    # Coupling frequency
    #-------------------
    export CPL_FREQ=21600
    
    # Restart files for OASIS
    #------------------------
    # If TRUE: create OASIS restart files from pre-existing atm/oce/wav outputs.
    # If FALSE: create OASIS restart files from calm conditions (need to read at least the grid for each model)
    export CPL_restart="FALSE"
    export oce_rst_file="${OCE_FILES_DIR}/croco_grd.nc"
    export oce_rst_timeind=-1 # time index (-1 is last) in the file to extract as restart
    export atm_rst_file="${ATM_FILES_DIR}/wrfinput_dXX_2005_01_01_00" # the domain dXX will be automatically replaced
    export atm_rst_timeind=-1 # time index (-1 is last) in the file to extract as restart
    export wav_rst_file="${WAV_FILES_DIR}/ww3.200501.nc"
    export wav_rst_timeind=-1 # time index (-1 is last) in the file to extract as restart
    

    You should check the coupling frequency, the restart flag and path towards model files to use to create oasis restart files.

    For CROCO settings, first indicate if you request CROCO compilation:

    #-------------------------------------------------------------------------------
    # OCE
    #-------------------------------------------------------------------------------
    
    # Where to find or put the croco execuatble
    export OCE_EXE_DIR=${CHOME}/CROCO_IN
    
    # Online Compilation
    #-------------------
    #!!!!!!! IMPORTANT NOTE !!!!!!!
    # If activated, creates croco executable depending on this namelist.
    #   - In param.h it modifies the grid size, the number of procs in x and y direction with those given in myjob.sh
    #   - In cppdefs.h it modifies the following options with informations given below
    #       MPI, OA_COUPLING, OW_COUPLING, MRL_WCI,
    #       XIOS, LOGFILE, MPI_NOLAND,
    #       AGRIF, AGRIF_2WAY,
    #       BULK_FLUX, ONLINE, AROME, ARPEGE, ERA_ECMWF
    #       FRC_BRY, CLIMATOLOGY
    #       TIDES, PSOURCE, PSOURCE_NCFILE, PSOURCE_NCFILE_TS
    # Other changes of parameterizations, numerical schemes, etc should be made "by hand" in CROCO_IN/cppdefs.h.base
    #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    export ONLINE_COMP=1
    

    Then, edit model time steps:

    # Time steps
    #-----------
    export DT_OCE=3600
    export NDTFAST=60
    

    Then several options for zooms, wave coupling are provided. They are generally automatically set-up depending on the RUNtype defined at the beginning.

    Then, edit the options regarding the forcing files:

    # Forcings
    #---------
    export ini_ext='ini_SODA' # ini extension file (ini_SODA,...)
    export bdy_ext='bry_SODA' # bry extension file (clm_SODA,bry_SODA,...)
    
    # flag for surface forcing should be true except in the case of atm coupling
    if [[ $MOD =~ .*a.* ]] ; then
        export surfrc_flag="FALSE"
    else
        export surfrc_flag="TRUE"
    fi
    export interponline=0 # switch (1=on, 0=off) for online surface interpolation. Only works with MONTHLY input files!
    export frc_ext='blk_ERA5' # surface forcing extension(blk_ERA5, frc_ERA5,...). If interponline=1 precise the type (ERA_ECMWF or AROME,  [CFSR by default], names as cppkey name in croco)
    
    export tide_flag="FALSE" # the forcing extension must be blk_??? otherwise tide forcing overwrites it
    export river_flag="FALSE"
    

    Finally edit CROCO output settings:

    # Output settings
    #----------------
    #!!! WARNING: when XIOS is activated the following values (for the model) are not taken into account
    export oce_his_sec=86400     # history output interval (in number of second)
    export oce_avg_sec=86400     # average output interval (in number of second)
    

    Then go to the WAVE model section:

    #-------------------------------------------------------------------------------
    # WAV
    #-------------------------------------------------------------------------------
    
    # Where to find the wav executable
    if [[ $RUNtype =~ .*owa.* ]] ; then
        export WAV_EXE_DIR=${WAV}/exe_owa
    elif [[ $RUNtype =~ .*ow.* ]] ; then
        export WAV_EXE_DIR=${WAV}/exe_ow
    elif [[ $RUNtype =~ .*aw.* ]] ; then
        export WAV_EXE_DIR=${WAV}/exe_aw
    else
        export WAV_EXE_DIR=${WAV}/exe_frc
    fi
    
    # Namelist
    #---------
    # Chosing the ww3_shel.inp.base.SHELL_EXT (see options in WW3_IN)
    if [[ $RUNtype =~ .*toy.* ]] ; then
        export SHELL_EXT=$MOD
    else
        export SHELL_EXT=$RUNtype
    fi
    

    First the paths for the executable are defined, and the ww3_shel namelist to use is defined. Those are generally set-up with the RUNtupe but can be changed if necessary.

    Then edit the ww3 time steps and grid settings (these will be updated in ww3_grid.inp):

    # Time steps
    #-----------
    export DT_WAV=3600     # TMAX = 3*TCFL
    export DT_WW_PRO=1200  # TCFL = 0.8 x dx/(g/fmin4pi) with fmin=0.0373 => 3-4 % of dx
    export DT_WW_REF=1800  # TMAX / 2
    export DT_WW_SRC=10    # TSRC = usually 10s  (could be between 5s and 60s)
    
    # Grid size
    #----------
    export wavnx=41 ; export wavny=42
    
    # Parameter
    #----------
    export hmin=75; # e.g. minimum water depth in CROCO (will be used to delimit coastline in WW3)
    

    Then, choose the forcing to use for ww3:

    # Forcing files
    #--------------
    # forcin: forcing file(s) PREFIX list (input file are supposed to be in the form: PREFIX_Y????M??.nc)
    # forcww3: name of ww3_prnc.inp extension, e.g current or wind/era5, see in WW3_IN directory
    if [[ $RUNtype =~ .*owa.* ]]; then
        export forcin=()
        export forcww3=()
    elif [[ $RUNtype =~ .*Afrc.* || $RUNtype =~ .*ow.* ]] ; then
        export forcin=(ERA5_wind)
        export forcww3=(wind.era5)
    elif [[ $RUNtype =~ .*Ofrc.* ]] ; then
        export forcin=(CROCO_current CROCO_level)
        export forcww3=(current level)
    elif [[ $RUNtype =~ .*frc.* ]] ; then
        export forcin=(ERA5_wind CROCO_current CROCO_level)
        export forcww3=(wind.era5 current level)
    fi
    

    As well as the boundary data (can be left empty):

    # Boundary files
    #---------------
    # prefix for boundary files (leave empty is none), there are supposed to be in WAV_FILES_DIR
    export bouncin=
    

    Finally, set output settings:

    # Output settings
    #----------------
    export wav_int=21600  # output interval (s)
    export wav_pnt=0  # point output interval. Put 0 if no point output
    export point_output_list=${WAV_FILES_DIR}/my_point_output_test.txt # file where to find list of point (format: lon lat name) to output spectrum
    export wav_trck=0     # track output interval. Put 0 if no track output
    export flagout="TRUE" # Keep (TRUE) or not (FALSE) the ww3 output binary files (e.g. out_grd.ww3)
    
  • Now that you have completed the necessary files, you are ready to run your simulation. To do so, simply do

    ./submitjob.sh