Add new workflows that leverage enroot+pyxis

Author: fluidnumericsJoe

.github/workflows/docker-publish.yml

@@ -12,6 +12,7 @@ CEOP
 
 #ifndef EXF_OPTIONS_H
 #define EXF_OPTIONS_H
+#endif
 #include "PACKAGES_CONFIG.h"
 #include "CPP_OPTIONS.h"
 

Dockerfile

@@ -43,17 +43,14 @@
  nIter0=0,
  baseTime=0.0,
 # End at end of day 365
-# endTime=31536000.0,
- endTime=120.0
+ endTime=31536000.0,
  deltaT=120.0,
  abEps=0.01,
  pChkptFreq=31536000.0, 
  chkptFreq=3153600.0,
  dumpFreq=120.0,
  monitorSelect=2,
- monitorFreq=120.0,
-# monitorFreq=3600.0,
- periodicExternalForcing=.FALSE.,
+ monitorFreq=3600.0,
 # Daily forced
 # externForcingPeriod=86400.0,
 # 365 day year

simulations/glorysv12-curvilinear/code/EXF_OPTIONS.h

@@ -1,10 +1,7 @@
 &CAL_NML
   TheCalendar = 'gregorian',
   calendarDumps = .TRUE.,
-  startDate_1 = 20020601,
+  startDate_1 = 20020701,
   startDate_2 = 000000,
-# TODO : Update end date for simulation
-  endDate_1   = 20101231,
-  endDate_2   = 235959,
 /
 

simulations/glorysv12-curvilinear/input/data

@@ -1,5 +1,5 @@
 ! -----------------------------------------------------------------------------
-! data.exf -- example for ERA-Interim 6h, 0.75deg, and a curvilinear model grid
+! data.exf -- example for ERA-Interim 3h, 0.25deg, and a curvilinear model grid
 ! Edit file names, dates, periods, and the lon/lat metadata to match your files.
 ! -----------------------------------------------------------------------------
 
@@ -25,111 +25,116 @@
 
 &EXF_NML_02
   ! --- filenames, startdates and periods (seconds) ---
+
   ! Wind components (10 m)
-  uwindfile       = 'uwind.bin',
-  uwindstartdate1 = 20020701,
-  uwindstartdate2 = 000000,
-  uwindperiod     = 10800.0,     ! 3 hours
-  uwind_lon0 = -90.0 
-  uwind_lon_inc = 0.25
-  uwind_lat0 = 20.0
-  uwind_lat_inc = 0.25
-  uwind_nlon = 321
-  uwind_nlat = 161
-
-  vwindfile       = 'vwind.bin',
-  vwindstartdate1 = 20020701,
-  vwindstartdate2 = 000000,
-  vwindperiod     = 10800.0,
-  vwind_lon0 = -90.0 
-  vwind_lon_inc = 0.25
-  vwind_lat0 = 20.0
-  vwind_lat_inc = 0.25
-  vwind_nlon = 321
-  vwind_nlat = 161
+  uwindfile        = 'uwind.bin',
+  uwindstartdate1  = 20020701,
+  uwindstartdate2  = 000000,
+  uwindperiod      = 10800.0,
+
+  vwindfile        = 'vwind.bin',
+  vwindstartdate1  = 20020701,
+  vwindstartdate2  = 000000,
+  vwindperiod      = 10800.0,
 
   ! 2-m air temperature (K)
-  atempfile       = 'atemp.bin',
-  atempstartdate1 = 20020701,
-  atempstartdate2 = 000000,
-  atempperiod     = 10800.0,
-  atemp_lon0 = -90.0 
-  atemp_lon_inc = 0.25
-  atemp_lat0 = 20.0
-  atemp_lat_inc = 0.25
-  atemp_nlon = 321
-  atemp_nlat = 161
+  atempfile        = 'atemp.bin',
+  atempstartdate1  = 20020701,
+  atempstartdate2  = 000000,
+  atempperiod      = 10800.0,
 
   ! 2-m specific humidity (kg/kg)
-  aqhfile         = 'aqh.bin',
-  aqhstartdate1   = 20020701,
-  aqhstartdate2   = 000000,
-  aqhperiod       = 10800.0,
-  aqh_lon0 = -90.0 
-  aqh_lon_inc = 0.25
-  aqh_lat0 = 20.0
-  aqh_lat_inc = 0.25
-  aqh_nlon = 321
-  aqh_nlat = 161
-
-  ! shortwave radiation
-  swdownfile      = 'swdown.bin',
-  swdownstartdate1= 20020701,
-  swdownstartdate2= 000000,
-  swdownperiod    = 10800.0,
-  swdown_lon0 = -90.0 
-  swdown_lon_inc = 0.25
-  swdown_lat0 = 20.0
-  swdown_lat_inc = 0.25
-  swdown_nlon = 321
-  swdown_nlat = 161
+  aqhfile          = 'aqh.bin',
+  aqhstartdate1    = 20020701,
+  aqhstartdate2    = 000000,
+  aqhperiod        = 10800.0,
+
+  ! Shortwave radiation
+  swdownfile       = 'swdown.bin',
+  swdownstartdate1 = 20020701,
+  swdownstartdate2 = 000000,
+  swdownperiod     = 10800.0,
 
   ! Longwave radiation
