...
casper% cp -r Gen_Data_NEWGRID Gen_Data_ne0np4.SAM01.ne30x4 casper% cd Gen_Data_ne0np4.SAM01.ne30x4 casper% mv Gen_MERRA2_NEWGRID.csh Gen_MERRA2_ne0np4.SAM01.ne30x4.csh |
|---|
...
- Edit your new script, and set the PBS commands:
#!/bin/csh module load ncl |
|---|
After making the changes below you can submit this job on casper:
casper% qsub Gen_MERRA2_ne0np4.SAM01.ne30x4.csh
2. In the Configuration section, set the reference date corresponding to the first day of data you desire, then number of days of data to process from that date, and the path where you wish to have the data stored:
#============================================================= # CONFIGURATION SECTION: #============================================================= # Set a REFERENCE (Starting) Date and the number of days to process #-------------------------------------------------------------------------------------------------------- set RUNNUM=01 set REF_DATE='20121201' set NUM_DAYS=400 # Set INPUT/OUTPUT/TMP directories #-------------------------------------------------------- set NAMELIST='./Config/Config_makeIC-'$RUNNUM'.nl' set MYLOGDIR='./LOG/LOG_002.'$RUNNUM'/' set MYTMPDIR='./TMP/TMP_002.'$RUNNUM'/' set MYOUTDIR='/path/to/my/repo/ne0np4.SAM01.ne30x4/nudging/MERRA2_ne0np4.SAM01.ne30x4_L32/' set INPUTDIR='/glade/collections/rda/data/ds313.3/orig_res/' # Set ESMF options #--------------------------- set ESMF_interp='conserve' set ESMF_pole='none' set ESMF_clean='False' set ESMF_clean='True' |
|---|
3. For the processing options, set the CASE name. This is the root filename for your nudging data files. Note that some reanalysis datasets store winds in the form of vorticity and divergence values rather then U,V. It is important that the VORT_DIV_TO_UV flag is set to True for these datasets, this is a common source of processing errors. Finally, set the fname_grid_info value to point to a file containing the desired ouput grid, and set fname_phis_output to point to a file containing the model topography.
# Set Processing options #------------------------------------- set CASE = 'MERRA2_ne0np4.SAM01.ne30x4_L32' set DYCORE = 'se' set PRECISION = 'float' set VORT_DIV_TO_UV = 'False' set SST_MASK = 'False' set ICE_MASK = 'False' set OUTPUT_PHIS = 'True' set REGRID_ALL = 'False' set ADJUST_STATE_FROM_TOPO = 'True' set MASS_FIX = 'True' # Set files containing OUTPUT Grid structure and topography #-------------------------------------------------------------------------------------------- set fname_grid_info = '/path/to/my/repo/ne0np4.SAM01.ne30x4/inic/cami-mam4_0000-01-01_ne0np4.SAM01.ne30x4_L32_c200309.nc' set fname_phis_output = '/path/to/my/repo/ne0np4.SAM01.ne30x4/topo/topo_ne30np4.SAM01.ne30x4_blin_200309.nc' set ftype_phis_output = 'SE_TOPOGRAPHY' |
|---|
4. Update the path to WRAPIT: /glade/u/apps/opt/ncl/6.6.2/bin/WRAPIT MAKEIC.stub MAKEIC.f90
The last modification that is needed for this example is to edit 5. Edit the file makeIC_se_002.ncl. At about line 430, the dstGridName variable has to be set to use the SCRIP file for your new grid:
dstGridName="/path/to/my/repo/ne0np4.SAM01.ne30x4/grids/SAM01_ne30x4_SCRIP.nc" |
|---|
After making the changes above you can submit this job on casper:
casper% qsub Gen_MERRA2_ne0np4.SAM01.ne30x4.csh
With these changes, the script can be submitted to generate the desired 400 days of data. The user must alway check the Log and output file to verify that the dataset was processed correctly. Selecting some files and browsing variables for valid values using ncdump can save a lot of time and effort if there were problems during the processing.
...