Regridding for the new Landice Component: Difference between revisions

Initial add up to regrid
 
Add the regrid from MERRA section
Line 109: Line 109:
  turb_import_rst.e19901128_21z.Ganymed-1_0_p2.Ganymed-1_0_M_CF0090x6C_DE0360xPE0180
  turb_import_rst.e19901128_21z.Ganymed-1_0_p2.Ganymed-1_0_M_CF0090x6C_DE0360xPE0180
  turb_internal_rst.e19901128_21z.Ganymed-1_0_p2.Ganymed-1_0_M_CF0090x6C_DE0360xPE0180
  turb_internal_rst.e19901128_21z.Ganymed-1_0_p2.Ganymed-1_0_M_CF0090x6C_DE0360xPE0180
== Step 2: Regrid from MERRA ==
The next step is to use Joe Stassi's regrid.pl script to create some restarts at 21z on 1990-NOV-28. We'll first run two sets of runs at 144x91 and c48 to get a command to use for automated generation.
=== Run regrid.pl for 144x91 ===
$ cd $ESMADIR/Linux/bin
$ ./regrid.pl
<nowiki>Regrid from MERRA restarts (y/n)? [n] y
Enter ymd (yyyymmdd): 19901128
Enter hour (hh): 21
----------------
Atmosphere Grids
----------------
Lat/Lon Grids      Cubed-Sphere Grids
-------------      ------------------
b = 2 deg          C48    C360    C1000
c = 1 deg          C90    C500    C1440
d = 1/2 deg        C180    C720    C2880
e = 1/4 deg
Getting INPUT atmosphere grid resolution from fvcore file ... Done
INPUT atmosphere grid: D
Enter OUTPUT atmosphere grid: b
Ocean Grids
-----------
c =  360x180  (Reynolds)
f = 2880x1440  (OSTIA)
INPUT ocean grid: c
Enter OUTPUT ocean grid: [c]
GCM BCS tags
------------
Fortuna-1_4
Fortuna-1_5
Fortuna-2_0
Fortuna-2_1
Ganymed-1_0
Ganymed-1_0_m1
Ganymed-1_0_m2
Ganymed-1_0_M
Ganymed-1_0_D
Enter OUTPUT GCM BCS tag: [Ganymed-1_0_M]
Regrid upperair restarts (y/n)? [y]
Regrid surface restarts (y/n)? [y]
Enter OUTPUT directory for new restarts: /gpfsm/dnb31/mathomp4/Restarts-G20/144x91
Enter Experiment ID for OUTPUT restarts [b_d5_merra_jan89]
Write rst.lcv restart (y/n)? [y] n
Append tag/grid label to names of final restarts (y/n)? [n]
===============
REGRID RESTARTS
===============
ESMABIN: /gpfsm/dswdev/mathomp4/Models/Ganymed-2_0_UNSTABLE/GEOSagcm/Linux/bin
ESMATAG: Ganymed-2_0_UNSTABLE
logfile: /gpfsm/dnb31/mathomp4/Restarts-G20/144x91/b_d5_merra_jan89.19901128_21z.log
upper-air: yes
surface:  yes
bkg/satb:  no
rst.lcv:  no
namelabel: no
# inputs
#-------
. expid:      d5_merra_jan89
. date:        19901128
. hour:        21
. atmos grid:  540x361 (D)
. ocean grid:  360x180 (c)
. GCM tag:    Fortuna-1_4
. rstdir:      /archive/g_proj5/production/GEOSdas-2_1_4/d5_merra_jan89/rs/Y1990/M11
# intermediate-1 (surface only)
#---------------
. expid:      b_d5_merra_jan89-mid1
. date:        19901128
. hour:        21
. atmos grid:  576x361 (d)
. ocean grid:  360x180 (c)
. GCM tag:    Fortuna-2_1
. rescale:    no
. outdir:      /gpfsm/dnb31/mathomp4/Restarts-G20/144x91
. workdir:    /gpfsm/dnb31/mathomp4/Restarts-G20/144x91/19901128_7086
# intermediate-2 (surface only)
#---------------
. expid:      b_d5_merra_jan89-mid2
. date:        19901128
. hour:        21
. atmos grid:  576x361 (d)
. ocean grid:  360x180 (c)
. GCM tag:    Ganymed-1_0_m1
. rescale:    yes
. outdir:      /gpfsm/dnb31/mathomp4/Restarts-G20/144x91
. workdir:    /gpfsm/dnb31/mathomp4/Restarts-G20/144x91/19901128_7086
# outputs
#--------
. expid:      b_d5_merra_jan89
. date:        19901128
. hour:        21
. atmos grid:  144x91 (b)
. ocean grid:  360x180 (c)
. GCM tag:    Ganymed-1_0_M
. rescale:    no
. outdir:      /gpfsm/dnb31/mathomp4/Restarts-G20/144x91
. workdir:    /gpfsm/dnb31/mathomp4/Restarts-G20/144x91/19901128_7086
===============
Continue (y/n)? [y]
...<snip></nowiki>
When this succeeds you'll have a file, b_d5_merra_jan89.19901128_21z.CMD, which details how the regrid.pl ran:
$ cat b_d5_merra_jan89.19901128_21z.CMD
#!/usr/local/bin/csh
#===============================================================================
# The following command was used to create the b_d5_merra_jan89.19901128_21z restarts
#===============================================================================
set ESMABIN = /gpfsm/dswdev/mathomp4/Models/Ganymed-2_0_UNSTABLE/GEOSagcm/Linux/bin
$ESMABIN/regrid.pl -ymd 19901128 -hr 21 -grout b -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/144x91 -nobkg -nolbl -nolcv \
                    -tagout Ganymed-1_0_M -newid b_d5_merra_jan89 -rs 3 -oceanout c