-  lwdownfile      = 'lwdown.bin',
-  lwdownstartdate1= 20020701,
-  lwdownstartdate2= 000000,
-  lwdownperiod    = 10800.0,
-  lwdown_lon0 = -90.0 
-  lwdown_lon_inc = 0.25
-  lwdown_lat0 = 20.0
-  lwdown_lat_inc = 0.25
-  lwdown_nlon = 321
-  lwdown_nlat = 161
+  lwdownfile       = 'lwdown.bin',
+  lwdownstartdate1 = 20020701,
+  lwdownstartdate2 = 000000,
+  lwdownperiod     = 10800.0,
 
   ! Precipitation
-  precipfile      = 'precip.bin',
-  precipstartdate1= 20020701,
-  precipstartdate2= 000000,
-  precipperiod    = 10800.0,
-  precip_lon0 = -90.0 
-  precip_lon_inc = 0.25
-  precip_lat0 = 20.0
-  precip_lat_inc = 0.25
-  precip_nlon = 321
-  precip_nlat = 161
+  precipfile       = 'precip.bin',
+  precipstartdate1 = 20020701,
+  precipstartdate2 = 000000,
+  precipperiod     = 10800.0,
 
   ! Evaporation
-  evapfile      = 'evap.bin',
-  evapstartdate1= 20020701,
-  evapstartdate2= 000000,
-  evapperiod    = 10800.0,
-  evap_lon0 = -90.0 
-  evap_lon_inc = 0.25
-  evap_lat0 = 20.0
-  evap_lat_inc = 0.25
-  evap_nlon = 321
-  evap_nlat = 161
+  evapfile         = 'evap.bin',
+  evapstartdate1   = 20020701,
+  evapstartdate2   = 000000,
+  evapperiod       = 10800.0,
 
   ! Runoff
-  runofffile      = 'runoff.bin',
-  runoffstartdate1= 20020701,
-  runoffstartdate2= 000000,
-  runoffperiod    = 10800.0,
-  runoff_lon0 = -90.0 
-  runoff_lon_inc = 0.25
-  runoff_lat0 = 20.0
-  runoff_lat_inc = 0.25
-  runoff_nlon = 321
-  runoff_nlat = 161
+  !runofffile       = 'runoff.bin',
+  !runoffstartdate1 = 20020701,
+  !runoffstartdate2 = 000000,
+  !runoffperiod     = 10800.0,
+/
+
+&EXF_NML_03
+  ! climatological relaxation fields (climsst, climsss, etc.) -- unused here
+/
+
+&EXF_NML_04
+  ! --- input grid metadata for EXF interpolation (requires ALLOW_EXF_INTERPOLATION) ---
+
+  uwind_lon0   = -90.0,  uwind_lon_inc  = 0.25,
+  uwind_lat0   =  20.0,  uwind_lat_inc  = 0.25,
+  uwind_nlon   = 321,    uwind_nlat     = 161,
+
+  vwind_lon0   = -90.0,  vwind_lon_inc  = 0.25,
+  vwind_lat0   =  20.0,  vwind_lat_inc  = 0.25,
+  vwind_nlon   = 321,    vwind_nlat     = 161,
+
+  atemp_lon0   = -90.0,  atemp_lon_inc  = 0.25,
+  atemp_lat0   =  20.0,  atemp_lat_inc  = 0.25,
+  atemp_nlon   = 321,    atemp_nlat     = 161,
+
+  aqh_lon0     = -90.0,  aqh_lon_inc    = 0.25,
+  aqh_lat0     =  20.0,  aqh_lat_inc    = 0.25,
+  aqh_nlon     = 321,    aqh_nlat       = 161,
+
+  swdown_lon0  = -90.0,  swdown_lon_inc = 0.25,
+  swdown_lat0  =  20.0,  swdown_lat_inc = 0.25,
+  swdown_nlon  = 321,    swdown_nlat    = 161,
+
+  lwdown_lon0  = -90.0,  lwdown_lon_inc = 0.25,
+  lwdown_lat0  =  20.0,  lwdown_lat_inc = 0.25,
+  lwdown_nlon  = 321,    lwdown_nlat    = 161,
+
+  precip_lon0  = -90.0,  precip_lon_inc = 0.25,
+  precip_lat0  =  20.0,  precip_lat_inc = 0.25,
+  precip_nlon  = 321,    precip_nlat    = 161,
+
+  evap_lon0    = -90.0,  evap_lon_inc   = 0.25,
+  evap_lat0    =  20.0,  evap_lat_inc   = 0.25,
+  evap_nlon    = 321,    evap_nlat      = 161,
+
+  runoff_lon0  = -90.0,  runoff_lon_inc = 0.25,
+  runoff_lat0  =  20.0,  runoff_lat_inc = 0.25,
+  runoff_nlon  = 321,    runoff_nlat    = 161,
+/
+&EXF_NML_OBCS
+  useOBCSYearlyFields  = .FALSE.,
+
+  obcsNstartdate1 = 20020701,  obcsNstartdate2 = 000000,
+  obcsNperiod     = 86400.0,   obcsNrepCycle   = 0.0,
+
+  obcsSstartdate1 = 20020701,  obcsSstartdate2 = 000000,
+  obcsSperiod     = 86400.0,   obcsSrepCycle   = 0.0,
+
+  obcsEstartdate1 = 20020701,  obcsEstartdate2 = 000000,
+  obcsEperiod     = 86400.0,   obcsErepCycle   = 0.0,
+
+  obcsWstartdate1 = 20020701,  obcsWstartdate2 = 000000,
+  obcsWperiod     = 86400.0,   obcsWrepCycle   = 0.0,
 /
