Skip to content

Conversation

@sanAkel
Copy link
Contributor

@sanAkel sanAkel commented Jan 2, 2026

This PR provides what is needed to be done:

That is the following:

  1. As in RTOFS v2.5, use the NCODA output file: ${typec}_${dtg}_0000_analfld where typec=icecov_sfc_1o${SIZN}. and SIZN="${IDM}x${JDM}"
  2. Replace these lines that convert: icecov_${dtg}_analfld ➡️ ssmi.$dtg.r; note CICE4 reads file: ssmi.r with something else that will work with CICE6.

Proof of the ⬆️ statement:

santha.akella@dlogin02:~/my_stmp/test> pwd -P
/lfs/h2/emc/stmp/santha.akella/test
santha.akella@dlogin02:~/my_stmp/test> 
santha.akella@dlogin02:~/my_stmp/test> ls -atlh1
total 33M
drwxr-sr-x 2 santha.akella emc 4.0K Jan  2 02:19 .
-rw-r--r-- 1 santha.akella emc  57M Jan  2 02:19 icecov_20251215_analfld.nc
-rw-r--r-- 1 santha.akella emc  201 Jan  2 02:19 output.log
lrwxrwxrwx 1 santha.akella emc   78 Jan  2 02:19 depth_GLBb0.08_09m11ob2_mom6.nc -> /lfs/h2/emc/ptmp/santha.akella/RTOFS_GLO//fix//depth_GLBb0.08_09m11ob2_mom6.nc
lrwxrwxrwx 1 santha.akella emc  156 Jan  2 02:19 icecov_20251215_analfld -> /lfs/h2/emc/couple/noscrub/zulema.garraffo/COMDIR/prod/com/rtofs/v2.5/rtofs.20251216//ncoda/hycom_var/restart/icecov_sfc_1o4500x3297_2025121500_0000_analfld
-rwxr-xr-x 1 santha.akella emc  924 Jan  2 02:19 test.sh
drwxr-sr-x 3 santha.akella emc 4.0K Jan  1 20:45 ..
santha.akella@dlogin02:~/my_stmp/test> 
santha.akella@dlogin02:~/my_stmp/test> cat output.log 
STATUS: All inputs validated. Beginning conversion for icecov...
Summary for sic: Min=0.0000, Max=1.0000
SUCCESS: Created NetCDF [sic] at /lfs/h2/emc/stmp/santha.akella/test/icecov_20251215_analfld.nc
santha.akella@dlogin02:~/my_stmp/test> 
santha.akella@dlogin02:~/my_stmp/test> ncdump -h icecov_20251215_analfld.nc 
netcdf icecov_20251215_analfld {
dimensions:
        ny = 3297 ;
        nx = 4500 ;
        time = 1 ;
variables:
        float sic(ny, nx) ;
                sic:_FillValue = -9.99e+33f ;
                sic:units = "fraction" ;
                sic:long_name = "Ice Concentration" ;
                sic:scale_factor_applied = 0.01 ;
        int64 time(time) ;
                time:units = "days since 2025-12-15 00:00:00" ;
                time:calendar = "proleptic_gregorian" ;
}

cat test.sh:

#!/bin/sh

export DATAROOT=/lfs/h2/emc/stmp/santha.akella/test/
export jobid=icecov_sic
export HOMErtofs=/lfs/h2/emc/ptmp/santha.akella/RTOFS_GLO/
export USHrtofs=${HOMErtofs}/ush/
export FIXrtofs=${HOMErtofs}/fix/
export PDY=20251216
export dtg=20251215

# py module-- JOB card needs to have it -- Dan --
source $HOMErtofs/scripts/load_py_modules.sh

CONMIN=/lfs/h2/emc/couple/noscrub/zulema.garraffo/COMDIR/prod/com/rtofs/v2.5/rtofs.${PDY}/
IDM=4500
JDM=3297
var_name="icecov"
#var_name="icethk"
#var_name="icetmp"
#var_name="mixlyr"

SIZN="${IDM}x${JDM}"
typec=${var_name}_sfc_1o${SIZN}

cd ${DATAROOT}

# link to topography
ln -f -s ${FIXrtofs}/depth_GLBb0.08_09m11ob2_mom6.nc depth_GLBb0.08_09m11ob2_mom6.nc

# input file
ln -sf $CONMIN/ncoda/hycom_var/restart/${typec}_${dtg}00_0000_analfld ./${var_name}_${dtg}_analfld

# run
$HOMErtofs/ush/rtofs_glo2d_ice.sh depth_GLBb0.08_09m11ob2_mom6.nc ${var_name}_${dtg}_analfld ${var_name} ${DATAROOT}
err=$? 
echo " error from ${HOMErtofs}/ush/rtofs_glo2d_ice.sh =" $err

hence it resolves #99

@sanAkel sanAkel added the Ready to merge Tested and works label Jan 7, 2026
@sanAkel
Copy link
Contributor Author

sanAkel commented Jan 13, 2026

@DanIredell-NOAA Can you please DO NOT merge yet!

As you said, merging @zgarraffo's incremental update needs to be reconciled in a top level script: scripts/exrtofs_glo_ncoda_inc.sh (GitHub has identified.)

@sanAkel sanAkel removed the Ready to merge Tested and works label Jan 13, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Handle sea ice analysis field: NCODA ➡️ CICE6

2 participants