'''Note:''' The above last command is actually on one line, but I've split it for display's sake.
=== Run regrid.pl for c48 ===
Do a similar run for c48 and get out this command:
$ cat C48_d5_merra_jan89.19901128_21z.CMD
#!/usr/local/bin/csh
#===============================================================================
# The following command was used to create the C48_d5_merra_jan89.19901128_21z restarts
#===============================================================================
set ESMABIN = /gpfsm/dswdev/mathomp4/Models/Ganymed-2_0_UNSTABLE/GEOSagcm/Linux/bin
$ESMABIN/regrid.pl -ymd 19901128 -hr 21 -grout C48 -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/c48 -nobkg -nolbl -nolcv \
                    -tagout Ganymed-1_0_M -newid C48_d5_merra_jan89 -rs 3 -oceanout c
=== Regrid all the other resolutions ===
We now have commands we can use to regrid to other resolutions. Create a simple script to do so:
<nowiki>$ cat regrid.latlon.csh
#!/usr/local/bin/csh
set ESMABIN = /gpfsm/dswdev/mathomp4/Models/Ganymed-2_0_UNSTABLE/GEOSagcm/Linux/bin
# 144x91
# ------
$ESMABIN/regrid.pl -merra -ymd 19901128 -hr 21 -grout b -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/144x91 -nobkg -nolbl -nolcv \
                    -tagout Ganymed-1_0_M -newid b_d5_merra_jan89 -rs 3 -oceanout c
# 288x181
# -------
$ESMABIN/regrid.pl -merra -ymd 19901128 -hr 21 -grout c -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/288x181 -nobkg -nolbl -nolcv \
                  -tagout Ganymed-1_0_M -newid c_d5_merra_jan89 -rs 3 -oceanout c
# 576x361
# -------
$ESMABIN/regrid.pl -merra -ymd 19901128 -hr 21 -grout d -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/576x361 -nobkg -nolbl -nolcv \
                  -tagout Ganymed-1_0_M -newid d_d5_merra_jan89 -rs 3 -oceanout c
# 1152x721
# --------
$ESMABIN/regrid.pl -merra -ymd 19901128 -hr 21 -grout e -outdir /gpfsm/dnb31/mathomp4/Restarts-G20/1152x721 -nobkg -nolbl -nolcv \
                  -tagout Ganymed-1_0_M -newid e_d5_merra_jan89 -rs 3 -oceanout c</nowiki>
and something similarly ugly for the cubed-sphere.
== Step 3: Regrid the good Landice restarts ==
Now we need to work at getting a good LANDICE restart. Again we'll use Joe's script.