================================================================================
E3SM Initial Condition Remapping Workflow
================================================================================
Case name: v3.LR.historical_0321
Timestamp: 2015-01-01-00000
Source grid: ne30np4
Target grid: northamericax4v1np4
Vertical levels: L80
Keep intermediate files: yes
================================================================================

Step 1: Horizontal remapping...
  Input: v3.LR.historical_0321.eam.i.2015-01-01-00000.nc
  Output: v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc
  Weights: map_ne30np4_to_northamericax4v1np4_esmfbilin.20250813.nc
  ✓ Horizontal remapping completed

Step 2: Creating vertical file with topo-adjusted PS...
  Step 2.1: Copying vertical coordinate template...
  ✓ Vertical coordinate file copied

  Step 2.2: Remapping PHIS...
    Source: /global/cfs/cdirs/e3sm/inputdata/atm/cam/topo/USGS-gtopo30_ne30np4pg2_x6t-SGH.c20210614.nc
    Output: PHIS_d.ne30np4_mapped_northamericax4v1np4.nc
  ✓ PHIS remapping completed

  Step 2.3: Adjusting PS for topography...
    IC file: v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc
    Target topo: USGS-gtopo30_northamericax4v1np4pg2_oroshp_x6t.c20250813.nc
    Mapped PHIS: PHIS_d.ne30np4_mapped_northamericax4v1np4.nc
Vert file: E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc
Raw IC file: v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc
Topo file: USGS-gtopo30_northamericax4v1np4pg2_oroshp_x6t.c20250813.nc
Mapped PHIS file: PHIS_d.ne30np4_mapped_northamericax4v1np4.nc
PHIS varname: PHIS_d

Reading topography file...
    Reading PHIS_d: shape=(130088,), dims=('ncol_d',)
    Reading lat: shape=(57816,), dims=('ncol',)
    Reading lon: shape=(57816,), dims=('ncol',)
  PHIS shape: (130088,)

Reading initial condition file...
    Reading PS: shape=(1, 130088), dims=('time', 'ncol_d')
    Reading hyam: shape=(80,), dims=('lev',)
    Reading hybm: shape=(80,), dims=('lev',)
    Reading hyai: shape=(81,), dims=('ilev',)
    Reading hybi: shape=(81,), dims=('ilev',)
    Reading T: shape=(1, 80, 130088), dims=('time', 'lev', 'ncol_d')
    Reading lat: shape=(130088,), dims=('ncol_d',)
    Reading lon: shape=(130088,), dims=('ncol_d',)
  PS shape: (1, 130088)
  T shape: (1, 80, 130088)
  hybm shape: (80,)

Reading mapped PHIS file...
    Reading PHIS_d: shape=(130088,), dims=('ncol_d',)
  Mapped PHIS shape: (130088,)

Number of vertical levels: 80
Original PS shape: (1, 130088)
Processed PS shape: (130088,)
Original T shape: (1, 80, 130088)
Processed T shape: (130088,)

Final data shapes:
  PS: (130088,)
  T: (130088,)
  PHIS_src: (130088,)
  PHIS_tgt: (130088,)

Physical constants:
  gravit: 9.80616
  rair: 287.04
  xlapse: 0.0065
  alpha: 0.1902640789055043
  Bottom sigma (hybm[-1]): 0.998161220860788

Bottom level temperature (tbot) stats:
  Min: 226.867411 K
  Max: 311.325950 K
  Mean: 278.622193 K

Original PS stats:
  Min: 52670.216724 Pa
  Max: 104625.128447 Pa
  Mean: 98256.024108 Pa

Bottom level pressure (pmid_bot) stats:
  Min: 52573.367828 Pa
  Max: 104432.745943 Pa
  Mean: 98075.352981 Pa

Derived surface air temperature (tsair) stats:
  Min: 226.946928 K
  Max: 311.435069 K
  Mean: 278.719849 K

PHIS_src (mapped) stats:
  Min: -258.905748 m²/s²
  Max: 51189.316409 m²/s²
  Mean: 2720.446500 m²/s²

PHIS_tgt (topo) stats:
  Min: -262.130223 m²/s²
  Max: 51213.600855 m²/s²
  Mean: 2720.820523 m²/s²

PHIS difference (phis1 - phis2) stats:
  Min: -13815.897614 m²/s²
  Max: 7844.537865 m²/s²
  Mean: -0.374023 m²/s²

Adjusted PS stats:
  Min: 52652.024757 Pa
  Max: 104625.128447 Pa
  Mean: 98273.585864 Pa

PS change stats:
  Min: -13808.091929 Pa
  Max: 9684.281871 Pa
  Mean: 17.561756 Pa
  Std: 1127.252868 Pa

nlev and ncol: 80, 130088

Existing dimensions in E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc:
  lev: 80 
  ilev: 81 
Creating ncol dimension with size 130088
Creating time dimension (unlimited)
Creating new PS variable with dimensions: ('time', 'ncol')
Writing PS with shape (1, 130088) to dimensions (time, ncol)
Successfully wrote adjusted PS to E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc
  ✓ PS adjustment completed

Step 3: Vertical remapping...
  Using vertical file: E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc
  Output: v3.LR.historical_0321_mapped_northamericax4v1np4-topoadj.eam.i.2015-01-01-00000.nc
Input #00: /global/cfs/cdirs/e3sm/www/tang30/remap_atmIC_m2m/v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc
Vertical : E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc
ncks: INFO nco_ntp_vrt() temperature extrapolated toward/into surface assuming constant moist adiabatic lapse rate = 6.5 K/(100 hPa) for variable T
  ✓ Vertical remapping completed

Step 4: Copying adjusted PS to final file...
  ✓ PS copied to final file

Step 5: Keeping intermediate files (as requested)
  Intermediate files preserved:
    - v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc (horizontally remapped)
    - E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc (vertical coordinate with adjusted PS)
    - PHIS_d.ne30np4_mapped_northamericax4v1np4.nc (remapped PHIS)

================================================================================
SUCCESS: All steps completed
================================================================================
Final output file: v3.LR.historical_0321_mapped_northamericax4v1np4-topoadj.eam.i.2015-01-01-00000.nc

Intermediate files:
  - v3.LR.historical_0321_mapped_northamericax4v1np4.eam.i.2015-01-01-00000.nc (horizontally remapped)
  - E3SM_vert_coor_L80_northamericax4v1pg2.2015-01-01-00000.nc (vertical coordinate with adjusted PS)
  - PHIS_d.ne30np4_mapped_northamericax4v1np4.nc (remapped PHIS)
================================================================================