-  

simulations/glorysv12-curvilinear/input/data.cal

@@ -4,13 +4,12 @@
 # mnc_indir_str='start/',
  mnc_use_outdir=.TRUE.,
  mnc_outdir_str='mnc_',
- mnc_outdir_date=.FALSE.,
+ mnc_outdir_date=.TRUE.,
  monitor_mnc=.FALSE.,
 # snapshot_mnc=.FALSE.,
 # timeave_mnc=.TRUE.,
  pickup_write_mnc=.FALSE.,
  pickup_read_mnc=.FALSE.,
- mnc_filefreq=2635200.0,
-# mnc_filefreq=31622400.0,
+ mnc_filefreq=86400.0,
  &end
 

simulations/glorysv12-curvilinear/input/data.cheapaml

@@ -0,0 +1,41 @@
+# Open-boundaries
+ &OBCS_PARM01
+ OB_Iwest=424*1,
+ OB_Ieast=424*-1,
+ OB_Jsouth=768*1,
+ OB_Jnorth=768*-1,
+ OBCS_monitorFreq=7200.0,
+ useOBCSsponge=.TRUE.,
+ useOBCSprescribe=.TRUE.,
+ useOBCSbalance=.TRUE.,
+ OBCS_balanceFacE=1.0,
+ OBCS_balanceFacS=1.0,
+ OBCS_balanceFacN=1.0,
+ OBNtFile='T.north.bin',
+ OBNsFile='S.north.bin',
+ OBNuFile='U.north.bin',
+ OBNvFile='V.north.bin',
+ OBStFile='T.south.bin',
+ OBSsFile='S.south.bin',
+ OBSuFile='U.south.bin',
+ OBSvFile='V.south.bin',
+ OBEtFile='T.east.bin',
+ OBEsFile='S.east.bin',
+ OBEuFile='U.east.bin',
+ OBEvFile='V.east.bin',
+ OBWtFile='T.west.bin',
+ OBWsFile='S.west.bin',
+ OBWuFile='U.west.bin',
+ OBWvFile='V.west.bin',
+ &end
+# Orlanski parameters
+ &OBCS_PARM02
+ &end
+# Sponge-layer parameters 
+ &OBCS_PARM03
+  spongeThickness=8,
+  Urelaxobcsinner=1800.0,
+  Urelaxobcsbound=600.0,
+  Vrelaxobcsinner=1800.0,
+  Vrelaxobcsbound=600.0,  
+ &end

simulations/glorysv12-curvilinear/input/data.exf

@@ -1,8 +1,10 @@
 # Packages
  &PACKAGES
- useOBCS=.FALSE.,
- useKPP=.FALSE.,
+ useOBCS=.TRUE.,
+ useKPP=.TRUE.,
  useCHEAPAML=.FALSE.,
- useMNC=.FALSE.,
+ useMNC=.TRUE.,
  useDIAGNOSTICS=.FALSE.,
+ useCAL=.TRUE.,
+ useEXF=.TRUE,
  &end